一种升降轨dem融合方法的应用.docxVIP

一种升降轨dem融合方法的应用 摘要：针对InSAR技术在地形复杂区域生成的DEM可能存在错误或空洞等无效点导致可用性降低的问题，本文提出了一种升降轨DEM融合方法，利用Sentinel-1 SARVV单视复数升降轨数据干涉处理得到的DEM与已有TanDEM-X升降轨DEM为例，分析升降轨DEM存在问题并将外部辅助DEM的方差作为参考进行升降轨DEM异常点判别，依据干涉对的相干系数与垂直基线高度加权融合得到最终DEM。结果表明，两种产品在地形复杂区域融合得到的DEM异常点占比较升降轨DEM明显减少，Sentinel-1选取较平坦区域融合后DEM异常点占比由升轨、降轨的5.9%和3.7%降至0.1%，融合后地形较完整，TanDEM-X选取包含水域区域与山地区域融合后异常点占比由升轨、降轨的1.2%和1.8%、2.1%和0.7%减少到0.2%和0.3%，融合后DEM较好地消除了异常点与修复了错误区域，可用性提高。 0 引言数字高程模型(Digital Elevation Model,DEM)是重要的地理空间数据，广泛地应用在国民经济、国防建设、人文和自然科学等各个领域升降轨DEM融合已有方法大多采用加权融合本文依据外部辅助DEM进行升降轨DEM异常点判别，利用升降轨干涉对相干性和垂直基线高度进行加权融合以减少DEM中的异常区域，融合后的DEM地形较完整，可用性较高。1 In SAR技术生成DEM及升降轨融合方法In SAR技术由于其主动成像、可穿透云雾等特性，逐渐成为获取大范围地表形态的重要技术手段1.1 In SAR生成DEM基本流程In SAR生成DEM即利用主辅影像干涉相位值和卫星轨道参数求出目标区域高程值式中，ΔφIn SAR生成DEM流程如图1所示。首先，将主辅影像进行配准以减小同名像点的坐标偏移误差，共轭相乘得到干涉图，平地相位，由于SAR侧视成像，接收地面点回波的顺序为由近及远，在高程变化不大的区域仍会出现较密集的干涉条纹，可利用外部参考DEM去除平地相位，相干系数是衡量干涉测量精度的重要指标，在升降轨数据融合时可作为判别是否用来融合的参考指标，将干涉图滤波以减小相位噪声，相位解缠即将去平干涉相位恢复整周模糊数至真实相位，解缠后经过相高转换和地理编码得到DEM。1.2 升降轨DEM融合方法In SAR生成的DEM质量受侧视成像导致的数据缺失首先，分别进行升降轨DEM错误与空洞区域的判别，该类区域一般为高程突变，且升降轨DEM同一区域高程值变化趋势一般不相同，采用滑动窗口方式计算方差，累计计算单个像素所经过的窗口的方差之和，如下式:式中，x式中，H为融合后高程值;h由于利用异常点方差来判别，方差阈值与所选区域地形直接相关，故在DEM融合时将升降轨重合区域选择需要处理区域裁剪小块进行处理。2 实验及分析2.1 实验地区及数据集实验采用的数据包含欧空局(European Space Agency,ESA)发布的4景Sentinel-1 SAR卫星干涉宽幅模式双极化数据，影像分辨率为5 m×20 m，选取VV极化方式C波段单视复数影像，重访周期一般为6 d或12 d，由于Sentinel-1 SAR卫星的设计主要是为了服务合成孔径雷达差分干涉测量(Differential Synthetic Aperture Radar Interferometry,D-In SAR)技术，为选择到较合适的时间、空间基线，将实验区选于土耳其经纬度40°附近，地形起伏较大，地面高程变化范围约为从1 100—3 500 m，所选影像具体参数见表1。已有Tan DEM-X升降轨DEM位于菲律宾某山地地区，包含高山与大面积水域，地形起伏较大，地面高程变化范围约1 800 m，具体参数见表2。2.2 Sentinel-1升降轨DEM生成及融合处理首先，按照目标区域裁剪主辅影像，Sentinel-1 SAR卫星采用TOPSAR成像方式，可按照Burst条带分割，使用精密定轨星历数据(Precise Orbit Ephemerides,POD)以减小轨道误差保证影像配准精度;然后，采用增强谱分集校正相邻Burst间系统偏差，单视复数影像复共轭运算得到干涉相位影像如图4(a)所示(上图为升轨，下图为降轨)，裁剪目标区域，利用1弧秒SRTM去除平地相位，将去除平地相位的干涉图进行Goldstein滤波，得到升降轨干涉图滤波如图4(b)所示，得到相干系数如4(c)所示，相位解缠采用统计费用网络流相位解缠(SNAPHU)方法，初始化流方法为最小费用流，分块大小为1×1，解缠结果如图4(d)所示。由于升轨垂直基线只有123 m，得到的干涉结果中部分干涉条纹较模糊，可能