Generic Atmospheric Correction Online Service for InSAR (GACOS) 是纽卡斯尔大学Li Zhenhong 教授团队提供大气改正数据,数据可以在网站 http://ceg-research.ncl.ac.uk/v2/gacos/上进行申请,时间分辨率可以精确到分钟,空间分辨率为90m。
申请到数据之后可以利用该数据对InSAR的干涉图进行大气改正,网站上也有相应的matlab示例代码。写在这里给自己做个备忘。
1. 要利用这个大气改正数据需要将干涉图进行地理编码,
2.将干涉图的参数文件(示例代码里的phs.rsc)进行指定.
WIDTH 4940
FILE_LENGTH 2111
XMIN 1
XMAX 4940
YMIN 1
YMAX 2111
X_FIRST -2.5700000
Y_FIRST 52.9066667
X_STEP 8.3333333e-04
Y_STEP -8.3333333e-04
X_UNIT degres
Y_UNIT degres
Z_OFFSET 0
Z_SCALE 1
PROJECTION LATLON
DATUM WGS84
3.剩下就是调用make_correction相应的操作了
注意:
1. 个人认为这个流程比较适合InSAR产品级别的大气相位改正,对于需要在雷达坐标系下进行大气相位去除的小伙伴们就需要自己去对这个大气相位进行地理编码
2. 网站上申请数据那个时间应该是UTC时间,我觉得这个应该说明一下。