CT投影视野受限图像重建.pdfVIP

CTCT 投影投影视野受限视野受限的图像重建的图像重建 CTCT 投影投影视野受限视野受限的的图像重建图像重建 赵双任 多伦多 DTI 成像行 shrzhao@126.com 扬新铁 西北工业大学 yangx@ 摘要：CT(Computer Tomography)设备的探头宽度有限可引起探头视野 受限 (Limited Field of View, LFOV)， 由此产生的投影数据被截 断。由被截断的投影数据进行图像重建产生截断膺像 (Truncation Artifacts)。外插算法 (Extrapolation)可以减小截断膺像。但此种算 法常常会矫枉过正或校正不足。传统的重建重投影 (Reconstruction Reprojection)迭代算法适合视角 (Limited Angle of View) 受限的情 况，但对视野受限情况并没有明显的优势。本文为 CT 投影视野受限的 图像重建建立一个理论基础。由此可以对出一个新的迭代算法。 1.引言1.引言 1.1.引言引言 常用 CT 探头的宽度为 40cm。探头视野大约为 26cm 取决于 CT 系统的摄 像的缩小率。常用 CT 系统采用扇形和锥形光束。探头视野是 CT 系统 旋转中心处探头的宽度的对应值。此处缩小率指光源到 CT 系统旋转中 心的距离与光源到探头的距离的比值。病人的身体一般大于 26cm，因 此此情况下探头视野小于病人的身体。此种情况称为探头视野受限。 探头视野受限会引起由探头测量所得到的投影数据在探头边缘处被截 断。此处截断是指投影数据在探头边缘处以外突然为零。 FPB 算法 [1]常常用于用扇形和锥形光束 CT 投影数据的图像重 建。若将 FPB 算法直接用于截断的投影数据，重建的图像将包含截断 膺像。截断膺像是在感兴趣区域边缘处的亮环。此处感兴趣区域即要 进行图像重建的区域，此区域为圆形，正好落在所有 CT 投影的视野 内。CT 视野以外区域重建的图像质量很差，一般不予以考虑。感兴趣 区域的直径应小于 CT 视野的宽度，即感兴趣区域加上一个环状边缘区 域所构成区域的直径正好为 CT 视野宽度相同。文献 [2]采用的边缘宽 度为 10 个像素点。但截断膺像主要集中在感兴趣区域的边缘处，增加 环状区域的宽度只是选择更容易的问题去求解，而没有解决问题本 身。因此本文假定上述环状边缘区域的宽度为零，即感兴趣区域的直 1 径都与 CT 视野的宽度正好相同。 消除截断膺像的常用算法为外插算 法[3,4]。外插算法估计投影视野外丢失了的投影数据。该算法的结果 依赖于用于外插的函数，例如，常数外插，线性函数外插，多项式韩 数外插，COS 函数外插，指数函数等等。因此该算法重建的图形常常被 矫枉过正，或矫正不足。如果矫枉过正，感兴趣区域边缘处的亮环变 成暗环。如果矫正不足此亮环有所减弱，但仍不消失。小波 (Wavelet) 算法也常被用来减轻截断膺像。但是小波算法常常接合外插算法一起 使用 [2,3]。文献[3]承认在所谓小波算法中对截断膺像起作用的还是 外插算法。已知先决条件的方法被用来解决探头视野受限的问题[5]。 如果病人的身体落在感兴趣区域外只是一小部分比如两个胳膊， Hibert 变换算法可给出解析解 [6]。局部断层扫描(LocalTomography) 算法[8,9]可完全消除截断膺像。但该算法在消除截断膺像的同时图像 凹陷。 图像重建重投影的迭代算法 [9,10]给出了对视角受限(Limited angle of view)问题的计算结果并提到该算法对探头视野受限问题也 有效，但没有给出这方面的计算结果。利用此算法至今也没有看到比 外插算法更好的计算结果。文献 [11]将该算法用于解决偏移探头 (Offset Detector)的投影数据截断问题。此种情况下，对应投影一个 扫描角度，探头有两个位置。探头的视野由两