[发明专利]一种SVD-RANSAC亚像素相位相关匹配方法

说明书

技术领域

本发明涉及一种图像区域匹配算法，尤其是涉及一种SVD-RANSAC亚像素相位相关匹配方法。

背景技术

影像的亚像素精确匹配是摄影测量与遥感领域的研究热点与难题之一。通过影像间的亚像素匹配，能够精确的获取立体视差、地物位移、地表形变场等数据，基于亚像素匹配的图像配准，在DEM/DSM生成、图像镶嵌、影像融合、信息提取和变形监测等领域有着极其重要的应用。

图像匹配一般分为特征匹配与区域匹配两大类。相对于特征匹配，区域匹配具有精度高，更方便的粗差剔除，分布均匀等优势。相位相关是一种频域下的区域匹配算法，相位相关的理论基础是傅里叶变换的平移特性，即影像间平移在频域下对应为线性相位差，理论上，影像的平移只会引起傅里叶系数相位角的变化而不会改变其幅值。与传统的灰度相关系数相比，具有匹配精度高、速度快、抗噪性强和受辐射差异影响小等特点，近年来得到了广泛的应用与关注。

亚像素的相位相关方法研究主要分为两类，第一类方法通过确定互功率谱的逆傅里叶变换后精确的峰值位置来获取偏移量，第二类方法通过直接估计互功率谱的线性相位差来获取偏移量。然而，现有的相位相关亚像素匹配方法依然存在问题：通过插值确定峰值的方法精度相对较低；利用线性相位差的方法，普遍具有计算复杂度高，抗粗差性差等缺点；另外，亚像素匹配精度受像素锁现象(pixel locking effect)影响。上述的问题直接影响到基于相位相关的亚像素匹配精度和匹配结果的稳健性，并对运算效率提出了较高的要求。

文献“A subspace identification extension to the phase correlation method[MRI application][J].Medical Imaging，IEEE Transactions on，2003，22(2)：277-280，利用最小二乘估计来计算相位角向量直线的斜率，线性拟合的质量依赖于解缠后相位角的线性度。即使奇异值分解获取主奇异值向量本身是一个提高SNR的过程，但相关过程中的混叠，噪声和其他误差仍然会影响两幅影像间的相位信息，从而影响归一化互功率谱矩阵的部分量值，导致奇异值分解后的相位角向量出现偏差，并不呈现严格的直线关系，由于最小二乘的抗差性弱，估计的偏移值结果会受到偏差的严重影响，从而导致亚像素级估计结果的精度和稳定性下降。

发明内容

本发明的目的就是为了克服上述现有技术存在的缺陷而提供一种结果精度和稳定性高、有效抑制pixel locking现象的SVD-RANSAC亚像素相位相关匹配方法。

本发明的目的可以通过以下技术方案来实现：

一种SVD-RANSAC亚像素相位相关匹配方法，该方法包括以下步骤：

1)获取存在水平方向偏移a和垂直方向偏移b的两幅影像g(x，y)和h(x，y)，对两幅影像分别进行Hanning窗函数加权后进行离散傅里叶变换，得到G(u，v)和H(u，v)，则H(u，v)＝G(u，v)exp{-i(au+bv)}，计算两幅影像的归一化互功率谱矩阵Q(u，v)：