乘代数重建法和最大熵重建法比较研究-光子学报.PDF

第27 卷第5 期 光 子 学 报 Vo127 No. 5 1998 年 5 月 ACTA PHOTONICA SINICA May 1998 乘代数重建法和最大熵重建法比较研究 江少恩 刘忠礼 郑志坚 唐道源 ( 中国工程物理研究院核物理 化学研究所, 高温高密度等离子体物理重点实验室, 成都市525 信箱77 分箱, 610003) 本文从最大熵原理引入拉格朗日因子推导出 代数重建法( MART ) 和最大熵重建法 ( MENT) , 同时他们是等价的, 只是迭代方法不同而已. 而且用实例进行模拟计算, 计算结果也 表明这两种方法收敛情况和重建精度完全一致, 从而进 一步证明MART 和MENT 的等价性. 最大熵原理; 代数重建法( MART) ; 最大熵重建法( MENT) 0 CT 技术( Com uted Tomogra hy) 原则上可以通过多个不同方位的二维投影重建出物体某种量的 三维分布. 这使得CT 技术在物理研究、医学和工业无损检测等领域中得到广泛的应用. 但在物理研究 中的CT 重建问题 医学中的重建问题相比, 它有不同的特点: 物理研究中经常遇到的问题是 非完全 数据, 这指的是待测物体中的投影方位角受限, 或者投影数目有限. 这是由于物理过程一般来说时间 短, 不象医学CT 那样旋转扫描器, 只能固定多个方位的探头( 探测器) 在物理过程发生时同时测量. 因 为考虑到成本因素, 只能用少数几个探头; 另外物理研究中可能待测的区域很小, 使观测方位受限, 这也 使得只能在实验中用少量的探头. 因此在物理研究中的重建数据常常为 非完全数据. 如在惯性约 1~ 3 束聚变( ICF) 物理实验过程就是这种情况 , 它的时间过程短( 纳秒级) , 待测区域小( 靶球为百微米级) . 对于 非完全数据的重建问题, 迭代方法比傅里叶变换法( 包括卷积反投影法) 的重建效果好得 4 5 6 7 多迭代重建法主要有ART ( 代数重建法) 、MART ( 乘代数重建法) 和 MENT ( 最大熵重建法) 三 种. 对ART, 有多种变化形式, MART 和MENT 的变化形式则较少. 本文通过最大熵原理对MART 和 MENT 进行比较研究. 在第一部分, 由最大熵原理推导MART 和MENT 的迭代形式解, 从而说明它们 的解的等价性. 在第一部分, 由最大熵原理推导MART 和 MENT 的迭代形式解, 从而说明它们的解的 等价性. 第二部分, 通过实例计算说明这两种方法的等价性. 在文的最后作一总结. 1 CT 技术就是从多幅二维投影P, ( u, v ) 重建出三维分布 f ( x , y, z ) 的技术. CT 技术分为两大类重建算法: 一类是傅里 4 8 叶变换重建法 , 另一类为代数重建法 , 代数法一般为迭代过 程. 变换重建法适用于投影数目较多的情况, 如医学CT 扫描 仪, 可达100 多个投影. 而代数重建法则适用于投影数目较少 的情况. 在 ICF 实验中, 靶尺寸很小, 观测方位角有限, 且时间 过程短, 不可能在激光 靶互作用过程中将探头( 针孔相机) 作 旋转扫描, 因此需要同时安装多个探头对X 光 成象, 因此, 在ICF 实验中, 一般采用代数重建法. CT 技术一 般来说, 探头越多, 重建精度越高. 但由于观测方位角有限, 图1 投影坐标示意图 收稿日期: 1997 07 02 Fig. 1 Geometry of the coordinate system 5 期 江少恩等. 乘代数重建法和最大熵重建法比较研究 407 且要涉及成本因素, 要求尽量少用探头