1、双能 CT 图像域基材料分解算法的研究进展摘要:能谱 CT 可产生不同X 射线能量下的基材料图像,所产生的基材料图像可用于组织成分和造影剂 分布的定性与定量评价,且对成像物质分离、鉴别的能力明显优于传统单能 CT。能谱 CT 中双能谱技术 是最常用的模式之一,在临床应用中发挥了重大作用。本文就双能谱 CT 图像域基材料分解的两物质分 解、多物质分解方法进行总结,最后展望未来可能的发展方向。关键词:能谱CT;双能CT;图像域;基材料分解计算机断层成像(computed tomography,CT)技术可以提供人体的断层解剖结构和形态学特 征1 。在很多情况下,不同组织有着非常相似的CT 值,单能
2、 CT 无法对其进行区分,因此能谱 CT 应 运而生,它拓展了 CT 的应用范围,也使 CT 成为临床应用中不可或缺的一种影像诊断设备2。当前,能谱 CT 在临床应用最广泛的模式是双能 CT(dual-energy computed tomography,DECT), 也是当前研究最多的能谱技术,其将物质分解为组成元素的能力源于 X 射线衰减的能量和元素的依 赖性3 。此时,X 射线的衰减主要为光电效应和康普顿效应4 。早期,Maass 等5采用图像域分解算法 成功得到两物质分解图像,从而能够有效地替代线性衰减系数进行组织区分,对物质成分进行判断, 进而实现物质识别。但往往由于人体组织成分并不
3、只是由两种物质构成,故在肝脂肪定量检测等多 组织分析时,就需要用到多物质分解(multi-material decomposition,MMD),对物质成分进行精确定性定量判断6-7。目前,基材料分解方法分为两大类:一类是将图像重建和基材料分解两个过程结合在一起的直 接迭代基材料分解法,也叫“一步法”,此方法由于计算过程复杂,求解过程耗时较长,因此应用较 少8;另一类为“两步法”,即把图像重建和基材料分解分成了两步进行;“两步法 ”根据高低能投 影数据集是否沿着相同的射线路径扫描可分为投影域和图像域基材料分解法9-10 。基于投影域的分解 法理论上能有效消除射束硬化伪影的影响,但大多数 DEC
4、T 无法满足高低能投影数据集沿着相同的射 线路径进行扫描;基于图像域的分解法往往更加灵活,适用于大多数 DECT,在计算上也更为简便11。本文主要针对“两步法 ”中的第二步,即对两幅高、低能 CT 图像分解成两物质分解图像或多物 质分解图像所运用到的主要的图像域基材料分解算法展开综述。1 两物质分解方法1.1 直接矩阵求逆分解根据基于图像的分解理论,CT 图像的线性衰减系数近似为两个基材料图像的加权和12 ,即:其中, XH/L 表示以mm-1 为单位的在不同能级下重建的CT 图像; d1/2 表示基材料图像或称为分解图像, 是无量纲的; i; j 是图像像素指数;成分矩阵由 xkH/L (
5、k = 1 或 2)元素组成, xkH/L 是在高、低能量 下测量的基材料线性衰减系数或质量衰减系数。此时,求解基材料图像只需将矩阵求逆然后计算即 可,过程如(2)式所示:直接矩阵求逆分解是图像域基材料分解的基础方法,其他方法都是在此法上衍变出来的。该方法的 优点是计算简便,缺点是分解过程中放大了图像的噪声,分解图像质量较差。为解决通过直接矩阵求逆生成分解图像的噪声问题,人们提出了分解后去噪和分解前去噪两种 解决方案(表 1),最初的方案是分解后去噪,这类方法在分解之后直接对分解图像进行噪声抑制, 以此提高分解图像的信噪比。 自迭代技术成功应用于 CT 图像重建之后,分解前去噪的方案得到广泛
6、应用。通过对 CT 图像重建过程进行优化,从而获得去噪后的 CT 图像,以尽可能降低噪声在后续分 解过程中的干扰。两种解决方法均可在一定程度上抑制噪声,提高图像的信噪比,但这两种方法都 是独立于分解过程的去噪方法,没有考虑分解过程的噪声对图像的影响,所以分解图像还是会存在 噪声干扰。表 1 两种解决方案的主要研究Table 1 Main research of two solutions方法文献时间优点缺点分解后去噪131976方法实现简单,计算效率高由于图像分辨率损失很大,效果有限141984实现简便效果有限151985实现简单效果有限161988缓解了空间分辨率损失的问题有边缘伪影1719
7、95算法可以在不考虑噪声相关性的情况下实现 噪声抑制分解后的图像中高频噪声被过度抑制,导致图 像纹理的改变182003算法利用CT 或分解图像的冗余结构或统计信 息进行噪声抑制,可以更好地抑制噪声没有完全描述 DECT 图像和分解图像之间的映 射关系分解前去噪192014使重建的两幅 CT 图像噪声变得强烈相关,进 而使得分解图像的噪声得到显著抑制CT 重建和图像分解的结合增加了计算的复杂 性,并且算法需要大量迭代才能收敛202015可以在保留定量测量和高频边缘信息的同 时显著降低噪声在心肌成像中仍会存在边缘效应212018在抑制噪声的同时可以保持图像边缘细节没有考虑分解过程的噪声对图像的影响
8、222019可获得高质量的重建 CT 图像以便后续分解没有考虑分解过程的噪声对图像的影响1.2 迭代分解迭代分解在直接矩阵求逆分解的基础上,充分挖掘分解过程的噪声统计特性,将物质分解和抑 制噪声两步骤结合在一起,在分解过程中逐步抑制噪声,可产生较好的图像质量。结果表明,迭代 分解可以在很大程度上抑制图像的噪声,明显提高分解图像的质量,但同时也增加了计算量。比较经典的图像域迭代分解是 Niu 等23在2014 年提出的,该算法采用光滑正则化的最小二乘估 计模型作为分解模型框架,在最小二乘项中包括分解图像的估计变量协方差矩阵的逆作为惩罚权重, 以获得更好的噪声抑制性能;与传统惩罚加权最小二乘模型去
9、噪不同,此模型结合分解和去噪,采 用共轭梯度迭代方法,边去噪边分解,最后得到高质量的分解图像。然而,该算法并没有完全描述 CT 图像和分解图像之间的映射关系,不能在 CT 图像的对象边界分布和噪声行为上实现明显的区分。为解决上述问题,Tang 等24充分探讨了CT 图像边界和噪声在频域的分布差异,在 2015 年提出 了一种基于多尺度惩罚加权最小二乘(penalized weighted least square,PWLS)的 DECT 图像域 分解算法,利用基于各向同性扩散的多尺度分解算法将高能量和低能量 CT 图像分解到尺度空间,该 方法充分利用了信号处理在每个尺度上的灵活性,以抑制放大的
10、噪声,同时保持物质特定图像的高 空间分辨率。Li 等 25对 分解模 型框架 中 的 正 则化项进行创 新 ,提 出 了 DECT-MULTRA( mixed union of learned transforms,MULTRA)算法,运用机器学习的方法对正则化项进行改进,正则化项使用的 是基于混合学习稀疏变换的混合联合模型,加了正则化项的分解模型可以加强分解后图像的平滑度、 同时保留边缘锐度。结果显示,该方法显示出了巨大的优越性,无论是在增加低对比度边缘的清晰 度还是减少不同软组织图像边界的伪影上,DECT-MULTRA 方法分解的图像质量最好。此外,与其他 方法相比较,DECT-MULTR
11、A 迭代收敛速度快,但是需要仔细调整合适的参数,计算量很大。王冲旭 等26使用了一种深度迭代残差网络(IR-Net)对正则化项进行改进,实验结果显示,相较于 DECT- MULTRA 调参过程复杂,运行时间过长,使用 IR-Net可以在很短的时间内完成分解,提高了基材料 分解效率。根据文献27对直接矩阵求逆分解法和Niu 提出的迭代分解法的模拟复现结果如图 1 所示,从图 像质量上可见迭代分解优于直接矩阵求逆分解。(a) (b)(c) (d)注:(a)和(b)是直接矩阵求逆法分解的骨和软组织的基材料图像;(c)和(d)是迭代分 解法分解的骨和软组织的基材料图像。图 1 两种分解方法的模拟复现结
12、果Fig.1 Simulation results obtained with the two decomposition methods1.3 基于深度学习方法的分解迭代分解显著提高了图像的信噪比,但因其计算量过大等问题限制了算法的发展。近年来,随着深度学习在医学领域的广泛应用,图像域基材料分解问题也可以用深度学习的方法来解决(表 2)28-30。 结果显示,使用基于深度学习方法的分解算法可以在极大程度上抑制基材料图像的噪声,甚至一些 神经网络分解的基材料图像不存在噪声干扰的问题,此外,还提升了基材料分解的精度和效率,但 是基于深度学习的分解算法需要大量的数据集训练网络模型,训练样本的数量制
13、约着算法的发展。表 2 基于深度学习方法分解图像的主要研究Table 2 Main research of image decomposition based on deep learning method文献时间网络模型优点缺点312018.03U-Net分解图像有较低的噪声水平方法没有降噪能力,只具备分解能力272018.09FCN + FCL具有高分解精度和噪声鲁棒性在边缘保持方面没有明显改进322019.04蝴蝶网极大程度上抑制了图像噪声,提高了分解 图像的质量如果改变 X 射线源设置,则需要重新 训练网络332021.03DIWGAN分解图像与真实分解图像很接近,并且在 噪声和伪影抑
14、制方面效果较好在分解的软组织图像中有一些软组织 结构丢失2 多物质分解方法2.1 直接矩阵求逆分解假设 DECT 重建图像中每个体素都是由各种基材料组成的,那么每个体素所反映的吸收系数可以 理解为是由各种基材料吸收系数共同组合而成34 ,如公式(3)所示:式中 (E) 表示在能量级 E 下各种基材料衰减系数的总和,N 表示基材料的数目, i (E) 表示在能量 级 E 下第 i 种基材料的衰减系数, i 表示第 i 种基材料体积占所有基材料总体积的体积分数。利用 DECT 的能谱数据可重建出两幅单能图像,如公式(4): 0 i 1; i = 1; 2; 3。(4)求解 DECT 测量的多物质分
15、解是一个不适定问题,因为要从两幅图像中得出多幅图像,为了解决 这个问题,2009 年,Liu 等35提出了一种使用质量分数守恒为约束条件的三物质分解方法。该算法 成功求出基材料的质量分数,但前提是在混合物的组成已知的情况下,所以该方法适用性不高。 2014 年,Mendona 等36独创的提出了一种图像域MMD 方法,成功将一次扫描得到的低能量和高能量 的滤波反投影图像分解为多个基材料图像。该方法分解的 MMD 图像已成功地应用于虚拟非对比增强 图像、脂肪肝和肝纤维化的应用,也成为了图像域 MMD 直接矩阵求逆的基本方法,为其他方法的拓 展奠定了基础。然而,直接矩阵求逆分解的图像容易受噪声的干
16、扰。Jiang 等37提出了一种用于DECT 的 MMD 图像域噪声抑制方法,算法推导了 MMD 噪声的传播规律, 采用多维旋转策略来寻找噪声干扰最小的主轴,沿主轴方向抑制噪声。该方法在不改变噪声分布和 空间结构的前提下,有效地抑制了分解后图像的噪声,提高了图像质量,但没有考虑分解前噪声对 图像的干扰。Lee 等38提出了一个用于DECT 的 MMD 框架,该框架包括分解前和分解后阶段,将全变 分去噪方法作为一种噪声抑制手段应用于分解前和分解后阶段,结果表明,该方法能在提高分解精 度的同时更清晰地区分密度差较小的基材料,但没有考虑分解图像的噪声特性。为此,降俊汝等39 利用字典学习充分发掘分解
17、图像的稀疏性,进一步提高了分解图像的准确性,但处理步骤繁多,参 数难调,增加了计算问题的复杂度。2.2 迭代分解MMD 直接矩阵求逆分解方法中,高低能量的线性衰减系数只逐像素估计基材料的体积分数,而 不考虑测量的噪声统计特性,为此,采用统计模型建模,迭代求解来抑制噪声。Long 等40提出了一 种基于边缘保持正则化的惩罚似然(penalized-likelihood,PL)方法,PL 方法显著降低了分解图 像中的噪声、条纹和交叉污染伪影。然而,这种方法计算量大,主要是由于多个基材料图像之间的 正向和反向投影以及对低能量和高能量的正弦图像的映射,而且没有考虑每个像素中基材料类型的 稀疏性。随后,
18、Xue 等41建立带有负对数似然项和边缘保持( edge-preserving,EP)正则化的 PWLS 估 计模型,后称为 PWLS-EP-LOOP 方法。与 Mendona 提出的图像域直接矩阵求逆相比,PWLS-EP-LOOP 方法在保持基材料图像解剖结构和分解精度的同时,能显著抑制噪声。由于不需要正向或反向投影, 因此在计算上比 PL 方法更实用。同 PL 方法一样,没有考虑每个像素中基材料类型的稀疏性,为此, Ding 等42考虑了分解图像梯度的低秩性、基材料成分的稀疏性以及质量和体积守恒提出了一种带有 3 个正则化项的 PWLS 模型,该方法减少了交叉污染,提高了分解图像的精度,但
19、图像中高频噪声被 过度抑制,易导致图像纹理的改变。3 总结和展望本文对双能 CT 图像域基材料分解的两物质分解和多物质分解法进行了总结。对于两物质分解法, 深度学习方法在图像去噪和保持图像细微结构方面好于直接矩阵求逆、迭代分解法33;对于多物质 图像分解方法中,有报道42迭代分解法在基材料分解精度方面明显优于直接矩阵求逆法,但未见深 度学习分解法的应用。因此基于深度学习的基材料图像分解方法将是下一步研究的重要内容。尽管基材料分解算法理论上可以帮助临床医生更好的分析病灶,但低 X 射线能量图像噪声会影 响物质分解的准确性,鲁棒性好的基材料分解算法可能会克服这一限制,因此基材料分解算法还需 要进一
20、步提升43。图像域基材料分解有 3 个重要的发展方向, 减小噪声对基材料分解的影响至关重要; 由 于图像重建过程中产生的光束硬化伪影无法消除,消除伪影的影响可以极大程度的提高分解图像的 精确度; 以光子计数探测器 CT(photon-counting detector CT,PCD CT)为代表的能谱 CT 成像 技术一经问世便引起了广泛关注,PCD CT 成像较 DECT 成像的优势在于重建过程受能谱硬化效应的 影响小,可以有效抑制射束硬化伪影的产生,但其尚未在临床上应用。能谱 CT 的图像域基材料分解研究将是一个新的发展方向,为完全消除射束硬化伪影的影响带来 了曙光。参考文献1 高海英.
21、能谱 CT 成像关键参数检测技术研究D. 广州: 南方医科大学, 2015.GAO H Y. Study on testing techniques of spectral CT imaging key parametersD. Guangzhou: Southern Medical University, 2015. (in Chinese).2 韩文艳. CT 能谱成像的基本原理与临床应用优势J. 中国医疗设备, 2015,30(12): 90-91.HAN W Y. Basic principle and advantages of clinical application of CT
22、energy spectrum imagingJ. China Medical Devices, 2015, 30(12): 90-91. (in Chinese).3 ALVAREZ R E, MACOVSKI A. Energy-selective reconstructions in X-ray computerized tomographyJ. Physics in Medicine & Biology, 1976, 21(5): 733-744.4 王丽新, 孙丰荣, 仲海, 等. 双能 CT 成像的数值仿真J. 航天医学与医学工程, 2015,28(5): 350-357.WANG
23、 L X, SUN F R, ZHONG H, et al. Numerical simulation of dual-energy CT imagingJ. Space Medicine & Medical Engineering, 2015, 28(5): 350-357. (in Chinese).5 MAASS C, BAER M, KACHELRIESS M. Image-based dual energy CT using optimized precorrection functions: A practical new approach of material decomposition in image domainJ. Medical Physics,