记录一下StamPS+PS的过程

1. 下载数据：alaska.edu

2. 利用SNAP处理数据

流程：split->apply orbit->back-geocoding - >deburst->subset->interferogram->stamps export

3. Split，选择IW1，VV极化

4. apply orbit，应用精密轨道

5. back-geocoding，用stack overview确定主影像，以其为准,每次配准5副影像，共做5次，注意每次配准主影像放置最高

6. deburst，去除条带

7. subset，打开deburst后的一个文件，右键subset即可，注意需保存save

8. interferogram，对sub文件生成干涉图

9. stamps格式导出，subset在上，干涉图fig在下

10. 每个文件包含dem、diff0、geo、rslc四个文件夹，整理成一个文件

11. 在Ubuntu18中配置好stamps、matlab等环境

12. 在matlab中运行脚本1：

% amplitude difference dispersion (0.4 for PS and 0.6 for SBAS)

ADD_PS=0.4;

% number of patches in range and azimuth

NPatchesRange=5;

NPatchesAzimuth=5;

% paramters for StaMPS

density_rand=1;

weed_standard_dev=1;

weed_time_win=365;

merge_resample_size=100;

unwrap_grid_size=100;

unwrap_time_win=365;

masterdate= '20210825 ';

% prepare data

cmd=['mt_prep_snap ' masterdate pwd '/Export ' num2str(ADD_PS) ' ' num2str(NPatchesRange) ' ' num2str(NPatchesAzimuth) ' 50 50' ];

system(cmd);

stamps(1,1);

% set some parameters

setparm('density_rand',density_rand);

setparm('weed_standard_dev',weed_standard_dev);

setparm('weed_time_win',weed_time_win);

setparm('merge_resample_size',merge_resample_size);

setparm('unwrap_grid_size',unwrap_grid_size);

setparm('unwrap_time_win',unwrap_time_win);

setparm('insar_processor','snap');

13. 进入每个生成的patch文件夹，运行：stamps(1,5)

14. 返回上一级目录，运行：stamps(5,5)

15. 去除大气相位：GACOS，设置经纬度、时间等，提交后邮箱收到文件，整理成APS文件夹：

16. 将APS放置patch同一级目录下，运行脚本2：

stamps(5,5)

% remove atmospheric phase by gacos

getparm_aps;

load('parms.mat');

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

setparm_aps('gacos_datapath','./APS');

setparm_aps('UTC_sat','10:29');

aps_weather_model('gacos',1,3);

% phase unwrapping and time series analysis

setparm('tropo','a_gacos');

setparm('subtr_tropo','y');

stamps(6,6);

stamps(7,7);

%stamps(6,7);

% re-do phase unwrapping

unwrap_grid_size=getparm('unwrap_grid_size');

setparm('unwrap_method','3D');

setparm('unwrap_hold_good_values','y');

setparm('unwrap_grid_size',unwrap_grid_size+10);

stamps(6,6)

setparm('scla_deramp','y');

stamps(7,7)

stamps(6,7);

17. 绘制时间序列形变图：ps_plot('V-dao','a_gacos',1,0,0,'ts')

主要参考up主尼莫：

https://www.bilibili.com/video/BV1vv411i7se?p=2&vd_source=5897ef2245bbf96656b07e9ef7ba73e6