复杂交会条件下区域网平差连接点自适应定权方法及系统技术方案

本发明专利技术提出一种复杂交会条件下区域网平差连接点自适应定权方法及系统，根据卫星影像初始RPC模型，以影像仿射变换参数和像点物方坐标作为参数构建区域网平差模型；解算区域网平差模型，得到连接点物方空间坐标的平差值；由连接点所对应的各影像上同名像点，计算连接点对应同名光线之间的最大交会角，并以该交会角作为观测误差方程的权矩阵中所有同名像点的权，返回由新的权矩阵重新解算平差模型，直到平差结果收敛。根据不同连接点的光线交会角确定不同权，能自适应不同平差模型和情形，充分体现了不同贡献的平差量应取定不同权的原则，提高了平差结果的可靠性。本发明专利技术可行、有效，平差解算结果稳定、可靠。

【技术实现步骤摘要】
复杂交会条件下区域网平差连接点自适应定权方法及系统
本专利技术属于遥感影像处理领域，涉及光学卫星区域网平差过程中在复杂交会条件下连接点自适应定权的方法。
技术介绍
连接点是进行区域网平差的重要载体。由于光学卫星影像区域网覆盖范围广，区域网中影像数量大，相互重叠，交会关系复杂，这些复杂的影像间相对定向关系都要通过连接点反映。在普通的平差解算方法中，由于连接点是匹配得到，没有经过人工量测，因此认为各连接点像点观测值是等权观测值，这样的设定在区域网范围小、网中连接点数量较少的情况下对平差结果影响不大。然而光学卫星影像具有覆盖范围广、数量多、交会复杂的特点，经自动匹配出来的连接点对区域网平差贡献差别很大，不再适用等权观测的假定，通过改变区域网平差模型中连接点像点观测值的权矩阵可以解决这一问题。但尚未有相关技术方案出现。
技术实现思路
本专利技术所要解决的问题是一种光学卫星区域网平差过程中连接点观测值自适应定权技术方案。本专利技术的技术方案提供一种复杂交会条件下区域网平差连接点自适应定权方法，包括以下步骤：步骤1，根据卫星影像初始RPC模型，以影像仿射变换参数和像点物方坐标作为参数构建区域网平差模型；设像点坐标为(s,l)，构建的平差模型公式如下，其中，Fx(Lat,Lon,Height)、Fy(Lat,Lon,Height)分别为沿轨方向和垂轨方向的像点坐标函数模型，Δl、Δs分别为沿轨方向和垂轨方向的像方附加参数模型，a0、a1、a2、b0、b1、b2为仿射变换参数；步骤2，解算区域网平差模型，得到连接点物方空间坐标的平差值，包括对于连接点物方坐标其赋予初值并进行线性化处理得到观测误差方程，进行一次平差解算；步骤3，判断是否收敛，是则输出平差结果，否则进入步骤4；步骤4，由连接点所对应的各影像上同名像点，计算连接点对应同名光线之间的最大交会角，并以该交会角作为观测误差方程的权矩阵中所有同名像点的权，返回步骤2由新的权矩阵重新解算平差模型，直到平差结果收敛。而且，权矩阵初始给定为单位阵I。本专利技术相应提供一种复杂交会条件下区域网平差连接点自适应定权系统，包括以下模块：第一模块，用于根据卫星影像初始RPC模型，以影像仿射变换参数和像点物方坐标作为参数构建区域网平差模型；设像点坐标为(s,l)，构建的平差模型公式如下，其中，Fx(Lat,Lon,Height)、Fy(Lat,Lon,Height)分别为沿轨方向和垂轨方向的像点坐标函数模型，Δl、Δs分别为沿轨方向和垂轨方向的像方附加参数模型，a0、a1、a2、b0、b1、b2为仿射变换参数；第二模块，用于解算区域网平差模型，得到连接点物方空间坐标的平差值，包括对于连接点物方坐标其赋予初值并进行线性化处理得到观测误差方程，进行一次平差解算；第三模块，用于判断是否收敛，是则输出平差结果，否则命令第四模块工作；第四模块，用于由连接点所对应的各影像上同名像点，计算连接点对应同名光线之间的最大交会角，并以该交会角作为观测误差方程的权矩阵中所有同名像点的权，命令第二模块工作，由新的权矩阵重新解算平差模型，直到平差结果收敛。而且，权矩阵初始给定为单位阵I。本专利技术的优点在于：根据不同连接点的光线交会角确定不同权，能自适应不同平差模型和情形，充分体现了不同贡献的平差量应取定不同权的原则，提高了平差结果的可靠性。附图说明图1为本专利技术实施例的流程示意图；图2为本专利技术中交会角大小与高程误差关系示意图；图3为本专利技术中交会角计算示意图。具体实施方式以下结合附图和实施例详细说明本专利技术具体实施方式。本专利技术提出一种根据连接点反映的光线交会角大小来更新同名像点观测值权矩阵的方法。根据不同贡献的平差量取不同权的原则，在一次平差解算完成后，以物方点作为标准将所有像点归类，即把所有对应同一物方点的像点归为一类，然后结合平差得到的物方点三维坐标，计算这类像点中每组同名像点两两之间的交会角，取最大交会角作为下次平差中这类像点观测值的权。这种方法能在区域网平差解算过程中自适应计算权矩阵，可应对不同区域网情况，提高了平差结果的可靠性。参见图1，实施例的流程具体实现如下：1.构建区域网平差模型并解算影像序列中有多副影像1、2、…m，首先构建自由网平差模型，然后进行连接点物方坐标解算，接着可以判断是否收敛，是则输出平差结果，否则进入第2部分。对于单景卫星影像，根据初始RPC模型建立以影像仿射变换参数和连接点物方空间坐标为参数的连接点像点观测值的误差方程。平差过程中应用的RPC模型为：式中，Snαβγ,SdαβγLnαβγ,Ldαβγ表示有理多项式系数(rationalpolynomialcoefficients,RPC)，(Uα,Vβ,Wγ)简化记为(U,V,W)，(U,V,W)和(x,y)分别表示正则化的地面点大地坐标和像点坐标，它们与非正则化的地面点大地坐标(Lat,Lon,Height)和像点坐标(s,l)之间的关系如公式(2)，其中，LAT_OFF、LAT_SCALE、LON_OFF、LON_SCALE、H_OFF和H_SCALE为地面坐标的正则化参数；S_OFF、S_SCALE、L_OFF和L_SCALE为影像坐标的正则化参数。附加影像仿射变换参数模型后，则平差模型公式为：式中，Fx(Lat,Lon,Height)、Fy(Lat,Lon,Height)分别为沿轨方向和垂轨方向的像点坐标函数模型，Δl、Δs分别为沿轨方向和垂轨方向的像方附加参数模型，a0、a1、a2、b0、b1、b2为仿射变换参数。对于连接点物方坐标(Lat,Lon,Height)而言，式为一非线性方程，需要对其赋予初值(Lat,Lon,Height)0并进行线性化处理得到观测误差方程：式中，vl、vs分别为影像上沿轨和垂轨方向的改正数，写成矩阵形式即：V＝Ax+Bt-LP(7)其中，V代表像点坐标观测值残差向量；x＝[X1…Xi…Xm]T(i＝1,2…m)代表各景影像RPC像方附加模型参数向量，Xi＝(a0,a1,a2,b0,b1,b2)i代表影像Imgi的RPC模型像方附加参数向量，m代表待平差影像数；t＝[T1…Tj…Tn]T(j＝1,2…n)代表各连接点物方坐标改正值向量，Tj＝d(Lat,Lon,Height)j代表连接点TPj的物方坐标改正数，n代表连接点个数；A、B则分别为对应未知数的偏导数系数矩阵，L和P分别为相应的常向量和权矩阵。具体实施时，本领域技术人员可自行预先设定初值(Lat,Lon,Height)0，或者设定初值求取规则，可采用现有技术，例如使用初始定向参数，采用前方交会的方式确定初值。其中连接点像点观测值的权矩阵为P，初始给定为单位阵I，平差过程中根据平差结果更新变化。观测误差方程建立完成后即可进行一次平差解算，计算待解参数的法方程为：依照矩阵运算法则即可求出待解参数改正数。2.计算各连接点对应同名光线间的最大交会角，更新连接点像点观测值权矩阵，然后可重新迭代进行连接点物方坐标解算连接点在多张影像上有同名像点，多条同名光线两两之间相交。如图2所示，L为左片上一像点坐标，该点与右片上的点R1，R2为可能为同名像点，同名光线的交会角分别为θ1,θ2，左片的像点L有一像方误差dx，该误差在物方引起的误差分别为dh1，dh2本文档来自技高网...

【技术保护点】
一种复杂交会条件下区域网平差连接点自适应定权方法，其特征在于，包括以下步骤：步骤1，根据卫星影像初始RPC模型，以影像仿射变换参数和像点物方坐标作为参数构建区域网平差模型；设像点坐标为(s,l)，构建的平差模型公式如下，l+Δl=Fx(Lat,Lon,Height)s+Δs=Fy(Lat,Lon,Height)]]>Δl=a0+a1l+a2sΔs=b0+b1s+b2s]]>其中，Fx(Lat,Lon,Height)、Fy(Lat,Lon,Height)分别为沿轨方向和垂轨方向的像点坐标函数模型，Δl、Δs分别为沿轨方向和垂轨方向的像方附加参数模型，a0、a1、a2、b0、b1、b2为仿射变换参数；步骤2，解算区域网平差模型，得到连接点物方空间坐标的平差值，包括对于连接点物方坐标其赋予初值并进行线性化处理得到观测误差方程，进行一次平差解算；步骤3，判断是否收敛，是则输出平差结果，否则进入步骤4；步骤4，由连接点所对应的各影像上同名像点，计算连接点对应同名光线之间的最大交会角，并以该交会角作为观测误差方程的权矩阵中所有同名像点的权，返回步骤2由新的权矩阵重新解算平差模型，直到平差结果收敛。...

【技术特征摘要】
1.一种复杂交会条件下区域网平差连接点自适应定权方法，其特征在于，包括以下步骤：步骤1，根据卫星影像初始RPC模型，以影像仿射变换参数和像点物方坐标作为参数构建区域网平差模型；设像点坐标为(s,l)，构建的平差模型公式如下，其中，Fx(Lat,Lon,Height)、Fy(Lat,Lon,Height)分别为沿轨方向和垂轨方向的像点坐标函数模型，Δl、Δs分别为沿轨方向和垂轨方向的像方附加参数模型，a0、a1、a2、b0、b1、b2为仿射变换参数；步骤2，解算区域网平差模型，得到连接点物方空间坐标的平差值，包括对于连接点物方坐标其赋予初值并进行线性化处理得到观测误差方程，进行一次平差解算；步骤3，判断是否收敛，是则输出平差结果，否则进入步骤4；步骤4，由连接点所对应的各影像上同名像点，计算连接点对应同名光线之间的最大交会角，并以该交会角取定观测误差方程的权矩阵中所有同名像点的权，包括取定权值为同名光线之间的最大交会角正相关的dr为高差H导致的像点位移，返回步骤2由新的权矩阵重新解算平差模型，直到平差结果收敛。2.如权利要求1所述复杂交会条件下区域网平差连接点自适应定权方法，其特征在于：权矩阵初始给定为单位阵I。3.一种复杂交会条...

【专利技术属性】
技术研发人员：杨博，王密，李德仁，
申请(专利权)人：武汉大学，
类型：发明
国别省市：湖北;42