SARscape数据处理SAR数据笔记——DINSAR微小形变提取

转发来源：https://blog.csdn.net/u011322358/article/details/105436686/

在InSAR技术的基础上，如果重复进行干涉成像或结合己有的精细DEM数据来消除干涉图中地形因素的影响，可以检测出地表的微小形变，这是。DInSAR的技术基础。
数据类型：震间两景SentinellA IW slc VV极化方式影像。

• 1.设置系统参数

/SARsca pe/Preferences，设置Load Preferences为SENTINEL TOPSAR，这套参数是专门针对哨兵1A数据的丁OPSAR (IW)模式做InSAR处理时的系统参数。

• 2.数据导

入
网站下载的数据解压，打开数据导入工具//SARscape/Import Data/SAR Spaceborne/SENTINEL 1，在Input File List分别输入两景哨兵1A数据的元数据文件 manifest.safe。 Optional Input Orbit File List，轨道文件，该文件是可选文件，在这里没有轨道文件，就 不输入。
输出的数据文件包括：整景图像的强度图数据(_pwr), slc索引文件(.slc_list), 带有地理坐标的外边框矢量文件（.shp）,GoogIeEarth文件(.kml)。打开矢量文件，view->Reference Map Link,可以自动链接Arcgis Online的在线底图，查看数据的地理范围。

• 3.研究区裁剪

在Google earth上根据所要分析的范围进行范围的选取，生成.kml文件，然后在global mapper上将kml文件转换成.shp文件。 点击/SARscape/General Tools/Sample Selections/Sample Selection SAR Geometry Data工具，打开Sample Selection SAR Geometry Dat面板。 选择需要裁剪的两景数据的VV _slc _list文件。

• 4.基线估算

打开具/SARscape/Interferometry/InterferometricTools/Baseline Estimation,输入主从影像的_slc _list文件，点击Exec执行，输出基线估算的结果。

• 5.DI nSAR工作流

打开/SARscape/Interferometry/DInSAR Displacement Workflow工具。
干涉图生成
这一步生成的结果有: INTERF_ out_ dint:去平后的干涉图 INTERF_ out_ int:干涉图 INTERF_ out_master_pwr:主影像强度图 INTERF_out_ slave_pwr:从影像强度图 INTERF_out_par:文本文件保存的配准时的偏移量数据， INTERF out sint:合成的相位 INTERF_ out_ srdem:斜距几何下的参考DEM INTERF_out_par orbit_ off:估算轨道偏移时用到的点矢量数据，包括在方位向和距离 向的点的位置坐标、测量到的偏移量、计算出的线性偏移量。 INTERF_out_par winCC_off:从强度数据的相差上估算交叉相关偏移量的点矢量数据。 包含每个点上交叉相关的值(CC)，范围是0-10 INTERF_out_par winCoh_off:从相位数据的相差上估算相干性的点矢量数据，包含信 噪比(SNR)和相干性，相干性值的范围是0-1
滤波和相干性计算(Adaptive Filter and Coherenace Generation）
点击Next按钮，进行干涉图滤波和相干性生成处理，处理完成之后，自动加载滤波后的干涉图_fi nt和相干性系数图_cc。这一步处理之后生成的结果有: INTERF_ out_ fint:滤波后的干涉图 INTERF out cc:相干性系数图
相位解缠(Phase Unwrapping)##
点击N ext按钮，进行干涉图滤波和相干性生成处理，处理完成之后，自动加载滤波后的相 位解缠结果图_upha。这一步处理之后生成的结果有: INTERF_out upha:相位解缠结果
控制点选择(GCP Selection )
在控制点生成面板上，点击Next，打开控制点选择工具，鼠标变为选点状态，单击鼠标 左键就可以选点。在图像周边远离形变区域的地方、相位好的地方选择若干控制点。
轨道精炼和重去平(Refinement and Reflattening)
除可能的斜坡相位，对卫星轨道和相位偏移进行纠正。进行轨道精炼和相位偏移的计算，消这一步对解缠后的相位是否能正确转化为形变值很关键。
一步得到的结果有: INTERF out_ reflat_ fint:重去平后的干涉图 INTERF out reflat sint:重去平后的合成相位图 INTERF_out_reflat_upha:重去平后的解缠结果 INTERF out reflat srdem:重去平后的斜距DEM INTERF_ out_ reflat.txt:轨道精炼所用的轨道修正参数 INTERF_out_reflat_refinement.shp:轨道精炼用到的有效的控制点文件，斜距坐标系。 INTERF_out_reflat_refinement_geo.shp:轨道精炼用到的有效控制点文件，地理坐标系。
相位转形变以及地理编码(Phase to Displacement Conversion and Geocoding)
将经过绝对校准和解缠的相位，结合合成相位，转换为形变数据以及地理编码到制图坐标系统，默认得到的是LOS方向的形变。
地理编码参数(Geocoding
点击Next按钮，进行相位转形变和地理编码处理。地理编码的坐标系是以参考DEM的坐标 系为准。这一步得到的结果有: _ slc_out_disp_dem:重采样到制图输出分辨率上的参考DEM数据，经过地理编码，范 围和输出的SA R产品一致。 _ slc_out_disp_ cc_geo:地理编码的相干性系数图 _ slc_out_disp_precision:数据质量的估算结果图，代表形变的精度。 _ slc_out_disp_ ALOS:视线方位角，正值是正北的顺时针方向，负值是正北的逆时针方向。 _ slc_out_d isp_I LOS:视线入射角，视线和水平面垂线的夹角。 _ slc_out_disp:传感器观测方向的形变，即LO S方向上的形变。

• 结果输出(Output )

结果默认输出在ENVI的默认输出路径下，文件名命名为sentinell_ 69_20160109_100020274_1W_ SIW1_ A_ VV_slc_list_ out_disp。若想保留中间结果便于查看点击Finish,将Delete Temporary Files不勾选。输出结果。结束DInSAR Diaplcement处理的工作流。生产的形变数据自动进行 密度分割配色展示。

• 结果分析

在结果密度分割图层上Slices右键选择Export Color Slices->Class Image，将密度分割结果
1
保存为分类文件。 在工具栏中，Annotations下拉框中提供很多的注记要素，包括图例、指北针、比例尺、 公里网，以及文字、图形、图像等元素，对结果进行制图输出。如下为简单的分类结果制图 效果。