GMTSAR软件InSAR时序处理流程

菲林可乐

已于 2024-01-06 11:27:45 修改

阅读量1.1k

点赞数 4

文章标签： linux

于 2023-11-07 13:04:34 首次发布

本文链接：https://blog.csdn.net/m0_60086135/article/details/134264835

版权

1.Download InSAR data

 nohup python3 ./download-all-2023-11-01_12-08-11.py > T85_download.log &

2.generated SAFE filelist

ls -d /home/wx/gmtsar_sbas_test/data/*SAFE > filelist 

sed 's|^/data2/xinw/HMF_Sentinel1_data/T85/T85_SAFE/||' filelist | sort -t '_' -k 7,7 -k 8,8n -k 9,9n | sed 's|^|/data2/xinw/HMF_Sentinel1_data/T85/T85_SAFE/|'

sed 's|^/data2/xinw/HMF_Sentinel1_data/T85/T85_SAFE/||' filelist | sort -t '_' -k 7,7 -k 8,8n -k 9,9n | sed 's|^|/data2/xinw/HMF_Sentinel1_data/T85/T85_SAFE/|' > filelist_sort

sort -t '_' -k 7,7 -k 8,8n -k 9,9n filelist > filelist_sort

3.Download orbit

download_sentinel_orbits_linux.csh filelist 1

filelist:
absolutepathto/filename1.SAFE
absolutepathto/filename2.SAFE

     mode 1 = precise orbits (POEORB)
            (most users should choose precise orbits)
    mode 2 = temporary (restituted) orbits (RESORB)

4.Download DEM

make_dem.csh -156.1 -154.8 18.8 19.7

5.拼接burst

organize_files_tops_linux.csh filelist pins.ll 1 看看多少数据覆盖

organize_files_tops_linux.csh filelist pins.ll 2

 ls -lht F0301_F0304/*/measurement

ls *SSV_20151006*/preview 

scp -r  xinw@114.214:/data2/xinw/HMF_Sentinel1_data/T85/T85_SAFE/S1A_IW_SLC__1SSV_20151006T122221_20151006T122248_008032_00B3EF_006D.SAFE/preview ./preview3

6. 超链接

~/gmtsar_sbas_test/F2/raw$ ln -s  ../../organized/F0301_F0304/*/*/*iw2*vv*tiff .
~/gmtsar_sbas_test/F2/raw$ ln -s  ../../organized/F0301_F0304/*/*/*iw2*vv*xml

ln -s ../topo/dem.grd  .
ln -s ../../organized/*EOF .

7. data.in

prep_data_linux.csh

8.baselines

preproc_batch_tops.csh  data.in dem.grd  1
gmt psconvert baseline.ps -Tf -A

preproc_batch_tops.csh  data.in dem.grd  2

9. intf.in

 cp /usr/local/GMTSAR/gmtsar/csh/batch_tops.config .
 
vi batch_tops.config

master_image = S1_20180514_ALL_F2
threshold_geocode = 0


cp raw/baseline_table.dat  .

ls *ALL*PRM > prmlist
get_baseline_table.csh prmlist S1_20180514_ALL_F2.PRM
cp baseline_table.dat ../


select_pairs.csh baseline_table.dat 90 150

10.

ls ../F2/intf_all/ > intflist

create_merge_input.csh intflist .. 2 > mergelist

merge_batch.csh  mergelist  batch_tops.config

 ls 20*/corr.grd > grid_list

nohup stack_coherence_mask.csh grid_list  0.075 >& stack_coherence_mask.log &

gmt grdinfo 2021069_2021117/phasefilt.grd -I-
landmask.csh 0/68376/0/12660
gmt grdsample landmask_ra.grd -R2020003_2020039/phasefilt.grd -Gtmp.grd
mv tmp.grd landmask_ra.grd

find . -type d -exec sh -c 'test -e "$0/unwrap.grd" || basename $0' {} \;

sudo cp unwrap_intf.csh /usr/local/GMTSAR/bin/

rep_sbas.csh intf.in baseline_table.dat ../merge unwrap.grd corr.grd

sbas intf.tab scene.tab 3 3 5963 2067 -smooth 1.0 -wavelength 0.0554658 -incidence 40 -range 892000 -rms -dem

proj_ra2ll.csh trans.dat vel.grd vel_ll.grd 
gmt makecpt -Cpolar -T-7000/7000/1000 -Z -D > vel.cpt
grd2kml.csh vel_ll vel.cpt

gmt makecpt -Cjet -T-40/40/5 -Z -D > vel_ll.cpt

gmt grdmask mask.txt -Rlos_ll.grd -N1/NaN/NaN -Gmask.grd

gmt grdmath los_ll.grd  mask.grd MUL = los_ll_mask.grd