[发明专利]分布式光源CT图像重建方法与系统

权利要求书

1.一种分布式光源CT图像重建方法，其特征在于，包括：

获取N段扫描投影数据；以及

根据获取N段扫描投影数据，利用加权滤波反投影重建算法，获取重建图像，其中，N≥2，

所述获取N段扫描投影数据，包括：对物体进行N段扫描，获取N段扫描投影数据p_I(l，t)～p_N(l,t)，其中p_N(l,t)表示进行第N段扫描过程中，直线光源阵列中坐标为l的光源点闪烁时，探测器上坐标为t处的单元采集到的投影数据值，

对所述N段扫描投影数据p_I(l，t)～p_N(l，t)进行平移变换，得到N段平移后的投影数据q_I(l，t)～q_N(l，t)，

对所述N段平移后的投影数据进行加权滤波处理，得到N段加权滤波后的数据Q_I(l，t′_I)～Q_N(l,t′_N)，t′_N表示进行第N段扫描过程中，在直线光源阵列中坐标为l的光源点闪烁时，经过点(x,y)的射线在等效探测器上的采样位置，

对所述N段加权滤波后的数据进行反投影，得到所述N段扫描投影数据的重建图像f_I(x,y)～f_N(x,y)，(x,y)表示以物体中心为原点所建立的坐标系坐标，

对所述N段扫描投影数据的重建图像进行加和，得到最终重建图像f(x,y)，

在N为2时，所述根据获取N段扫描投影数据，利用加权滤波反投影重建算法，获取重建图像的重建公式为：

f(x,y)＝f_I(x,y)+f_II(x,y)

其中，

D和D′分别表示光源阵列和探测器阵列到物体中心的距离，h为卷积函数核，理论值为W_I(l,t)和W_II(l,t)是为了解决两段投影数据冗余重叠问题而提出的权重因子，W_I(l,t)对应于第一段投影数据、W_II(l,t)对应于第二段投影数据。

2.根据权利要求1所述的方法，其特征在于，所述获取N段扫描投影数据，包括：利用N段直线分布式光源及N段直线分布式探测器采用静态扫描模式获取N段扫描投影数据。

3.根据权利要求2所述的方法，其特征在于，在获取N段扫描投影数据之前，还包括：对直线分布式光源之间的间距及夹角，直线分布式光源与直线分布式探测器之间的间距，直线分布式探测器之间的夹角进行标定。

4.根据权利要求1所述的方法，其特征在于，对所述N段平移后的投影数据进行加权滤波处理，得到N段加权滤波后的数据，包括：

分别将所述N段平移后的投影数据q_I(l,t)～q_N(l,t)乘上权重因子，用卷积核h沿着直线分布式探测器所在的直线方向做一维卷积得到N段加权滤波后的数据Q_I(l,t′_I)～Q_N(l,t′_N)。

5.根据权利要求4所述的方法，其特征在于，对所述N段平移后的投影数据进行加权处理，以消除多段投影数据冗余重叠所带来的重建伪影。

6.根据权利要求4所述的方法，其特征在于，对所述N段加权滤波后的数据进行反投影，得到所述N段扫描投影数据的重建图像f_I(x,y)～f_N(x,y)，包括：

沿直线分布式光源所在的直线方向反投影Q_I(l,t′_I)～Q_N(l,t′_N)，得到所述N段扫描投影数据的重建图像f_I(x,y)～f_N(x,y)。