CT图像振动伪影的计算机仿真.pdf

中国体视学学会cT理论与应用分会、中国工程物理研究院 2008年7月 应用电子学研究所、1(CT理论与应用研究》编辑部论文集 CT图像振动伪影的计算机仿真 李英辅1，于新波2，赵文仓1’孔4，杨殿才3 1．青岛科技大学认知与信息科学研究所，山东青岛266042 2．青岛大学医学院附属医院口腔科，山东青岛266003 3．青岛高校软控股份有限公司青岛市工业信息化重点实验室，山东青岛266042 4．西安交通大学电子与信息工程学院，西安710049 摘要：计算机断层成像技术(cT)在实际应用中，由于软、硬件的特性，会受到很多干扰影响， 重建图像中往往会出现伪影，而机械振动导致的条状伪影是比较常见的一种形式。为能从根源 上削弱或去除这种伪影对结果的影响，本文在分析cT技术和图像重建算法等原理的基础上，针 对机械振动伪影进行仿真，得到机械振动伪影的特征和变化趋势，利用投影、图像的特征有效 地帮助判断该种伪影的成因、影响等信息。 关键词：条状伪影；机械振动；仿真 1 引言 在重建算法的研究中，模拟投影数据是必不可少的一部分，一方面可以避免实测投影 数据的误差，另一方面也可以根据需要有选择地模拟失真的投影数据。黄建林、吕东辉Il】等 人通过对图像重建中x射线投影模拟的常用方法的研究一离散法、矢量法、解析法，得出 结论：矢量法和离散法一般用于三维状态下，解析法和矢量法只能对规则形状物体及密度 分布均匀的物体产生精确的投影。刘泽[21等人研究的“基于3D Shepp-Logan头部模型的三 维医学图像重建仿真”给出一种以3DShepp—Logan头部模型作为三位医学图像重建领域进 行仿真实验和算法性能评测的基本参考模型，并以此模型进行断层扫描仿真、二维以及三 维重建，最后得到三维图像。然而目前类似的研究的模拟没有考虑到射线的散射和硬化等 现实情况下的干扰现象，而本文所做的仿真与重建实验是关于现实情况下常见干扰因素之 一一机械振动对投影的影响，仿真和重建研究是针对机架在旋转平面内的振动导致的CT 图像伪影，可用于投影仿真过程干扰影响的研究。本文结合CT原理与算法等理论基础，基 的成因等特性进行CT图像重建，以期得到振动伪影明显削弱或消除的CT图像。由于伪影 成因的常见性，该种cT图像伪影相应的仿真、重建结果可以广泛的应用到CT图像优化中， 在一定程度上可以有效的改善CT图像伪影，同时可以辅助判断图像伪影成因等。有效地模 拟干扰条件作用于模型上得到的投影能够描绘出该干扰条件对物体在实际情况中的影响， 降低成本，提高效率。 基金项目：山东省自然科学基金(Y2007G32)和山东省博士后创新项目专项资金资助项目资助． 通讯作者；赵文仓，理学博士。研究领域：图像处理与模式识别、认知信息学、cT理论与应用研究等， zhaocenter一,journal@yahoo．COm．cn· 李英辅等：cT图像振动伪影的计算机仿真 2基础理论知识 计算机断层重建依赖于不同角度X射线束流测量结果，以形成图像。假设输入X光子 是单能的，在均匀材料入射和出射两边测量的X射线强度，可以通过线积分来描述，如式 I=／oe一肛 (1) 式中厶是入射X射线强度，，是出射X射线的强度，工是厚度，∥是材料的线性衰减系 数。对于现实中非均匀物体(即物体由多个具有不同衰减系数的材料组成)，其总体的衰减 系数可以通过将物体分割成小单元来计算。当单元足够小时可以看作是均匀的物体。被测 物体在线积分测量中，射线透射过待测物体，设物体各段线性衰减系数为“，鲍，鸬，……， 相应长度为毛，乇，毛，……，可得 以=m(Io／／) (2) 物体在矽平面内部都不均匀，即衰减系数∥=∥似力，则在某一方向Z沿吕东辉Ⅲ、刘 泽啪某一路径三的总衰减为式 I∥比=In(厶／D (3) 其中f却为投影。可得到投影p的公式