使用GDAL对静止卫星圆盘数据进行校正（以FY2为例子）

使用GDAL对静止卫星圆盘数据进行校正（以FY2为例子）

 

前言

使用GDAL对静止卫星数据的校正的方式与之前写的校正MODIS数据一样，只不过对于静止卫星的圆盘数据而言，经纬度查找是固定的，不会发生变化。

经纬度查找表数据

对于FY2气象卫星而言，经纬度查找表文件（NOM_ITG_2288_2288(0E0N)_LE.zip）可以从网上下载。数据解压之后，里面有三个文件，分别是：

序号	名称	说明
1	NOM_ITG_2288_2288(0E0N)_LE.dat	经纬度查找表数据
2	NOM_ITG_2288_2288(0E0N)_LE.hdr	经纬度查找表数据头文件
3	Readme for FY latlon.txt	数据说明文件

在数据说明文件中，详细说明了经纬度查找表数据的使用方法，同时针对FY2系列卫星的经度进行了说明，由于上述经纬度数据是以中心点为0度经度生成的，所以对于FY2的数据，需要在经度数据上加上卫星所在的经度。

序号	卫星名称	星下点经度
1	FY-2F	112°E
2	FY-2E	104.5°E
3	FY-2D	86.5°E
4	FY-2C	104.5°E

影像校正

将上面的经纬度查找表数据按照博客的内容构造一个Geolocation信息，然后调用gdalwarp就可以处理。 
需要注意的是，需要事先对经纬度数据中的经度全部加上卫星的星下点经度，然后再构造GeoLocation信息。也可以不用这么处理，直接校正，校正结束之后，修改图像的仿射变换参数，先通过GDALDataset::GetGeotransform函数获取仿射变换参数，然后给第一个值加上卫星的星下点经度之后，再使用GDALDataset::SetGeotransform函数更新到校正结果数据中即可。

使用GDAL工具校正

具体的处理方式参考《使用GDAL工具对FY3系列卫星数据进行校正》以及《使用GDAL对HDF数据进行校正》两篇博客。

自己编写代码校正

具体处理步骤参考《使用GDAL对HDF数据进行geoloc校正》的相关内容。

示例数据