GMTSAR的D-InSAR基本处理过程（少走弯路）

ytkz2024-08-202024-09-09

之前成功地在linux安装了GMTSAR，我这个linux是指win10系统的WSL。

安装了GMTSAR，就要用起来，现在来讲讲如何使用GMTSAR进行D-Insar处理。

以下均为实际操作。

首先在wsl创建新文件夹insar_test，记住这文件夹，后面一切的数据数据都在这个文件夹，或该文件夹的子文件夹下处理。

mkdir /home/insar_test

进入该文件夹

cd /home/insar_test

1.raw文件夹

在/home/insar_test创建raw文件夹

mkdir /home/insar_test/raw

raw文件夹存放SENTINEL1影像和对应的精准轨道文件。

在ASF官网下载哨兵数据和哥白尼官网下载轨道数据。

ASF下载哨兵数据官网

esa下载数据官网

我的影像数据存放如下：

精准轨道文件的下载可以参考这篇文章：

一定要看一下上面的文章，不然分不清该下载哪个轨道数据。

貌似https://s1qc.asf.alaska.edu/aux_poeorb/的轨道数据没有2023年和2024年，所以我去欧空局下载对应的哨兵1轨道数据。

tips：其实GMTSAR提供了自动下载轨道数据的方法，但是我在测试的时候发现一直不成功，这可能是以前欧空局不需要登陆账号可以下载数据，但现在下载数据需要登陆账号。所以为了稳定使用GMTSAR，建议先手动下载轨道数据，等后面有批量下载的需求再研究如何自动下载。

2.topo文件

在/home/insar_test创建raw文件夹

mkdir /home/insar_test/topo

进入该文件夹

cd /home/insar_test

topo文件存放对应区域的dem数据。已知区域为 1

有两种方法下载。

dem数据下载方法1：

格式: make_dem.csh W E S N [mode]

例如我要下载经度为110 -111，纬度为22-23的区域的SRTM-1s格式的dem数据。

则在wsl中输入：

make_dem.csh 110 111 22 23 1

dem数据下载方法1：

在gmtsar官网下载对应的dem数据（4度以内）

gmtsar官网dem

我的dem数据存放如下：

3.准备配置文件

进入/home/insar_test文件夹

cd /home/insar_test

再输入：

pop_config.csh S1_TOPS > config.s1a.txt

此时在/home/insar_test文件夹中生成一个config.s1a.txt文件，编辑该文件。至于怎么用vim编辑文件大家自己去网上搜，这里不做介绍。

vim config.s1a.txt

修改config.s1a.txt中的threshold_snaphu = 0，把0改为0.1

4.运行p2p_S1_TOPS_Frame.csh命令

此时说明一下数据准备的情况。

在/home/insar_test有一个文件config.s1a.txt

在/home/insar_test有两个文件夹raw、topo

在/home/insar_test/raw文件夹有两景S1影像和对应的轨道文件

在/home/insar_test/topo文件夹有一景dem影像

在终端进入/home/insar_test

cd /home/insar_test

在终端运行这个命令：

p2p_S1_TOPS_Frame.csh S1A_IW_SLC__1SDV_20231005T105411_20231005T105433_050629_06196B_87AB.SAFE S1A_OPER_AUX_POEORB_OPOD_20231025T080724_V20231004T225942_20231006T005942.EOF S1A_IW_SLC__1SDV_20240109T105408_20240109T105430_052029_06499F_4C45.SAFE S1A_OPER_AUX_POEORB_OPOD_20240129T070738_V20240108T225942_20240110T005942.EOF config.s1a.txt vv 1

上面命令的格式说明如下：

Usage: p2p_S1_TOPS_Frame.csh Master.SAFE Master.EOF Aligned.SAFE Aligned.EOF config.s1a.txt polarization parallel

Example: p2p_S1_TOPS_Frame.csh S1A_IW_SLC__1SDV_20150607T014936_20150607T015003_006261_00832E_3626.SAFE S1A_OPER_AUX_POEORB_OPOD_20150615T155109_V20150525T225944_20150527T005944.EOF S1A_IW_SLC__1SSV_20150526T014935_20150526T015002_006086_007E23_679A.SAFE S1A_OPER_AUX_POEORB_OPOD_20150627T155155_V20150606T225944_20150608T005944.EOF config.s1a.txt vv 1

    Place the .SAFE file in the raw folder, DEM in the topo folder
    During processing, F1, F2, F3 and merge folder will be generated
    Final results will be placed in the merge folder, with phase
    corr [unwrapped phase].
    polarization = vv vh hh or hv
    parallel = 0-sequential  1-parallel

Reference: Xu, X., Sandwell, D.T., Tymofyeyeva, E., González-Ortega, A. and Tong, X.,
    2017. Tectonic and Anthropogenic Deformation at the Cerro Prieto Geothermal
    Step-Over Revealed by Sentinel-1A InSAR. IEEE Transactions on Geoscience and Remote Sensing.