一种图像配准方法和装置制造方法及图纸

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

【技术实现步骤摘要】
一种图像配准方法和装置
本申请涉及图像处理，尤其涉及一种图像配准方法和装置。
技术介绍
近年来，一种新的超声容积导航技术得到发展，其原理为在超声查找病变的过程中，使实时或静态声像图与术前采集的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所示，本申请的图像配准方法，其一种实施方式，包括以下步骤：步本文档来自技高网...

【技术保护点】
1.一种图像配准方法，其特征在于，包括：/n输入待配准图像和参考图像并分别进行小波变换，得到对应的待配准图像子图像和参考图像子图像；/n对所述待配准图像子图像进行配准，根据所述参考图像子图像和配准后的所述待配准图像子图像计算所述待配准子图像的CC函数值；/n根据所述CC函数值获取修正后的配准图像。/n

【技术特征摘要】
1.一种图像配准方法，其特征在于，包括：
输入待配准图像和参考图像并分别进行小波变换，得到对应的待配准图像子图像和参考图像子图像；
对所述待配准图像子图像进行配准，根据所述参考图像子图像和配准后的所述待配准图像子图像计算所述待配准子图像的CC函数值；
根据所述CC函数值获取修正后的配准图像。


2.如权利要求1所述的方法，其特征在于，所述待配准图像子图像和参考图像子图像，分别包括低频分量、水平高频分量、垂直高频分量和对角高频分量；
所述对待配准图像子图像进行配准，包括：
分别对所述待配准图像子图像的所述低频分量、所述水平高频分量、所述垂直高频分量和所述对角高频分量中的一个或多个进行拉普拉斯金字塔修正配准。


3.如权利要求2所述的方法，其特征在于，所述分别对所述待配准图像子图像的所述低频分量、所述水平高频分量、所述垂直高频分量和所述对角高频分量中的一个或多个进行拉普拉斯金字塔修正配准，包括：
分别对所述待配准图像子图像的所述低频分量、所述水平高频分量、所述垂直高频分量和所述对角高频分量的多层图像从底层图像开始，通过平移量和旋转角估计该层图像的配准参数，并计算每一分量的CC函数值；
通过该层子图像的配准参数调整下一层子图像分量直至最高层子图像；
记录最高层子图像的所述低频分量、所述水平高频分量、所述垂直高频分量和所述对角高频分量的CC函数值。


4.如权利要求3所述的方法，其特征在于，所述根据CC函数值获取修正后的配准图像，包括：
分别累加所述最高层子图像的所述低频分量、所述水平高频分量、所述垂直高频分量和所述对角高频分量的CC函数值；
输出累加和最大且超过预定阈值的图像作为配准图像。


5.一种图像配准装置，其特征在于，包括：
转换模块，用于输入待配准图像和参考图像并分别进行小波变换，得到对应的...

【专利技术属性】
技术研发人员：刘林泉，杨加成，张强，席奎，
申请(专利权)人：深圳市恩普电子技术有限公司，
类型：发明
国别省市：广东;44