使用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校正》的相关内容。
示例数据
