基于GMTSAR软件的富士山时序形变探测（二）

本篇从第三章操作内容的第二步开始写，为了更好地衔接上一部分，本篇直接从“2”开始。

2、选取主影像及建立时空关系对

2.1、链接影像

这里需要知道Sentinel-1A数据是由三个IW条带组成，每个条带都会有9个左右的burst，这也是文件夹分成F1、F2、F3的原因，其中数字1、2、3分别对应条带1、2、3，然后我们需要知道自己的研究区落在哪个条带，然后开始链接数据，此次实验主要研究富士山附近，预览图：

通过预览图可以看出最右侧的条带（黄色圆圈）包含富士山，因为是降轨数据，所以最右面为IW1条带，因此我们开始在F1/raw下链接，这里需要链接xml、tiff、EOF、dem四个文件。命令行为ln -s +文件的路径：
ln -s ../../data/*.SAFE/*/*iw1*vv*xml
ln -s ../../data/*.SAFE/*/*iw1*vv*tiff
ln -s ../../orbit/*EOF
ln -s ../topo/dem.grd
以此次实验为例，链接tiff为：

其余的文件和上面类似，这里不做赘述，当然你可以写个小脚本自动链接，小编这里没有实现，大神可以自己编写，全部链接完成后如图所示。这是F1文件下的，也就是IW1条带的，如果研究区覆盖两个条带或者三个条带，那么需要重复此操作在F2、F3中依次进行，后续干涉完成后需要合并。

2.2、建立影像精轨对应文件data.in

链接完文件后，我们需要建立一个“data.in”文件，里面放置数据名和精密轨道名对应起来，原说明书里也没有给自动化的脚本，希望有后续大神可以补充一下这方便的便携性。在原说明书中说有这样一个命令prep_data.csh可以帮助做到这些，但是这个命令只能生成数据目录和精密轨道目录，但是还是得自己手动对应起来，但是需要联网才能生成精密轨道目录，而且是好多年的都生成，并不是对应你的Sentinel-1A数据进行生成，我也不太理解为啥，运行完prep_data.cs后如图所示。

其中text.dat里是你的影像数据，orbit.list里是所有的精密轨道数据，我建议还是复制精密轨道文件的名字对应上，对应完成后如图所示：

2.3、选取主影像

这里需要在raw文件下输入命令行preproc_batch_tops.csh data.in dem.grd 1，生成一个时空极限列表和图像，所有生成文件如下图所示，然后运行mv baseline_table.dat ../将baseline_table.dat文件复制到上一目录F1下。

时空基线图像如下图所示，我们需要选取一个主影像master，主影像的位置大约在时空基线的中间，此次实验S1_20210427_ALL_F1就比较符合要求。

然后你需要将你所选的主影像，此次实验为放S1_20210427_ALL_F1在data.in文件的最上方，因为系统会默认为第一个为主影像，如下图。

然后输入命令行preproc_batch_tops.csh data.in dem.grd 2，再跑一次模式2，该模式为配置工作，如果有F2或者F3需要，要进行同样的操作。
到这里基本上选取主影像和配准工作就完成了，但是要注意在进行2.3步时一定要基于足够大的DEM可以覆盖你一个条带，不然很容易报错，尤其是有几景影像覆盖区差异偏大的情况下。

2.4、建立时空干涉对

这里需要返回F1文件夹，输入时空参数建立干涉对，考虑此次实验情况，设置时间基线200，空间基线100，参数没有固定要求，可按照需求构建，干涉对数量不宜过少，大致为全部数据时间跨度的一半就可以。本文输入命令行select_pairs.csh baseline_table.dat 100 360时空基线图如下：

干涉对列表在文件intf.in里，此次实验生成17干涉对如下：