[发明专利]基于图形处理器的乳腺CT图像重建方法

摘要

一种基于图形处理器的乳腺CT图像重建方法，依次包括下述步骤：（1）数据接收与传输：将由扫描检测装置检测到的模拟信号经模/数转换器转换而成的数字信号传输到系统内存，并把系统内存中的数字信号传输到图形处理器的显存；（2）按并行的有限角度重建成像算法，在图形处理器中并行进行数据运算；在图形处理器中并行进行的数据运算为迭代计算；（3）判定步骤（2）的计算结果是否达到预期目标，如达到预期目标则执行步骤（4），否则重新回到步骤（2）；（4）图形处理器将迭代计算的结果传送至系统内存，进行图像后处理。本发明在图形处理器上进行重建运算，可非常大的提升迭代计算速度，快速准确地实现在有限角度下的乳腺CT图像重建。

主权项

1.一种基于图形处理器的乳腺CT图像重建方法，其特征在于依次包括下述步骤：（1）数据接收与传输：将由乳腺CT设备的扫描检测装置检测到的模拟信号经模/数转换器转换而成的数字信号经过系统总线传输到系统内存，然后根据数字信号的数据量大小申请显存，并把系统内存中的数字信号传输到图形处理器的显存中；所述数字信号是在有限角度重建成像方式下获得的信号；（2）按并行的有限角度重建成像算法，在图形处理器中并行进行数据运算；定义算子B和C，有限角度重建成像算法中GP迭代的过程为：B=T_F,C=FT_IF¹①⁰=_k②ⁱ⁺¹=C_k+(I-CB)ⁱ③其中Fourier空间的已知数据定义为_k，全部数据表示为，F表示Fourier变换，F¹表示Fourier逆变换，T_I和T_F是图像空间和频率空间的二值函数矩阵，I是单位矩阵；GP迭代最终以(1-λ_i)ⁿ的速度收敛到，其中，{λ_i}是CB的特征值，并且0<λ_i<1；在图形处理器中并行进行的数据运算为迭代计算，迭代计算包括下述步骤：（2－1）Fourier变换与逆变换；（2－2）求空域与频域的特征值；（2－3）求算子B、C的特征值；（3）利用预定条件判定步骤（2）的计算结果是否达到预期目标，如达到预期目标则执行步骤（4），否则重新回到步骤（2）进行迭代计算；（4）图形处理器将迭代计算的结果传送至系统内存，进行图像后处理。