[发明专利]基于复共线性分析的RPC参数优选方法

说明书

技术领域

本发明涉及RPC(Rational Polynomial Coefficients)参数的优选方法，属于摄影测量与遥感领域。

背景技术

1999年9月24日，美国空间成像公司(Space Imaging)将IKONOS卫星成功送入预定轨道并接受卫星遥感影像，标志着高空间分辨率卫星遥感时代的到来。自从高空间分辨率遥感卫星IKONOS采用有理函数模型(Rational FunctionModel，RFM)替代以共线条件方程为基础的严格几何处理模型进行影像几何处理以来，RFM受到了遥感学界的广泛关注，并被学者们深入研究，取得了很好的应用效果。Grodecki和Tao等证实了在单线阵推扫式卫星遥感影像几何处理中可用RFM取代严格几何处理模型，可以用于影像纠正、正射影像制作和目标三维重建等。Fraser和Hu等研究了如何利用地面控制点提高RFM精度的方法。Grodecki和刘军等研究了基于RFM的区域网平差方法。

就以上应用而言，RFM实质是对严格几何处理模型的数学逼近。其关键在于精确求解有理函数模型各多项式的系数(Rational Polynomial Coefficients，RPC)，简称为RPC参数。研究表明，对IKONOS和QuickBird高分辨率卫星遥感影像，采用地形无关的解算方案建立RFM可以达到极高的拟合精度。由于RFM形式简单，适用于各种类型的遥感传感器，包括新型的航空/航天遥感传感器，而且无需使用成像过程中的各种几何参数，如卫星星历、传感器姿态角以及物理特性参数和成像方式等；加之RPC参数不具有明确的物理意义，能够很好地隐藏传感器的核心信息。因此，RFM被广泛应用于高分辨率卫星遥感影像的几何处理中，已成为可以替代严格几何处理模型的一种独立于传感器的广义几何处理模型。

然而，为了RFM应用上的便利和保证遥感影像的几何处理精度，常常需要采用地形相关方案求解RPC参数。地面控制点分布不均匀或者RFM过度参数化，常常导致求解RPC参数的误差方程设计矩阵的列向量之间存在近似的线性关系(数学上称之为复共线性)，即RPC参数之间存在相关性，导致按照最小二乘平差原理建立的法方程严重病态，使得解不稳定或者出现较大的偏差。目前采用的岭估计法虽然能在一定程度上改善法方程的病态，但所能达到的精度依然十分有限。Tao和Hu对一景SPOT影像进行过试验，人工量测了71个地面控制点，根据其中的40个控制点按岭估计法求解RPC参数，所得到的控制点和检查点的坐标残差中误差分别为±6.42×10^-3pixels和±6.13pixels。由此不难看出，虽然RFM对控制点的拟合精度很高，但RFM用于高分辨率卫星遥感影像几何处理的精度还是非常有限的。

发明内容

本发明的目的就在于克服上述现有技术的不足而试图从分析求解RPC参数的误差方程设计矩阵列向量间的复共线性着手，提出一种去相关的RPC参数优选方法。根据一定准则优选出适量的显著性RPC参数，剔除部分相关的RPC参数，使RFM中的参数保持相对独立，以消除法方程病态对RPC参数求解精度的影响，实现在缺少地面控制点情况下，能够采用地形相关方案精确求解RPC参数，并保证基于此RPC参数构建RFM的遥感影像几何处理精度。

实现本发明目的采用的技术方案是：首先建立求解RPC参数的严密误差方程；然后分析误差方程设计矩阵列向量间的复共性线，并根据一定的原则优选RPC参数；最后根据最小二乘平差原理求解优选出的RPC参数。

上述严密误差方程的建立是将有理函数模型按照泰勒级数展开得到的，具体过程如下：

有理函数模型将像点坐标(l，s)表示为含地面点坐标(P，L，H)的多项式的比值，即：