使用python的h5py包读取hdf文件时出错

Question

我想从MODIS网站下载的hdf文件中提取数据。链接中提供了一个示例文件。我正在使用以下代码行读取hdf文件:

>>> import h5py
>>> f = h5py.File( 'MYD08_M3.A2002182.051.2008334061251.psgscs_000500751197.hdf', 'r' )

得到的错误:

Traceback (most recent call last):
    File "<pyshell#3>", line 1, in <module>
f = h5py.File( 'MYD08_M3.A2002182.051.2008334061251.psgscs_000500751197.hdf', 'r' )
    File "C:Python27libsite-packagesh5py_hlfiles.py", line 165, in __init__
fid = make_fid(name, mode, userblock_size, fapl)
    File "C:Python27libsite-packagesh5py_hlfiles.py", line 57, in make_fid
fid = h5f.open(name, h5f.ACC_RDONLY, fapl=fapl)
    File "h5f.pyx", line 70, in h5py.h5f.open (h5pyh5f.c:1640)
IOError: unable to open file (File accessability: Unable to open file)

我已经尝试了来自不同来源的其他几个hdf文件，但我得到了同样的错误。这里有什么问题吗?

Answer 1

我认为可能有两个问题:

1)由于文件扩展名是"hdf"，可能这是一个HDF4文件。HDF5文件通常有"。Hdf5"或"。h5·"扩展。我不确定h5py是否能够读取HDF4文件。

2)也许你必须改变文件本身的权限。如果您在linux机器上，请尝试:chmod +r file.hdf

你可以尝试用HDFView打开你的文件。该软件可在多个平台上使用。你可以很容易地检查文件的属性

Answer 2

这听起来像是文件权限错误，甚至是文件存在。也许可以添加一些检查，如

import os
hdf_file = 'MYD08_M3.A2002182.051.2008334061251.psgscs_000500751197.hdf'
if not os.path.isfile(hdf_file):
    print 'file %s not found' % hdf_file
if not os.access(hdf_file, os.R_OK):
    print 'file %s not readable' % hdf_file
f = h5py.File(hdf_file, 'r')

Answer 3

我遇到了同样的问题，后来发现我的文件只有"只读"，由于某种原因，它阻止了h5py读取它。将权限修改为"写"后，我就可以读了。不知道为什么要这样设置