CT图像重建实验中伪影的分析与消除.pdfVIP

修改稿 CT 图像重建实验中伪影的分析与消除 吴晓锋 ,刘力 （中国科学院高能物理研究所，北京100039） 摘要：本文通过对扇束 CT （Computed Tomography ）实验数据分析，用 FBP-Ramp 重建法对不同采集条件下投影 数据进行重建，分析了各种伪影与采集条件及重建参数的关系，指出进一步克服这些伪影的方法。此外对 图像质量与采集条件的关系也进行了分析。 关键词：CT 断层扫描；图像重建；伪影 中图分类号： R445.2 TP751 文献标识码：A 文章编号： 前言：CT技术 （Computed Tomography ）已广泛应用于人体断层扫描和工业物品的断层扫描[1] 。 图 1 为一个典型的第三代扇束CT工作原理示意图，其中X光源S发射一个扇面X光束扫过被测物体 所在断层区域（即FOV ，视野），弧线形等角间距的探测器阵列D接收并记录经过物体衰减后的X 光计数，通常由光源和探测器围绕FOV 中心旋转测量 360o构成对该物体断层的一次完整测量。由 这些投影数据可重建出该物体断层中的质量密度分布 μ（x,y ）。 D S 2γm FOV 图 1 扇束 CT 工件原理示意图， 2γ 为扇束张角 m I ( , ) β γ ln ( , ) 成像方程为 ∫µ x y dl ( , ) I o β γ L 其中 β为扇束位置（角度）指标， γ为扇束中等角间距射线束指标，对应于不同探测器小条 （Bin ）;I( β, γ)为实际测量的射线强度，I ( β, γ) 为空扫描得到的射线强度，积分表示沿该射线 0 的质量密度和。 图像重建方法有滤波反投影法（Filtered Back Projection, FBP ）、OSEM法（Ordered-subset Expectation Maximization ）等多种算法。目前关于各种重建算法的研究已开展大量工作[2] 。对重建 中出现的各种伪影的分析及解决方法的研究对提高重建图像质量也是至关重要的[3-4] ，因为伪影可 导致临床诊断图象出现假阳性和假阴性[5] 。本文通过对一个扇束CT原理性实验的数据进行分析， 用FBP-Ramp重建法对不同采集条件下投影数据进行重建，并分析了相应伪影及其与采集条件及 重建参数的关系，指出进一步克服这些伪影的方法。此外对图像质量与采集条件的关系也进行了 分析。 收稿日期： 吴晓锋（1978--），男（汉族），河南汝阳人，硕士研究生，现从事图象重建算法研究。 1 扇束 CT 实验及数据处理方法 实验用被测物体为一个中等质量密度的圆盘，其中偏心放置一高密度的圆形螺拴，见示意图 2 。用一平面扇束X 光源进行 完全覆盖扫描（见图 1），其中紧密 放置 1024 （即n γ＝1024）个探 测器单元（bin ）,形成对物体的 等角间距投影扫描，整个系统 围绕物体旋转3600进行n β次 等角距扫描，每扇束最大张角为: 2 γ =0.008 弧度， m 重建图像矩阵为 512 ＊512， 用 FBP-Ramp 进行重建， 平均计数约为 3054 光子/bin 。 图2 被测物模型 参数n n 构成了本实验的测量条件；n 反映了径向采样密度，n