高能X光照相中密度重建的共轭梯度算法.pdfVIP

维普资讯 计 算 物 理 第23卷 第2期 CHINESEJOURNALOFCOMPUTATIONALPHYSICS Vo1．23，No．2 2OO6年 3月 Mar．．2006 [文章编号】 1001．246X(2006)02．0144．07 高能X光照相中密度重建的共轭梯度算法 许海波 ， 魏素花 (北京应用物理与计算数学研究所，北京 100088) [摘 要】 通过对高能X光照相过程的物理分析，得到了辐射照相的点扩展函数 ，给出了估价函数的表达式． 并以典型的法国试验客体模型为例，将共轭梯度算法应用到高能x光照相的密度重建中，得到了较满意的结果． 这种方法从模拟辐射照相过程出发，将辐射照相的模拟结果与实际测量结果做比较并找出最大可能解 ．可在一 定程度上克服其它密度重建算法通过反卷积过程消除模糊效应所引起 的不确定性 ． [关键词】 X光照相；点扩展函数；估价函数 ；共轭梯度算法；密度重建 [中图分类号】 0443 [文献标识码】 A O 引言 用物体 的横截面投影的数集来重构物体的三维图像是一种有用的处理过程 ，这种横截面投影是在不损 伤物体本身的条件下获得的，例如 x光照相 ．因此这种利用横截面投影图集来重建三维图像的技术广泛地 应用于放射医学、电子显微镜、射电天文学等领域 ．投影的信息收集模式有透射、发射及反射 ．不同的信息收 集模式具有不同的重建模式 ． 高能X光照相的信息收集模式属于透射模式，是利用X射线的强穿透能力，将物体投影到平面上，也就 是将三维物体压缩成二维图像 ．由于照相是物体各层次投影的重叠(或者说投影时丢失了沿射线方向的分辨 能力)，难以从投影图像中看出物体 内部的结构 ，而且，照相过程会引起图像的退化 ．图像处理就是将退化 的 图像恢复，包括图像 的去噪音 、去模糊 、边界探测等技术 ．图像重建就是利用一组物体的投影重建三维物体图 像 ，重建过程由计算机完成，通常称为计算机层析成像技术 (简称 CT)．主要的重建方法有傅立 叶变换重建 、 滤波反投影重建 、代数法重建等… ． 实际的图像包含噪音、模糊、散射的影响．因此，在利用传统的密度重建方法求密度之前，必须消除这些 因素的影响．整个过程如图 1所示 ．由于高能X照相中的图像模糊很严重 ，消模糊过程容易引起振荡现象，复 原难度很大 ． 在重建过程中尽可能考虑到闪光照相的实际物理过程 ，以提高层析重建的精度 ，而且为了克服消模糊过 程引起的振荡现象，美国洛斯-阿拉莫斯国家实验室 (u L)开发了称为HARDES的射线径迹程序，并与约束 共轭梯度优化算法耦合，进行闪光照相的密度重建 ．本文将共轭梯度算法应用到高能X光照相的密度重 建 中，整个过程如图2所示．我们采用了标准的共轭梯度程序并做了修改，可以根据实际需求加入一定的约 由r·--·--—一 ———【———————————————．—．——一 I度【算法I ．．．．． ．．．．． ． 图l 传统的密度重建算法流程图 图2 基于数值模拟的密度重建算法流程图 Fig．1 Flow char