[发明专利]一种低剂量CT图像重建方法

说明书

技术领域

本发明涉及一种图像重建方法，尤其涉及一种低剂量CT图像重建方法。 

背景技术

计算机断层成像（Computer Tomography,CT）是一种有效的无损三维检测技术，它能从样品的二维投影中解出断层切片，实现三维可视化，在医学、材料学、地球物理学、考古学及生物学等科学领域以及工业无损检测领域发挥着重要的作用。参见图1和2，典型的CT成像系统包括具有一定穿透能力的辐射源（平行束110、扇形束210及锥束）和探测器（平面探测器130和非平面探测器230），两者之间是样品台120，通过样品与辐射线和探测器之间的相对旋转（旋转中心O），获取样品不同角度的投影图像140，然后将投影图像输入计算机150，通过图像重建算法计算出断层切片，即重建图像160。其中辐射源包括X射线、伽玛射线、中子、质子、电子、声波和离子束等等。 

对于传统的医用X射线CT及荧光造影，控制辐射剂量是必须考虑的一个关键问题，病人受到的辐射剂量随着获取投影数呈近似线性关系，长时间接受电离辐射增加了对病人的伤害。另外，对于电子显微技术来说，由于样品必须被放置在衬底上，导致某些角度的投影无法获取到，最终导致重建结果分辨率的降低。 

传统CT重建方法中，投影数据是按照等角度增量的方式进行获取，即不同的投影构成一个极坐标系，而待重建的断层切片是笛卡尔坐标，因此在数据重建过程中不可避免地使用到插值，从而增加了重建结果的误差，降低了重建结果的分辨率和信噪比。 

传统CT重建方法中使用最为广泛的是滤波反投影（FBP），该方法较容易实现且重建速度快，但是当投影角度不完全时，滤波反投影往往不能给出满意的结果。同时由于该方法在反投影过程中使用插值算法，重建结果中会不 可避免地出现伪影，最终导致图像质量的降低，更重要的是，FBP算法需要大量的投影角度，使得数据获取过程中样品接收较大的辐射剂量，在医用CT应用中，将给受检者带来辐射伤害。其他的CT重建算法，如代数迭代算法（ART）及期望最大化算法（EM）等，虽然可以有效解决投影角度不完全的问题，并能利用少量投影实现图像重建，有效降低样品受到的辐射剂量，但是其重建过程只涉及笛卡尔坐标系和极坐标系，因而需要大量的插值，降低了重建结果的精确性，同时重建速度慢，限制了其使用范围。 

发明内容

针对上述现有技术中存在的缺陷，本发明的目的在于提供一种低剂量CT图像重建方法，在保证重建图像质量的同时，还能有效降低样品所受到的辐射剂量。 

为了实现上述目的，本发明采用如下技术方案： 

一种低剂量CT图像重建方法，该方法包括以下步骤： 

步骤（1），获取不同投影角度θ的等斜率增量的投影数据； 

步骤（2），对所述投影数据进行修正； 

步骤（3），对修正后的所述投影数据执行傅里叶变换，以使其转换成极坐标系下的频域空间数据； 

步骤（4），将所述极坐标系下的频域空间数据转换成伪极坐标系下的频域空间数据； 

步骤（5），通过迭代法使所述步骤（4）中的所述频域空间数据在伪极坐标系下的频域空间与时域空间之间重复转换，直至满足预设的最终条件； 

步骤（6），输出满足最终条件时的时域图像。 

进一步地，所述步骤（1）中的所述投影数据为通过等斜率增量方式获取的数据。 

进一步地，所述步骤（1）中的所述投影数据为通过等角度增量方式获取并从获取的数据中选取出的近似等斜率增量的数据。 

优选地，所述步骤（1）中的所述投影角度θ满足以下公式：