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

说明书

本发明提供了一种CT‑CBCT影像形变配准方法。属于医学图像处理技术领域，具体步骤：对CBCT图像与CT图像进行预处理；建立图像A与B局部控制点并确定控制点的领域与尺寸；比较图像A第i个控制点Q_Ai邻域灰度均值E_A与图像B第i个控制点Q_Bi邻域灰度均值E_B的大小；完成一组CT和一组CBCT图像的形变配准。本发明对控制点邻域灰度校正配准采用两级评价判断出暗区与伪影区，从而根据是否暗区以及是否伪影区有针对性的对两组图像进行局部灰度校正配准；另外，本发明基于参与配准的两组图像利用梯度变化特征并采用网格划分机制建立了局部控制点，基于控制点邻域进行局部灰度校正配准避免了全局盲配。

技术领域

本发明涉及医学图像处理技术领域，具体涉及一种CT-CBCT影像形变配准方法，特别是涉及一种基于控制点邻域灰度校正的CT-CBCT图像形变配准方法。

背景技术

医学图像配准是精准放射治疗的关键技术之一，并为病情诊断、靶区勾画、精准摆位和精准治疗提供依据。多模态图像形变配准能够为医生和物理师提供更加丰富的医学图像信息。

传统的形变配准方法在处理CT-CBCT形变配准问题时往往由于图像对比度不足以及灰度不连续性等问题导致CT-CBCT图像误配或错配现象的发生，为解决这些问题时，国内外研究学者提出全局强度校正方法来对CBCT图像进行预处理，目的是将CBCT-CT多模态的形变配准转换成同模态的形变配准问题，但是在实际的CT-CBCT形变配准研究中发现伪影总是局部存在的，这会破坏原本的强度或灰度分布对应关系，因此需要对存在伪影的CBCT图像进行局部灰度校正，边对CBCT图像进行校正边对CT-CBCT图像进行匹配从而完成CT-CBCT的最终配准。

本发明克服了CBCT伪影带来的强度分布或灰度分布不对应致使CT-CBCT图像误配或错配的不足，不是盲目地对CT-CBCT图像进行全局配准，而是对CT-CBCT图像进行局部灰度校正配准，首先基于CBCT图像的灰度特征并利用网格划分机制构建了局部控制点，然后对控制点邻域灰度校正配准采用两级评价判断出暗区与伪影区，从而根据是否暗区以及是否伪影区有针对性的对两组图像进行局部灰度校正配准，最后对所有待配准层的CBCT图像进行局部灰度校正配准得到所有层CBCT图像的形变矩阵，完成了整组CT与CBCT图像的形变配准。

发明内容

本发明要解决的技术问题为：克服CBCT伪影带来的强度分布或灰度分布不对应致使CT-CBCT图像误配或错配的不足，提出了一种基于控制点邻域灰度校正的CT-CBCT图像形变配准方法，对控制点邻域灰度采用两级评价判断出暗区与伪影区，从而根据是否暗区以及是否伪影区进行相应的局部灰度校正，得到CBCT图像的形变矩阵，完成CT-CBCT图像形变配准。

本发明的技术方案是：一种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一致；