本发明涉及图像处理技术领域，具体为一种图像透视变换的优化方法。
背景技术：
图像透视变换是利用透视中心、像点、目标点三点共线的条件将图片投影到一个新的视平面，也称作投影映射。透视变换的变换公式为x＝h·x'，其中x是像点，即参考图像的坐标，x是目标点，即目标图像的坐标，h是变换矩阵。通过参考图像和目标图像的部分坐标得到h，根据h对参考图像进行透视变换，就可以根据参考图像得到目标图像的拟合，此方法可以应用于个人相册压缩和云端图像存储等领域中。由参考图像和目标图像的部分坐标得到h的过程中，为使得h更加精确，需要使用最优化方法对h进行修正。目前，比较常用的最优化方法有：梯度下降法、牛顿法和levenberg-marquardt(简称为l-m)法等。梯度下降法也称为最速下降法，用负梯度方向为搜索方向，越接近目标值，步长越小，前进越慢，存在靠近极小值时收敛速度减慢和直线搜索时得不到理想结果等一些问题。牛顿法收敛速度快，但对初始点要求严格，方向构造困难，计算复杂且占用内存较大。l-m方法是利用梯度求最大或最小值的方法，形象的说，属于“爬山”法的一种。l-m方法同时具有梯度法和牛顿法的优点。在l-m算法中，每次迭代会寻找一个合适的阻尼因子λ。当λ很小时，步长等于牛顿法步长，当λ很大时，步长约等于梯度下降法的步长。l-m方法对过参数化问题不敏感，能有效处理冗余参数问题，使代价函数陷入局部极小值的机会大大减小，这些特性使得l-m方法在计算机视觉等领域得到广泛应用。所以，在对h进行修正时都采用l-m方法。但是在透视变换中采用最初的l-m方法存在计算复杂度高，运行时间长的问题。
所以，如何设计一种图像透视变换的优化方法，成为我们当前要解决的问题。
技术实现要素：
本发明的目的在于提供一种图像透视变换的优化方法，以解决上述背景技术中提出的问题。
为实现上述目的，本发明提供如下技术方案：一种图像透视变换的优化方法，包括如下步骤：
1)从粗匹配的图像中随机抽取4对匹配点并判断是否有3点共线的情况，如果有，则舍去这一组匹配点重新随机抽取，直到4对匹配点中没有3点共线的情况；
2)对选取的4对匹配点进行数据归一化，根据h元素的三个方程计算出单应矩阵h；
3)利用单应矩阵h，采用对称转移误差方法，来估算一组对应的误差，如果d2transfer大于设定的阈值，则作为外点舍去，否者认定其对应的匹配点为内点；
4)在估计次数n内重复采样与验证，记录内点数目，选取具有较多内点数目的数据集进行下一步操作；
5)设定估计次数n的阈值为1000，选取内点数目与角点总数的比值ω(ω≥0.5)作为数据是真实数字模型内点的概率，那么一次估计中所有4个点都是内点的概率为ω4，n次采样后，记录内点数目最大时的一致集；
6)利用数目最大时数据集的所有内点从新估计单应矩阵hr，并用hr实现图像的配准，即ir＝hri,其中i代表配准前的图像，ir代表配准后的图像。
根据上述技术方案，所述步骤3中对称转移误差方法的计算公式为：
d2transfer＝d(pi,h-1pi’)2+d(pi’,hpi)2
采用glt算法估算变换矩阵，该方法可以得到更准确的参量，由4对匹配点对应pi→p’i，确定变换方程为p’i＝hpi。
根据上述技术方案，所述步骤5)中如果经过n次至少有一次估计中的所有数据点都是内点的概率是p(p＝0.99)，则n需要满足：
根据上述技术方案，所述步骤2)中h元素的方程为：
这些方程都有aih＝0的形式，其中ai是3×9的矩阵，h是由矩阵h的元素组成的9维矢量，h＝[h1h2h3]t,hi是h的i个元素，每组点对应给出关于h元素的两个独立方程，利用4组点或得方程组ah＝0，h是h未知元素的矢量，对a进行svd分解，v的值就是h的解；应用上述公式与pi→p'i，求得单应矩阵h'；令h＝t-1h't，求得单应矩阵h。
根据上述技术方案，所述步骤1)中匹配点抽取采用harris算法，即在两幅图像中提取的角点数目不一定相等，并且有相当多的冗余点，这时就需要去除可能导致匹配出现较大误差或导致匹配失败的冗余点，得出一一对应的精确匹配点；首先利用ncc进行特征点的粗匹配，然后利用ransac剔除错误的匹配点，实现图像的精确配准。
根据上述技术方案，所述ncc算法通过计算特征点邻域像素灰度的互相关系数，并以该系数为匹配准则进行匹配，其原理为：分别在参考图像和待配准图像中以每一个角点为中心，取一个(2n+1)×(2n+1)大小的相关窗，然后以参考帧中的每一个特征点为参考点在持配准图像中搜索对应的匹配点，互相关最大时的搜索窗口位置决定了模板图像在待匹配图像中的位置。
与现有技术相比，本发明的有益效果是：本发明针对图像末制异过程中图像配准的精确性和实时性的具休要求，提出了一种透视变换模型下基于harris算法的图像配准方法，该方法增强了传统harria算法的旋转不变性，能够更好地适用于发生大角度旋转的图像，结合ncc算法和ransac算法并引人dlt算法计算稳定的透视变换单应知矩阵，利用该单应矩阵实现图像配准，能有效提升该图像透视的变换的效果，有效提升使用的方便性。
附图说明
图1是本发明的图像透视变换的优化方法方框图；
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
请参阅图1，本发明提供一种图像透视变换的优化方法，包括如下步骤：
1)从粗匹配的图像中随机抽取4对匹配点并判断是否有3点共线的情况，如果有，则舍去这一组匹配点重新随机抽取，直到4对匹配点中没有3点共线的情况；
2)对选取的4对匹配点进行数据归一化，根据h元素的三个方程计算出单应矩阵h；
3)利用单应矩阵h，采用对称转移误差方法，来估算一组对应的误差，如果d2transfer大于设定的阈值，则作为外点舍去，否者认定其对应的匹配点为内点；
4)在估计次数n内重复采样与验证，记录内点数目，选取具有较多内点数目的数据集进行下一步操作；
5)设定估计次数n的阈值为1000，选取内点数目与角点总数的比值ω(ω≥0.5)作为数据是真实数字模型内点的概率，那么一次估计中所有4个点都是内点的概率为ω4，n次采样后，记录内点数目最大时的一致集；
6)利用数目最大时数据集的所有内点从新估计单应矩阵hr，并用hr实现图像的配准，即ir＝hri,其中i代表配准前的图像，ir代表配准后的图像。
根据上述技术方案，步骤3中对称转移误差方法的计算公式为：
d2transfer＝d(pi,h-1pi')2+d(pi',hpi)2
采用glt算法估算变换矩阵，该方法可以得到更准确的参量，由4对匹配点对应pi→p'i，确定变换方程为p'i＝hpi。
根据上述技术方案，步骤5)中如果经过n次至少有一次估计中的所有数据点都是内点的概率是p(p＝0.99)，则n需要满足：
根据上述技术方案，步骤2)中h元素的方程为：
这些方程都有aih＝0的形式，其中ai是3×9的矩阵，h是由矩阵h的元素组成的9维矢量，h＝[h1h2h3]t,hi是h的i个元素，每组点对应给出关于h元素的两个独立方程，利用4组点或得方程组ah＝0，h是h未知元素的矢量，对a进行svd分解，v的值就是h的解；应用上述公式与pi→p'i，求得单应矩阵h'；令h＝t-1h't，求得单应矩阵h。
根据上述技术方案，步骤1)中匹配点抽取采用harris算法，即在两幅图像中提取的角点数目不一定相等，并且有相当多的冗余点，这时就需要去除可能导致匹配出现较大误差或导致匹配失败的冗余点，得出一一对应的精确匹配点；首先利用ncc进行特征点的粗匹配，然后利用ransac剔除错误的匹配点，实现图像的精确配准。
根据上述技术方案，ncc算法通过计算特征点邻域像素灰度的互相关系数，并以该系数为匹配准则进行匹配，其原理为：分别在参考图像和待配准图像中以每一个角点为中心，取一个(2n+1)×(2n+1)大小的相关窗，然后以参考帧中的每一个特征点为参考点在持配准图像中搜索对应的匹配点，互相关最大时的搜索窗口位置决定了模板图像在待匹配图像中的位置。
基于上述，本发明的优点在于，本发明针对图像末制异过程中图像配准的精确性和实时性的具休要求，提出了一种透视变换模型下基于harris算法的图像配准方法，该方法增强了传统harria算法的旋转不变性，能够更好地适用于发生大角度旋转的图像，结合ncc算法和ransac算法并引人dlt算法计算稳定的透视变换单应知矩阵，利用该单应矩阵实现图像配准，能有效提升该图像透视的变换的效果，有效提升使用的方便性。
尽管已经示出和描述了本发明的实施例，对于本领域的普通技术人员而言，可以理解在不脱离本发明的原理和精神的情况下可以对这些实施例进行多种变化、修改、替换和变型，本发明的范围由所附权利要求及其等同物限定。