GMTSAR时序InSAR处理流程（Sentinel-1A）（后续一直更新）

##预处理

cd data

ls -d $PWD/*.SAFE > SAFE_filelist

less SAFE_filtelist   #检查文件是否为SAFE，路径是否正确

ls -d *EOF > orbits.list

download_sentinel_orbits.csh SAFE_filelist 1>& dso_mode1.log &  #下载轨道数据

make_dem.csh   #下载DEM

#创建pins.ll文件,选择研究区所在的burst。IW:升轨：从左到右IW1、IW2、IW3,降轨反之

pins.ll格式

升轨：E  S

          W  N

降轨：E  N

           W  S

organize_files_tops_linux.csh SAFE_filelist ../reframed/pins.ll 1 >& oft_mode1.log &  

organize_files_tops_linux.csh SAFE_filelist ../reframed/pins.ll 2 >& oft_mode2.log &

#查看log屏幕输出

tail -f oft_mode2.log

#创建超链接文件

ln -s ../../topo/dem.grd .

ln -s ../../data/F*_F*/*.SAFE/*/*iw1*vv*xml .

ln -s ../../data/F*_F*/*.SAFE/*/*iw1*vv*tiff .

ln -s ../../data/*.EOF .

ln -s ../topo/dem.grd .

#建data.in文件

ls *.EOF > orbits.list

prep_data_linux.csh

##粗配准，在F1/raw目录下输入命令

preproc_batch_tops.csh data.in dem.grd 1

##干涉

#在F1路径下，生成基线图

select_pairs.csh baseline_table.dat 60 200

#查看多少个基线对

wc -l intf.in

cp intf.in ../F2(F3)

vi intf.in

#在F2路径下修改intf.in

:%s/F1/F2/g (将F1替换为F2)

#设置batch_tops.config参数以及主影像

#测试、查看一个干涉对

head -1 intf.in > one.in

intf_tops.csh one.in batch_tops.config

# 将batch_tops.config中的Pro_stage修改为2运行所有干涉对

intf_tops_parallel.csh intf.in batch_tops.config 6 >& ifp.log &

##合并子条带 （merge_batch.csh）

#准备文件

ls -d 20* >intflist(在F*/intf_all下)

cp ../F1/batch_tops.config .

ln -s ../topo/dem.grd .

#创建merge_list

create_merge_input.csh intflist .. 2 > mere_list

#合并

merge_batch.csh merge_list batch_tops.config（注：再次运行merge_batch.csh生成的trans.dat不会覆盖，需手动删除）

##相位解缠

#掩膜低相干地区: landmask.csh

cd 到合并的干涉图目录之一,查看minX/maxX/minY/maxY

gmt grdinfo phasefilt.grd

#merge下创建landmask_ra.grd文件

landmask.csh minX/maxX/minY/maxY

 (初次写，记录实验过程，后续更新并完善！！！)