[发明专利]一种CT-CBCT影像形变配准方法

权利要求书

1.一种CT-CBCT影像形变配准方法，其特征在于，具体操作步骤如下：

步骤(1.1)、对CBCT图像与CT图像进行预处理；

步骤(1.2)、建立图像A与B局部控制点并确定控制点的邻域与尺寸；

步骤(1.3)、对图像A的局部控制点邻域进行第一级评价，即暗区评价，比较图像A第i个控制点Q_Ai邻域灰度均值E_A与图像B第i个控制点Q_Bi邻域灰度均值E_B的大小；

若E_A≥m_EE_B，其中，m_E为常数，取值是：0.5≤m_E≤1.0，则判定图像A第i个控制点Q_Ai邻域为非暗区，对图像A第i个控制点邻域参照图像B第i个控制点邻域进行直方图校正配准，输出图像A局部控制点邻域形变矩阵M_1i，矩阵大小与D_Ai一致；

最后，汇总图像A的所有控制点邻域的形变矩阵映射到图像A中对应位置，得到图像A的形变矩阵同时完成了图像A与图像B的形变配准；

若E_A＜m_EE_B，其中，m_E为常数，取值是：0.5≤m_E≤1.0，则判定图像A第i个控制点Q_Ai邻域为暗区，此时需对图像A的控制点Q_Ai邻域进行第二级评价，即伪影评价，再比较图像A第i个控制点Q_Ai邻域灰度均方差与图像B第i个控制点Q_Bi邻域灰度均方差平方根的大小，

若其中m_σ为常数，取值是：0.1≤m_σ≤0.5，则判定图像A第i个控制点Q_Ai邻域为非伪影区，对图像A第i个控制点邻域参照图像B第i个控制点邻域进行线性灰度值映射，输出图像A局部控制点邻域形变矩阵M_2i，矩阵大小与D_Ai一致；

若其中m_σ为常数，0.1≤m_σ≤0.5，则判定图像A第i个控制点Q_Ai邻域为伪影区，则先对Q_Ai邻域灰度进行局部放大，放大系数为图像B邻域灰度均值与图像A邻域灰度均值的比值，再对Q_Ai邻域参照Q_Bi邻域进行直方图校正配准，输出图像A局部控制点邻域形变矩阵M_3i，矩阵大小与D_Ai一致；

最后，汇总图像A的所有控制点邻域的形变矩阵映射到图像A中对应位置，得到图像A的形变矩阵同时完成图像A与图像B的形变配准；

步骤(1.4)、通过将图像A的所有控制点邻域的形变矩阵映射到图像A中对应位置，得到图像A的形变矩阵后，再对所有待配准层的同时完成了图像CBCT图像进行局部校正配准，最终完成一组CT和一组CBCT图像的形变配准。

2.根据权利要求1所述的一种CT-CBCT影像形变配准方法，其特征在于，

在步骤(1.1)中，所述对CBCT图像与CT图像进行预处理具体操作方法如下：取CBCT图像与CT图像的交集区域并基于该交集区域分别裁剪CBCT图像与CT图像，同时将CBCT图像参照CT图像进行采样，得到同分辨率同尺寸CBCT图像与CT图像，得到待配准层的CBCT图像A与CT图像B的所有灰度信息作为输入，记图像A的宽度与高度均为W、H，则图像B的宽度与高度均为W、H，其中，1≤i≤H、1≤j≤W。