SARscape5.2中哨兵1A IW数据InSAR处理
本文发表于2016-01-11
哨兵1A数据免费分发已经有一段时间了，SARscape5.2全面支持哨兵1A数据的自动下载、读取、强度数据处理、干涉测量处理。本文以一对哨兵1A IW数据对为数据源，进行InSAR数据处理。
-
数据源：哨兵1A IW SLC数据
-
时相：20150325、20150418
-
入射角：39.58
-
极化方式：VV
-
覆盖区域：拉萨市及北部区域
-
覆盖范围： 183*160平方公里
1、数据准备
（1）在ENVI下设置默认的输出路径，打开File->Preferences，在Directories面板设置数据默认输入输出路径。
ENVI系统参数设置默认的输入输出路径
（2）数据导入：打开工具/SARscape/Import Data/SAR Spaceborne/SENTINEL 1，输入两景数据的元数据文件.safe和轨道参数文件，在本例中数据用于做InSAR，所以下载了相应的轨道文件，在Optional Input Orbit File List输入。
数据导入面板
数据导入后，查看数据覆盖范围是西藏的拉萨市及北部区域，加载arcgis online的底图显示如下：
数据覆盖范围
（3）下载参考DEM文件：打开工具/SARscape/General Tools/Digital Elevation Model Extraction/SRTM-3 Version 2，在Input File输入两个时相的slc_list文件，DEM/Cartographic System设置GEO GLOBLE WGS84，其他参数按照默认。
注：1）这一步也可以在InSAR界面下进行。
2）如果网速许可，最好下载SRTM Version4版本的DEM数据。SRTM Version2版本的数据有空洞，会对最终的结果有所影响。
DEM下载面板
下载到的SRTM DEM数据
（4）设置SARscape系统参数：打开/SARscape/Preferences，选择Load Preferences->SENTINEL_TOPSAR，点击OK。
SARscape系统参数设置
2、InSAR工作流
（1）基线估算：打开工具/SARscape/Interferometry/Interferometric Tools/Baseline Estimation，输入主从影像的_slc_list文件，点击Exec执行，输出基线估算的结果。
基线估算面板
基线估算结果
基线估算的结果显示，这两景数据的空间基线是34.166米，时间基线是24天，一个相位的2π变化代表的模糊高程变化是455.452米。
（2）InSAR工作流：打开/SARscape/Interferometry/InSAR DEM Workflow工具
在Input面板，输入主从影像，参考DEM文件，Parameters按照默认的20。切换到干涉图生成、干涉图滤波和相干性生成、相位解缠面板，参数都按照默认，点击Next>>>按钮，程序自动执行这几步。处理时间约为2小时30分。
数据输入面板
干涉图生成面板
干涉滤波和相干性计算面板
相位解缠面板
经过这几步处理之后，生成干涉图、去平后的干涉图、滤波后的干涉图、相干性图以及相位解缠图。
干涉图_int
去平后的干涉图_dint
滤波后的干涉图_fint
相干性图_cc
（3）选择控制点：干涉工作流的几步处理完成后，界面停留在轨道精炼和重去平这一步，点击Great GCP File，弹出控制点选择的流程化工具，按照默认，点击Next，到选控制点界面，选择相干性高、相位质量好的位置作为控制点，本例中选择的控制点位置如图所示：
选择的控制点
（4）轨道精炼和重去平：用选择的相位好的点作为控制点，进行轨道精炼和重去平，将轨道精炼方法设置为自动Automatic Refinement，其他按照默认，点击Next按钮。轨道精炼和重去平后，输出轨道精炼的结果面板。
轨道精炼和重去平
轨道精炼结果报告
（5）相位转高程：切换到Geocoding选项卡，把Relax Interpolation选项和Dummy Removal选项都改为True，结果会对相干性低的区域做内插处理。
相位转高程面板默认参数
相位转高程过程的地理编码参数
在Output面板，按照默认路径和文件名输出DEM结果，如果要保留中间结果，不选择Delete Temporary Files。点击Next。最终输出的DEM数据如下图所示：
DEM结果
相干性低的区域，主要是水体和雪盖，这些区域生成的DEM不平滑。
注：Sentinel1数据空间基线短，适合用于InSAR地表形变测量，提取DEM效果一般，本文仅说明操作流程。
