StaMPS中使用GACOS数据进行大气校正，解缠、计算形变速率

1 GACOS数据下载

GACOS官网

• Time of insterest（in UTC）时间设置
  这个在数据名称里面可以看到，如果数据的条带一致，那么每日拍摄时间应该也是一样的，如下图红框所示，T后面的数字就是HHMMSS，即10时18分43秒
  

• 选择二进制文件Binary grid
  

• 提交之后，如果没有出错，会弹出这个界面
  

• 然后等待邮件即可，邮箱里面有下载链接，一定要及时下载，数据只会在服务器上保存72小时，一般要几分钟至几个小时，取决于提交的数据空间范围和时间跨度

• 下载好了之后解压可以看到如下内容
  

2 查看GACOS用法

打开压缩包里面的readme.pdf，可以看到GACOS的使用步骤，在stamps中，一般使用前三个

3 在StaMPS中使用GACOS进行大气校正

• 在StaMPS中，做完stamps（1,5）之后，开始做大气校正

1. 查看当前大气校准参数

getparm_aps

2. 加载parms.mat参数配置文件，保存配置参数

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

3. 设置gacos数据路径

setparm_aps('gacos_datapath','/APS')

4. 设置SAR数据成像时间

setparm_aps('UTC_sat','10:18')

5. 再次运行getparm_aps，检查参数是否设置成功

6. 运行aps_weather_model('gacos',1,3)，进行大气校正

大气校正完成后会生成tca2.mat文件

4 stamps6 相位解缠 stamps7时间形变速率

1. 在参数中设置去除大气噪声的算法为gacos，并将去除大气噪声的参数改成y，默认是n

setparm('tropo','a_gacos')
setparm('subtr_tropo','y')

2. 然后做stamps（6,7）

3. 因为上面的stamps（6,7）之后的结果是一次相位解缠和时间形变速率计算的结果，所以去除了一些误差，现在基于上面stamps7的结果重做相位解缠stamps(6,6)，误差就更小了。
4. 设置setparm('scla_deramp','y') 去除deramp相位
5. stamps(7,7)时间形变速率
6. stamps(6,7)再做一遍，效果更平滑

5 制图

ps_plot('v-dao','a_gacos',1,0,0,'ts')
v表示形变速率
d表示dem误差
a表示大气atmosphere误差
o表示轨道orbit误差
ts表示time series时间序列