基于光谱估计和图像重建的双能量CT迭代算法-CT理论与应用.PDFVIP

第27 卷 第1 期 CT 理论与应用研究 Vol.27, No.1 2018 年2 月 （45-54） CT Theory and Applications Feb., 2018 常少杰, 牟轩沁. 基于光谱估计和图像重建的双能量 CT 迭代算法[J]. CT 理论与应用研究, 2018, 27 (1): 45-54. doi:10.15953/j.1004-4140.2018.27.01.06. Chang SJ, Mou XQ. An iterative algorithm with simultaneously updating the spectrum and the image of dual energy computed tomography [J]. CT Theory and Applications, 2018, 27 (1): 45-54. (in Chinese). doi:10.15953/j. 1004-4140.2018.27.01.06. 基于光谱估计和图像重建的双能量CT 迭代算法 常少杰，牟轩沁 （西安交通大学电信学院图像处理与识别研究所，西安710049） 摘要：已知的光谱信息对双能量CT （DECT）重建具有重要的意义，但在实际应用中，光谱信息 很难直接获得，此外，双能量CT 的重建算法对噪声敏感，因此，在重建算法中需要设计对应的 噪声抑制方法来提升重建图像质量。本文提出了一种同时更新双能量CT 估计光谱和重建图像的 迭代方法，该方法不需要提前获得光谱信息，估计光谱和重建图像可以通过提出的算法同时获 得。两组实验用来验证算法的有效性。实验结果表明，提出方法不仅能够准确估计光谱，也可 以同时提升DECT 重建图像质量。 关键词：计算机断层成像 （CT）；双能CT；X 射线光谱估计；迭代重建 doi:10.15953/j.1004-4140.2018.27.01.06 中图分类号：TP 391.41 文献标志码：A 双能量CT （dual energy computed tomography，DECT）在材料分解、骨密度检测和无 [1－4] 损检测方面具有重要的作用 。在 DECT 扫描下，可以得到高、低能谱下的投影数据。利 用这一组投影数据，可以重建得到电子密度图像和有效原子序数图像，或者基材料密度图 像。 [5] 近10 年来，已有的DECT 重建方法可以分为3 大类，即基于投影域重建方法 ，基于图 [6－7] [8－10] 像域重建方法 以及迭代方法 。对于基于投影域重建方法，投影数据可以通过非线性 [5] 模型分成两种基材料的投影数据 。随后，利用一般的重建方法将基材料密度图像从分配好 的投影数据中重建出来。对基于图像域重建方法，两种基材料密度图像可以通过高低能量 [6] 下的重建图像直接利用线性组合得到 。考虑到投影数据中噪声的统计特性和材料分解模 型，一些迭代方法也被提出，材料分解和图像重建可以在迭代过程中同时进行。例如， [9] [10] Fessler 等 研究了关于DECT 的最大似然重建方法。随后，Xu 等 提出了一种基于似然正 则的宽能谱DECT 迭代重建方法。此外，在DECT 迭代方法中，可以加入去噪约束项来降低 [11] [12] 噪声对重建图像的影响，例如全变分（Total Variation，TV）正则 ，字典正则 等。在 医用X 射线球管范围内，X 射线满足宽能谱特性，在不同电压下产生相应的光子数量，因此，