静止轨道卫星几何校正——GOES-16几何校正

前言

事先声明，GOES-16的几何校正是由标称投影转换为等经纬度投影，我没找到经纬度查找表，因此，采用了python GDAL尝试校正，没想到可以成功转换，也能满足应用类研究的需求，如果你需要精确校正，请谨慎使用！！！

认识GOES-16数据

这个数据是按照波段分开的，每一种波段都有各自的空间分辨率，代码运行的效率也随着波段不同而不同。如下图所示：

校正方法

1.使用下面的代码，python3 GDAL

# _*_ coding: utf-8 _*_
#Author:QHWGIS1997
#Date:2020/11/26

import gdal,numpy,time,os
start=time.clock()
GOES_16_DATA=gdal.Open(r"E:\AIRCAS\GOES\JANUARY\01\data\ABI-L1b-RadF_2019_001_16_OR_ABI-L1b-RadF-M3C01_G16_s20190011600366_e20190011611133_c20190011611180.nc")
41175.nc")#打开.nc文件

datasets=GOES_16_DATA.GetSubDatasets()


gdal.Warp('E:/AIRCAS/GOES/JANUARY/01/band1/201901011610.tif',datasets[0][0],dstSRS='EPSG:4326')#进行等经纬度投影

#关闭数据集
root_ds=None
end = time.clock()
print('Running time: %s Seconds'%(end-start))

2.校正后结果