[发明专利]一种卫星光学遥感相机内方元素在轨检校方法

权利要求书

1.一种卫星光学遥感相机内方元素在轨检校方法，其特征在于包括如下步骤：

建立OXYZ三维坐标系，其中相机在坐标系中的位置为O₁，地球中心在坐标系中的位置为O₂，坐标系中X向为线阵推扫方向，Y向为垂轨方向，Z向为遥感相机视轴，设地面某控制点P，向量O₁P与地球相交于P点，O₁P与焦平面相交于P₂点，由于光学部件设计与加工缺陷的存在，视向量O₁P在焦平面的实际成像位置为P₁点，设P₂点的Y坐标为y₂，P₁点的Y坐标为y₁，则根据径向畸变模型公式得出y₂与y₁的关系式：

y₂＝c₀+c₁y₁+c₂y₁²+c₃y₁³

计算模型参数c₀、c₁、c₂和c₃，将计算得到的模型参数c₀、c₁、c₂和c₃代入上式，得到P₂点修正后的Y坐标y₂′，进一步得到拟合残差Δy₂＝y′₂-y₂，调整相机焦距f，使得拟合残差Δy₂最小，完成在轨检校，其中计算模型参数c₀、c₁、c₂和c₃的步骤如下：

(1)选择几何检校区，卫星在几何检校区成像后得到1级图像，同时选择与1级图像具有相同分辨率或更高分辨率的正射影像和DEM高程图作为参考图，在1级图像上选取k个控制点，同时在所述参考图上选取k个同名点，记录1级图像上k个控制点的坐标(m，n)，m表示列号，n表示行号，记录参考图上k个同名点的经纬度坐标和高度值(Lon，Lat，h)，Lon表示经度，Lat表示纬度，h表示高度，其中k为正整数；

(2)根据列号m，计算OXYZ三维坐标系下控制点的y坐标，

y＝c(m-M/2)，其中M表示一行的像元个数，c表示像元尺寸，

得到控制点线阵向量w₀(0，y，f)，进一步归一化得到控制点线阵向量w_i；

(3)根据行号n，计算控制点所在行的成像时刻秒计数t，

t＝n·d-t₀，其中d为积分时间，t₀为1级图像第一行秒计数，

根据秒计数t，利用卫星下传的GPS数据，计算成像时刻卫星在地固系的位置坐标向量P_t，并假设卫星和相机的坐标重合，根据参考图上同名点的经纬度坐标和高度值(Lon，Lat，h)，计算得到地固系下控制点的坐标向量P_j，控制点坐标向量P_j减去卫星位置坐标向量P_t得到控制点视向量v₀，进一步归一化得到控制点视向量v_i；

(4)在k个控制点的线阵向量w_i中选择最接近相机视轴的向量w_ij，计算向量w_ij与视轴的夹角α：

y₀表示离相机视轴最近的控制点的y坐标；

分别以最接近相机视轴的向量w_ij、v_ij为基准，计算其它k-1个向量与所述向量w_ij、v_ij的夹角，得到k-1个控制点线阵向量间夹角Aw_i和k-1个控制点视向量间夹角Av_i，每个夹角都加上一个夹角α，得到Aw’_i和Av’_l；

(5)计算线阵坐标Bw_i和Bv_i，计算公式为：

Bw_i＝f*tan(Aw′_i)

Bv_i＝f*tan(Av′_i)，f为相机焦距；

线阵坐标Bw_i和Bv_i分别对应公式y₂＝c₀+c₁y₁+c₂y₁²+c₃y₁³中的y₁和y₂，将Bw_i和Bv_i代入公式y₂＝c₀+c₁y₁+c₂y₁²+c₃y₁³，用最小二乘法拟合得到模型参数c₀、c₁、c₂和c₃。