一种头颈部多模态医学图像自动配准方法与流程

本发明涉及医学图像处理的技术领域，尤其涉及到一种头颈部多模态医学图像自动配准方法。

背景技术：

ct(computedtomography，电子计算机断层扫描)成像可以更好显示骨骼信息，mri(magneticresonanceimaging,磁共振成像)对软组织成像明显，两种图像融合可以精准定位病灶位置，提高放射治疗的效果。

目前，放射治疗计划基本上都是基于ct图像的，但是由于ct图像对于软组织成像不明显，导致物理师在制定放射治疗计划时，无法准确勾勒出靶区。为提高放射治疗的效果，需要结合mr图像，将mr图像配准到ct图像，使mr图像在空间上与ct图像一致，使得物理师在制定放疗计划时，更加准确勾画靶区，减少放疗时对健康组织的副作用。

然而，目前医院所使用的配准软件，都需要物理师手动勾画出感兴趣区域(见图9)，再基于感兴趣区域进行配准。这种人工干预的配准方法，操作更加复杂，并会影响配准结果。

技术实现要素：

本发明的目的在于克服现有技术的不足，提出一种能够自动勾画感兴趣区域，精准度高的头颈部多模态医学图像自动配准方法。

为实现上述目的，本发明所提供的技术方案为：

一种头颈部多模态医学图像自动配准方法，包括以下步骤：

s1：导入固定图像和浮动图像；其中，ct图像作为配准的固定图像，mr图像作为配准的浮动图像；

s2：将步骤s1导入的两种模态的三维体数据下采样至低尺度，对低尺度下的体数据进行滤波去噪预处理；

s3：从ct和mr图像中提取精准的表面轮廓点集；

s4：通过icp算法对提取的ct和mr图像的表面轮廓点集进行表面匹配，输出匹配的误差和两点集的变换矩阵；

s5：根据变换矩阵对mr图像进行刚性变换，得到经变换的mr图像，利用变换矩阵对mr图像的表面轮廓点集进行重采样，得到经变换的mr图像的表面轮廓点集；

s6：求取ct图像与经变换的mr图像的表面轮廓点集的重叠区域点集；

s7：根据ct图像与经变换的mr图像的表面轮廓点集的重叠区域点集，得到三维的感兴趣区域；

s8：利用获取的感兴趣区域，进行刚性配准和非刚性配准。

进一步地，所述步骤s2中的预处理还包括对mr图像进行偏移场矫正，以及灰度归一化操作，以解决偏移场的存在导致图像灰度不均匀的问题，减少配准时的干扰。

进一步地，所述步骤s3中从ct和mr图像中提取精准的表面轮廓点集的具体步骤如下：

s3-1：利用mc算法提取预处理后的ct图像和mr图像的表面轮廓；

s3-2：步骤s31提取的表面轮廓比较粗糙，对粗提取的表面轮廓进行最大连通域计算，得到ct图像和mr图像精准的表面轮廓，即从ct、mr图像中提取精准的表面轮廓点集。

进一步地，所述步骤s6求取ct图像与经变换的mr图像的表面轮廓点集的重叠区域点集的具体步骤为：根据匹配的误差，定义一个以k倍误差大小为搜索阈值，遍历ct图像表面轮廓点集中的每个点，如果一点在阈值范围内能找到经变换的mr表面轮廓点集中的点，则该点为重叠区域的ct图像表面轮廓点，找到的mr表面轮廓点也为重叠区域的点，并保存下来；否则，该点不为重叠区域的点，予以剔除；遍历所有ct表面轮廓点，即可同时得到重叠区域的ct和mr轮廓点集。

进一步地，在步骤s7中，所述获取三维感兴趣区域的具体步骤如下：

s7-1：对ct图像的表面轮廓点集的重叠区域点集，以及经变换的mr图像的表面轮廓点集的重叠区域点集进行表面重建；

s7-2：根据下采样后的ct图像，经变换的mr图像的大小、像素间距，创建一个与ct图像一样大小的平面，与经变换的mr图像一样大小的平面；

s7-3：使用创建的平面与表面重建模型相交，得到交线，将z方向的像素间距，作为平面的上移间距，得到每张切片上的交点；

s7-4：由于图像是经过下采样的，所以有的切片上的交点是不连续的，需要对每张切片上的交点按照顺序进行首尾相连，形成封闭的轮廓；

s7-5：对封闭轮廓进行填洞处理；

s7-6：使用数学形态学中的闭运算对填洞图像进行处理，得到感兴趣区域的二值图像；

s7-7：对感兴趣区域的二值图像进行上采样，使得二值图像尺度与原体数据相同，然后使用上采样后的二值图像和原始ct、mr图像进行按位与操作，得到ct图像的感兴趣区域和经变换的mr图像的感兴趣区域。

进一步地，在步骤s8中，对ct图像的感兴趣区域和经变换的mr图像的感兴趣区域进行下采样，利用下采样后的感兴趣区域进行刚性配准，得到刚性配准的变换参数，利用变换参数，对未采样的mr图像进行重采样，得到经过刚性配准后的mr图像；然后通过demons算法对ct图像以及经过刚性配准后的mr图像进行非刚性配准进行局部形变，输出经过icp匹配、刚性配准后的mr图像配准到ct图像的形变场，使用形变场对mr图像进行形变，得到非刚性配准后的mr图像。

与现有技术相比，本方案原理和优点如下：

1.通过提取固定图像和浮动图像的表面轮廓点集进行表面匹配，获取表面匹配时的变换矩阵、以及固定图像和浮动图像的重叠区域点集，缩小配准寻优范围，提高精准度。

2.对重叠区域点集，进行表面重建。通过表面重建模型与平面相交的方法，提取重叠区域点集对应的三维图像区域,作为配准的感兴趣区域，进而对感兴趣区域进行自动配准，减少人工参与，简化配准流程。

附图说明

图1为本发明一种头颈部多模态医学图像自动配准方法的流程图；

图2为本发明中ct图像表面轮廓点集与mr图像表面轮廓点集进行表面匹配的示意图；

图3为本发明中ct图像表面轮廓点集与匹配后的mr图像表面轮廓点集的重叠区域示意图；

图4为本发明中ct图像表面轮廓点集与平面相交所得交线的示意图；

图5为本发明中mr图像表面轮廓点集与平面相交所得交线的示意图；

图6为为本发明中获取感兴趣区域的二值图像；

图7为为本发明中ct图像与感兴趣区域的二值图像相与后得出的ct图像的感兴趣区域的示意图；

图8为本发明中的配准结果示意图；

图9为现有技术中的瓦里安配准界面图。

具体实施方式

下面结合具体实施例对本发明作进一步说明：

参见图1所示，本实施例所述的一种头颈部多模态医学图像自动配准方法，包括以下步骤：

s1：导入固定图像和浮动图像；其中，ct图像作为配准的固定图像，mr图像作为配准的浮动图像；

s2：将步骤s1导入的两种模态的三维体数据下采样至低尺度，对低尺度下的体数据进行滤波去噪，对mr图像进行偏移场矫正，以及灰度归一化操作，以解决偏移场的存在导致图像灰度不均匀的问题，减少配准时的干扰；

s3：从ct和mr图像中提取精准的表面轮廓点集，具体步骤如下：

s3-1：利用mc算法提取预处理后的ct图像和mr图像的表面轮廓；

s3-2：步骤s31提取的表面轮廓比较粗糙，对粗提取的表面轮廓进行最大连通域计算，得到ct图像和mr图像精准的表面轮廓，即从ct、mr图像中提取精准的表面轮廓点集；

提取到ct、mr图像精准的表面轮廓点集后，进入步骤s4；

s4：通过icp算法对提取的ct和mr图像的表面轮廓点集进行表面匹配(如图2所示)，输出匹配的误差和两点集的变换矩阵；

s5：根据变换矩阵对mr图像进行刚性变换，得到经变换的mr图像，利用变换矩阵对mr图像的表面轮廓点集进行重采样，得到经变换的mr图像的表面轮廓点集；

s6：如图3所示，求取ct图像与经变换的mr图像的表面轮廓点集的重叠区域点集，具体为：

根据匹配的误差，定义一个以k倍误差大小为搜索阈值，遍历ct图像表面轮廓点集中的每个点，如果一点在阈值范围内能找到经变换的mr表面轮廓点集中的点，则该点为重叠区域的ct图像表面轮廓点，找到的mr表面轮廓点也为重叠区域的点，并保存下来；否则，该点不为重叠区域的点，予以剔除；遍历所有ct表面轮廓点，即可同时得到重叠区域的ct和mr轮廓点集；

s7：根据ct图像与经变换的mr图像的表面轮廓点集的重叠区域点集，得到三维的感兴趣区域，具体步骤如下：

s7-1：对ct图像的表面轮廓点集的重叠区域点集，以及经变换的mr图像的表面轮廓点集的重叠区域点集进行表面重建；

s7-2：根据下采样后的ct图像，经变换的mr图像的大小、像素间距，创建一个与ct图像一样大小的平面，与经变换的mr图像一样大小的平面；

s7-3：使用创建的平面与表面重建模型相交，得到交线(如图4和5所示)，将z方向的像素间距，作为平面的上移间距，得到每张切片上的交点；

s7-4：由于图像是经过下采样的，所以有的切片上的交点是不连续的，需要对每张切片上的交点按照顺序进行首尾相连，形成封闭的轮廓；

s7-5：对封闭轮廓进行填洞处理；

s7-6：使用数学形态学中的闭运算对填洞图像进行处理，得到感兴趣区域的二值图像，如图6所示；

s7-7：对感兴趣区域的二值图像进行上采样，使得二值图像尺度与原体数据相同，然后使用上采样后的二值图像和原始ct、mr图像进行按位与操作，得到ct图像的感兴趣区域和经变换的mr图像的感兴趣区域，如图7所示；

s8：对ct图像的感兴趣区域和经变换的mr图像的感兴趣区域进行下采样，利用下采样后的感兴趣区域进行刚性配准，得到刚性配准的变换参数，利用变换参数，对未采样的mr图像进行重采样，得到经过刚性配准后的mr图像；然后通过demons算法对ct图像以及经过刚性配准后的mr图像进行非刚性配准进行局部形变，输出经过icp匹配、刚性配准后的mr图像配准到ct图像的形变场，使用形变场对mr图像进行形变，得到非刚性配准后的mr图像，如图8所示。

本实施例通过提取固定图像和浮动图像的表面轮廓点集进行表面匹配，获取表面匹配时的变换矩阵、以及固定图像和浮动图像的重叠区域点集，缩小配准寻优范围，提高精准度。另外，对重叠区域点集，进行表面重建。通过表面重建模型与平面相交的方法，提取重叠区域点集对应的三维图像区域,作为配准的感兴趣区域，进而对感兴趣区域进行自动配准，减少人工参与，简化配准流程。

以上所述之实施例子只为本发明之较佳实施例，并非以此限制本发明的实施范围，故凡依本发明之形状、原理所作的变化，均应涵盖在本发明的保护范围内。