04-InSAR与DInSAR处理技术论述.pptx

InSAR与DInSAR处理技术;InSAR处理技术;InSAR技术是利用雷达系统获取同一地区两幅SAR影像所提供的相位信息进行干涉处理，来获取地表的三维信息，可以建立目标地区的数字高程模型。 InSAR技术的另外一个重要应用是获取地表的变形;InSAR技术流程;COSMO SkyMed SLC 像对 成像时间：9/7/2009 和 9/8/2009 极化方式：HH 极化 HIMAGE mode 入射角：43.95° 覆盖地区：横断山脉 范围: 4060 x 6130pixels (约为12 x 18 km2) 像元间隔: 3 x 3 m 辅助数据：90米分辨率的SRTM DEM。 ;选取一套适合于COSMO SkyMed HIMAGE 3米数据的系统参数 /SARscape/Preferences Load Preferences—VHR（better than 10m）;从像对轨道参数中得到基线估算结果 /SARscape/Interferometry/Interferometric Tools/Baseline Estimation;Normal Baseline (m) = -180.532 空间基线 Critical Baseline min-max (m) = [-6400.098] - [6400.098] 临界基线 2 PI Ambiguity height (m) = 58.916 2π模糊高程 2 PI Ambiguity displacement (m) = 0.016 2π模糊位移 Range Shift (pixels) = -0.015 Azimuth Shift (pixels) = 0.662 Doppler Centroid diff. (Hz) = 719.488 多普勒图心偏差 Critical min-max (Hz) = [-3363.229] - [3363.229] 多普勒图心偏差临界值 Time Baseline (Days) = 1.000 时间基线 PlotSampleNumber=100 用于基线估算的采样点个数;InSAR处理流程化工具 /SARscape/Interferometry/InSAR DEM Workflow ;InSAR处理工作流; 生成干涉图和使用参考DEM对干涉图去平 ;;对上一步生成的去平地效应之后的干涉图进行滤波处理，降低了相干噪声 ;Interferogram_fint;;区域增长法、最小费用流算法、改进的最小费用流法 设置一个不太大的相关性阈值（0.15和0.2之间） ;;用于轨道精炼;当轨道参数不够精确，影响从干涉相位转变为地形高度 使用GCP重新定义基线参数 计算相位偏移(如获取绝对相位值) 重新修改解缠图像的头文件中的轨道参数 ;;;;D-InSAR处理技术;重复轨道InSAR测量DEM时实际上假设地表没有变化。实际上在发生地震、火山活动或者地壳运动的情况下，地表会有或大或小的形变。 在InSAR技术的基础上，如果重复进行干涉成像或结合已有的精细DEM数据来消除干涉图中地形因素的影响，可以检测出地表的微小形变，这是D-InSAR的技术基础 三种技术方法：双过差分、三过差分、四过差分。从可靠性上讲，双过差分干涉最可靠，而且目前全球大部分地区都有免费的SRTM3的DEM，可以满足很多应用需求。 ;地形信息来自已有的高精度DEM 假设两幅SAR图像获取的时间段中存在SAR观测斜距方向的形变δr ，干涉位相φ1可表示为;需要三幅SAR图像，其中两幅应该是在事件发生前，其中一幅作为公共主图像，其干涉结果得到本底DEM信息；公共主图像与事件后的图像进行干涉，再与得到DEM的干涉图进行差分，得到形变。其原理与二过差分干涉基本相同，但是处理流程有一些不同。;四过差分干涉的数据要求较高，时空基线对适合的情况比较少。4幅SAR图像中3幅是事件发生前的，1幅是事件发生后的，但是无法找到公共主图像，只能分二组，一组生产DEM，一组包含事件形变信息。差分干涉流程类似于二过差分干涉，不同的是DEM是由另一组干涉生成。;DInSAR技术流程;两景ENVISAT ASAR SLC数据作为例子 ENVISAT ASAR SLC像对 IS2模式，VV极化 入射角：22.74° 成像时间：12/3/2003（震前）和2/11/2004（震后） 地区：伊朗巴姆地区 覆盖范围：56*40km，含有卫星轨道信息保证精度； 像元间隔：25*25米 辅助数据：25米分辨率的DEM数据，作为参考DEM。;选择适用于ASAR的一套系统参数 /SARscape/Preferences Load Preferences—