[发明专利]一种环境一号卫星CCD影像正射校正方法

摘要

一种环境一号卫星CCD图像正射校正方法，首先导入HJ星CCD一级影像以及对应的高程数据与参考数据，然后在每景HJ星CCD一级影像上均匀采集地面控制点，将采集的控制点坐标导入三次项有理函数模型，计算出有理多项式系数(RPC)，即获得其对应的有理函数模型,再将环境一号卫星CCD一级影像所有点的像素坐标输入到有理函数模型中获得对应正射校正后图像各点的坐标，然后采用双线性插值方法将校正后图像各点的灰度值插值到像素坐标为整行整列的点上，获得正射校正图像。本发明方法相对传统多项式方法能获得更高的精度。

主权项

一种环境一号卫星CCD影像正射校正方法，其特征在于步骤如下：(1)导入环境一号卫星CCD的一级影像；导入格网间距为M1米的SRTM高程数据与分辨率为M2米的LandSat全色参考数据；所述M1的范围为[30,150]，M2的范围为[7.5,30]；(2)在步骤(1)导入的一级影像中均匀采集N个地面控制点，获得该控制点在参考影像的空间坐标(P，L，H)和控制点在一级影像中的像素坐标(X，Y)；所述N大于等于19；(3)将步骤(2)采集的地面控制点在参考影像的空间坐标(P，L，H)和控制点在一级影像中的像素坐标(X，Y)导入有理函数模型，建立4N个等式，计算出有理多项式系数，即RPC，确定有理函数模型，具体计算公式如下：P=P1(X,Y,H)P2(X,Y,H)L=R3(X,Y,H)R4(X,Y,H)]]>其中，P、L和H分别为采集的控制点在参考影像的横坐标、纵坐标和高程值；X和Y分别为采集的控制点在一级影像中的像素横坐标和像素纵坐标；P1、P2、P3和P4均为形式相同的多项式，展开后分别为：P1＝a1+a2X+a3Y+a4H+a5XY+a6YH+a7XH+a8Y2+a9X2+a10H2+a11XYH+a12Y3+a13YX2+a14YH2+a15Y2X+a16X3+a17XH2+a18Y2H+a19X2H+a20H3P2＝b1+b2X+b3Y+b4H+b5XY+b6YH+b7XH+b8Y2+b9X2+b10H2+b11XYH+b12Y3+b13YX2+b14YH2+b15Y2X+b16X3+b17XH2+b18Y2H+b19X2H+b20H3P3＝c1+c2X+c3Y+c4H+c5XY+c6YH+c7XH+c8Y2+c9X2+c10H2+c11XYH+c12Y3+c13YX2+c14YH2+c15Y2X+c16X3+c17XH2+c18Y2H+c19X2H+c20H3P4＝d1+d2X+d3Y+d4H+d5XY+d6YH+d7XH+d8Y2+d9X2+d10H2+d11XYH+d12Y3+d13YX2+d14YH2+d15Y2X+d16X3+d17XH2+d18Y2H+d19X2H+d20H3式中，aj,bj,cj,dj，j＝1,2,...,20为有理多项式系数；(4)将环境一号卫星CCD的一级影像所有像素坐标(P，L)以及对应高程值H代入步骤(3)中确定的有理函数模型，得到校正后图像的坐标值(X，Y)；(5)采用内插值方法对步骤(4)中得到的校正后图像的坐标值(X，Y)进行处理，得到整数坐标值像素的灰度值，从而得到校正后的环境一号卫星CCD影像。