一种面向多能谱X射线CT成像的金属伪影校正方法技术

本发明专利技术公开了一种面向多能谱X射线CT成像的金属伪影校正方法，首先对含有金属伪影的多能谱X射线CT投影图进行加权图像重建及材料分解，得到加权重建图像I0、金属材料及其他基材料的投影图；然后，根据金属材料的投影图进行图像重建得到金属材料图像Im，作为后续迭代重建过程中的先验知识；根据分解后所有基材料的投影图及其在指定能量下的X射线衰减系数生成虚拟单能投影图Ps；最后，以I0作为迭代图像重建的初始图像，Im作为迭代图像重建的先验知识，Ps作为迭代过程中的目标投影图，运用统计迭代重建算法进行图像重建，得到最终的金属伪影校正CT图像。

A metal artifact correction method for multi-spectral X-ray CT imaging

The invention discloses a metal artifact correction method for multi-energy spectrum X-ray CT imaging. Firstly, the multi-energy spectrum X-ray CT projection images containing metal artifacts are reconstructed by weighted image reconstruction and material decomposition, and the projection images of weighted reconstruction images I0, metal materials and other base materials are obtained. Then, the metal material images are reconstructed by image reconstruction according to the projection images of metal materials. Prior knowledge in the subsequent iterative reconstruction process; Ps are generated from the projection maps of all base materials after decomposition and their X-ray attenuation coefficients at specified energies; finally, I0 is used as the initial image of iterative image reconstruction, Im is used as the prior knowledge of iterative image reconstruction, Ps is used as the target projection map in the iterative process, and statistical iterative reconstruction algorithm is used to map. The final metal artifact correction CT image is obtained by image reconstruction.

【技术实现步骤摘要】
一种面向多能谱X射线CT成像的金属伪影校正方法
本专利技术涉及计算机断层成像
，尤其涉及一种面向多能谱X射线CT成像的金属伪影校正方法。
技术介绍
X射线CT(ComputedTomography)具有成像速度快、空间分辨率高等优势在医学诊断、工业检测等领域具有广泛应用。但是，CT成像时由于被检查物体中存在金属材料（如金属材质的假牙）将产生严重的金属伪影，表现为条纹状亮纹、暗带或阴影，使得CT图像失真，严重降低CT图像的质量，掩盖图像细节，不利于医学精确诊断及工业高精度检测等应用，因此有必要对金属伪影进行校正。现有技术一般采用插值法和迭代法进行金属伪影的校正。插值法是在投影域利用金属区域周围的投影值对金属伪影区域进行插值或者根据某种相似或对称原理对金属伪影区域进行成片插值修补，该方法会导致插值区域与附近投影区域的边界模糊，产生次级伪影。迭代法通过对金属投影区域之外的投影区域迭代反投影重建CT图像，该方法虽能在一定程度上校正金属伪影，但计算效率较低。进一步地，由于CT图像重建理论的前提是X射线是单能的，而在实际的成像系统中，X射线是连续能谱射线，因此上述方法难以对金属伪影进行有效的校正。在基于迭代重建的金属伪影校正中，重点是先利用反投影重建算法得到初始CT图像，在对该图像进行图像分割得到金属材料的图像，将其作为后续迭代重建的先验知识。例如中国专利申请号201710868075.5和201710084349.1。但由于重建图像中存在严重的金属伪影，难以准确分割出金属图像，只能对金属伪影起到部分抑制作用；进一步地，由于常规CT图像重建理论的基本前提是X射线是单能的，而实际的成像系统中，X射线是连续能谱射线，因此上述方法难以对金属伪影进行有效的校正。中国专利申请号201710868075.5提出“用于基于分割和光谱的金属伪影降低的计算机程序和CT系统”。在该专利中，基于在每个投影角度的金属的厚度和在每个投影角度的X-射线衰减，从分割的高-Z体积数据集中生成高-Z投影集，将这些集合归一化，并且然后使用低-Z投影和高-Z投影集创建校正的断层摄影体积数据集。中国专利申请号201710084349.1提出“一种CT图像金属伪影校正方法及装置”。在该专利中，通过在先验投影的约束下，根据原始投影和金属投影，对原始投影的金属区域迭代进行高斯能量下降处理，得到满足条件的初步未校正投影并进行反投影重建得到初步校正图像，将初步校正图像平滑去噪后用于更新先验图像，通过迭代更新先验图像，并将每次更新后的先验图像进行投影得到更新后的先验投影，在更新后的先验投影的约束下，获得满足条件的初步未校正投影并进一步得到初步校正图像，直至获得满足条件的初步校正图像，与金属图像融合，得到最终校正图像。其中，金属图像、先验图像和先验投影通过以下方法获得：利用CT系统中的原始投影反投影重建未校正CT图像，对所述未校正CT图像分别进行第一分割和第二分割，对应得到金属图像和先验图像，对所述金属图像进行投影得到金属投影；对所述先验图像进行投影，得到先验投影。近年来，随着探测器技术的发展，基于光子计数X射线探测器的多能谱X射线CT成像系统将逐步获得应用。多能谱X射线CT成像系统通过能量选通阈值对不同能窗的X射线光子进行数据采集和图像重建，可有效改善重建CT图像的质量，提供物质的电子密度、等效原子序数的信息，具有独特的材料分解与识别能力，可为CT图像伪影校正带来新的途径。
技术实现思路
本专利技术所要解决的技术问题是针对
技术介绍
中所涉及到的缺陷，提供一种面向多能谱X射线CT成像的金属伪影校正方法，能够有效提高CT图像质量，有利于临床诊断及工业无损检测。本专利技术为解决上述技术问题采用以下技术方案：一种面向多能谱X射线CT成像的金属伪影校正方法，所述金属伪影校正的步骤包括：步骤1：采用能谱CT探测器对含有金属的被检测物体进行多能谱CT投影，得到存在金属伪影的多能谱CT投影图，根据被检测物体中材料的特性以及所包含金属的特性，确定包含金属材料在内的各种基材料的类别数和材质；步骤2：根据多能谱CT投影图，利用加权重建算法得到重建图像I0；步骤3：根据多能谱CT投影图，运用材料分解算法得到包含金属材料在内的各种基材料的投影图；步骤4：根据包含金属材料在内的各种基材料的投影图及其在预设的能量阈值下的X射线衰减系数进行计算，得到虚拟单能投影图Ps；步骤5：根据金属材料的投影图，利用反投影重建算法得到金属材料的CT图像Im；步骤6：运用统计迭代重建算法进行金属伪影校正和CT图像重建，其中，迭代重建的初始图像是加权重建图像I0；迭代重建的先验知识是金属材料的CT图像Im，在迭代重建过程中，当前重建图像中与金属材料坐标相同处的像素值设置为金属材料的像素值；迭代重建的目标投影图是虚拟单能投影图Ps。作为本专利技术一种面向多能谱X射线CT成像的金属伪影校正方法进一步的优化方案，所述被检测物品中基材料的类别数不大于能谱CT探测器的能量分段数。所述预设的能量阈值和金属材料的材质相关。本专利技术采用以上技术方案与现有技术相比，具有以下技术效果：本专利技术根据多能谱X射线CT成像所具有的材料识别能力，运用材料分解算法获得金属材料图像及虚拟单能投影图，进而利用统计迭代重建算法进行图像重建，可有效克服金属伪影，获得高质量的CT图像，有利于临床诊断及工业无损检测。附图说明图1为一种面向多能谱X射线CT成像的金属伪影校正方法的流程图。具体实施方式下面结合附图对本专利技术的技术方案做进一步的详细说明：如图1所示，本专利技术公开了一种面向多能谱X射线CT成像的金属伪影校正方法，下面以患者前列腺内植入金标记物情况下多能谱CT图像金属伪影校正为例进行说明，校正的具体步骤如下：步骤1：采用能谱CT探测器（能量分段数不小于3）进行多能谱CT投影，得到存在金属伪影的多能谱CT投影图，并进一步确定基材料为软组织、骨组织和金；步骤2：根据多能谱CT投影图，利用投影加权重建算法或图像加权重建算法得到重建图像I0；步骤3：根据多能谱CT投影图，运用材料分解算法得到金、软组织和骨组织的投影图；步骤4：根据金、软组织和骨组织的投影图及其在50KeV下的X射线衰减系数计算得到虚拟单能投影图Ps；步骤5：根据金的投影图，利用反投影重建算法得到金的CT图像Im；步骤6：利用统计迭代重建算法进行金属伪影校正和CT图像重建，其中，迭代重建的初始图像是加权重建图像I0；迭代重建的先验知识是金的CT图像Im，在迭代重建过程中，当前重建图像中与金的坐标相同处的像素值设置为金的像素值；迭代重建的目标投影图是虚拟单能投影图Ps。本次探测仅包含软组织、骨组织和金三种基材料，很明显不大于能谱CT探测器的能量分段数。本
技术人员可以理解的是，除非另外定义，这里使用的所有术语（包括技术术语和科学术语）具有与本专利技术所属领域中的普通技术人员的一般理解相同的意义。还应该理解的是，诸如通用字典中定义的那些术语应该被理解为具有与现有技术的上下文中的意义一致的意义，并且除非像这里一样定义，不会用理想化或过于正式的含义来解释。以上所述的具体实施方式，对本专利技术的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本专利技术的具体实施方式而已，并不用于限制本专利技术，凡本文档来自技高网...

【技术保护点】
1.一种面向多能谱X射线CT成像的金属伪影校正方法，其特征在于，所述金属伪影校正的步骤包括：步骤1：采用能谱CT探测器对含有金属的被检测物体进行多能谱CT投影，得到存在金属伪影的多能谱CT投影图，根据被检测物体中材料的特性以及所包含金属的特性，确定包含金属材料在内的各种基材料的类别数和材质；步骤2：根据多能谱CT投影图，利用加权重建算法得到重建图像I0；步骤3：根据多能谱CT投影图，运用材料分解算法得到包含金属材料在内的各种基材料的投影图；步骤4：根据包含金属材料在内的各种基材料的投影图及其在预设的能量阈值下的X射线衰减系数进行计算，得到虚拟单能投影图Ps；步骤5：根据金属材料的投影图，利用反投影重建算法得到金属材料的CT图像Im；步骤6：运用统计迭代重建算法进行金属伪影校正和CT图像重建，其中，迭代重建的初始图像是加权重建图像I0；迭代重建的先验知识是金属材料的CT图像Im，在迭代重建过程中，当前重建图像中与金属材料坐标相同处的像素值设置为金属材料的像素值；迭代重建的目标投影图是虚拟单能投影图Ps。

【技术特征摘要】
1.一种面向多能谱X射线CT成像的金属伪影校正方法，其特征在于，所述金属伪影校正的步骤包括：步骤1：采用能谱CT探测器对含有金属的被检测物体进行多能谱CT投影，得到存在金属伪影的多能谱CT投影图，根据被检测物体中材料的特性以及所包含金属的特性，确定包含金属材料在内的各种基材料的类别数和材质；步骤2：根据多能谱CT投影图，利用加权重建算法得到重建图像I0；步骤3：根据多能谱CT投影图，运用材料分解算法得到包含金属材料在内的各种基材料的投影图；步骤4：根据包含金属材料在内的各种基材料的投影图及其在预设的能量阈值下的X射线衰减系数...

【专利技术属性】
技术研发人员：周正东，
申请(专利权)人：南京航空航天大学，
类型：发明
国别省市：江苏,32