一种非刚性医学图像配准方法与流程

本发明属于图像处理领域，具体涉及应用于医学图像的非刚性配准方法。

背景技术：

图像配准(imageregistration)技术是指将两幅或者多幅在不同取景设备、不同取景时间或者不同取景环境等不同条件下获得的图像进行匹配的过程。图像配准被应用于许多方面，例如运动分析，视频压缩与编码，目标跟踪以及图像稳定等。现阶段非刚性配准的方法主要分为基于特征的配准方法和基于灰度的配准方法。基于特征的配准算法速度较快，但是对于特征不明显的图像有局限性，典型的算法有sift(scale-invariantfeaturetransform)方法；基于灰度的配准方法可以对小形变图像进行高精度的配准，但配准时间较长，典型的算法有demons方法，activedemons方法，引入了多分辨率策略的改进activedemons(miad)方法，基于b样条的自由形变方法等。

b样条(basicspline)方法最早由d.rueckert等应用于图像配准中，他们在乳房mr图像的非刚性配准中应用并发展了seungyonglee等提出的用于离散数据插值的多层级均匀b样条算法，提出了基于b样条的自由形变(free-formdeformation，ffd)配准方法，获得了较好的结果，后续的一些方法主要也都基于此方法上的一些改进，如：j.a.schnabel和d.rueckert等则在此基础上继续改善算法，提出了非均匀多层级ffd配准方法，提高了配准效率；张红颖等提出了一种基于区域相似性的自动寻找对应标记点的方法，利用多层级b样条进行插值，能得到比较快速精确的配准结果。marcomora等修改了原有的b样条参数模型，约束了控制点组成的网格，使得形变更加平滑。但是经典的基于多层级b样条的自由形变配准方法对于形变较大或局部扭曲较大的一些图像，无法得到较精确的配准结果。

技术实现要素：

针对经典的基于多层级b样条的自由形变(free-formdeformation，ffd)配准方法对于形变较大或局部扭曲较大的一些图像，无法很好地模拟出其形变场，无法得到较为准确的结果的问题，本发明提出一种非刚性医学图像配准方法。

本发明的一种非刚性医学图像配准方法，包括以下步骤：

步骤1：读取参考图像s和浮动图像m，并将s和m转化为double类型；

步骤2：对两幅图像进行高斯平滑处理，并对其进行仿射变换，返回仿射变换参数矩阵；

步骤3：在经典多层b样条自由形变配准方法的代价函数上引入控制点参数矩阵的frobenius范数项作为l2规范正则项，作为新的配准模型(l2-normffd)，在经典的基于多层级b样条的自由形变配准方法中，其代价函数由相似性测度和平滑约束两项组成，表达式为：

式中，ψ表示控制点参数矩阵，t表示配准的空间变换，i(xo，yo)表示参考图像，i(x，y)表示浮动图像，t(i(x，y))表示配准后的结果，λ表示平滑约束系数，csimilarity和csmooth(t)分别表示相似性测度和空间变换的平滑约束，在式(1)的基础上引入了控制点参数矩阵的frobenius范数项作为l2规范正则项，其表达式为：

式中，||ψ||f表示引入的控制点参数矩阵的frobenius范数项，λl为控制该项权重的系数。

步骤4：利用参考图像和配准结果计算客观评价指标，为方法做出评价。

所述步骤3具体包括以下步骤：

(3-1)设置l2-normffd模型的参数，包括选用的相似性测度、l2-norm正则项系数、平滑约束项系数、初始间距、b样条层数；

(3-2)设b样条层数为m(m>1)，开始非刚性变换：

(3-2-1)制作第1层初始形变网格，计算相似性测度，进行第1层b样条网格形变；

(3-2-2)将控制点参数矩阵代入能量函数，采用l-bfgs算法对能量函数进行优化求解，直到找到最小值；

(3-2-3)返回控制点参数矩阵；

(3-3)加密第1层b样条网格，开始第2层形变，其实现步骤同(3-2)；

(3-4)加密第m-1层b样条网格，开始第m层形变，其实现步骤同(3-2)，结束形变；

(3-5)将找到的最优控制点参数矩阵作用到浮动图像上，进行线性插值或者双三次插值还原图像，得到配准结果。

本发明相比于经典的基于多层b样条的自由形变配准方法，在提出的新方法的代价函数中引入了控制点参数矩阵的frobenius范数项作为l2规范正则项，通过将l2规范正则项的系数与平滑约束项系数进行协调一致性调整，既能保持图像拓扑结构，又能解决因平滑约束太强而导致精度较低的问题。通过引入l2规范正则项，可以对控制点网格的形变场进行空间限制，从而在一定程度上避免了过配准，能够对局部大扭曲形变进行配准。除此之外，l2范数还具有使优化求解变得快速且稳定的良好性质，增强了图像配准的鲁棒性。

附图说明

图1为肝脏mr图像：(a)浮动图像(b)参考图像；(c)原差值图像。

图2为肝脏mr图像的传统方法配准结果：(a)demons配准；(b)ad配准；(c)miad配准；(d)sift配准；(e)demons差值；(f)ad差值；(g)miad差值；(h)sift差值。

图3为肝脏mr图像配准结果和差值；(a)ffd配准；(b)ffd差值；(c)l2-normffd配准(d)l2-normffd差值。

图4为dir-lab数据库肺部4dctcase1的肺部示意图，

图5为配准前后最大呼气相位与最大吸气相位之间300组对应标记点的空间偏移对比图：

(a)为配准前最大呼气相位与最大吸气相位之间300组对应标记点的空间偏移；

(b)为配准后最大呼气相位与最大吸气相位之间300组对应标记点的空间偏移。

具体实施方式

下面结合附图及具体实施例对本发明作进一步详细的说明。

本发明实现医学图像配准依次经过以下步骤：

步骤1：读取参考图像s和浮动图像m，并将s和m转化为double类型；

步骤2：对两幅图像进行高斯平滑处理，并对其进行仿射变换，返回仿射变换参数矩阵；

步骤3：在经典多层b样条自由形变配准方法的代价函数上引入控制点参数矩阵的frobenius范数项作为l2规范正则项，作为新的配准模型(l2-normffd)。设置l2-normffd模型的参数，包括选用的相似性测度、l2-norm正则项系数、平滑约束项系数、初始间距、b样条层数。设b样条层数为m，开始非刚性变换：

将找到的最佳形变网格映射作用到浮动图像上；

进行线性插值或者双三次插值还原图像，得到配准结果；

步骤4：利用参考图像和配准结果计算客观评价指标，为方法做出评价。

本发明的优点主要体现在以下两方面：

1.提高图像配准的精确度

本发明能够相对精确地实现医学图像的配准。图1至图3为对肝脏mr图像的配准实验结果，表1为该实验下本发明在峰值信噪比(psnr,peaksignaltonoiseratio)，均方差(mse,meansquareerror)，相关系数(cc,correlationcoefficient)，互信息(mi,mutualinformation)和结构相似度(ssim,structuralsimilarity)5个指标上与几种传统非刚性配准方法的客观评价比较。图4为dir-lab数据库肺部4dctcase1的肺部示意图；图5为配准前后300组标记点空间偏移量对比，图5(a)为配准前最大呼气相位与最大吸气相位之间300对对应标记点的空间偏移，图5(b)为配准后最大呼气相位与最大吸气相位之间300对对应标记点的空间偏移。表2为10组4dct数据配准实验下本发明与其他部分算法的实验结果对比，其中第2列tre为配准前的对应标记点平均偏移量(括号内为偏移量的标准差)，第3列至第5列分别为经过mls(movingleast-squares)、ffd和l2-normffd三种方法配准后的对应标记点空间偏移量均值，括号内为空间偏移量的标准差。通过配准结果对比可以看出，本发明对图像配准的精确度有着相对明显的提高。

表1肝脏mr图像配准客观评价

表2dir-lab肺部4dct数据配准客观评价

2.加强图像配准的鲁棒性

图像配准方法的鲁棒性直接影响着其实际应用价值。本发明相对于经典的基于b样条的自由形变配准方法，能够使形变场避免发生折叠，更好保持住图像的拓扑结构，相对增强了算法的鲁棒性。

综上所述，本发明设计了一种应用于医学图像的非刚性配准方法(l2-normffd)。该方法基于多层级b样条理论和l2规范正则化理论，能够更精确更鲁棒地实现对医学图像的配准。

最后应说明的是，以上实施实例仅用以说明本发明的技术方案而非限制，尽管参照较佳的实施实例对本发明进行了详细的说明，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到的修改或等同替换，而不脱离本发明技术方案的精神和范围，都应涵盖在本发明的保护范围之内。