面向InSAR的初始DSM的区域网平差方法技术

本发明专利技术公开了一种面向InSAR的初始DSM的区域网平差方法，该方法包括以下步骤：步骤S1，生成初始DSM及附加数据；步骤S2，计算控制点在主影像中的像素坐标；步骤S3，确定各个初始DSM的邻接关系矩阵；步骤S4，均匀抽取各个初始DSM重叠处的连接点，并计算所述连接点在所述初始DSM中的像素坐标以及与所述初始DSM相邻的影像中的像素坐标；步骤S5，确定所述连接点的经相位高程转换之后的高程值；步骤S6，构建最小二乘平差方程，对所述初始DSM进行区域网平差。本发明专利技术能够支持业务化InSAR地形测绘过程中，流程化生产的初始DSM产品的精度的进一步提升。升。升。

【技术实现步骤摘要】
面向InSAR的初始DSM的区域网平差方法


[0001]本专利技术涉及图像处理
，特别涉及一种面向InSAR的初始DSM的区域网平差方法。

技术介绍

[0002]合成孔径雷达干涉(InSAR)技术已经成为了最为有效的全球测图手段之一。2000年，SRTM提供了全球56
°
S～60
°
N范围内的数字表面模型(DSM)数据，其标称精度达到了16m，格网大小30米，在全球尺度的科研以及工程应用中发挥了不可估量的作用。在2016年，德宇航宣布使用TanDEM
‑
X双星绕飞星座完成了全球DSM的测绘，90％地区的精度为3.5米，如果不考虑极低精度区域，全球90％地区的精度则达到了0.8米，数据已经可以用于精细的军事化定位和军事目标识别。InSAR地形测绘的巨大潜力，使得InSAR技术在工程应用上迎来了一个新的时代。然而我国SAR卫星技术发展缓慢，现阶段尚没有民用的业务化InSAR地形测绘卫星，这导致了我国的应用技术发展迟缓。InSAR地形测绘技术，特别是面向大规模全球尺度的地形测绘技术，是我国将本文档来自技高网...

【技术保护点】

【技术特征摘要】
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)在经度方向是否有重叠，具体包括：如果lons2>lons1，且lons2<lone1，则判定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。5.根据权利要求4所述的方法，其特征在于，确定与所述最西北角的初始DSM的南向邻接的初始DSM具体包括以下子步骤：(1)判断f(i,j)与f(0,0)在纬度方向是否有重叠，具体包括...

【专利技术属性】
技术研发人员：李涛，唐新明，陈乾福，周晓青，高小明，张祥，李参海，王怀，禄競，张文君，
申请(专利权)人：自然资源部国土卫星遥感应用中心，
类型：发明
国别省市：