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

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

在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，将密度分割结果

保存为分类文件。
在工具栏中，Annotations下拉框中提供很多的注记要素，包括图例、指北针、比例尺、
公里网，以及文字、图形、图像等元素，对结果进行制图输出。如下为简单的分类结果制图
效果。