Python读取hdf文件并转化为tiff格式输出

2024-04-02 19:55

短信预约 -IT技能 免费直播动态提醒

1.Python3读取hdf文件

最开始使用Python导入pyhdf包的时候是可以的，但是当导入pyhdf.SD的时候就出现了以下问题：

我查询了很多的相关问题，但是都不能解决我的问题。
网络上解释，这种原因是由于 numpy 和 pyhdf 的版本号不对应所造成的，而我用的是Python3.9版本，确实找不到哪一个版本的numpy与之对应。
但是网上有人说是因为pyhdf只能在Python2中才能运行，于是我又下载了Python2.7，但其实pyhdf在Python2、3上都能运行：

2.使用Python2读取hdf文件

我在官网上下载了最新的Python2.7.18版本（目前已停止更新），并进行安装。
安装好Python2.7之后，需要重新配置环境变量：此电脑（右击）——属性——高级系统设置——环境变量（高级）——在系统变量（s）的path中添加Python2.7的路径：

配置好之后，可以在CMD命令窗口中输入python查看，是否成功。

配置成功！

3. 如何使Python2和3兼容

因为我电脑里面原本装有Python3.9，此时又装入了Python2.7，某些命令难免会冲突，比如，pip下载whl时，不知道是下载到那个版本中，所以这里需要使两个版本区别开来。

3.1 修改Python2.7的名字

修改 python.exe 名字为 python2.exe 和 python3.exe：

之后再在CMD窗口中运行Python：

python版本号查看：

配置成功！

3.2 pip更新

Python 安装包需要用到包管理工具pip，但是当同时安装python2和python3的时候，pip只是其中一个版本，以下将提供一个修改方式，即重新安装两个版本的pip，使得两个python版本的pip能够共存。

在DOS命令框输入命令：

1. pip3：

python3 -m pip install --upgrade pip --force-reinstall

2. pip2

python2 -m pip install --upgrade pip --force-reinstall

以上配置完成！

3.2 使用pip2下载python2.7对应的包

这里简单下载了numpy、GDAL、hypdf，其中GDAL和pyhdf只在官网下载的：

pip install GDAL-2.2.4-cp27-cp27m-win_amd64.whl

**注意：**要切换到文件所在文件夹，才能完成安装！

pyhdf.SD成功导入，问题解决：

4. 使用python2.7.18读取hdf文件并输出为TIFF格式

# -*- coding : utf-8 -*-
# C:/python2.7
from osgeo import gdal, ogr, osr
import warnings
warnings.simplefilter ("ignore")
from osgeo import gdal
from pyhdf.SD import SD, SDC, SDim
'''File path'''
filepath = r'E:\08frost\2008\NIEER_CGF-MODIS_SCE_20080101_DAILY_500m_V02.hdf'
'''Data'''
filename = filepath.split('\\')[-1].split('.')[0]
hdf = SD(filepath)
snow = hdf.select('Day_Snow_Cover_Area').get()
'''Spatial Reference'''
sr = osr.SpatialReference()
sr.ImportFromEPSG(4326)
s=sr.ExportToWkt()
driver=gdal.GetDriverByName("GTiff")
dataset=driver.Create(filename+'.tif',14000,8000,1,gdal.GDT_Int16, ["TILED=YES", "COMPRESS=LZW"])
im_geotrans=(72,0.005,0.0,56,0.0,-0.005)
dataset.SetGeoTransform(im_geotrans)
dataset.SetProjection(s)
dataset.GetRasterBand(1).WriteArray(snow)

文件展示：