SARscape5.6.2中GACOS数据的使用
SARscape5.6.2于2022年5月发布，其中一个新功能是支持干涉处理时使用GACOS数据进行大气校正，本文介绍在干涉处理时使用GACOS数据的方法。
GACOS利用对流层迭代分解(ITD)模型(Yu et al.， 2017)从对流层总延迟中分离分层和湍流信号，并生成高空间分辨率的天顶总延迟图，用于校正InSAR测量和其他应用。
GACOS数据可在http://www.gacos.net/网站免费下载。
在使用SARscape5.6.2进行干涉处理时，可用到GACOS数据的处理有三处，分别是干涉处理通用工具中的大气相位延迟校正工具：/SARscape/Interferometry/Interferometric Tools/Atmospheric Phase Delay Correction、PS处理的干涉工作流步骤：/SARscape/Interferometric Stacking/PS/2 - Interferometric Process、以及SBAS处理的干涉工作流步骤：/SARscape/Interferometric Stacking/SBAS/2 - Interferometric Process。
干涉处理模块中的大气相位延迟校正工具
PS干涉工作流处理时可选GACOS校正
SBAS干涉工作流处理时可选GACOS校正
SARscape中干涉处理时使用GACOS数据，步骤如下：（1）下载相应位置和时相的GACOS文件；（2）在SARscape中导入GACOS文件；（3）在干涉处理时输入GACOS文件。
下文以SBAS干涉工作流处理中使用GACOS数据为例，介绍其使用方法。
所用数据集为Sentinel1A数据，共20时相，具体见下表：
|
数据时相 |
数据时相 |
其他 |
|
20170321 20170414 20170520 20170625 20170719 20170824 20170917 20171023 20171116 20171222 |
20180115 20180220 20180328 20180421 20180515 20180620 20180714 20180819 20180912 20181006 |
入射角43.7° 极化方式VV 参考DEM：SRTM30m 工程区范围：19km*22km 数据量：1.75G |
第一步：下载每一时相的GACOS数据
打开网站http://www.gacos.net/，绘制工程区、设置时间、所有时相、以及接收数据的邮箱，提交数据下载请求。
注：时间信息可以在SAR数据导入之后的.sml文件里<AcquisitionTime>字段获取。
20个时相下载到本地的GACOS数据文件如下：
所有时相的GACOS数据
第二步：在SARscape中导入GACOS数据
打开工具/SARscape/Import Data/Other Format/GACOS，导入所有时相的GACOS数据，各界面数据选择和参数设置如下。
GACOS数据导入工具
数据导入之后可以在ENVI中打开，数据带有地理坐标，可以叠加SAR数据的范围文件，确保GACOS数据范围包含SAR数据范围。
导入之后的GACOS数据在ENVI中打开显示
第三步：在干涉处理时加入GACOS文件
正常进行SBAS工作流处理，在第二步SBAS干涉工作流处理/SARscape/Interferometric Stacking/SBAS/2 - Interferometric Process时，在可选文件中Optional Water Vapour File Lists项手动选择上一步导入后生成的所有时相的GACOS数据，并在干涉处理主参数界面设置大气校正传感器为GACOS。
SBAS干涉工作流中选择GACOS传感器进行大气校正
SBAS干涉工作流处理之后，每一对像对的干涉图快视图存放在路径
_SBAS_processing\interferogram_stacking\interf_tiff下， _fint_ql.tif为滤波后差分干涉图快视图，_clean_fint_ql.tif为进行大气校正之后的滤波后差分干涉图快视图，可以打开进行快视图查看。如下是一对像对的大气校正前后的干涉图。
IS_20170917_m_6_20171023_s_7_fint
IS_20170917_m_6_20171023_s_7_clean_fint
分别对该组数据分别使用GACOS数据和不使用GACOS对做全自动SBAS-InSAR处理，其他参数相同，可以看出，在SBAS最终结果的形变速率的值，以及历史形变量上，有一定差别。
使用GACOS数据和未使用GACOS数据SBAS-InSAR得到的形变速率图对比
同一位置历史形变量对比
在InSAR/DInSAR处理时，也可以使用GACOS数据对干涉图进行大气水汽信息的去除：在做完差分干涉、干涉图滤波和解缠之后，打开工具：/SARscape/Interferometry/Interferometric Tools/Atmospheric Phase Delay Correction，根据要求输入滤波后的干涉图、斜距参考DEM文件、在可选文件中分别输入主辅影像的GACOS水汽数据、解缠结果，点击运行，得到大气相位延迟校正后的干涉图_corr_fint。
干涉处理模块大气相位延迟校正工具
该工具处理得到的结果是进行大气相位延迟校正之后的干涉图，后缀为_corr_fint
原始差分干涉图_fint
大气校正后的差分干涉图_corr_fint
注：本文介绍的工具至少需要SARscape5.6.2版本软件。
