局部扫描情况下多项式拟合的X射线硬化校正.pdfVIP

第 27 卷第 11 期 2015 年 11 月 强激光与粒子束 HIGH POWER LASER AND PARTICLE BEAMS Vo l. 27 , No.11 Nov. , 2015 局部扫描情况下多项式拟合的 X 射线硬化校正椅 陈 浩1 ， 2 ， 陈云斌 1 ， 2 ， 李寿涛1，2 (1.中国工程物理研究院应用电子学研究所，四川绵阳 621900; 2. 国家 X 射线数字化成像仪器中心，四川绵阳 621000) 摘 要: 在进行 X-CT 成像时，容易诱发射束硬化效应。通过多项式拟合建立起多色投影同透射距离的 数学模型，然后将多色投影校正为单色投影，是一种简单高效的射束硬化校正方法。对于局部扫描的情况，在 重建图像域通过正投影方法计算得到的透射距离不再代表真实的射线穿过长度，拟合数据不再有效。因此借 助一圆柱模体，根据圆柱自身的几何特性，通过解析的方法计算得到射线真实的透射距离，并用于射束硬化校 正。实验结果表明，该方法能够有效抑制射束硬化引起的图像伪影，提高 CT 重建图像质量。 关键词: 硬化校正 s 局部扫描 s 多项式拟合; 正投影 中图分类号: TP391 文献标志码: A doi:10. 11884/HPLPB201527. 114001 CT 是医学和工业领域的常用设备之-，它采用特定的 CT 重建算法获得切片图像，对疾病诊断和缺陷检 查起着重要的作用。然而，传统的 CT 滤波反投影重建算法存在一个前提假设 :X 射线能谱单一，射线在传播 过程中，同种材料的物质对射线具有同等的衰减效应。然而这种理想假设背离了实际的物理特性，X 射线能谱 具有一定的带宽，射线同物质相互作用时，不同能量的光子具有不同的衰减程度，射线在传播过程中，能谱会发 生变化。具体地讲，射线谱峰向高能方向偏移，高能光子所占比例越来越大，射线变得更难衰减，俗称射束硬 化。射束硬化是诱发 CT 重建图像诸多伪影的物理原因之一，比如杯状伪影、金属伪影。常用的射束硬化校正 方法是线性化方法，即通过数学函数拟合建立起多色投影同透射距离的映射关系，比如典型的多项式拟合，然 后将多色投影校正为线性单色投影[1-4J 。李岭等报道了一种通过临界点坐标，构造一种权重校正函数用于硬化 校正。但该方法局限于圆柱形等对称扫描对象[5J 。在实际工程应用中，基于多项式拟合的硬化校正方法其关 键点是需要得到每个多色投影值所对应的射线透射距离。一种办法是，对物体进行完整的 CT 扫描并得到重 建图像，然后在重建图像域采用正投影的方法计算射线穿过物体的长度[6-8J 。上述方法的局限性在于， CT 的 扫描视野(FOV)必须完全覆盖物体，不能存在投影数据截断的现象(即局部扫描)。而对于小视野 CT，比如 配备小尺寸平板探测器的锥束牙科 CT，局部扫描的情况将难以避免。局部扫描只能重建物体部分区域，在重 建图像域通过正投影计算获得的透射距离，不再代表准确真实的 X 射线穿过物体的长度，拟合数据将不再有 效。 针对牙科 CT 局部扫描情况下的 X 射线硬化校正问题，本文提出一种解决办法，设计了一个直径 145 mm，由聚甲基丙烯酸甲醋(PMMA)组成的圆柱形模体，结合圆柱自身的几何特征，通过解析计算的方法来得 到真实的射线透射距离，并用于多项式拟合进行硬化校正。由于牙科 CT 的扫描对象多为人体头部，人体的软 组织与 PMMA 校正模体具有相近的密度，因此通过 PMMA 模体建立的硬化校正模型可适用于牙科扫描的硬 化预校正。 1 基于正投影技术的 X 射线硬化校正 透射距离和多色投影的非线性关系可以用多项式表达为 N L=~αiPi(L) i=O (1) 式中 :L为透射距离 ， ai为第i 阶多项式的系数 ， P(L)表示透射距离为L时所对应的多色投影 ， N为多项式阶数。 普收稿日期 :2015-05-15; 修订日期: 2015-07-24 基金项目:中国工程物理研究院应用电子学研究所基金项目 (2012XZY02) -11:者简介:陈 浩(1 978一九男，硕士，从事 X 射线辐射成像技术研究; ultrablue@263. net 。 114001-1 强激光与粒子束 由于单色投影 M(L)满足 N M(L) =μ(Eo)L =μ (Eo) ~ ai pi (L) (2) 式中 :μ (Eo)表示射线能量为 E。下的衰减系数，一般取为 P(L)在 0 处的一阶导数，即 μ(Eo) =p(O) =l/L(O) =l/aj (3) 这样就由非线性多色投影校正得到了线性单色投影。函数拟合需要选取多色投影及其对应的透射距离作 为数据拟合点。多色投影即是探测器测量得到的实际投影，透射距离可通过 ray-driven 方法川，在重建图像