双能CT能谱误差传递分析-CT理论与应用研究.PDF

第22卷 第 4 期 CT理论与应用研究 Vol.22, No.4 2013年10月（571-578） CT Theory and Applications Oct., 2013 刘圆圆, 郑鹏, 邢宇翔, 等. 双能CT 能谱误差传递分析[J]. CT 理论与应用研究, 2013, 22(4): 571-578. Liu YY, Zheng P, Xing YX, et al. Propagation analysis for spectrum error in dual-energy CT[J]. CT Theory and Applications, 2013, 22(4): 571-578. 双能 CT 能谱误差传递分析  1 1 2，3 2，3 2，3 刘圆圆 ，郑鹏 ，邢宇翔 ，陈志强 ，张丽 1.环境保护部核与辐射安全中心，北京 100082 2.清华大学工程物理系，北京 100084 3.粒子技术与辐射成像教育部重点实验室（清华大学），北京 100084 摘要：双能CT 技术可以重建物质的原子序数和电子密度信息，是一种有效的物质辨别技术。进 行 X 射线能谱估计是双能重建的前提，能谱估计的准确性直接影响双能重建结果。但是，目前 还没有学者对能谱估计的误差如何影响双能图像质量进行系统的研究。本文从双效应分解双能 重建出发，量化分析能谱估计误差对双能重建结果的影响并给出误差传递的理论计算方法，论 证影响双能图像准确重建的关键因素，同时由此定义适用于双能成像分析的能谱误差和针对双 能成像的等效能谱概念。初步实验结果验证了这些理论的有效性。 关键词：双能CT；物质辨别；能谱估计；误差分析 文章编号：1004-4140（2013）04-0571-08 中图分类号：TP391.41 文献标志码：A 传统的 X 射线计算机断层成像技术（CT）可以重建物质的线性衰减系数，得到断层的 结构信息，与之相比，双能 CT（Dual-Energy CT）可以重建物质的有效原子序数 Z 与电子 [1] 密度ρ ，从而达到物质辨别的效果。1976 年Alvarez 等 最早提出了双能 CT 重建方法。近 e 年来，双能 CT 在安全检查[2] [3] [4－7] 、医疗诊断 等领域得到了广泛的研究和应用 。目前，双能 [8] [9] [10] CT 重建方法主要包括前处理 、后处理 和迭代重建三种方法 。具体的双能分解形式又可 分为双效应分解[1]与基材料分解[11]，在非医学 CT 中，特别是安全检查领域，最为常用的双 能重建方法为双效应分解-前处理重建方法，本文针对此种重建方法进行研究。 在实际应用中，双能重建得到的 Z 和 ρe 的图像结果较差，尤其是 Z 的图像。重建结 果的误差来源很多，如 X 射线的量子噪声、探测器中的电子学噪声、线性衰减系数表达 公式的误差等，除此之外，X 射线能谱估计的误差对最终重建结果的影响也是一个重要 因素。 工程实践中，双能 CT 选择两个 X 射线能谱进行投影与重建[12－13]。由于 X射线能谱是未 知的，因此需要利用数值模拟或实验数据对 X 射线进行能谱估计，从而获得较高的双能重 建精度[14]。实践中，常选用蒙卡结果代替真实的 X 射线能谱，然而蒙卡结果与真实 X 射线 能谱存在误差，目前还没有学者对能谱误差进行合理化的定义，此误差对最终双能重建结 果的影响也没有明确的认识。 本文基于双效应分解-前处理重建方法，首先，对能谱误差进行合理化定义，其次，对