python读取hdf文件 高效_python 读hdf4文件，再转写成一个tif文件

1.安装pyhdf包

(2)在dos命令行下切换到包所在的文件夹，并运行pip install 包名进行安装。

2.读hdf4文件

#导入包

from pyhdf.SD import *

from osgeo import osr

import numpy as np

(1)#读取文件

file="3B43.20100501.7.HDF"

hdf=SD(file)

#获取hdf文件的属性信息

attr=hdf.attributes()

#获取hdf文件的图层以及每个图层所对应的行与列等信息，datasets是一个字典

datasets=hdf.datasets()

#读某一图层(我这里的图层是降雨图层),图层名可通过datasets看到。并获取图层数据

rainfall_array=hdf.select("precipitation").get()

#图层的旋转与转置(这里我要做旋转转置是由于我的图层需求，可通过查看自己的图层是否需要做这一步)

#对矩阵进行左右翻转

data=np.fliplr(rainfall_array)

#对矩阵进行转置

data=np.transpose(data)

#从图层中读取到的数据是5月份逐小时数据，我们这里转换成逐月数据

data=data*24*31

#将数据转换成两个字节的整型数据

d=data.astype(np.uint16)

(2)转tif文件