本发明涉及摄影测量与遥感测绘技术领域，尤其涉及一种基于奇异值分解的影像匹配提纯方法。
背景技术：
在影像匹配过程中，经常会存在一些误匹配点对，这时可以通过匹配提纯的方法对误匹配点对进行剔除。一般的提纯思路是寻找一个恰当的匹配点对约束模型，估计出正确的模型参数，利用约束模型进行匹配点对提纯。在匹配点对约束模型中，通常使用变换约束矩阵作为估计模型，其一般包括平移变换、刚体变换、相似性变换、仿射变换、影射变换、极线几何基本矩阵变换等。现行的变换约束矩阵模型已经能够很好的进行匹配点对间的几何约束，因此大多数学者更加倾向于模型参数估计方法的研究。
在模型参数估计中，常用的方法有稳健回归估计和随机参数估计等。对于稳健回归估计方法，如m‐估计，其核心思想是采用迭代加权最小二乘估计回归系数，但其仅能适应误匹配率较小的情况。对于随机参数估计方法，如lmeds(leastmedianofsquares)算法、mlesac(maximumlikelihoodestimationsampleandconsensus)算法和ransac(randomsampleconsensus)算法等，其核心思想是随机选择样本子集，迭代挑选出最佳模型参数。其中，ransac算法能够在存在大量外点的数据中找到内点，因而得到广泛应用并派生出一系列改进算法。改进的ransac方法主要从两个方面进行了优化。一方面是速度优化，如p-ransac(preemptiveransac)和r-ransac(randomizedransac)等；一方面是提纯率的优化，如m-ransac(multiransac)等。然而这些算法都需要选择一个合适的初始内点集，提纯结果容易陷入局部最优解，造成部分正确匹配点对遗漏；在误匹配率增大时，成功选择内点集的试验次数也急剧增大。
技术实现要素：
本发明的目的在于提供一种基于奇异值分解的影像匹配提纯方法，能够避免初始内点集选择不当造成的匹配点对提纯不准，保证匹配点对提纯的准确率。
为实现上述目的，本发明采用的技术方案是：
基于奇异值分解的影像匹配提纯方法，包括以下步骤：
a、记c为迭代索引，c＝1，2，3，……，pnum，其中pnum为迭代次数上限；
b、假设点集φ＝{ρk|k＝0，1，2，……，n}中包含n对匹配点对，对于任意一对匹配点对，分别采用基本矩阵建立误差方程，即
m′kfmk＝0；
其中，mk＝(xk,yk,1)t，mk′＝(x′k,yk′,1)t分别为匹配点对齐次坐标，f为基本矩阵，令f＝(fij)，则上式可写为如下形式
(xk′xk,x′kyk,x′k,y′kxk,y′kyk,y′k,xk,yk,1)f＝0；
其中，f＝(f11,f12,f13,f21,f22,f23,f31,f32,f33)t，利用各对匹配点对的误差方程建立约束模型：
其中，a为n×9维的系数矩阵；
c、对系数矩阵a进行奇异值分解，利用a的前t个奇异值及其对应的左、右奇异向量构建a的近似矩阵a′，计算差分矩阵δa＝a-a′，利用δa计算每对匹配点对的近似误差δk以及近似均方误差mse，以近似均方误差mse作为阈值，对每对匹配点对的误差方程进行筛选，得到约束后的匹配点对误差方程；
d、利用约束后的匹配点对误差方程再次组成新的系数矩阵，得到降维矩阵b，利用降维矩阵b建立约束模型bf＝0，对降维矩阵b进行奇异值分解，求解f的最小二乘解，并进行秩2约束，得到模型参数解f；
e、利用基本矩阵及求解的模型参数f，采用误差方程m′kfmk＝0对各个匹配点对进行约束，若匹配点对的mk′fmk值大于经验值α，则将该点标记为外点，否则标记为内点；
f、判断c值与pnum的大小，当c＜pnum时，判断模型参数f是否达到预定精度，若是则将标记的内点作为提纯点对，完成提纯过程，若否则将c值加1，并将标记的内点作为点集φ中包含的匹配点对，然后返回步骤b，当c＝pnum时，将标记的内点作为提纯点对，完成提纯过程。
所述的步骤c中，假设系数矩阵a的秩为r，则对系数矩阵a进行奇异值分解有
a＝usvt；
其中，u为n×n维正交矩阵，v为9×9维正交矩阵，s为n×9维对角阵，u的列向量u1，u2，……，un称为系数矩阵a的左奇异向量，v的列向量v1，v2，……，v9称为系数矩阵a的右奇异向量，s的对角线元素σ1，σ2，……，σr称为系数矩阵a的奇异值，且以递减顺序排列，即σ1≥σ2≥……≥σr；
由系数矩阵a的前t个奇异值重新组成对角阵st×t＝diag[σ1,σ2,……,σt]，由奇异值对应的左、右奇异向量重新组成矩阵un×t＝[u1,u2,……,ut]和v9×t＝[v1,v2,……,vt]，则利用可重构的系数矩阵a，得到
a′＝um×tst×tvt×nt；
对上式进行化简可得
其中，a′为a的近似矩阵，由系数矩阵a的主奇异值重构而成，包含了系数矩阵a的主要信息；
计算a与a′间的差分矩阵δa为
计算各对匹配点对的近似误差为
其中，δaij为δa中的元素。
所述的步骤c中，以近似均方误差mse作为阈值，对每对匹配点对的误差方程进行筛选时，若δk<mse则保留第k对匹配点对的误差方程，若δk≥mse则舍弃第k对匹配点对对应的误差方程，由此得到约束后的匹配点对误差方程。
所述的步骤f中，判断模型参数f是否达到预定精度时，计算当前f值与上次迭代所得f值的差值，若差值小于设定值β，则判定f达到预定精度，若差值大于等于设定值β，则判定f没有达到预定精度。
本发明采用特定模型描述影像匹配点对之间的变换关系，建立模型误差方程，并对模型含有匹配点对信息的系数矩阵进行奇异值分解，重构大奇异值对应的矩阵，进而得到差分矩阵，然后利用差分矩阵剔除可能的误匹配点对，构造降阶矩阵，进行模型求解，最后利用求解的模型参数进行匹配点对约束，达到匹配点对提纯目的；
本发明分析了奇异值分解中奇异值的对应含义，并将其引用到匹配点对的提纯过程中，从而避免了初始内点集的选择，大大提高了匹配点对提纯的效率和准确率。
附图说明
图1为本发明的流程图。
具体实施方式
如图1所示，本发明所述的一种基于奇异值分解的影像匹配提纯方法，包括以下步骤：
a、记c为迭代索引，c＝1，2，3，……，pnum，其中pnum为迭代次数上限。
b、假设点集φ＝{ρk|k＝0，1，2，……，n}中包含n对匹配点对，对于任意一对匹配点对，分别采用基本矩阵建立误差方程，即
m′kfmk＝0；
其中，mk＝(xk,yk,1)t，mk′＝(x′k,yk′,1)t分别为匹配点对齐次坐标，f为基本矩阵，令f＝(fij)，则上式可写为如下形式
(xk′xk,x′kyk,xk′,y′kxk,y′kyk,y′k,xk,yk,1)f＝0；
其中，f＝(f11,f12,f13,f21,f22,f23,f31,f32,f33)t，利用各个匹配点对的误差方程建立约束模型：
其中，a为n×9维的系数矩阵。
c、对系数矩阵a进行奇异值分解，利用a的前t个奇异值及其对应的左、右奇异向量构建a的近似矩阵a′，即假设系数矩阵a的秩为r，则对系数矩阵a进行奇异值分解有
a＝usvt；
其中，u为n×n维正交矩阵，v为9×9维正交矩阵，s为n×9维对角阵，u的列向量u1，u2，……，un称为系数矩阵a的左奇异向量，v的列向量v1，v2，……，v9称为系数矩阵a的右奇异向量，s的对角线元素σ1，σ2，……，σr称为系数矩阵a的奇异值，且以递减顺序排列，即σ1≥σ2≥……≥σr；
由系数矩阵a的前t个奇异值重新组成对角阵st×t＝diag[σ1,σ2,……,σt]，由奇异值对应的左、右奇异向量重新组成矩阵un×t＝[u1,u2,……,ut]和v9×t＝[v1,v2,……,vt]，则利用可重构的系数矩阵a，得到
a′＝um×tst×tvt×nt；
对上式进行化简可得
其中，a′为a的近似矩阵，由系数矩阵a的主奇异值重构而成，包含了系数矩阵a的主要信息。
计算a与a′间的差分矩阵δa为
计算各对匹配点对的近似误差为
其中，δaij为δa中的元素。
以近似均方误差mse作为阈值，对每对匹配点对的误差方程进行筛选，若δk<mse则保留第k对匹配点对的误差方程，若δk≥mse则舍弃第k对匹配点对对应的误差方程，由此得到约束后的匹配点对误差方程。
d、利用约束后的匹配点对误差方程再次组成新的系数矩阵，得到降维矩阵b，利用降维矩阵b建立约束模型bf＝0，对降维矩阵b进行奇异值分解，求解f的最小二乘解，并进行秩2约束，得到模型参数解f。
e、利用基本矩阵及求解的模型参数f，采用误差方程m′kfmk＝0对各对匹配点对进行约束，若匹配点对的mk′fmk值大于经验值α，则将该点标记为外点，否则标记为内点。
f、判断c值与pnum的大小，当c＜pnum时，判断模型参数f是否达到预定精度，若是则将标记的内点作为提纯点对，完成提纯过程，若否则将c值加1，并将标记的内点作为点集φ中包含的匹配点对，然后返回步骤b，当c＝pnum时，将标记的内点作为提纯点对，完成提纯过程。
判断模型参数f是否达到预定精度时，计算当前f值与上次迭代所得f值的差值，若差值小于设定值β，则判定f达到预定精度，若差值大于等于设定值β，则判定f没有达到预定精度。
本发明将奇异值分解引用到影像匹配点对的提纯过程中，避免了现有提纯方法中初始内点集的选择，大大提高了匹配点对提纯的效率和准确率。