一种基于融合图像引导的多视差图像序列配准方法技术

本发明专利技术属于图像配准技术领域，具体涉及一种基于融合图像引导的多视差图像序列配准方法。技术方案包括以下步骤：首先，利用非降采样小波变换对待配准的多视差图像序列的高低频子图进行融合；然后选择融合图像中分布于四角的四个标准特征子图，依次在待配准的多视差图像序列中对应位置设置搜索范围，通过选择搜索范围内相关系数最大值位置得到特征坐标点对；根据特征坐标点对计算每幅图像的拉伸系数；最后将待配准的多视差图像序列中所有位置与对应的拉伸系数相乘得到配准后的多视差图像序列。本发明专利技术的方法能够有效克服基准图像中纹理稀疏的问题，可实现多视差图像序列的精确配准。准。

【技术实现步骤摘要】
一种基于融合图像引导的多视差图像序列配准方法


[0001]本专利技术属于图像配准
，具体涉及一种基于融合图像引导的多视差图像序列配准方法。

技术介绍

[0002]传感器技术的迅猛发展，使得图像采集设备的终端尺寸不断趋于微型化与小型化，以手机、智能终端摄像头为代表的图像采集设备在分辨率和图像处理性能方面出现了质的提升。这使得以前需要基于大型图像运算设备计算的三维重建问题得以通过小型计算终端完成，现阶段通常运用小型终端采集多幅视差图像序列实现相关应用场景的三维重建。其中，图像配准作为多视差图像序列三维重建问题的前期处理步骤，其配准精度直接决定重建结果的好坏。然而，多视差图像序列在采集过程中由于受到相机抖动等因素会导致图像序列出现移位，进而导致三维重建结果出现较大偏差。因此，准确高效的多视差图像序列的配准尤为重要。
[0003]现有的多视差图像配准的方法主要分为图像特征点变换和多平面分割两类。其中，图像特征点变换的方法通常对多幅视差图像利用统一的特征描述算子(如SIFT、LBP)进行匹配，但这类方法仅适用于视差范围较小的情况，对于距离目标物体较远的视差图像，由于这类图像普遍具有纹理稀疏与图像平滑的特点，这可能导致特征点选择出现困难；多平面分割的方法主要用于克服视差范围较大的情况，通过将视差图像分割为多个区域，视差图像序列间的各区域分别采用不同的单应性矩阵建立关联，但此类方法对分割方法的准确性有较高要求，除此之外，相邻视差图像之间的特征区域可能会随着相机与物体间的距离发生模糊变化，导致特征点匹配的结果出现错误。
[0004]综上可知，图像特征点变换和多平面分割均无法直接用于三维重建中多视差图像序列的配准问题，而且这两类匹配方法通常采用首张视差图像作为基准图像，然后对其他视差图像序列进行匹配，假使首张图像中特征稀疏，将导致后续图像序列的匹配出现较大的积累误差。因此，我们认为在多视差图像序列配准过程中，基准图像是保证所有图像序列配准的关键。本专利采用融合图像作为基准图像对所有的多视差图像序列进行配准，确保了多视差图像序列可以实现精确匹配。

技术实现思路

[0005]为克服现有解决方案中不足之处，本专利技术的目的是提供一种基于融合图像引导的多视差图像序列配准方法。
[0006]本专利技术所采用的技术方案是，一种基于融合图像引导的多视差图像序列配准方法，包括以下步骤：
[0007]步骤1，设置非降采样小波变换的最大分解层数为L
max
，滤波器组为
‘
maxflat
’
；
[0008]步骤2，对待配准的多视差图像序列通过式(1)进行非降采样小波变换得到n
幅低频子图与n
×
L
max
幅高频子图
[0009][0010]其中，I
i
为多视差图像序列中第i幅图像，多视差图像序列的总数为n，l表示分解尺度且其取值为1≤l≤L
max
，NSWT表示非降采样小波变换；
[0011]步骤3，对步骤2得到的n幅低频子图按照式(2)得到对应的低频融合子图
[0012][0013]其中max(
·
)表示取最大值函数；
[0014]步骤4，对步骤2得到的n
×
L
max
幅高频子图通过式(3)得到对应的高频融合子图集
[0015][0016]其中abs(
·
)表示取绝对值函数；
[0017]步骤5，对步骤3和步骤4得到的低频融合子图与高频融合子图集通过式(4)的非降采样小波变换的反变换得到多视差图像序列的融合图像I
F
，
[0018][0019]其中，Inverse NSWT表示非降采样小波变换反变换；
[0020]步骤6，选择融合图像I
F
中处于四角位置的四个矩形区域作为标准特征子图，并在对应的待配准的多视差图像序列中相同位置设定标准特征子图的搜索范围其中搜索范围与标准特征子图的关系满足式(5)；
[0021][0022]步骤7，设四个标准特征子图的中心坐标为然后将标准特征子图依次在待配准的多视差图像序列对应的搜索范围内遍历，并根据式(6)计算出搜索范围内所有的相关系数，选择相关系数最大值作为配准点
[0023][0024]最终，待配准的多视差图像序列中各图像均得到4个特征坐标点对
其中G(
·
)表示相关系数；
[0025]步骤8，将步骤7得到的待配准的多视差图像序列中各图像的4个特征坐标点对根据式(7)依次计算其对应的拉伸系数
[0026][0027]步骤9，待配准的多视差图像序列按照步骤8中得到的拉伸系数按照式(8)得到配准后的多视差图像序列
[0028][0029]与现有技术相比，本专利技术具有以下优点：
[0030](1)本专利技术采用多视差图像序列的融合图像作为引导配准的基准图像，可有效的克服基准图像中纹理稀疏和特征点难以选择的问题；
[0031](2)本专利技术提出的配准方法具有较高的运算效率，可实现多视差图像序列的快速精确匹配，保证了后续三维重建结果的准确性。
附图说明
[0032]图1是本专利技术一种基于融合图像引导的多视差图像序列配准方法的示意图。
具体实施方式
[0033]如图1所示，本实施例所述一种基于融合图像引导的多视差图像序列配准方法，包括以下步骤：
[0034]步骤1，设置非降采样小波变换的最大分解层数为L
max
，滤波器组为
‘
maxflat
’
；
[0035]步骤2，对待配准的多视差图像序列通过式(1)进行非降采样小波变换得到n幅低频子图与n
×
L
max
幅高频子图
[0036][0037]其中，I
i
为多视差图像序列中第i幅图像，多视差图像序列的总数为n，l表示分解尺度且其取值为1≤l≤L
max
，NSWT表示非降采样小波变换；
[0038]步骤3，对步骤2得到的n幅低频子图按照式(2)得到对应的低频融合子图
[0039][0040]其中max(
·
)表示取最大值函数；
[0041]步骤4，对步骤2得到的n
×
L
max
幅高频子图通过式(3)得到对应的高频融合子图集
[0042][0043]其中abs(
·
)表示取绝对值函数；
[0044]步骤5，对步骤3和步骤4得到的低频融合子图与高频融合子图集通过式(4)的非降采样小波变换的反变换得到多视差图像序列的融合图像I
F
，
[0045][0046]其中，Inverse NSWT表示非降采样小波变换反变换；
[0047]步骤6，选择融合图像I
F
中处于四本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于融合图像引导的多视差图像序列配准方法，其特征是，包括以下步骤：步骤1，设置非降采样小波变换的最大分解层数为L
max
，滤波器组为
‘
maxflat
’
；步骤2，对待配准的多视差图像序列通过式(1)进行非降采样小波变换得到n幅低频子图与n
×
L
max
幅高频子图幅高频子图其中，I
i
为多视差图像序列中第i幅图像，多视差图像序列的总数为n，l表示分解尺度且其取值为1≤l≤L
max
，NSWT表示非降采样小波变换；步骤3，对步骤2得到的n幅低频子图按照式(2)得到对应的低频融合子图按照式(2)得到对应的低频融合子图其中max(
·
)表示取最大值函数；步骤4，对步骤2得到的n
×
L
max
幅高频子图通过式(3)得到对应的高频融合子图集图集其中abs(
·
)表示取绝对值函数；步骤5，...

【专利技术属性】
技术研发人员：闫涛，钱宇华，杜亮，
申请(专利权)人：山西大学，
类型：发明
国别省市：