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

3、干涉

3.1、进行1个干涉对干涉处理

我们首先进入F1文件下，设置batch_tops.config文件，主要参数如下：
（1）master_image = S1_20210427_ALL_F1
（2）Proc_stage = 1
（3）shift_topo = 0
（4）filter_wavelength = 200
（5）range_dec = 8
（6）azimuth_dec = 2
（7）threshold_snaphu = 0
（8）threshold_geocode = 0
（9）region_cut = 0/10000/10000/13400
为了减少工作量，这里可以进行研究区裁剪，范围如（9）所示，这里用到的范围是雷达坐标系中的Range和Azimuth，因为我不知道如何直接看SLC数据的范围，我先进行了一个干涉对的全F1条带干涉，然后使用结果图的范围来确定的裁剪范围。
使用命令head -1 intf.in > one.in选取一个干涉对单独放置到one.in文件下，然后使用命令intf_tops.csh one.in batch_tops.config干涉。这一步会在F1文件下的topo文件下生成trans.dat文件，这个后续地理编码的时候会用到，进行一个干涉对大约在10-15分钟左右。

3.2、进行全部干涉对处理

进行完一个干涉对之后，将配置文件batch_tops.config里的proc_stage改为2，然后输入intf_tops.csh intf.in batch_tops.config。如果想要多个进程并行处理可以输入命令intf_tops_parallel.csh intf.in batch_tops.config 3这样就可以三个干涉对同时进行，取决于你的电脑配置，如果提示安装parallel工具，输入命令apt get install parallel就可以。全部跑之后，在intf_all文件下会有生成的干涉对，此次实验有17个干涉对。

每一个干涉对中包括滤波后的干涉图，相干系数图等。如果研究区包括多个条带的话，在F2和F3文件中进行同样的操作。这里进行干涉没有裁剪，可以从看出是整个条带的干涉图，富士山区域大约在0/10000/10000/13400位置处。