[发明专利]基于投影基准面的卫星立体影像像对近似核线生成方法

权利要求书

1.基于投影基准面的卫星立体影像像对近似核线生成方法，包括以下步骤：

一、以卫星立体影像像对在参考椭球上地面覆盖区域内的任一点为原点，建立局部直角坐标系；投影基准面PRP是定义在局部直角坐标系下的平均高程面，高程为H；

二、对于组成卫星立体影像像对的左右原始影像，令f₁^L为左影像传感器成像几何模型坐标正解函数，f₁^R为右影像传感器成像几何模型坐标正解函数，f₂^L为左影像传感器成像几何模型坐标反解函数，f₂^R为右影像传感器成像几何模型坐标反解函数；基于下列运算，确定卫星立体影像像对核曲线在投影基准面上投影点轨迹的近似直线方向：

(1)在左原始影像的中心取一像点A，在其摄影光线上取两点P₁和P₂，对应的高程分别为H+h和H-h，这里，h可以取任意小于最大高程值的正数，基于式(1)(2)计算P₁和P₂的坐标；

(X_P1，Y_P1)＝f₁^L(l_A，s_A，H+h)            (1)

(X_P2，Y_P2)＝f₁^L(l_A，s_A，H-h)            (2)

上面函数式中：

(l_A，s_A)为像点A的像素坐标；

(X_p1，Y_p1，H+h)为P₁在局部直角坐标系下的坐标；

(X_p2，Y_p2，H-h)为P₂在局部直角坐标系下的坐标；

(2)基于式(3)(4)计算P₁、P₂在右原始影像上对应的投影点B、C的像素坐标；

(l_B，s_B)＝f₂^R(X_P1，Y_P1，H+h)             (3)

(l_C，s_C)＝f₂^R(X_P2，Y_P2，H-h)             (4)

上面函数式中：

(l_B，s_B)为投影点B的像素坐标；

(l_C，s_C)为投影点C的像素坐标；

(X_p1，Y_p1，H+h)为P₁在局部直角坐标系下的坐标；

(X_p2，Y_p2，H-h)为P₂在局部直角坐标系下的坐标；

(3)基于式(5)(6)计算像点B、C在投影基准面上的投影点P₃、P₄的坐标：

(X_P3，Y_P3)＝f₁^R(l_B，s_B，H)             (5)

(X_P4，Y_P4)＝f₁^R(l_C，s_C，H)             (6)

上面函数式中：

(l_B，s_B)为投影点B的像素坐标；

(l_C，s_C)为投影点C的像素坐标；

(X_P3，Y_P3，H)为P₃在局部直角坐标系下的坐标；

(X_P4，Y_P4，H)为P₄在局部直角坐标系下的坐标；

投影点P₃、P₄的连线方向即为卫星立体影像像对核曲线在投影基准面上投影点轨迹的近似直线方向，即卫星立体影像像对近似核线在投影基准面上的排列方向；

三、建立左右核线影像与左右原始影像之间的像点严格坐标变换关系：

(1)基于传感器成像几何模型坐标正解函数，计算左右原始影像在投影基准面上的投影覆盖范围；

(2)确定左右核线影像在投影基准面上的覆盖范围：为了使左右核线影像具有相同行数和列数，且保证同名核线对应的行号相同，令左右核线影像在投影基准面上对应相同矩形覆盖范围，取投影基准面上投影覆盖区域在近似核线排列方向上的最小外廓矩形；

(3)令左右核线影像每个像元所对应的投影基准面采样间隔相同，且等于左右原始影像的地面采样间隔，进而可确定左右核线影像的行数和列数；

(4)于是，左右核线影像像点与其对应的PRP投影点之间的坐标平移、旋转和缩放变换关系表示如下：

上面公式中：

(r_E，c_E)为左右核线影像像点的像素坐标；

(X_p，Y_p)为(r_E，c_E)在投影基准面上投影点的平面坐标；

G为左右原始影像的地面采样间隔；

R是由卫星立体影像像对近似核线在投影基准面上的排列方向确定的二维坐标旋转矩阵；

(dX，dY)为坐标平移量；

(5)基于传感器成像几何模型的正反解函数，可确定PRP投影点与左右原始影像像点之间的坐标变换关系，如式(9)(10)所示；

(l_L，s_L)＝f₂^L(X_p，Y_p，H)

                                                             (9)

(l_R，s_R)＝f₂^R(X_p，Y_p，H)

(X_P，Y_P)＝f₁^L(l_L，s_L，H)

                                                            (10)

(X_P，Y_P)＝f₁^R(l_R，s_R，H)

上面公式中：

(X_p，Y_p，H)为PRP投影点坐标；

(l_L，s_L)为投影点对应的左原始影像像点像素坐标；

(l_R，s_R)为投影点对应的右原始影像像点像素坐标；

因此，组合(7)和(9)便构成了左右核线影像到左右原始影像的像点坐标变换模型；反之，(8)和(10)组合起来便构成了左右原始影像到左右核线影像的像点坐标变换模型；

四、根据上述所建立的左右核线影像像点与左右原始影像像点之间的严格对应关系，采用间接法影像数字纠正的方式进行卫星立体影像像对近似核线重排，生成左右核线影像。