[发明专利]一种CT迭代重建中伪影校正方法

说明书

本发明公开的一种CT迭代重建中伪影校正方法，具体包括获取数字层析成像设备原始投影数据y、采用滤波反投影算法对获取的投影数据y进行重建，获得校正图像I₁，对校正图像I₁进行预处理，得到图像校正模板I₂，对图像校正模板I₂进行前向投影，得到投影校正模板y₁，利用投影校正模板y₁对原始投影数据y进行校正，得到校正后投影以及对校正后投影进行迭代重建，得到最终所需重建的图像μ六个步骤，本发明的CT迭代重建中伪影校正方法，与现有技术相比，操作简单，能够有效去除CT迭代重建中由rFOV引起的伪影，并且重建结果没有引起图像CT值的偏离。

技术领域

本发明属于医学影像的图像处理方法技术领域，涉及一种CT影像校正方法，具体涉及一种CT迭代重建中伪影校正方法。

背景技术

滤波反投影(FBP)算法一直是商用CT的主要重建方法，但是FBP算法并不适用于重建带有严重噪声的投影(即低剂量CT数据)。因此迭代重建算法(IterativeReconstruction，IR)成为了近年来的研究热点，因为IR算法可以通过构建合理的统计模型，同时根据采集数据的不同可以引入不同的先验信息，辅助迭代重建，可以得到优秀的成像质量。在过去几年中，不同的CT厂商已经发布了几种商业IR技术(例如GE Healthcare的ASIR，Siemens Healthcare的SAFIRE)。

但是在CT成像设备(例如，医用CT设备)重建过程中，身体的一些部分(例如，肩部或骨盆)可能超出重建视野(reconstruction field of view，rFOV)，若直接对原始投影进行迭代重建(Iterative Reconstruction，IR)，可能会在迭代重建过程中，由于rFOV外部的组织导致当前投影和测量投影之间的存在不一致性，并因此导致了重建图像中的伪影，以及出现CT值的偏离。

发明内容

本发明提供一种CT迭代重建中伪影校正方法，具体包括如下步骤：

步骤1、获取数字层析成像设备原始投影数据y；

步骤2、采用滤波反投影算法对步骤1中获取的投影数据y进行重建，获得校正图像I₁；

步骤3、对步骤2中得到的校正图像I₁进行预处理，得到图像校正模板I₂；

步骤4、对步骤3中得到的图像校正模板I₂进行前向投影，得到投影校正模板y₁；

步骤5、利用步骤4中得到的投影校正模板y₁对原始投影数据y进行校正，得到校正后投影

步骤6、对步骤5中得到的校正后投影进行迭代重建，得到最终所需重建的图像μ。

本方法的特点还在于，

步骤2中所述的采用滤波反投影算法进行重建时，其重建视野与投影数据采集时的扫描视野相同。

步骤3具体包括如下步骤：

3.1、对校正图像I₁中位于迭代重建视野区域内的像素置零，并保留重建视野之外的像素值，得到初步校正图像

3.2、对所述初步校正图像进行阈值处理，得到图像校正模板I₂；

阈值处理的具体方式为公式(1)：

其中，t表示可调阈值参数。

步骤5采用公式(2)：