r - 如何从 R 中的 HE5 (HDF5) 数据集中读取和提取属性数据

Question

数据集包含以下具有不同属性数据的组。如何使用包'rhdf5'、'rgdal'和'raster'提取和绘制具有纬度经度的属性编号8？

8  /HDFEOS/GRIDS/ColumnAmountNO2/Data Fields              ColumnAmountNO2Trop H5I_DATASET  FLOAT 1440 x 720

h5ls("D:/NO2/NO2 Data/2008/1Jan/OMI-Aura_L3-OMNO2d_2008m0101_v003-2018m0626t175046.he5")

                                       group                             name       otype dclass        dim
0                                          /                           HDFEOS   H5I_GROUP                  
1                                    /HDFEOS                       ADDITIONAL   H5I_GROUP                  
2                         /HDFEOS/ADDITIONAL                  FILE_ATTRIBUTES   H5I_GROUP                  
3                                    /HDFEOS                            GRIDS   H5I_GROUP                  
4                              /HDFEOS/GRIDS                  ColumnAmountNO2   H5I_GROUP                  
5              /HDFEOS/GRIDS/ColumnAmountNO2                      Data Fields   H5I_GROUP                  
6  /HDFEOS/GRIDS/ColumnAmountNO2/Data Fields                  ColumnAmountNO2 H5I_DATASET  FLOAT 1440 x 720
7  /HDFEOS/GRIDS/ColumnAmountNO2/Data Fields     ColumnAmountNO2CloudScreened H5I_DATASET  FLOAT 1440 x 720
8  /HDFEOS/GRIDS/ColumnAmountNO2/Data Fields              ColumnAmountNO2Trop H5I_DATASET  FLOAT 1440 x 720
9  /HDFEOS/GRIDS/ColumnAmountNO2/Data Fields ColumnAmountNO2TropCloudScreened H5I_DATASET  FLOAT 1440 x 720
10 /HDFEOS/GRIDS/ColumnAmountNO2/Data Fields                           Weight H5I_DATASET  FLOAT 1440 x 720
11                                         /               HDFEOS INFORMATION   H5I_GROUP                  
12                       /HDFEOS INFORMATION                 StructMetadata.0 H5I_DATASET STRING      ( 0 )

我使用以下命令读取组内的属性，但显示错误

E = h5read("D:/NO2/NO2 Data/2008/1Jan/OMI-Aura_L3-OMNO2d_2008m0101_v003-2018m0626t175046.he5","HDFEOS/GRIDS/ColumnAmountNO2/Data Fields/ColumnAmountNO2Trop")
Warning message:
In h5checktypeOrOpenLoc(file, readonly = TRUE, native = native) :
  An open HDF5 file handle exists. If the file has changed on disk meanwhile, the function may not work properly. Run 'h5closeAll()' to close all open HDF5 object handles.

任何帮助，将不胜感激。谢谢

标签： rrasterhdf5rgdalrhdf5