以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)
版权声明：本文内容由互联网用户自发贡献，该文观点仅代表作者本人。本站仅提供信息存储空间服务，不拥有所有权，不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容， 请发送邮件至 举报，一经查实，本站将立刻删除。
文章由极客之音整理，本文链接：https://www.bmabk.com/index.php/post/73508.html
