[发明专利]一种异源光学遥感卫星影像立体像对预处理方法

权利要求书

1.一种异源光学遥感卫星影像立体像对预处理方法，其特征在于，包括如下步骤：

S1、确定两张异源卫星的原始立体影像的基准分辨率和物方水平投影面的高程；

S2、将两张原始立体影像同时投影到物方水平投影面上，然后根据确定的基准分辨率对两张原始立体影像在物方进行重新采样，纠正生成两张立体影像；

S3、利用原始立体影像的RPC参数，结合影像纠正过程的几何变换关系，对纠正生成两张新立体影像进行重新拟合生成新的RPC参数；

S4、根据纠正生成的立体影像，采用SIFT算子获取两张纠正生成的立体影像上密集的同名点，并根据同名点坐标，获取同名点对应的影像灰度；

S5、根据同名像点的灰度值取值，结合影像灰度范围，分成若干段灰度，采用三次样条曲线进行辐射畸变拟合，计算三次样条曲线参数；

S6、根据样条曲线参数对纠正后左像或者右像进行灰度的重新计算，选择辐射质量好的影像作为基准，改正辐射质量差的影像，获得新的左像或者右像；

步骤S1中物方水平投影面及基准分辨率的确定方法包括如下步骤：

a1、利用左右卫星影像RPC参数获取两景影像覆盖区域的地形高程平均值，取平均高程面作为物方水平投影面；

a2、将两张原始立体影像投影到物方水平投影面上，分别统计物方水平投影面积和像方影像面积，两者的比值可以作为计算卫星立体影像分辨率：

其中Area_g表示原始影像在平均高程面的物方水平投影面积；Area_i表示原始影像的面积，即影像宽乘以影像高；gsd表示对应原始立体影像的空间分辨率；

a3、选择两张原始立体影像中高的空间分辨率值作为基准分辨率；

步骤S5采用三次样条曲线进行辐射拟合，并利用段与段之间的梯度平滑约束，即可形成三次样条曲线：

I′(p′)＝a_n·I³(p)+b_n·I²(p)+c_n·I(p)+d_n

式中，p、p′表示同名点；I(p)、I′(p′)表示同名点灰度；a、b、c、d表示样条参数；m、n表示相邻段。