一种图像配准方法和装置

技术领域

本申请涉及图像处理，尤其涉及一种图像配准方法和装置。

背景技术

近年来，一种新的超声容积导航技术得到发展，其原理为在超声查找病变的过程中，使实时或静态声像图与术前采集的CT或MR的H维数据进行图像融合，得到融合图像后，使用电磁场跟踪系统，协助医生在实时超声图像中更快地找到病变部位。这种技术是电磁导航技术与图像融合技术的结合，需要超声系统具备电磁跟踪装置及图像配准、融合算法，后者能够在经皮介入治疗或穿刺活检等治疗中实现多模式的图像配准、融合、视觉导航。图像配准、融合算法能够将患者术前的CT/MRI图像与术中的超声声像图进行实时配准、融合，送样就能够实时在声像图中显示针尖的位置和轨迹。然而目前采用的图像配准技术易得到错误窗口对，导致配准出错。

发明内容

本申请提供一种图像配准方法和装置。

根据本申请的第一方面，本申请提供一种图像配准方法，包括：

输入待配准图像和参考图像并分别进行小波变换，得到对应的待配准图像子图像和参考图像子图像；

对所述待配准图像子图像进行配准，根据所述参考图像子图像和配准后的所述待配准图像子图像计算所述待配准子图像的CC函数值；

根据所述CC函数值获取修正后的配准子图像。

进一步地，所述待配准图像子图像和参考图像子图像，分别包括低频分量、水平高频分量、垂直高频分量和对角高频分量；

所述对待配准图像子图像进行配准，包括：

分别对所述待配准图像子图像的所述低频分量、所述水平高频分量、所述垂直高频分量和所述对角高频分量中的一个或多个进行拉普拉斯金字塔修正配准。

进一步地，所述分别对所述待配准图像子图像的所述低频分量、所述水平高频分量、所述垂直高频分量和所述对角高频分量中的一个或多个进行拉普拉斯金字塔修正配准，包括：

分别对所述待配准图像子图像的所述低频分量、所述水平高频分量、所述垂直高频分量和所述对角高频分量的多层图像从底层图像开始，通过平移量和旋转角估计该层图像的配准参数，并计算每一分量的CC函数值；

通过该层子图像的配准参数调整下一层子图像分量直至最高层子图像；

记录最高层子图像的所述低频分量、所述水平高频分量、所述垂直高频分量和所述对角高频分量的CC函数值。

进一步地，所述根据CC函数值获取修正后的配准子图像，包括：

分别累加所述最高层子图像的所述低频分量、所述水平高频分量、所述垂直高频分量和所述对角高频分量的CC函数值；

输出累加和最大且超过预定阈值的图像作为配准图像。

根据本申请的第二方面，本申请提供一种图像配准装置，包括：

转换模块，用于输入待配准图像和参考图像并分别进行小波变换，得到对应的待配准图像子图像和参考图像子图像；

配准模块，用于对所述待配准图像子图像进行配准，根据所述参考图像子图像和配准后的所述待配准图像子图像计算所述待配准子图像的CC函数值；

输出模块，用于根据所述CC函数值获取修正后的配准子图像。

进一步地，所述待配准图像子图像和参考图像子图像，分别包括低频分量、水平高频分量、垂直高频分量和对角高频分量；

所述配准模块，还用于分别对所述待配准图像子图像的所述低频分量、所述水平高频分量、所述垂直高频分量和所述对角高频分量中的一个或多个进行拉普拉斯金字塔修正配准。

进一步地，所述配准模块包括：

处理单元，用于分别对所述待配准图像子图像的所述低频分量、所述水平高频分量、所述垂直高频分量和所述对角高频分量的多层图像从底层图像开始，通过平移量和旋转角估计该层图像的配准参数，并计算每一分量的CC函数值；

调整单元，用于通过该层子图像的配准参数调整下一层子图像分量直至最高层子图像；

记录单元，用于记录最高层子图像的所述低频分量、所述水平高频分量、所述垂直高频分量和所述对角高频分量的CC函数值。

进一步地，所述输出模块包括：

累加单元，用于分别累加所述最高层子图像的所述低频分量、所述水平高频分量、所述垂直高频分量和所述对角高频分量的CC函数值；

输出单元，用于输出累加和最大且超过预定阈值的图像作为配准图像。

根据本申请的第三方面，本申请提供一种图像配准装置，包括：

存储器，用于存储程序；

处理器，用于通过执行所述存储器存储的程序以实现上述方法。

根据本申请的第四方面，本申请提供一种计算机可读存储介质，包括程序，所述程序能够被处理器执行以实现上述方法。

本申请的有益效果是(或不写这句，对于要申请美国的申请，可省去有益效果部分)：

由于采用了以上技术方案，使本申请具备的有益效果在于：

在本申请的实施例中，由于包括对待配准图像和参考图像分别进行小波变换，得到对应的待配准图像子图像和参考图像子图像；对待配准图像子图像进行配准，根据参考图像子图像和配准后的待配准图像子图像计算待配准子图像的CC函数值；根据CC函数值获取修正后的配准图像。本申请对小波分解后多个尺度、多个角度的子图像求解CC函数，并通过其多级特性，对多个多尺度分量子图像进行配准，在一定程度上对误差、精度、方向性、对位移容错性进行了适量改进。

附图说明

图1为本申请的方法在一种实施方式中的流程图；

图2为本申请的方法在另一种实施方式中的流程图；

图3为本申请装置在一种实施方式中的程序模块示意图；

图4为本申请装置在另一种实施方式中的程序模块示意图。

具体实施方式

下面通过具体实施方式结合附图对本发明作进一步详细说明。本申请可以以多种不同的形式来实现，并不限于本实施例所描述的实施方式。提供以下具体实施方式的目的是便于对本申请公开内容更清楚透彻的理解，其中上、下、左、右等指示方位的字词仅是针对所示结构在对应附图中位置而言。

然而，本领域的技术人员可能会意识到其中的一个或多个的具体细节描述可以被省略，或者还可以采用其他的方法、组件或材料。在一些例子中，一些实施方式并没有描述或没有详细的描述。

本文中为部件所编序号本身，例如“第一”、“第二”等，仅用于区分所描述的对象，不具有任何顺序或技术含义。

此外，本文中记载的技术特征、技术方案还可以在一个或多个实施例中以任意合适的方式组合。对于本领域的技术人员来说，易于理解与本文提供的实施例有关的方法的步骤或操作顺序还可以改变。因此，附图和实施例中的任何顺序仅仅用于说明用途，并不暗示要求按照一定的顺序，除非明确说明要求按照某一顺序。

对小波分解后多个尺度、多个角度的子图像求解CC函数，并通过其多级特性，对多个多尺度分量子图像进行类拉普拉斯金字塔配准方式进行配准，在子图像最高层皆有一个CC函数，共有低频、水平高频、垂直高频、对角高频等4个CC函数值，且每个CC函数值带有不同的配准信息，若以单一为准则易引进误差，累加4个CC函数值，超过阈值且最大为其配准结果，在一定程度上对误差、精度、方向性、对位移容错性进行了适量改进。

实施例一：

如图1所示，本申请的图像配准方法，其一种实施方式，包括以下步骤：

步骤102：输入待配准图像和参考图像并分别进行小波变换，得到对应的待配准图像子图像和参考图像子图像。

步骤104：对待配准图像子图像进行配准，根据参考图像子图像和配准后的待配准图像子图像计算待配准子图像的CC函数值。

CC函数值具体可通过以下公式进行计算：

其中W为原始图像窗口对图像，E表示均值，I为待配准图像所选窗口，(i,j)为窗口内的组成图像的点的坐标。

步骤106：根据CC函数值获取修正后的配准图像。

如图2所示，本申请的图像配准方法，其另一种实施方式，包括以下步骤：

步骤202：分别输入待配准图像和参考图像。

步骤204：分别对待配准图像和参考图像进行小波分解，得到多个尺度、多方向和多层级的子图像。

对待配准图像进行小波分解后，得到待配准图像的子图像，待配准图像的子图像包括低频分量、水平高频分量、垂直高频分量和对角高频分量。对参考图像进行小波分解后，得到参考图像的子图像，参考图像的子图像也包括低频分量、水平高频分量、垂直高频分量和对角高频分量。

步骤206：分别获取待配准图像子图像的低频分量、水平高频分量、垂直高频分量和对角高频分量的底层子图像，此时，底层子图像为当前子图像。

步骤208：计算该层子图像的配准参数，配准参数包括计算平移量Δx和旋转角

步骤210：判断该层子图像是否为最高层子图像；若是转步骤214，否则转步骤212。

步骤212：使用该层子图像的配准参数调整下层子图像，配准层级加1，下层子图像转为当前子图像，转步骤208。

步骤214：获得参考图像的子图像各个分量的最高层子图像的CC函数值。

进一步地，本申请的方法可以分别对待配准图像子图像的低频分量、水平高频分量、垂直高频分量和对角高频分量中的一个或多个进行拉普拉斯金字塔修正配准。在本实施方式中，可以对待配准图像子图像的低频分量、水平高频分量、垂直高频分量和对角高频分量这四种分量进行拉普拉斯金字塔修正配准。具体可以先对四种分量中的一种进行配准，得到该分量的最高层子图像的CC函数值，再对另一种分量进行配准，得到其的最高层子图像的CC函数值，直到得到全部分量的最高层子图像的CC函数值。或者，可以同时对四种分量的底层子图像进行配准，然后再对四种分量的下一层子图像进行配准，直到得到全部分量的最高层子图像的CC函数值。

步骤216：累加最高层子图像的低频分量、水平高频分量、垂直高频分量和对角高频分量的CC函数值。

步骤218：输出累加和最大且超过预定阈值的图像作为配准图像。

实施例二：

如图3、图4所示，本申请的图像配准装置300，其一种实施方式，包括转换模块310、配准模块320和输出模块330。

转换模块310，用于输入待配准图像和参考图像并分别进行小波变换，得到对应的待配准图像子图像和参考图像子图像；

配准模块320，用于对待配准图像子图像进行配准，根据参考图像子图像和配准后的待配准图像子图像计算待配准子图像的CC函数值；

输出模块330，用于根据CC函数值获取修正后的配准图像。

进一步地，待配准图像子图像和参考图像子图像，分别包括低频分量、水平高频分量、垂直高频分量和对角高频分量；

配准模320块，还用于分别对待配准图像子图像的低频分量、水平高频分量、垂直高频分量和对角高频分量中的一个或多个进行拉普拉斯金字塔修正配准。

进一步地，配准模块320可以包括处理单元321、调整单元322和记录单元323。

处理单元321，用于分别对待配准图像子图像的低频分量、水平高频分量、垂直高频分量和对角高频分量的多层子图像从底层子图像开始，通过平移量和旋转角估计该层图像的配准参数，并计算每一分量的CC函数值；

调整单元322，用于通过该层子图像的配准参数调整下一层子图像分量直至最高层子图像；

记录单元323，用于记录最高层子图像的低频分量、水平高频分量、垂直高频分量和对角高频分量的CC函数值。

进一步地，输出模块330可以包括累加单元331和输出单元332。

累加单元，用于分别累加最高层子图像的低频分量、水平高频分量、垂直高频分量和对角高频分量的CC函数值；

输出单元，用于输出累加和最大且超过预定阈值的图像作为配准图像。

实施例三：

本申请的图像配准装置，其一种实施方式，包括存储器和处理器。

存储器，用于存储程序；

处理器，用于通过执行所述存储器存储的程序以实现实施例一中的方法。

实施例四：

一种计算机可读存储介质，包括程序，程序能够被处理器执行以实现实施例一中的方法。

本领域技术人员可以理解，上述实施方式中各种方法的全部或部分步骤可以通过程序来指令相关硬件完成，该程序可以存储于一计算机可读存储介质中，存储介质可以包括：只读存储器、随机存储器、磁盘或光盘等。

以上内容是结合具体的实施方式对本申请所作的进一步详细说明，不能认定本申请的具体实施只局限于这些说明。对于本申请所属技术领域的普通技术人员来说，在不脱离本申请构思的前提下，还可以做出若干简单推演或替换。