多金属物体CT图像的金属伪影校正.pdf

物 理 学报 Acta Phys．Sin．Vo1．63，No．5(2014) 058702 多金属物体CT图像的金属伪影校正术 魏星 ) 闫镔 ) 张峰 ) 李永丽2) 席晓琦 ) 李磊 )十 1)(信息工程大学信息系统工程学院，郑州 450002) 2)(河南省人民医院，郑州 450002) (2013年 l0月8日收到；2013年 l1月13日收到修改稿) 针对多金属伪影的校正问题，本文通过仿真实验分析了多金属伪影的成因，并提出了一种基于投影校正 的多金属伪影校正方法．该方法首先直接从投影域分割出金属区域，然后建立对金属区域投影值的校正模型， 最后通过调整模型参数达到校正目的．模型以重建图像的灰度熵为目标函数，采用单纯形法迭代求解使熵最 小时的校正参数．仿真和实际数据的实验结果表明，本文算法对多金属伪影的校正起到了良好的效果，且校 正后的图像质量优于插值校正法． 关键词：计算机断层成像，金属伪影，投影校正，灰度熵 PACS：87．59．_e．87．59．Fin DOI：10．7498／aps．63．058702 1 引 言 计算机断层成像 (computed tomography，CT) 过程中，被扫描物中类似金属的高衰减物质会导致 重建图像中产生金属伪影【1]．对CT系统的成像质 量有严重的影响．当被成像物中仅有单个金属时， 金属伪影表现为其周围存在放射状伪影：当金属物 数量较多且金属区域面积较大时，金属伪影不仅包 含金属物周围的放射状伪影，而且在金属连线方向 上还存在带状伪影 [2]．研究金属伪影的校正方法， 特别是多金属情况下的校正方法一直是CT领域的 难点与热点问题之一． 现阶段的金属伪影校正方法主要分为三类：插 值校正法、迭代校正法和混合校正法 [3】．迭代校正 法 [4-6】和混合校正法[7-9】对单个金属产生的金属 伪影有较好的校正效果，但它们均基于迭代重建算 法，具有计算量大、实用性差等不足．相比之下，对 插值环节进行改进并加入图像滤波处理后的插值 校正法，对多金属条件下的校正伪影校正起到了 一 定效果[1o]，是校正多金属伪影的主流算法．其 中，Abdoli等 [11]提出对插值投影和原始投影进行 加权的插值方法，Jeong等[12]在此基础上加入了对 重建图像的TV滤波来平滑带状伪影．Zhao[13]和 Kratz[14】等分别提出基于小波变换和非均匀傅里 叶变换的插值运算法．插值校正法具有原理简单、 校正速度快等优点，但在多金属情况下仍存在以 下不足：1)插值法需要分割重建图像中的金属区域， 但金属周围的白色带状伪影会造成金属区域分割 的不准确，从而影响最终的校正结果；2)由于分割 出的金属区域被直接填补到插值处理后的重建图 像中，因此金属内部由硬化效应导致的杯状伪影无 法得到消除；31若非金属区域被金属区域包围，非 金属区域会被错误的归为被插值对象，致使其内部 的有用信息无法得到恢复． 针对多金属情况下的金属伪影校正问题，本文 首先对多金属情况下的金属伪影进行了仿真，进而 提出了一种基于投影校正的多金属伪影校正方法， 该方法对插值校正法进行了以下两点改进：1)通过 阈值法直接在投影域上分割金属区域，克服了重建 图像分割环节对校正结果产生的影响；2)对金属区 域的投影值进行迭代校正，避免了金属区域插值环 }国家高技术研究发展计划(批准号：2012AA011603)和国家自然科学基金(批准号资助的课题 十通讯作者．E—mail：leehotline@gmail．corn ◎2014中国物理学会Chinese Physical Society htp：／／wulixb．iphy．ac．ca 物 理 学 报 Acta Phys．Sin． Vo1．63，No．5(2014) 058702 2 多金属伪影的仿真 金属伪影的产生包含硬化效应、散射效应、局 部容积效应等I大]素 【1l5]，其对金属伪影的贡献比重 与具体的实验条件有关．在高分辨率工业CT系统 中，多金属伪影受硬化和散射效应的综合影响． 图像重建算法假设X射线具有单一能量， 单 能X射线穿过匀质物体时，其强度衰减情况遵循比 尔定律： 9= 1n Io ， (1) 其中9为投影数据， 为物质的衰减系数， 为射线 的初始强度， 穿过物体后的射线强度，此时，投影 数据与射线穿过物体的长度呈线性关系．而实际的 CT系统中受硬化和散射效应的影响，投影值随射 线穿过物体的长度增加偏离真值的程度逐渐变大， 总体趋势呈现非线性关系．当被扫描物质中含有多 金属时，金属连线方向相比其他方向的射线穿过金 属物质的路径差距较大，这是导致金属连线方向上 出现带状伪影的主要原闪l1(=i】． 本文对多金属条件下的金属伪影进行