图像重建的分块迭代算法.pdfVIP

图像重建的分块迭代算法 1 2 赵双任 ，扬新铁 1 多伦多 DTI 成像行 2 西北工业大学 E-mail：shrzhao@126.com 摘 要： 图像重建和重新投影 (Reconstruction and Re-projection)的迭代算法 已在投影视野受限 （Limited Field of View, LFOV ）的情下使用。该情况下投影数据被截断。此处投影采用的平行光 束，扇形光束或锥形光束。该迭代算法对消除由于投影数据某处被截断所产生的截断膺像(Truncation Artifacts)非常有效。最近发现该算法不仅对消除截断膺像有效而且对消除条纹状膺像也有一些效果。 条纹状膺像是在投影数据未被截断情况下由滤波反投影(Filterred Backprojection, FBP )算法所产生 的膺像。条纹状膺像与一个完整扫描所含的投影数及一个投影内的点数有关，该数目越小膺像就越严 重。此发现引起了如下猜想，即该迭代算法或许也适用于投影数据未被截断的情况。本文对该迭代算法 也做了进一步改进，使它更有利于减小条纹状膺像。主要的改进是把图像划分成较小的区域，以及剔除 原迭代算法中用于数据被截断时的外插(Extrapolation)技术。本文投影数据采用平行光束。迭代算法 包含两个或者更多个图像重建过程。第一个过程采用 FBP 算法。第二次图像重建过程中第一次重建的图 像被切成小块；每一小块图像都被重新投影并再次重建；最终合并成一个完整图像。分块迭代算法的计 算结果与 FBP 算法的以及传统的迭代法的结果做了比较。 关键词：膺像，断层扫描，迭代，图像重建，投影，视野，截断，平行光束，CT 一.简介 对平行光束扇形光束和锥形光束的CT（Computer Tomography）系统已经有许多图像重 建算法[1-7]。如果探头宽度足够大，物体全部落在CT 探头的视野之内，测量所得到的投影 数据是完整的。相反如果物体的一部分落在CT 探头的视野之外，测量所得到的投影数据在 探头的边界处被截断。由截断了的投影数据所得的重建图像包含截断膺像(Truncation Artifacts) 。 将 局 部 断 层 扫 描 (Local Tomography) 方 法 [8-12] ， 图 像 重 建 再 投 影 (Reconstruction and Reprojection)的方法[13－16]和外差(Extrapolation)方法[17－19]相结合, 文献[20]研究了探头的视野受限，投影数据被截断的情况下的图像重建问题,并建立了一个新 颖的迭代算法。该算法将物体分成两部份。一部分为感兴趣的区域，另一部分是感兴趣之 外的区域。感兴趣的区域刚好全部落在探头视野之内，该区域的图像重建是探头视野受限 情况下所关注的。该迭代算法由至少两次图像重建构成。第一次采用FBP方法包括对在探头 之外未知投影数据进行外差。用FBP算法得到感兴趣区域之外部份的图像。将该部分图像再 次投影。将原投影数据与再次投影数据相结合。使用由此所得的投影数据对感兴趣的区域 进行第二次图像重建。该算法对减少截断膺像非常有效。最近发现，其实该算法也对投影 数据未被截断的情况起作用。跟据这一发现，本文将该迭代算法进一步改造，使之更适应 投影数据未被截断的情况。原算法中适用投影数据被截断部份比如外差方法被剔除。原算 法中感兴趣的区域为圆形，改造后成成为方形。原算法中只有一块感兴趣的区域。改造后 图像所在区域被分成许多子区域。这些子区域都是感兴趣的区域。迭代算法被用于每一个 子区域。子区域内的图像最终被拼接成整块图像。由此构成所谓分块迭代算法。与FBP方法 比较，分块迭代算法的抑制膺像作用很显著。但因为其他迭代算法也对膺像有一定抑制作 用。本文还将分块迭代算法与一种传统迭代算法进行了比较。此处提到传统迭代方法是一 种常见的线性方程迭代解法。尽管传统迭代方法对抑制膺像也有贡献，但分块迭代算法的 抑制膺像作用要更上一层。 二.传统的图像重建迭代算法