多金属物体CT图像的金属伪影校正-物理学报.PDF

物理学报 Acta Phys. Sin. Vol. 63, No. 5 (2014) 058702 多金属物体CT 图像的金属伪影校正∗ 魏星 闫镔 张峰 李永丽 席晓琦 李磊 1)(信息工程大学信息系统工程学院, 郑州 450002) 2)(河南省人民医院, 郑州 450002) ( 2013 年10 月8 日收到; 2013 年11 月13 日收到修改稿) 针对多金属伪影的校正问题, 本文通过仿真实验分析了多金属伪影的成因, 并提出了一种基于投影校正 的多金属伪影校正方法. 该方法首先直接从投影域分割出金属区域, 然后建立对金属区域投影值的校正模型, 最后通过调整模型参数达到校正目的. 模型以重建图像的灰度熵为目标函数, 采用单纯形法迭代求解使熵最 小时的校正参数. 仿真和实际数据的实验结果表明, 本文算法对多金属伪影的校正起到了良好的效果, 且校 正后的图像质量优于插值校正法. 关键词: 计算机断层成像, 金属伪影, 投影校正, 灰度熵 PACS: 87.59.–e, 87.59.Fm DOI: 10.7498/aps.63.058702 [11] 中, Abdoli 等 提出对插值投影和原始投影进行 [12] 1 引 言 加权的插值方法, Jeong 等 在此基础上加入了对 重建图像的TV 滤波来平滑带状伪影. Zhao[13] 和 计算机断层成像(computed tomography, CT) Kratz[14] 等分别提出基于小波变换和非均匀傅里 过程中, 被扫描物中类似金属的高衰减物质会导致 叶变换的插值运算法. 插值校正法具有原理简单、 重建图像中产生金属伪影[1] , 对CT 系统的成像质 校正速度快等优点, 但在多金属情况下仍存在以 量有严重的影响. 当被成像物中仅有单个金属时, 下不足:1) 插值法需要分割重建图像中的金属区域, 金属伪影表现为其周围存在放射状伪影; 当金属物 但金属周围的白色带状伪影会造成金属区域分割 数量较多且金属区域面积较大时, 金属伪影不仅包 的不准确, 从而影响最终的校正结果; 2) 由于分割 含金属物周围的放射状伪影, 而且在金属连线方向 出的金属区域被直接填补到插值处理后的重建图 上还存在带状伪影[2] . 研究金属伪影的校正方法, 像中, 因此金属内部由硬化效应导致的杯状伪影无 特别是多金属情况下的校正方法一直是CT 领域的 法得到消除; 3) 若非金属区域被金属区域包围, 非 难点与热点问题之一. 金属区域会被错误的归为被插值对象, 致使其内部 现阶段的金属伪影校正方法主要分为三类: 插 的有用信息无法得到恢复. 值校正法、迭代校正法和混合校正法[3] . 迭代校正 针对多金属情况下的金属伪影校正问题, 本文 法[4−6] 和混合校正法[7−9] 对单个金属产生的金属 首先对多金属情况下的金属伪影进行了仿真, 进而 伪影有较好的校正效果, 但它们均基于迭代重建算 提出了一种基于投影校正的多金属伪影校正方法, 法, 具