Python 读取hdf数据

 以Modis数据为例，要读取其中'Cloud_Phase_Infrared'变量

from pyhdf.SD import SD

def readModis(filename):
    ds = SD(filename)
    attributes = ds.select('Cloud_Phase_Infrared').attributes()##'Cloud_Phase_Infrared'变量的属性
    scales = attributes['scale_factor']#乘积缩放系数  在'Cloud_Phase_Infrared'是1
    offset = attributes['add_offset']#加法缩放系数  在'Cloud_Phase_Infrared'是0
    cloud_phase = ds.select('Cloud_Phase_Infrared')[:]#(406,270)
    # 127: 缺测      0: 晴空(cloud free)      1: 水云(water cloud)       2:冰云(ice cloud)      3:混合云(mixed phase cloud)       6:不确定相态 (undetermined phase)
    cloud_phase=cloud_phase.astype('float32')#原来类型是int8 数据只能在-128 - 127之间 改成float32好处理
    cloud_phase = cloud_phase * scales + offset
    cloud_phase[cloud_phase==127]=-999
    lon = ds.select('Longitude')[:]#(406,270)
    lat = ds.select('Latitude')[:]#(406,270)

    return cloud_phase,lon,lat

if __name__=="__main__":
    raw_phase_modis,raw_lon_modis,raw_lat_modis=readModis('/MOD06_L2.A2016063.0200.061.2017325162309.hdf')#(406,270)

另外，想查看hdf变量信息使用：

ds.info()，以我这个数据为例，输出是(127,14) 其中127是127个变量  14是属性.info(): get the number of datasets stored in the file and the number of attributes attached to it

打印数据中变量名：

for idx, sds in enumerate(ds.datasets().keys()):#ds.datasets().keys() 变量名称
    print(idx, sds)