[发明专利]心肌SPECT灌注显像的心脏轴向移动校正方法

权利要求书

1.一种心肌SPECT灌注显像的心脏轴向移动校正方法，其特征在于，包括以下步骤：

(1)投影数据采集完成后，初步估计各个投影轴向位移的大小；

(2)根据各个投影的轴向位移初步估计值，对各个投影进行初步校正；

(3)使用经初步轴向位移校正的投影数据，重建三维图像；

(4)将重建图像根据实际投影的采集条件进行再投影；

(5)对原始投影和对应的再投影进行分析，计算出实际投影的轴向位移；

(6)根据上一步的实际投影的轴向位移值，对各个投影进行校正；

(7)使用经上一步轴向位移校正的投影数据，重建三维图像；

(8)重复第(4)-第(5)步，当对应校正值最大差别小于要求值时，校正完成；否则，重复第(6)步-第(8)步。

2.根据权利要求1所述的心肌SPECT灌注显像的心脏轴向移动校正方法，其特征在于：所述步骤(1)具体为：首先将各个角度的投影图像沿横向求和，然后对这个一维序列进行互相关分析；

一维序列数据P_n与P_n-1的互相关函数可表示为：CCn(d)=Σj=1MPn(j)*Pn-1(j+d)]]>

其中：d代表偏移，M代表每组一维数据中数据的总个数；-k≤d≤k，k为最大可能的位移，一般设为10个像素；当j+d＜1或j+d＞M时，令P_n-1(j+d)＝0；互相关函数CC_n(d)表示P_n与偏移了d的P_n-1的相互匹配程度，CC_n(d)越大，匹配得越好；因此对CC_n(d)进行抛物线拟合，求出抛物线顶点对应的d，即可计算出位移。

3.根据权利要求1所述的心肌SPECT灌注显像的心脏轴向移动校正方法，其特征在于：所述步骤(2)具体为：通过线性插值的方法，实现非像素大小整数倍的图像的轴向位移，插值后对图像的总计数进行归一，即每个像素的计数乘以原始投影的总计数后，除以插值后图像的总计数。

4.根据权利要求1所述的心肌SPECT灌注显像的心脏轴向移动校正方法，其特征在于：所述步骤(3)具体为：图像重建采用滤波反投影法，使用Butterworth滤波器重建图像，不加衰减校正。

5.根据权利要求1所述的心肌SPECT灌注显像的心脏轴向移动校正方法，其特征在于：所述步骤(4)具体为：通过拉登变换计算实际采集条件的、上一步重建图像的投影。

6.根据权利要求1所述的心肌SPECT灌注显像的心脏轴向移动校正方法，其特征在于：所述步骤(5)具体为：对每个角度投影的计算值与测量值根据互信息进行图像配准，计算出各个角度投影沿轴向位移的距离。

7.根据权利要求1所述的心肌SPECT灌注显像的心脏轴向移动校正方法，其特征在于：所述步骤(6)具体为：通过拉登变换计算实际采集条件的、上一步重建图像的投影。

8.根据权利要求1所述的心肌SPECT灌注显像的心脏轴向移动校正方法，其特征在于：所述步骤(7)具体为：图像重建采用滤波反投影法，使用Butterworth滤波器重建图像，不加衰减校正。

9.根据权利要求1所述的心肌SPECT灌注显像的心脏轴向移动校正方法，其特征在于：所述步骤(8)具体为：重复第(4)-第(5)步，比较两次校正值，当对应校正值最大差别小于0.5mm时，校正完成。