一、数据下载
1、哨兵1数据下载
2、轨道数据下载
SLC数据与轨道数据下载见博文:InSAR个人笔记(三)ISCE之DInSAR处理--topsApp.py_玊363的博客-CSDN博客
3、dem数据下载
Dem数据需转换为grd格式,名为dem.grd
GMTSAR官网下载:https://topex.ucsd.edu/gmtsar/demgen/
命令行下载:make_dem.csh W E S N [mode] mode 1:SRTM-1s 2:SRTM-3s)
二、预准备工作
1、raw文件夹
放置两个SAR数据与轨道数据
2、topo文件夹
放置dem数据
3、config.s1a.txt文件
可使用pop_config.csh生成卫星的参数文件,具体说明:
https://gmtsar.github.io/pop_config.csh.html
(1)功能:生成一个配置文件,其中包含每个卫星的默认参数。
(2)输入:卫星 ERS, ENVI, ALOS, ALOS_SLC, ALOS2, ALOS2_SCAN, S1_STRIP, S1_TOPS, CSK_RAW, CSK_SLC, TSX, RS2
(3)输出: 干涉图处理的配置文件 config.s1a.txt
三、p2p_S1_TOPS_Frame.csh执行
具体说明:p2p_S1_TOPS_Frame.csh — GMTSAR 6.2 documentation
(1)功能:用于从Sentinel-1 TOPS数据生成干涉图的处理脚本
注意:scansar-扫描模式是共享多个独立sub-swaths的操作时间,最后获取一个完整的图像覆盖区域。它能解决Stripmap模式较小的刈幅
(2)输入:主影像 轨道文件 待配准影像 轨道文件 config.s1a.txt vv 0
0: 将按顺序处理subswaths,有利于内存较少的计算机,但执行时间更长。
vv 指定数据处理的极化方式。4种极化方式vv,vh,hv,hh
1: subswaths并行处理
config.s1a.txt 干涉图处理的配置文件。可以使用生成默认文pop_config.csh
(3)输出:F1, F2, and F3文件夹,每一个文件夹中包含一个subswath处理文件
(4)示例:p2p_S1_TOPS_Frame.csh S1A_IW_SLC__1SDV_20180317T235700_20180317T235727_021062_0242D6_6F45.SAFE S1A_OPER_AUX_POEORB_OPOD_20180406T120838_V20180316T225942_20180318T005942.EOF S1A_IW_SLC__1SDV_20180329T235700_20180329T235727_021237_024860_7623.SAFE S1A_OPER_AUX_RESORB_OPOD_20180330T044439_V20180329T235400_20180330T031129.EOF config.s1a.txt vv 1
最终结果生成的文件夹如图所示:
其中merge文件中包含的文件有: