03-InSAR与DInSAR处理技术课件

InSAR与DInSAR处理技术技术支持邮箱：ENVI-IDL@技术支持热线：400-819-2881-7/enviidlInSAR处理技术InSAR技术流程SLC数据对数据聚焦SARRAWData聚焦模块基线估算参考椭球轨道精炼&重去平相位编辑相干系数基于DEM/参考椭球高度生成干涉图&干涉去平多视强度数据

相位解缠相位转高程&地理编码自适应滤波&相干性计算数字高程模型参考DEMCosmoSkyMedSLC像对成像时间：2009-09-07和2009-09-08HH极化HIMAGEMode入射角：43.95°横断山脉(China)范围:4061x6131pixels

(15x15km)空间分辨率:5x5m辅助数据：SRTMDEMInSAR练习——横断山脉InSAR提取DEM数据准备SLC数据读取参考DEM下载选择一套默认的处理参数/SARscape/DefaultValues第一步：基线估算从像对轨道参数中得到基线估算结果第二步：干涉图生成使用参考DEM生成干涉图和去平的干涉图Toolbox/SARscape/Interferometry/1B-InterferometricWorkflow(MultiSteps)/1-InterferogramGeneration生成的干涉图和去平后的干涉图Interferogram_dint干涉图Interferogram_int去平干涉图第三步：干涉图滤波&相干性计算经过滤波的干涉图，降低了相干噪声相干图像去除系统性空间相关去除主从相对的相关性在相位解缠之前确定干涉相位的质量Toolbox/SARscape/Interferometry/1B-InterferometricWorkflow(MultiSteps)/2-AdaptiveFilterandCoherenceGeneration_fint_cc滤波后的干涉图和相干系数图第四步：相位解缠绝对相位p-p未解缠的相位p-p第四步：相位解缠区域增长法、最小费用流算法、改进的最小费用流法设置一个不太大的相关性阈值（0.15和0.2之间）可以避免相位解缠结果中不连续的区域产生“相位孤岛”Toolbox/SARscape/Interferometry/1B-InterferometricWorkflow(MultiSteps)/3-PhaseUnwrapping解缠后的相位_upha第五步：轨道精炼&重去平当轨道参数不够精确，影响从干涉相位转变为地形高度使用GCP重新定义基线参数计算相位偏移(如获取绝对相位值)重新修改解缠图像的头文件中的轨道参数Toolbox/SARscape/Interferometry/2-RefinementandRe-flattening第五步：轨道精炼&重去平第五步：轨道精炼&重去平生成控制点第五步：轨道精炼&重去平重去平前后干涉图结果对比_fintreflat_fint重去平前后相位解缠结果对比_uphareflat_upha第六步：相位转高程Toolbox/SARscape/Interferometry/3A-PhasetoHeightConversionandGeocoding地理编码距离多普勒方法大地测量和制图转换干涉图生成到相位解缠一步设置Toolbox/SARscape/Interferometry/1A-InterferometricWorkflow(SingleStep)InSAR提取的DEM（5米）SARscapeInSAR应用：四川麻溪乡水坝GoogleEarthSRTMDEM聚束式Cosmo-SkymedInSARInSAR获取DEM©Cosmo-SkyMedASIDatadistributedbye-GEOS

withSARscape

TerraSAR-X-1data,3m分辨率,高程精度<1mData©Infoterra

BoliviaSARscapeInSAR应用：从高分星载SAR中获取DEMSARscapeInSAR应用：从机载SAR上获取的DEM

Data©AeroSensing,Processingsarmap

0.5m

采样间隔D-InSAR处理技术重复轨道InSAR测量DEM时实际上假设地表没有变化。实际上在发生地震、火山活动或者地壳运动的情况下，地表会有或大或小的形变。在InSAR技术的基础上，如果重复进行干涉成像或结合已有的精细DEM数据来消除干涉图中地形因素的影响，可以检测出地表的微小形变，这是D-InSAR的技术基础三种技术方法：双过差分、三过差分、四过差分。从可靠性上讲，双过差分干涉最可靠，而且目前全球大部分地区都有免费的SRTM3的DEM，可以满足很多应用需求。地表形变监测——DInSAR地形信息来自已有的高精度DEM假设两幅SAR图像获取的时间段中存在SAR观测斜距方向的形变δr，干涉位相φ1可表示为D-InSAR的方法：双过（2－pass）差分φ0是由DEM按照干涉基线和入射角模拟的位相差分干涉结果只与波长有关，与基线无关实际差分干涉处理时，干涉位相仍然经过去平，而模拟位相只需模拟地形位相即可。需要三幅SAR图像，其中两幅应该是在事件发生前，其中一幅作为公共主图像，其干涉结果得到本底DEM信息；公共主图像与事件后的图像进行干涉，再与得到DEM的干涉图进行差分，得到形变。其原理与二过差分干涉基本相同，但是处理流程有一些不同。D-InSAR的方法：三过差分干涉四过差分干涉的数据要求较高，时空基线对适合的情况比较少。4幅SAR图像中3幅是事件发生前的，1幅是事件发生后的，但是无法找到公共主图像，只能分二组，一组生产DEM，一组包含事件形变信息。差分干涉流程类似于二过差分干涉，不同的是DEM是由另一组干涉生成。D-InSAR的方法：四过差分干涉SLC1dataSLC2dataSLC3dataSLC4data干涉图干涉图差分干涉相干变化DInSAR技术流程SLC数据对数据聚焦SARRAWData聚焦模块基线估算参考椭球轨道精炼&重去平相位编辑相干系数基于DEM/参考椭球高度生成干涉图&干涉去平多视强度数据

相位解缠相位转形变&地理编码自适应滤波&相干性计算形变图参考DEM覆盖地区：伊朗巴布地区数据源：ENVISATASARSLC数据，IS2模式成像时间：震前的2003年12月3日；震后47天2004年1月11日入射角：23度覆盖区域：56*40km，含有卫星轨道信息保证精度；研究区的25米SRTM数据

数据路径：\Basedata\3.DInSARExamplesDInSAR练习(1):基于DInSAR的地表形变监测Toolbox/SARscape/Interferometry/1A-InterferometricWorkflow(SingleStep)干涉图生成到相位解缠一个面板中完成第二步：干涉图生成与相位解缠第三步：轨道精炼&重去平当轨道参数不够精确，影响从干涉相位转变为地形高度使用GCP重新定义基线参数计算相位偏移(如获取绝对相位值)重新修改解缠图像的头文件中的轨道参数Toolbox/SARsc