基于GLT方法的风云三号C（FY-3C）卫星影像几何校正

ENVI5.1版本开始直接支持HDF5格式的读取，可以直接处理风云三号系列气象卫星HDF5格式的影像处理，风云三号A卫星处理可参考：ENVI5.1下基于GLT方法的风云三号卫星影像几何校正。风云三号C卫星影像与风云三号A数据有点差别，本文以风云三号C的VIRR数据为例介绍操作过程。

软件版本：ENVI5.5.1，ENVI5.1及更高版本类似。

第一步、数据打开

风云三号C数据以HDF5格式存储，包括一个数据文件和一个几何定位文件。

（1） File > Open As >Scientific Formats >HDF5，选择风云3的数据文件：FY3C_VIRRX_GBAL_L1_20170101_0235_1000M_MS.HDF文件，进入Select HDF5 Datasets面板。

图1 HDF5 Datasets面板

（2）文件中包含很多信息，选择图像数据EV_RefSB，点击中间的加载箭头，将所选数据加载到右边的列表中，单击右下角Open Rasters，打开图像。

图2 FY-3 EV_RefSB数据（假彩色显示）

（3）同样的方法，点击File > Open As >Scientific Formats >HDF5，选择FY-3的几何定位文件：FY3C_VIRRX_GBAL_L1_20170101_0235_GEOXX_MS.HDF文件，进入Select HDF5 Datasets面板相同方法，同时将定位文件（Latitude和Longitude）加载并在ENVI中打开。

图3 通过HDF Datasets加载FY-3经纬度数据

图4 FY-3 经纬度数据

（4）部分FY-3几何定位数据，即经纬度数据存在无效值，如无效值为-999.900024表示。对于存在无效值的经纬度数据，在构建GLT时会出现异常，因此需要对经纬度数据的无效值进行修复。

注：对于无无效值的经纬度数据，跳过此步骤。

图5 FY-3 经纬度数据无效值查看

无效值修复方法为：

1) 在ToolBox中，选择/Terrain/Replace Bad Values

2) 在弹出的对话框中选择FY3C_VIRRX_GBAL_L1_20170101_0235_GEOXX_MS.HDF文件，点击OK。

3) 点击OK后弹出参数面板，Max填入-200，设置输出文件名称，其它参数保持默认，点击OK执行。

注：min不填自动读取图像最小值，Max大于经纬度最小-180即可

图6 FY-3 经纬度数据无效值修复

第二步、几何定位

（1）在Toolbox中，选择/Geometric Correction/Build GLT。

（2）在弹出的对话框中，Input X Geometry Band选择修复后的经度longitude波段， Input Y Geometry选择修复后的纬度latitude波段。在接下来弹出的对话框中填写输出GLT文件的投影信息，这里选择Geographic Lat/Lon，WGS-84，如下图所示。