使用pyhdf读取HDF4数据集失败：提示数据集不存在

我来帮你搞定这个问题！你遇到的pyhdf.error.HDF4Error: select: non-existent dataset错误，核心原因是你指定的数据集名称和HDF文件里实际存在的名称不匹配。

问题分析 #

你的代码里已经加了print(hdf.datasets())来打印所有可用数据集，但你可能没注意到：MOD13A2这类MODIS产品的HDF4文件中，NDVI数据集的名称是带空格的，而不是你用的下划线分隔格式。比如实际名称应该是"1 km 16 days NDVI"，而不是"1_km_16_days_NDVI"。

解决方案步骤 #

1. 确认正确的数据集名称
   先运行你现有的代码，查看print(hdf.datasets())输出的内容，找到NDVI对应的数据集名称。输出会类似这样（示例）：

   {'1 km 16 days NDVI': (('YDim: 1200', 'XDim: 1200'), 3, 0), ...}
   

   从这里复制准确的名称字符串。

2. 修改代码中的数据集名称
   将代码里的DATAFIELD_NAME替换为你找到的正确名称，比如：

   DATAFIELD_NAME = "1 km 16 days NDVI"
   

3. 修改后的完整代码

   import os
   import numpy as np
   from pyhdf.SD import SD, SDC
   
   # Open file.
   FILE_NAME = 'MOD13A2.A2016001.h23v05.006.2016029070140.hdf'
   hdf = SD(FILE_NAME, SDC.READ)
   
   # List available SDS datasets (确认数据集名称用)
   print(hdf.datasets())
   
   # Read dataset with correct name
   DATAFIELD_NAME = "1 km 16 days NDVI"
   data2D = hdf.select(DATAFIELD_NAME)
   data = data2D[:,:]
   
   # 可以加个打印验证数据形状
   print(f"读取的数据形状: {data.shape}")
   

额外注意事项 #

• 确保你的FILE_NAME路径正确，如果文件不在当前工作目录，要写完整的绝对路径。
• 确认pyhdf库已经正确安装，并且支持你读取的HDF4文件版本。

内容的提问来源于stack exchange，提问作者Yogesh Sathe