详细信息
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,如下图所示。

图7 GLT投影参数设置
(3) 在Build Geometry Lookup File Parameters面板中,填写GLT输出参数,像元大小选择默认,旋转角度(Rotation)为0,即正上方为北,选择保存路径和文件名输出。

图8 构建GLT参数设置

图9 GLT构建结果
(4) 在Toolbox中,选择/Geometric Correction/Georeference from GLT,在弹出的Input Geometry Lookup File对话框中选择GLT文件,在Input Data File对话框中选择待校正文件,点击OK。
(5) 设置输出文件名,其它参数保持默认,点击OK执行。

图10 利用GLT校正FY-3数据
(6) 得到的结果具有地理坐标,对于安装了Google Earth的用户,可以点击菜单栏File>Clip View to>Google Earth将校正结果叠加到Google Earth查看。

图11 FY-3几何校正结果叠加到Google Earth显示