一种RPC模型正变换方法

摘要

本发明涉及一种RPC模型正变换方法，该方法将反变换结果引入正变换过程进行迭代变换，其过程包括以下步骤：设置高程、和RPC模型参数，设定地面矩形区域；计算地面矩形区域4角点和中心点的影像像素坐标；计算仿射变换系数；计算量测像点坐标对应的近似经纬度坐标；计算理论像点坐标；计算残差e，如果残差e满足设定的差值范围，输出近似经纬度坐标作为量测像点坐标对应的经纬度坐标，否则调整地面矩形区域，然后返回迭代，直到残差e满足差值范围。本发明的技术方案能够确保在影像供应商提供的有限条件下获得精确的地面点位置，在遥感影像应用领域具有重要意义。

主权项

1.一种RPC模型正变换方法，其特征在于：将RPC模型反变换结果引入RPC模型正变换过程进行迭代变换，其过程包括以下步骤，步骤1，设置高程数值Dhei和RPC模型参数，设定中心点经纬度坐标为(Dlat_off，Dlon_off)、边长为Dlat_scale和Dlon_scale的地面矩形区域，其中4个角点的经纬度坐标为(Dlat_off+Dlat_scale/2.0，Dlon_off+Dlon_scale/2.0)；(Dlat_off-Dlat_scale/2.0，Dlon_off+Dlon_scale/2.0)；(Dlat_off+Dlat_scale/2.0，Dlon_off-Dlon_scale/2.0)；(Dlat_off-Dlat_scale/2.0，Dlon_off-Dlon_scale/2.0)；步骤2，利用RPC反变换模型计算地面矩形区域4个角点和中心点对应的影像像点坐标；步骤3，利用上述5个点的经纬度坐标和影像像点坐标，根据仿射变换关系式$\begin{array}{c}{D}_{\mathrm{lat}}=f_0+f_{1}x+f_{2}x\\ D_{\mathrm{lon}}=g_0+g_{1}x+g_{2}x\end{array},$ 利用最小二乘原理计算影像的近似仿射变换参数f0，f1，f2，g0，g1，g2，x、y分别表示正则化的影像像点坐标的列数值和行数值，Dlat、Dlon分别表示正则化的地面点坐标的纬度数值和经度数值；步骤4，根据影像的近似仿射变换参数f0，f1，f2，g0，g1，g2求解量测给定的量测像点坐标(xa，ya)对应的近似经纬度坐标(Dlata，Dlona)；步骤5，根据RPC模型的反变换模型计算地面点坐标(Dlata，Dlona，Dhei)对应的理论像点坐标(xb，yb)；步骤6，计算残差e，计算公式为e＝(|xa-xb|)2+(|ya-yb|)2，如果残差e满足设定的差值范围，直接输出(Dlata，Dlona)作为量测像点坐标(xa，ya)对应的经纬度坐标(Dlatb，Dlonb)，否则以(Dlata，Dlona)为中心点经纬度坐标缩小地面矩形区域范围，然后返回步骤2迭代运行，直到残差e满足设定的差值范围。