[发明专利]利用激光测高数据辅助提升遥感立体影像几何精度的方法

权利要求书

1.一种利用激光测高数据辅助提升遥感立体影像几何精度的方法，其特征在于，包括以下步骤：

S1，根据卫星下传的姿态、轨道数据，利用实验室标定的几何成像模型内外参数，构建立体影像几何成像模型；

S2，基于高分七号激光测高点的物方坐标、足印影像和激光测高点足印影像像点坐标，以及步骤S1构建的几何成像模型，经过影像匹配，获取激光测高点在立体影像上的像点坐标；

S3，对立体影像进行连接点匹配，获取分布均匀、数量足够的连接点，连接点的物方坐标由步骤S1建立的几何成像模型前方交会获得；

S4，根据连接点的像方坐标和物方坐标，分别对几何成像模型内外参数以及连接点物方平面坐标和高程坐标求偏导，构建连接点误差方程；

S5，根据激光测高点的像方坐标和物方坐标，分别对几何成像模型内外参数以及激光测高点物方平面坐标求偏导，构建激光测高点误差方程；

S6，根据最小二乘原理，对步骤S4和S5构建的误差方程进行法化，构建法方程，求解几何成像模型内外参数和连接点、激光测高点物方坐标改正数；

S7，根据步骤S6求得的改正数更新几何成像模型内外参数和连接点物方坐标以及激光测高点物方平面坐标，继续步骤S4-S6，直至迭代结束，得到最终几何成像模型内外参数；

S8，根据新求解的几何成像模型内外参数更新立体影像几何成像模型，提升立体影像几何精度；

步骤S2中激光测高点在立体影像上的像点坐标获取方法具体包括：

S21,将激光测高点物方坐标代入立体影像几何成像模型，计算得到当前激光测高点在立体影像上的初略位置，设其像点坐标为(x，y)，则该点与激光测高点足印影像像点坐标(x₀,y₀)构成共轭点对；

S22,以像点坐标(x,y)为中心，以激光光斑半径为最大搜索窗口，按公式(4)逐点计算与以激光测高点足印影像像点坐标(x₀,y₀)为中心的区域影像相关系数，取最大相关系数点位(x_1’,y_1’)作为像素级配准点位：

式中，g,g’分别为足印影像和立体影像灰度值，为匹配窗口的平均灰度值，w，h分别表示相关系数搜索窗口的宽和高；

S23,以坐标(x₀,y₀),(x_1’,y_1’)为初值，按公式(5)进行最小二乘匹配，获取子像素级配准点位(x₁,y₁)，即可获得激光测高点在当前影像上的像点坐标:

式中，h₀,h₁为影像辐射畸变参数，a_i,b_i为影像几何变换参数；

S24,重复步骤S21-S23，获取当前激光测高点在所有落点影像上的像点坐标；

S25,重复步骤S21-S24，获取所有激光测高点在立体影像中的像点坐标。

2.根据权利要求1所述的利用激光测高数据辅助提升遥感立体影像几何精度的方法，其特征在于，步骤S1中构建的几何成像模型表示为公式(1)：

其中，是影像点在WGS84坐标系下的物方坐标，是卫星GPS接收机相位中心在WGS84坐标系下的坐标，m是比例系数，是相机坐标系到卫星本体坐标系的转换矩阵，是卫星本体坐标系到J2000惯性坐标系的转换矩阵，是J2000惯性坐标系到WGS84坐标系的转换矩阵，(ψ_x,ψ_y)是影像探元的指向角，代表了影像内部精度；

式中，s是影像的列坐标，a_i,b_j(i,j≤5)是需要求解的几何成像模型内参数；

式中，ω,κ分别是围绕卫星本体坐标系Y,X,Z轴旋转角，是需要求解的几何成像模型外参数。