快捷導(dǎo)航

Python讀取hdf文件并轉(zhuǎn)化為tiff格式輸出

更新時(shí)間：2022年07月17日 15:43:09 作者：Amyniez

這篇文章主要介紹了Python讀取hdf文件并轉(zhuǎn)化為tiff格式輸出，文章圍繞主題展開詳細(xì)的內(nèi)容介紹，具有一定的參考價(jià)值，需要的小伙伴可以參考一下

1.Python3讀取hdf文件

最開始使用Python導(dǎo)入pyhdf包的時(shí)候是可以的，但是當(dāng)導(dǎo)入pyhdf.SD的時(shí)候就出現(xiàn)了以下問題：

我查詢了很多的相關(guān)問題，但是都不能解決我的問題。
網(wǎng)絡(luò)上解釋，這種原因是由于 numpy 和 pyhdf 的版本號(hào)不對(duì)應(yīng)所造成的，而我用的是Python3.9版本，確實(shí)找不到哪一個(gè)版本的numpy與之對(duì)應(yīng)。
但是網(wǎng)上有人說是因?yàn)閜yhdf只能在Python2中才能運(yùn)行，于是我又下載了Python2.7，但其實(shí)pyhdf在Python2、3上都能運(yùn)行：

2.使用Python2讀取hdf文件

我在官網(wǎng)上下載了最新的Python2.7.18版本（目前已停止更新），并進(jìn)行安裝。
安裝好Python2.7之后，需要重新配置環(huán)境變量：此電腦（右擊）——屬性——高級(jí)系統(tǒng)設(shè)置——環(huán)境變量（高級(jí)）——在系統(tǒng)變量（s）的path中添加Python2.7的路徑：

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

配置成功！

3. 如何使Python2和3兼容

因?yàn)槲译娔X里面原本裝有Python3.9，此時(shí)又裝入了Python2.7，某些命令難免會(huì)沖突，比如，pip下載whl時(shí)，不知道是下載到那個(gè)版本中，所以這里需要使兩個(gè)版本區(qū)別開來。

3.1 修改Python2.7的名字

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

之后再在CMD窗口中運(yùn)行Python：

python版本號(hào)查看：

配置成功！

3.2 pip更新

Python 安裝包需要用到包管理工具pip，但是當(dāng)同時(shí)安裝python2和python3的時(shí)候，pip只是其中一個(gè)版本，以下將提供一個(gè)修改方式，即重新安裝兩個(gè)版本的pip，使得兩個(gè)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對(duì)應(yīng)的包

這里簡(jiǎn)單下載了numpy、GDAL、hypdf，其中GDAL和pyhdf只在官網(wǎng)下載的：

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

**注意：**要切換到文件所在文件夾，才能完成安裝！

pyhdf.SD成功導(dǎo)入，問題解決：

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)

文件展示：