StaMPS结合TRIAN 做GACOS大气校正

StaMPS使用GACOS做大气校正

1. 软件安装

1.1 StaMPS的下载及安装.
可参考以前发表的文章 ubuntu16.04下安装ROI_PAC，ubuntu 16.04 安装StaMPS 及常见问题解析，StaMPS的step_read_whole_ERS，step_read_whole_RSAT，step_read_whole_TSX等出现DORIS_SCR: Undefined variable等。

1.2 TRAIN的下载及安装。
软件下载地址：https://github.com/dbekaert/TRAIN。
下载完成解压后可直接在matlab里设置函数调用地址，即可使用。
具体操作如下图，在matlab里点击SetPath，然后点击Add folder添加路径，将TRAIN软件的matlab函数路径~TRAIN-master/matlab复制粘贴即可。

2. GACOS数据下载及使用

2.1 GACOS数据下载网址：http://www.gacos.net。在网页上填写数据四至的坐标、UTC时间以及每景影像获取的时间（包括主影像）。Output format一定要选择Binary grid！Output format一定要选择Binary grid！Output format一定要选择Binary grid！ 否则GACOS无法被找到并正常使用。上述参数可以在matlab中运行aps_weather_model('gacos',0,0)获取（一般是在执行完StaMPS第五步stamps(1,5)后开始去大气，此获取gacos下载参数的命令也是要在第五步后，否则参数读取错误），正确示例结果如下图。输入完参数后填写自己的邮件地址，等几分钟后就会收到数据下载链接。

2.2 GACOS数据及其存放位置

下载完成后的数据如下图所示，直接复制到自己指定的新建文件夹即可。

2.3 参数设置。
在matlab中输入getparm_aps即可查看当前去除大气的参数设置。

设置参数（‘heading’, ‘lambda’, ‘UTC_sat’, 'gacos_datapath’等），先用getparm_aps检查参数设置是否正确，然后用setparm_aps将参数进行修改。以我处理的TerraSAR为例，只需要修改’gacos_datapath’的路径即可，其它参数默认都是正确的。

2.4 进行大气校正
大气校正要在StaMPS的第五步之后，即stamps(1,5)之后，使用命令 aps_weather_model('gacos',0,0) 可获得gaocos数据下载的相关参数。使用命令 aps_weather_model('gacos',1,3) 即可在2.3中参数设置正确的情况下进行大气校正。

2.5 后续StaMPS处理
首先要将StaMPS中大气相关的参数改正。执行以下两个命令，setparm('tropo','a_gacos');，及setparm('subtr_tropo','y');,然后进行相位解缠等后续操作，即stamps(6,7);。可在此基础上重复一次相位解缠stamps(6,6);,然后setparm('scla_deramp','y'); stamps(7,7); stamps(6,7);即可。

2.6 绘制结果图
ps_plot('w'); %缠绕相位图
ps_plot('w-a','a_gacos'); %PS缠绕相位图减去大气相位
ps_plot('u'); % 解缠相位图
ps_plot('u-a','a_gacos'); %PS解缠相位图减去大气相位
ps_plot('a','a_gacos'); %PS大气相位图
ps_plot('w-asb','a_gacos'); %SBAS缠绕相位图减去大气相位
ps_plot('u-asb','a_gacos'); %SBAS解缠相位图减去大气相位
ps_plot('asb','a_gacos'); %SBAS大气相位图