一种CT-CBCT影像形变配准方法

摘要

本发明提供了一种CT‑CBCT影像形变配准方法。属于医学图像处理技术领域，具体步骤：对CBCT图像与CT图像进行预处理；建立图像A与B局部控制点并确定控制点的领域与尺寸；比较图像A第i个控制点QAi邻域灰度均值EA与图像B第i个控制点QBi邻域灰度均值EB的大小；完成一组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

若E

最后，汇总图像A的所有控制点邻域的形变矩阵映射到图像A中对应位置，得到图像A的形变矩阵

若E

若

若

最后，汇总图像A的所有控制点邻域的形变矩阵映射到图像A中对应位置，得到图像A的形变矩阵

步骤(1.4)、通过将图像A的所有控制点邻域的形变矩阵映射到图像A中对应位置，得到图像A的形变矩阵

进一步的，在步骤(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。

进一步的，在步骤(1.2)中，所述建立图像A与B局部控制点并确定控制点的领域与尺寸的具体操作步骤如下：

(1.2.1)、基于图像A构建局部控制点，得到图像A的所有控制点集Q

(1.2.2)、选取图像A局部控制点Q

进一步的，

对图像A建立R×R均匀网格，其中，3≤R≤10，R为正整数，在每个网格中寻找梯度变化明显的点作为图像A的控制点，比较R分别取3、4、5、6、7、8、9、10这八种得到的形变矩阵

进一步的，在步骤(1.3)中，所述的m

进一步的，在步骤(1.4)中，得到对当前层CBCT图像即图像A的形变矩阵P，再对所有待配准层的CBCT图像进行局部灰度校正配准得到所有层CBCT图像的形变矩阵

本发明的有益效果是：(1)、克服了CBCT伪影带来的强度分布或灰度分布不对应致使CT-CBCT图像误配或错配的不足，对控制点邻域灰度校正配准采用两级评价判断出暗区与伪影区，从而根据是否暗区以及是否伪影区有针对性的对两组图像进行局部灰度校正配准，避免了对CBCT图像全局盲目校正；(2)、基于参与配准的两组图像利用梯度变化特征并采用网格划分机制建立了局部控制点，基于控制点邻域进行局部灰度校正配准避免了全局盲配。

附图说明

图1是本发明的结构流程图；

图2中，(a)为CT与CBCT图像；(b)为本发明的形变配准输出图像。

具体实施方式

为了更清楚地说明本发明的技术方案，下面结合附图对本发明的技术方案做进一步的详细说明：

如图1所述，一种CT-CBCT影像形变配准方法，具体操作步骤如下：

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

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

步骤(1.3)、对图像A的局部控制点邻域进行第一级评价，即暗区评价，比较图像A第i个控制点Q

若E

最后，汇总图像A的所有控制点邻域的形变矩阵映射到图像A中对应位置，得到图像A的形变矩阵

若E

若

若

最后，汇总图像A的所有控制点邻域的形变矩阵映射到图像A中对应位置，得到图像A的形变矩阵

步骤(1.4)、通过将图像A的所有控制点邻域的形变矩阵映射到图像A中对应位置，得到图像A的形变矩阵

进一步的，在步骤(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。

进一步的，在步骤(1.2)中，所述建立图像A与B局部控制点并确定控制点的领域与尺寸的具体操作步骤如下：

(1.2.1)、基于图像A构建局部控制点，得到图像A的所有控制点集Q

(1.2.2)、选取图像A局部控制点Q

进一步的，

对图像A建立R×R均匀网格，其中，3≤R≤10，R为正整数，在每个网格中寻找梯度变化明显的点作为图像A的控制点，比较R分别取3、4、5、6、7、8、9、10这八种得到的形变矩阵

进一步的，在步骤(1.3)中，所述的m

进一步的，在步骤(1.4)中，得到对当前层CBCT图像即图像A的形变矩阵P，再对所有待配准层的CBCT图像进行局部灰度校正配准得到所有层CBCT图像的形变矩阵

如图2所述，(a)为CT与CBCT图像(左图为CT，右图为CBCT)；(b)为本发明的形变配准输出图像；通过图2所述，本发明对控制点邻域灰度校正配准采用两级评价判断出暗区与伪影区，从而根据是否暗区以及是否伪影区有针对性的对两组图像进行局部灰度校正配准；另外，本发明基于参与配准的两组图像利用梯度变化特征并采用网格划分机制建立了局部控制点，基于控制点邻域进行局部灰度校正配准避免了全局盲配。

最后，应当理解的是，本发明中所述实施例仅用以说明本发明实施例的原则；其他的变形也可能属于本发明的范围；因此，作为示例而非限制，本发明实施例的替代配置可视为与本发明的教导一致；相应地，本发明的实施例不限于本发明明确介绍和描述的实施例。