一种基于GPU并行计算的不同分辨率医学图像配准系统的制作方法

本发明涉及图像处理技术领域，特别是涉及一种基于gpu并行计算的不同分辨率医学图像配准系统。

背景技术：

医学图像配准是图像处理研究领域中的一项重要且复杂的任务，由于医学图像会来自不同的采集设备、不同的扫描时间、不同的拍摄角度、不同的研究对象等。因此，往往在进行图像融合、图像分割等后续处理前需要先进行图像配准已达到针对同一解剖结构在不同的采集图像中对应于空间中的同一位置。

传统的图像配准方法主要使用cpu，并使用包括基于特征的方法，该方法往往需要利用边缘、表面、像素等信息进行特征提取，配准的成功与否很大程度上依赖于提取的特征。这就导致其存在注入适用性差、处理速度慢、配准评估不统一等问题。

现有专利文献cn111091567a公开了医学图像配准方法，该方法包括获取目标部位的参考医学图像序列以及待配准医学图像序列；确定待配准医学图像序列的各个图像层面中的血管边界，以得到各个图像层面的边界子图；分别对参考医学图像序列的至少一个图像层面以及待配准医学图像序列的至少一个图像层面进行血管边界的采样，得到第一边界点集合以及第二边界点集合；基于第一边界点集合以及第二边界点集合中对应边界点的偏移量，确定最优变换；根据最优变换对各个边界子图进行调整，以对所述待配准医学图像序列进行配准。由此可见，该方法需要提取血管边界作为特征点，当信息量较少时，且配准精度将受到较大的限制，另外，该方法中待配准图像和参考图像的大小一致，因此其不适用于复杂工作场景，如穿刺导航和不同成像设备间配准。

技术实现要素：

本发明所要解决的技术问题是提供一种基于gpu并行计算的不同分辨率医学图像配准系统，能够在不提取特征点的前提下实现医学图像的配准。

本发明解决其技术问题所采用的技术方案是：提供一种基于gpu并行计算的不同分辨率医学图像配准系统，包括：第一图像获取模块，用于获取参考图像；第二图像获取模块，用于获取待配准图像；控制点生成模块，用于在空间中生成多个初始变形参数为零的控制节点；配准模块，用于进行不同尺度下全局配准并通过曲线b样条拟合的方式拟合所述控制节点上的变形参数。

所述带配准图像的分辨率和所述参考图像的分辨率不同。

所述配准模块包括仿射配准和自由变形配准，所述仿射配准是指通过优化所述参考图像的放射变形参数，最小化所述参考图像和待配准图像差异以消除图像全局仿射变形，所述自由变形配准是指通过b样条曲线拟合所述控制节点上的变形参数。

所述仿射配准具体为：给定初始仿射变换参数x，将所述仿射变换参数x由cpu复制到gpu，对所述待配准图像使用仿射变换参数x进行仿射变换，记为imoving，在gpu上计算惩罚函数f，f为参考图像iref和仿射变换后的待配准图像imoving的差异，表示为：f＝0.5∑(iref-imoving)²；在gpu上通过链式求导法则计算惩罚函数f对仿射配准参数x的偏导数，记为j，将其传递到cpu上，再用levenberg-marquardt算法求迭代过程仿射配准参数x的增量dh，具体为：(j^t*j+mu*i)*dh＝j^t*f，其中，mu为迭代中根据惩罚函数f和j进行调整的参数，i为单位矩阵，dh通过cholesky分解求得，以x＝x+dh更新仿射配准参数x，当满足迭代终止条件，迭代终止。

根据仿射配准计算求得仿射变换参数时，还对待配准图像imoving进行缩放并通过插值计算，使待配准图像imoving的分辨率和参考图像iref一致。

所述自由变形配准具体为：给定常数α0为初始步长；将控制节点上的变形参数β初始化为0，并由cpu复制到gpu，在gpu中对图像进行b样条自由变形，得到变形后的图像，记为imoving'，在gpu上计算惩罚函数f'，f'为参考图像iref和变形后图像imoving'的差异，表示为：f＝0.5∑(iref-imoving')²；在gpu上通过链式求导法则计算惩罚函数f'对变形参数β偏导数，记为j'，通过线搜索方法找到最佳步长，记为α，以x'＝x'-αj更新步长x'，当满足迭代终止条件，迭代终止。

有益效果

由于采用了上述的技术方案，本发明与现有技术相比，具有以下的优点和积极效果：本发明通过控制节点的变形参数直接对原始二维或三维图像配准，不用提取特征点，无信息损失，配准精度高。本发明使用高分辨率术前图像做为基准参考图像，并通过插值配准与术中低分辨率图像进行配准，实现不同分辨率医学图像的配准，可以有效降低了术中成像扫描的ct辐射量和扫描时间。本发明还使用gpu进行并行计算，使得配准效率高，适用于更广泛医学场景。

附图说明

图1是本发明的结构方框图。

具体实施方式

下面结合具体实施例，进一步阐述本发明。应理解，这些实施例仅用于说明本发明而不用于限制本发明的范围。此外应理解，在阅读了本发明讲授的内容之后，本领域技术人员可以对本发明作各种改动或修改，这些等价形式同样落于本申请所附权利要求书所限定的范围。

本发明的实施方式涉及一种基于gpu并行计算的不同分辨率医学图像配准系统，如图1所示，包括：第一图像获取模块、第二图像获取模块、控制点生成模块、配准模块和输出模块。

本实施方式中，第一图像获取模块，用于获取参考图像。该第一图像获取模块通常在对病人术前进行全面扫描时使用，例如可以使用ct，mri或超声波成像等方式使用较高的分辨率和较小的成像层间距对病人进行全面三维术前扫描，获得二维或三维的参考图像，记为iref。

第二图像获取模块，用于获取待配准图像。该第二图像获取模块通常在手术中使用，例如在消融/穿刺手术中对病人进行快速扫描或使用较大的成像层间距进行扫描，获得低分辨率的二维或三维图像作为待配准图像，记为imoving。

所述控制点生成模块前还包括输入模块，所述输入模块用于将所述参考图像和待配准图像输入至gpu中，通过gpu进行并行计算。

控制点生成模块，用于在空间中生成多个初始变形参数为零的控制节点，通过改变这些控制节点的变形参数，能够使得带配准图像和参考图像之间的偏差缩小，从而实现配准的目的。

配准模块，用于进行不同尺度下全局配准并通过曲线拟合的方式拟合所述控制节点上的变形参数。由于本实施方式中，参考图像的分辨率高于待配准图像，因此在配准模块中还需要对参考图像进行插值处理，具体处理方式为：根据仿射变换参数，对待配准图像imoving进行缩放和插值，使得所述待配准图像和参考图像具有相同的分辨率和矩阵大小。当两者处于相同的分辨率时，使用b样条曲线拟合控制节点上的变形参数，并使用最速下降法和线搜索方法对变形参数进行非线性优化以最小化待配准图像imoving和参考图像iref之间的偏差，从而实现配准的目的。

本实施方式中包括仿射配准和自由变形配准，所述仿射配准是指通过优化所述参考图像的放射变形参数，最小化所述参考图像和待配准图像差异以消除图像全局仿射变形，所述自由变形配准是指通过b样条曲线拟合所述控制节点上的变形参数。

所述仿射配准具体为：给定仿射配准参数x，将所述仿射配准参数x由cpu复制到gpu，对所述参考图像和待配准图像采用不同降采样系数进行降采样；对所述待配准图像使用仿射配准参数x进行仿射变换，记为imoving，在gpu上计算惩罚函数f，f为参考图像iref和仿射变换后的待配准图像imoving的差异，表示为：f＝0.5∑(iref-imoving)²；在gpu上通过链式求导法则计算惩罚函数f对仿射配准参数x的偏导数，记为j，将其传递到cpu上，再用levenberg-marquardt算法求迭代过程仿射配准参数x的增量dh，具体为：(j^t*j+mu*i)*dh＝j^t*f，其中，mu为迭代中根据惩罚函数f和j进行调整的参数，i为单位矩阵，dh通过cholesky分解求得，更新仿射配准参数x＝x+dh，当满足迭代终止条件，迭代终止。

所述自由变形配准具体为：给定常数α0为步长x'；将控制节点上的变形参数β初始化为0，并由cpu复制到gpu，在gpu中对图像进行b样条自由变形，得到变形后的图像，记为imoving'，在gpu上计算惩罚函数f'，f'为参考图像iref和变形后图像imoving'的差异，表示为：f＝0.5∑(iref-imoving')²；在gpu上通过链式求导法则计算惩罚函数f'对变形参数β偏导数，记为j'，通过黄线搜索方法(包括但不限于wolfe-powell方法，黄金分割搜索方法)最佳步长，记为α，更新步长x'＝x'-αj，当满足迭代终止条件，迭代终止。

值得一提的是，本实施方式还可以使用其他的非线性优化方法实现对变形参数进行非线性优化，比如lm算法和高斯-牛顿法等。

输出模块，用于将控制节点上的变形参数的值和配准后的图像从gpu复制到cpu并进行输出。

不难发现，本发明通过控制节点的变形参数直接对原始二维或三维图像配准，不用提取特征点，无信息损失，配准精度高。本发明使用高分辨率术前图像做为基准参考图像，并通过插值配准与术中低分辨率图像进行配准，实现不同分辨率医学图像的配准，并且降低了术中成像扫描ct辐射量和扫描时间。本发明还使用gpu进行并行计算，使得配准效率高，适用于更广泛医学场景。