[发明专利]层析成像的重建方法

说明书

技术领域

本发明是关于层析图像的重建(reconstruct)方法，特别是关于使用图形处理单元(GPU，GraphicProcessingUnit)的层析成像的重建方法。

背景技术

放射线层析成像，如数字X线层析成像是医学图像领域中的重要组成部分。通过在有限的角度内获取多幅二维投影，放射线层析成像可基于这些二维投影重建出容积(volume)图像，从而提供三维的图像信息用于医学或生物学等领域的图像分析。然而，尽管层析成像的重建效果被广泛认可，但重建过程中庞大的计算量制约了重建速度，限制了层析成像的实际应用。

常见的层析成像在重建中主要涉及解析和/或迭代两类高计算量的运算，这些运算主要由中央处理器单元(CPU)承担而图形处理单元仅分担一小部分。由于CPU、GPU的根本设计目的不同，两者在计算中所扮演的角色也不尽相同：例如CPU为了更为通用，不断的增加缓存和流程控制单元，而GPU则是为了完成复杂的图形渲染工作，需要强大的并行计算能力，因此具有更多的计算单元。此外，CPU主要用于不规则的数据结构，其存取模式不可预测，需要不断增加缓存来解决存储器的延迟问题；而GPU则使用单程序多数据执行模式、规则的数据结构，使得存储可预测，从而可以有更多的空间用于增加计算单元。

因此，如何有效利用CPU与GPU的性能差异，使两者优缺互补从而改善层析成像的重建速度成为业界的研究热点。

发明内容

本发明的目的之一在于提供一层析成像的重建方法，其在不增加现有层析成像装置硬件结构的前提下，大幅度加快重建的速度。

根据本发明的一实施例，一层析成像的重建方法包含：获取投影数据；图形处理单元使用开放计算语言对该投影数据进行滤波；及该图形处理单元使用该开放计算语言对所滤波的投影数据进行反投影，对经该反投影而重建的容积图像进行归一化处理。

在本发明的部分实施例中，对该投影数据进行滤波进一步包含建立斜坡滤波器。对该投影数据进行滤波进一步包含对该投影数据分别按照投影采集的顺序进行正序和逆序扩充。例如，由该图形处理单元对该投影数据进行正序扩充，或由中央处理单元对该投影数据进行逆序扩充，或由中央处理单元与该图形处理单元协作对该投影数据进行正序或逆序扩充。使快速傅立叶变换后的滤波器数据与快速傅立叶变换后的扩充投影数据相乘而得到滤波后的频域结果。在本发明的一些实施例中，对该投影数据进行滤波进一步包含对投影数据进行加权，并对加权后的投影数据进行滤波。根据本发明的部分实施例，对所滤波的投影数据进行反投影进一步包含由中央处理器单元根据预设几何关系生成投影矩阵，基于该投影矩阵生成一体素在二维投影的投影点并获得投影值，其中射线方向由射线源和该体素的位置决定。获得投影值可由该图形处理单元通过双线性插值获得。根据本发明的一些实施例，反投影进一步包含加权所获得的投影值并累加加权的投影值。此外，在本发明的部分实施例中，归一化处理是根据射线穿过各体素的次数进行归一化。

本发明实施例所提供的层析成像的重建方法可充分利用现有硬件中中央处理器与图形处理单元的结构特点和性能优势，并结合开放计算语言的操作特点，加快数据处理的速度，节约中央处理器与图像处理单元间数据传输的时间，从而在不增加任何硬件成本的前提下大幅度提高层析成像的重建速度。

附图说明

图1所示是根据本发明一实施例的层析成像中的扫描投影示意图

图2所示是根据图1中的扫描投影示意图得到的相应几何模型

图3所示是根据本发明一实施例的组合使用CPU与GPU的层析图像的重建方法的流程示意图

图4所示是根据本发明一实施例的反投影方法的流程示意图

具体实施方式

为更好的理解本发明的精神，以下结合本发明的部分优选实施例对其作进一步说明。

在目前的层析图像重建过程中，庞杂的数据计算，如滤波所涉及的运算主要是在中央处理器单元上完成的。相应的，层析图像重建的速度瓶颈也反映在中央处理器单元的运算速度上。为提高重建速度当然可以选择具有更高运算能力的中央处理单元，但这必将导致成本的大幅增加。

本发明的实施例提供了解决上述问题的另外一种途径，其通过中央处理单元与图形处理单元的优势互补合作，并巧妙的利用开放计算语言(OpenCL，opencomputinglanguage)的特点，可在不改变现有层析图像装置硬件配置的基础上大幅度提高重建速度。