[发明专利]一种投影后影像的区域网平差方法

权利要求书

1.一种投影后影像的区域网平差方法，其特征在于包括以下内容：

获取遥感影像及其地图投影信息；

对所有遥感影像进行匹配，获取遥感影像间的同名点；

建立误差方程，求解遥感影像的精化参数，具体过程为：

利用影像的投影信息，将同名点影像坐标转换为投影坐标：

其中，为影像左上角点的投影坐标，dx,dy分别为影像列和行方向的分辨率；(x,y)为影像的行列坐标，(Dx,Dy)为(x,y)对应的投影坐标；

根据区域网平差数学模型建立误差方程，采用最小二乘法求解误差方程，得到每张遥感影像的精化参数，具体为：

当区域网平差的数学模型为仿射模型时，其精化坐标表示为：

其中，j为影像编号，i为点编号，为平差前的坐标，为精化后的坐标，为仿射变换系数；

对于同名点，其真实坐标值由未知数(Dx_i,Dy_i)表示；

误差方程为：

其中，分别为列方向和行方向上的坐标残差，当影像数为m，点数为n时，未知数总数为6*m+2*n；

误差方程可以写为：

v＝Ax＝0

当在求解过程中引入控制点时，控制点的投影坐标(Dx_ci,Dy_ci)为已知值，控制点误差方程为：

v₁＝A₁x-L₁＝0

L₁＝[Dx_c1 Dy_c1 ... Dx_ck Dy_ck]^T

其中，k为控制点的数量，控制点分布在l张影像上，p和q为控制点所在影像的编号；

综合同名点和控制点，上式利用分块矩阵的方式表示为：

其中，A₁为控制点误差方程系数，A₂，B为同名点误差方程系数，L₁为控制点坐标，X为第一类未知数，即精化系数，Y为第二类未知数，即同名点坐标值；在平差过程中，先消除第二类未知数，单独求解第一类未知数，方程为：

当在求解过程中无控制点时，A₁＝0，方程为：

根据最小二乘法求出第一类未知数X，求解出每张遥感影像的精化参数；

利用精化参数对遥感影像进行重采样；

影像镶嵌，生产影像地图。

2.根据权利要求1所述的投影后影像的区域网平差方法，其特征在于，对所有遥感影像进行匹配，获取遥感影像间的同名点，具体过程为：

根据遥感影像投影信息左上角点投影坐标值、格网大小及遥感影像的行数和列数，计算遥感影像的范围：

其中，为左上角点投影坐标值，(dx,dy)为行列方向的分辨率，sample，line为遥感影像的大小，为右下角点的投影坐标值；

根据影像范围，计算影像的重叠区，并在重叠区利用自动匹配方式计算影像同名点，其中，影像的重叠区为地理坐标矩形范围的交集，影像的同名点是指影像对同一地区摄影时，同一地物在不同影像上对应的像点。

3.根据权利要求1所述的投影后影像的区域网平差方法，其特征在于，区域网平差的数学模型为平移、仿射或多项式模型。

4.据权利要求1或2所述的投影后影像的区域网平差方法，其特征在于，利用精化参数对遥感影像进行重采样，具体过程为：获取每幅影像的精化参数，利用精化参数重采样影像，生成DOM产品，DOM产品的生成采用正解法或反解法。