gamma预处理的StaMPS-PSI

具体可参考B站尼莫大佬的视频StaMPS_PS时间序列分析_哔哩哔哩_bilibili

一、新建文件夹INSAR_20200317

./rslc/yyyymmdd.rslc

./rslc/yyyymmdd.rslc.par

./diff0/yyyymmdd_yyyymmdd.diff

./diff0/yyyymmdd_yyyymmdd.base

./geo/yyyymmdd_dem.rdc

./geo/yyyymmdd.lat

./geo/yyyymmdd.lon

二、mt_prep_gamma 20200317 '/media/liuan/LADATA/INSAR_20200317' 0.4 1 1 50 50

三、stamps(1,1)

help ps_plot

ps_plot('w') the wrapped phase of the selected initial PS candidate pixels

ps_plot('w',5,0,0,1) choose the second ifg with the mean amplitude of all the images as background (value 5).

四、设置参数

setparm('density_rand',1)

setparm('weed_standard_dev',0.6)

setparm('weed_time_win',365)

setparm('merge_resample_size',100)

setparm('unwrap_grid_size',100)

setparm('unwrap_time_win',365)

setparm('weed_neighbours','y')

五、cd PATCH_1

stamps(2,2)

stamps(3,3)

stamps(4,4)

stamps(5,5)

六、cd ..

从patch list中删除patch_1

stamps(1,4)

在patch list中添加patch_1

stamps(5,5)

ps_plot(‘w’) plot the phase of the selected pixels for the full study area 干涉图

ps_plot(‘w’,5,0,0,1)干涉图

七、去大气相位

GACOS Binary grid格式 一次最多输入20个日期

将文件.ztd/.ztd.jpg/.ztd.rsc放入./APS/yyyymmdd文件夹中

getparm_aps;

load('parms.mat')

save('parms_aps.mat','heading','lambda','-append');

setparm_aps('UTC_sat','22:12');

setparm_aps('stamps_processed','y');

setparm_aps('sounding_dir','/media/liuan/LADATA/INSAR_20200317/sounding_data');

setparm_aps('ll_matfile','/media/liuan/LADATA/INSAR_20200317/ps2.mat');

setparm_aps('gacos_datapath','/media/liuan/LADATA/INSAR_20200317/APS');

setparm_aps('demfile','/media/liuan/LADATA/INSAR_20200317/dummy.dem');

setparm_aps('bperp_matfile','/media/liuan/LADATA/INSAR_20200317/bp2.mat');

setparm_aps('hgt_matfile','/media/liuan/LADATA/INSAR_20200317/hgt2.mat');

setparm_aps('ifgday_matfile','/media/liuan/LADATA/INSAR_20200317/ps2.mat');

setparm_aps('phuw_matfile','/media/liuan/LADATA/INSAR_20200317/phuw2.mat');

setparm_aps('look_angle','/media/liuan/DATA/ds_20201217/la2.mat');

aps_weather_model('gacos',1,3)

八、解缠

setparm('tropo','a_gacos');

setparm('subtr_tropo','y');

stamps(6,7);

stamps(6,6);

setparm('scla_deramp','y');

stamps(7,7);

stamps(6,7);

九、绘制形变速率图

ps_plot('v-dao','a_gacos',1,0,0,'ts')