[发明专利]面向InSAR的初始DSM的区域网平差方法

权利要求书

1.一种面向InSAR的初始DSM的区域网平差方法，其特征在于，所述方包括以下步骤：

步骤S1，利用合成孔径雷达干涉(InSAR)生成初始DSM及附加数据，所述附加数据包括主影像头文件、从影像头文件、控制点质量评估文件、初始DSM对应的SRTM和初始DSM对应的SRTM误差；

步骤S2，构建斜距-多普勒方程和椭球方程，并利用构建的所述斜距-多普勒方程和椭球方程，计算所述控制点质量评估文件中的控制点在主影像中的像素坐标；

步骤S3，根据各个所述初始DSM的经纬度范围，确定各个所述初始DSM的邻接关系矩阵；

步骤S4，基于所述邻接关系矩阵，均匀抽取各个所述初始DSM重叠处的连接点，并计算所述连接点在所述初始DSM中的像素坐标以及与所述初始DSM相邻的影像中的像素坐标；

步骤S5，基于所述连接点在所述初始DSM中的像素坐标，确定所述连接点的经相位高程转换之后的高程值；

步骤S6，基于所述控制点的像素坐标和所述连接点在所述初始DSM中的像素坐标以及所述连接点的高程，构建最小二乘平差方程，对所述初始DSM进行区域网平差，确定所述初始DSM的误差多项式系数。

2.根据权利要求1所述的方法，其特征在于，步骤S3具体包括以下子步骤：

步骤S3.1，构建各个所述初始DSM的初始邻接关系矩阵，所述初始邻接关系矩阵F表达式为：

其中，f(i,j)为第i列第j行的初始DSM，g为东西向的初始DSM数量，k为南北向的初始DSM数量；

步骤S3.2，对步骤S3.1中构建的初始邻接关系矩阵中的所有初始DSM逐一进行迭代分析，确定各个初始DSM的邻接关系矩阵。

3.根据权利要求2所述的方法，其特征在于，步骤S3.2具体包括以下子步骤：

步骤S3.2.1，确定位于最西北角的初始DSM，并将该DSM记为f(0,0)；

步骤S3.2.2，分别确定与所述最西北角的初始DSM的东向邻接的初始DSM和南向邻接的初始DSM；

步骤S3.2.3，依次判断除最西北角初始DSM的其他初始DSM的邻接初始DSM，以确定所有初始DSM的邻接初始DSM，并最终确定初始DSM的邻接关系矩阵。

4.根据权利要求3所述的方法，其特征在于，将所述最西北角的初始DSM f(0,0)的经度范围记为[lons1,lone1]，纬度范围记为[lats1,late1]；将任一其他初始DSM f(i,j)的经度范围记为[lons2,lone2]，纬度范围记为[lats2,late2]；

其中，确定与所述最西北角的初始DSM的东向邻接的初始DSM具体包括以下子步骤：

(1)判断任一其他初始DSMf(i,j)与所述最西北角的初始DSMf(0,0)在经度方向是否有重叠，具体包括：

如果lons2lons1，且lons2lone1，则判定f(i,j)与f(0,0)在经度方向有重叠，并继续执行步骤(2)，否则判定无重叠，即，判定f(i,j)与f(0,0)在东西向无邻接关系，不需要执行下述步骤；

(2)计算f(i,j)与f(0,0)的纬度重叠范围；

纬度重叠范围表达为

lat_com＝min(late2,late1)-max(lats2,lats1)

其中，lat_com为纬度方向的重叠范围，min()为求取最小值，max()为求取最大值；

(3)判断f(i,j)与f(0,0)是否具备东向邻接关系；

具体的，如果lat_com(late1-lats1)*0.7或者lat_com(late2-lats2)*0.7，则判定f(i,j)与f(0,0)具备东向邻接关系，并赋值i＝0+1，j＝0，以寻找f(1,0)的东向和南向邻接DSM；否则，判定f(i,j)与f(0,0)不具备东向邻接关系，并确定与所述最西北角的初始DSM的南向邻接的初始DSM。