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

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

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值