前言
事先声明,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.校正后结果