[发明专利]高分辨率卫星遥感定位的组合岭估计光束法平差方法

说明书

技术领域

本方法涉及一种高分辨率卫星遥感定位的光束法平差法，尤其是涉及一种高分辨率卫星遥感定位的组合岭估计光束法平差方法。

背景技术

在空中三角测量中，光束法平差在理论上是最严密的方法，因为这种方法直接使用了像点坐标观测值，但是实验结果表明光束法平差并不是精度最高的，主要是由于还存在残余的系统误差，解决的办法主要是自检校法。李德仁等学者将光束法平差应用于高分辨率卫星影像定位中，但是由于法方程的病态性导致计算的结果不甚理想，因为未知外方位线角元素强相关，现有的解决方法主要有四种：增加虚拟观测方程、合并强相关项、线角分求法、岭估计。在处理高分辨率遥感影像定位法方程病态问题时岭估计特别是广义岭估计近来受到越来越多的重视。但是广义岭估计首先需要获取最小二乘解来解算岭参数，而当法方程系数阵严重病态时，可能无法获得最小二乘解，因为法方程系数阵无法求逆。本方法针对该问题结合文献的方法提出一种组合岭估计方法进行单线阵高分辨率影像的光束法平差，取得了良好的效果，可达0.5m的平面定位精度和0.7m的高程定位精度。

方法内容

本方法的目的就是为了克服上述现有技术存在的缺陷而提供的一种能提高定位精度的高分辨率卫星遥感定位的组合岭估计光束法平差方法

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

一种高分辨率卫星遥感定位的组合岭估计光束法平差方法，其特征在于，包括以下步骤：

1)读取像对数据、控制数据、传感数据，并设定限差与迭代次数；

2)数据预处理；

3)通过L曲线法计算岭参数；

4)将岭参数通过组合岭估计计算残差

5)将残差作为广义岭估计的岭参数矩阵，应用广义岭估计计算平差

6)判断是否中的元素值是否超过限差或迭代次数未超过设定值，若为是，返回步骤3)，若为否，执行步骤7)；

7)打印计算结果并进行精度分析。

所述的步骤2)数据预处理包括空间后方交会、空间前方交会。

所述的步骤3)中的岭参数计算过程如下：

通过L曲线法计算岭参数K值，即选择L曲线上曲率最大的点相对应的K值作为所求的参数，根据岭估计的估计准则可得和都是岭参数K的函数，选择不同的K值并以为横坐标、为纵坐标画图，得到L曲线，则L曲线上的点的曲率q的计算公式为：

q=ρ^′η^′′-ρ^′′η^′[(ρ^′)2+(η^′)2]3/2]]>

对上式求最大值，就可以得到最大曲率所对应的K值即为岭参数。

所述的4)中的残差计算过程如下：将步骤3)计算得到岭参数K代入可计算得到估计的残差

所述的5)中的平差计算过程如下：将步骤4中的残差作为广义岭估计的岭参数矩阵并将其带入可求得平差