[发明专利]基于混合高斯/泊松最大似然函数的CBCT三维重建方法

说明书

本发明公开了一种基于混合高斯/泊松最大似然函数的CBCT三维重建方法，属于统计迭代重建算法，是通过对投影数据的统计特性进行分析建立目标函数，利用一定的估计方法重建出高质量的图像。本发明首先通过anscombe变换将服从高斯‑泊松混合分布的实测投影数据变换为近似服从高斯分布，达到稳定噪声方差的目的，再通过梯度下降法求解混合高斯/泊松最大似然函数，结合三维全变分正则化方法(3D‑Total Variation,3D‑TV)进行迭代重建，从而达到改善重建图像质量的目的。

技术领域

本发明属于图像重建技术领域，主要涉及到混合高斯/泊松噪声下的CBCT三维重建方法。

背景技术

计算机断层成像(Computed Tomography，CT)技术通过无损方式获取物体内部结构信息、材质、缺陷等信息，广泛用于医学辅助诊断、工业无损检测、安全检查等领域。随着技术的快速发展出现了平板探测器，锥束CT(Cone Beam Computed Tomographyct，CBCT)系统通过它得到的投影数据是二维的，可以直接利用三维重建算法得到三维图像，是一种真正的三维重建。且CBCT具有体积小，重量轻，移动灵活，扫描速度快，射线利用率高、重建分辨率高等显著优点，广泛用于介入手术治疗、无损检测等。CBCT重建图像各向同性，提高了图像准确性和重建的实时性，能够更精确快速的分析病灶，提高患者的生活质量。CBCT价格低廉，可紧凑安装在口腔外科诊所和私人诊所，更方便为患者提供服务。为了降低CBCT成像中射线的辐射问题，低辐射剂量的CBCT成像方法应运而生，但是这种情况下信号的信噪比将会随之降低。因此，有效地去除低剂量条件下的噪声问题，对于CBCT重建方法而言是非常重要的。在低剂量条件下，系统中最常见的噪声包括高斯噪声和泊松噪声。高斯噪声主要来源于图像采集期间，往往由于不良照明或者高温引起的传感器噪声。泊松噪声是是由图像传感器中的光电转换过程引起的，在低剂量的CBCT重建过程中，影响更为严重。

1984年，Feldkamp、Davis和Kress针对圆形扫描轨迹提出了最著名的锥束三维图像重建算法--FDK算法。它基于圆轨迹扫描，是一种易于实现、重建速度快的算法，并且它也是近似算法中最为经典的算法。至今FDK算法仍然是商业CT应用的主流，在CT重建算法中有着无可替代的地位。该算法是一种近似的三维图像重建算法，FDK算法在小锥角时具有良好的重建效果，但是当锥角增大时，与中心层距离越远重建效果越差。针对FDK算法的缺点，出现了多种FDK算法的衍生算法，改善了FDK算法重建图像质量，然而，这些重建算法不足之处在于低剂量条件下复杂噪声模型对重建图像质量的影响，且随着噪声增加重建图像质量退化严重。以FDK为主的解析算法相当于是对一个积分方程的求解，而后出现迭代算法相当于是对一个大型线性方程组的求解。由于线性方程组对成像几何结构和成像物理效应的刻画比积分方程更加的真实，还可以加入先验知识，约束条件等，所以迭代算法比解析算法对扫描模式的依赖性小，更适合于实际成像问题，且迭代算法得到的结果更加准确，对噪声的抑制效果更好。由于迭代算法的重建速度慢是其致命缺点，很多学者对迭代算法的加速进行了研究，通过设计提高计算机的并行处理能力，加速或者利用专门的图像处理器进行加速或是对算法进行改进加速。随着硬件技术的不断更新和各种专用服务器的开发使得计算机的运行速度越来越快，迭代算法的优势也将更好的发挥。迭代重建算法分为代数重建和统计重建，统计迭代重建算法(SIR，Statistical Iterative Reconstruction)是通过对投影数据的统计特性进行分析建立目标函数，利用一定的估计方法重建出高质量的图像。统计迭代算法的研究，国内主要是针对放射型CT进行的。

本发明将利用统计迭代重建算法进行锥形束CT的重建，提出一种基于混合高斯/泊松最大似然函数的CBCT三维重建算法，从而重建出高质量的图像。

发明内容

本发明提出一种基于混合高斯/泊松最大似然函数的CBCT三维重建算法，属于统计迭代重建算法，是通过对投影数据的统计特性进行分析建立目标函数，利用一定的估计方法重建出高质量的图像。主要包含以下步骤：

步骤1，设置系统参数和重建图像估计值，并加载真实投影数据；