[发明专利]通过微修正RPC参数实现卫星遥感影像高精度几何纠正方法

权利要求书

1.一种通过微修正RPC参数实现卫星遥感影像高精度几何纠正方法，包括以下步骤：

(1)以RFM模型为基础，建立像点行列坐标参数R，C显函数形式的RFM方程；

(2)在影像控制点的支持下，采用传统的像方仿射变换定向方法，求解得到影像的像方仿射变换模型系数m₀,m₁,m₂,n₀,n₁,n₂，并获取影像控制点像方行坐标的定向精度m_R和像方列坐标的定向精度m_C；

(3)以RFM模型为基础，按照函数的求导算法，利用影像像点行坐标参数R对RPC参数中有理多项式系数a_i,b_j的导数以及影像像点列坐标参数C对RPC参数中有理多项式系数c_i,d_j的导数计算式，其中i＝1…20；j＝2…20，计算所有影像控制点物方坐标点处的78个参数的导数值，并记录各导数绝对值在不同影像控制点处的各自最大值，分别记为其中，i＝1…20；j＝2…20；

(4)以像点行列坐标参数R，C为观测值，以RPC参数中的有理多项式系数a_i，b_j，c_i，d_j的增量da_i,db_j,dc_i,dd_j作为未知数依次构建所有影像控制点的误差方程组，其中，i＝1…20；j＝2…20，以及将da_i,db_j,dc_i,dd_j作为虚拟观测值构建相应的误差方程组；

(5)根据步骤(2)得到的像方定向精度m_R,m_C和步骤(3)得到的自适应构造法方程的正则化矩阵P_I和正则化参数s，并求法方程的正则化解；

正则化矩阵P_I为对角线元素非0而其他元素为0的矩阵，对角线元素p_kk按下式计算，其中p_kk表示第k行k列元素，k＝1,2……78；

其中分别为步骤(3)中计算得到的RPC参数的导数在各控制点物方坐标位置导数绝对值的最大值

正则化参数s按下式计算：

s＝k₁·N_P/39

其中N_p为影像控制点的数目，39为行列坐标误差方程中未知数的个数，k₁为平衡系数，取值1.0或按经验值给定；

法方程的正则化解按下式计算：

X＝(A^TP_pA+sP_I)^-1(A^TP_pL_p+sP_IL_I)

其中A^TP_pA为像点误差方程的法方程，A为像点误差方程系数矩阵，P_p为像点量测坐标的权矩阵，P_I为步骤(5)中得到的正则化矩阵，L_p为像点坐标误差方程的常数项向量，L_I为步骤(4)中78个RPC参数增量未知数da_i，db_j，dc_i，dd_j误差方程的常数项向量其中，i＝1…20；j＝2…20。