本发明是一种基于权重相位相关峰值拟合的高精度图像配准方法，尤其涉及一种权重的相位相关峰值拟合方法，属于计算机视觉、遥感图像处理领域。
技术背景
基于相位相关方法的配准参数估计是计算机视觉及遥感图像处理领域的热门研究方向之一，在三维目标识别、大地测量、机器人视觉导航等方向有着广泛的应用。基于相位相关的图像配准参数估计是利用傅里叶变换平移性质对待配准图像进行傅里叶变换，构建归一化互功率谱，再对归一化互功率谱进行傅里叶逆变换，通过高效、高精度的峰值拟合方法来求解配准参数。
由于相位相关方法对光照强度不敏感、配准速度快，且具有亚像素配准精度等特性，使得它在图像处理领域得到来广泛应用。在基于相位相关方法的图像配准中，峰值拟合是相位相关方法的关键步骤之一，受旁瓣能量、系统累积误差的影响，传统拟合方法，常常不能高效、稳定地求解亚像素级图像配准参数。
技术实现要素：
本发明主要解决的技术问题是提供一种基于权重的相位相关峰值拟合方法，以解决图像配准中，传统基于相位相关的图像配准方法受旁瓣能量以及系统累积误差的影响，导致其不能稳定、高效地对图像配准参数进行高精度估计的情况。
本发明的解决方案是：1)构建一个反比例函数，并证明其满足kroneckerdelta函数的平移性质；2)构建一种基于权重反比例函数的相位相关拟合方法，基于三个拟合点实现图像配准参数的高效、高精度估计。
本发明为实现上述解决方案，步骤如下：1、图像预处理，求解归一化互功率谱及其傅里叶反变换；2、构造一个反比例函数，并证明它满足kroneckerdelta函数的平移性质；3、构建一种基于权重反比例函数的相位相关拟合方法，基于三个拟合点实现图像配准参数的高效、高精度估计，并根配准参数融合待配准图像；
1.图像预处理，求解归一化互功率谱及其傅里叶反变换：
1)通过高斯滤波函数，对待配准图像进行平滑处理；
2)待配准图像傅里叶变换求解：
假设存在两幅图像fi(x)，x＝[x,y]t∈r2,i＝1,2，且图像f1(x)、f2(x)之间存在一个平移量d＝[δx,δy]t∈r2，那么f1(x)、f2(x)存在如下关系：f1(x)＝f2(x+d)。根据傅里叶变换理论，若f1(x)、f2(x)之间存在平移关系，则在傅里叶频域空间下，他们的傅里叶变换关系可表示为
3)待配准图像归一化互功率谱求解：
为了获取相位差，我们计算f1(x)、f2(x)经傅里叶变换后的归一化互功率谱c(u)：
其中，u表示频率域坐标，j表示复数，w表示图像尺寸，t表示转置。
4)对归一化互功率谱c(u)进行傅里叶反变换求解：
其中，
2.构造一个反比例函数，并证明它满足kroneckerdelta函数的平移性质：
1)将归一化互功率谱c(u)的傅里叶反变换可以被看作kroneckerdelta函数
i(x)≈δ(x+d)，x＝[x,y]t，d＝[δx,δy]t
2)根据kroneckerdelta函数的多维性质，可将δ(x+d)表示为
δ(x+d)＝δ(x+δx)δ(y+δy)，x＝[x,y]t，d＝[δx,δy]t
3)根据kroneckerdelta函数的平移性质(以x方向为例)
4)构建反比例函数
利用极限定理证明所构造的反比例函数满足kroneckerdelta函数的定义，在a→0条件下，根据洛必达法则得：
5)当a→0时，拟通过积分公式证明所构造的反比例函数满足kroneckerdelta函数的平移性质，积分公式定义为：
3.提出一种基于权重反比例函数的相位相关峰值拟合方法，实现基于三个拟合点的图像配准参数高效、高精度估计，并根配准参数融合待配准图像：
1)根据步骤2的证明，构建的比例函数满足kroneckerdelta函数的平移性质，可将其近似为kroneckerdelta函数：
并将反比例拟合函数定义为：
并将其用于i(x)的亚像素级峰值位置估计。
2)根据峰值点(x,y)＝argmaxx,yi(x)，在x方向上选择包括峰值点在内的三个点，作为拟合点，
其中，(x1，p(x1))表示在x方向上峰值点的左边领域点，(x1，p(x1))表示在x方向上的峰值点，(x3，p(x3))表示在x方向上峰值点的右边领域点。
通过p(x1)除以p(x2)、p(x2)除以p(x3)，得
为了简化复杂公式的运算过程，设a＝0，得
3)考虑
第一，如果峰值位置位于待配准图像的左边缘，表示x1不存在，则使用位置x2、x3及其峰值p(x2)，p(x2)来拟合亚像素级峰值位置，定义如下：
选择x2和x3(
第二，当峰值位置位于待配准图像的右边缘，将导致x3消失，那么位置x1、x2将被用于求解峰值点的位置，定义如下：
选择[x1,x2]中包含的解为最佳峰值位置；
第三，如果峰值位置不位于待配准图像的边缘，那么峰值位置将会出现在x2的两侧，这里我们构建一个基于权重的反比例函数拟合算法，该算法通过一个合理的加权函数提供一个更稳定、更高效的峰值位置求解，定义如下：
其中，ω1和ω2表示加权函数，定义为
其中，k表示经验参数，p(x1)是x1坐标的峰值，p(x3)是x3坐标的峰值。在y方向上，可根据权利要求书中权利要求3和4的步骤对y方向上的平移参数进行求解。
同理，针对待配准图像中出现的旋转、放缩变化，在傅里叶梅林变换的辅助下，将旋转、放缩参数转换成符合权利要求书中权利要求3和4中的平移参数求解形式，并通过权利要求书中权利要求3和4的步骤对旋转、放缩参数进行高效、高精度求解。最后，根据基于权重相位相关峰值拟合方法求得的高精度图像配准参数，通过线性变换矩阵来实现待配准图像的融合。