[发明专利]一种具有不可压缩器官的医学图像配准方法及系统

权利要求书

1.一种具有不可压缩器官的医学图像配准方法，其特征在于，包括：

S1，接收固定图像和待配准图像，将所述固定图像和待配准图像根据上一次迭代获得的形变场进行预配准后，基于高斯牛顿算法和高效二阶最小化算法计算获得所述待配准图像向所述固定图像配准后的本次迭代的第一速度场；

其中，所述固定图像和待配准图像为包含相同的不可压缩的器官的医学图像；

S2，根据所述本次迭代的第一速度场和上一次迭代获得的速度场，基于Baker-Campbell-Hausdorff公式计算获得本次迭代的第二速度场；

S3，将所述本次迭代的第二速度场通过霍其-亥姆霍兹分解，获得所述本次迭代的第二速度场的无源场和调和场；同时计算所述调合场的自适应权重，通过所述调合场的自适应权重对所述调合场进行放大，将放大后的调合场与所述无源场叠加，获得本次迭代的速度场；

其中，所述将所述本次迭代的第二速度场通过霍其-亥姆霍兹分解，获得所述本次迭代的第二速度场的无源场和调和场，包括：

将所述第二速度场分解为无旋场无源场和调合场所述无旋场为标量势函数的梯度，所述无源场为向量势函数的旋度；

基于所述无旋场确定所述标量势函数；

基于所述无源场确定所述向量势函数；

最小化所述标量势函数和所述第二速度场的欧式距离，确定最终的无旋场；

最小化所述向量势函数和所述第二速度场的欧氏距离，确定最终的无源场；

基于所述最终的无旋场、所述最终的无源场和所述第二速度场，确定调合场。

S4，判断获知所述本次迭代的速度场使得配准能量场的能量小于预设阈值，将所述本次迭代的速度场通过指数变换进行逆变换，获得本次迭代的形变场，将本次迭代获得的形变场作为令所述待配准图像向所述固定图像配准的形变场。

2.根据权利要求1所述的方法，其特征在于，所述步骤S1之前还包括：若本次迭代为第一次迭代，则初始化速度场，使得初始的速度场为0，同时初始化所述形变场，使得初始的形变场为1。

3.根据权利要求1所述的方法，其特征在于，所述步骤S2之前还包括：对所述本次迭代的第一速度场进行高斯平滑，获取高斯平滑后的本次迭代的第一速度场。

4.根据权利要求3所述的方法，其特征在于，所述步骤S3之前还包括：对所述本次迭代的第二速度场进行高斯平滑，获取高斯平滑后的本次迭代的第二速度场。

5.根据权利要求1所述的方法，其特征在于，所述步骤S4中还包括，判断获知所述本次迭代的速度场使得配准能量场的能量大于预设阈值，则提取本次迭代的速度场和本次迭代的形变场，返回执行S1。

6.根据权利要求1所述的方法，其特征在于，所述步骤S4中还包括计算所述本次迭代的速度场的能量函数，具体为：

根据所述本次迭代的速度场对所述固定图像和待配准图像进行相似性测度计算；

计算本次迭代的速度场和所述上一次迭代获得的速度场的差异项，将所述相似性测度和所述差异项进行叠加，获得所述本次迭代的速度场的能量函数。

7.根据权利要求1所述的方法，其特征在于，所述步骤S3中，所述自适应权重通过以下公式计算获得；

式中，k为自适应权重，a为常数，Jac(.)为雅克比行列式，std(.)为标准差函数，表示第i次迭代中分辨率为j的第二速度场。