HDF文件转Tif 温室气体 python代码

原创

于 2022-02-14 19:53:29 发布 · 1.1k 阅读

这篇博客分享了两种使用Python将HDF文件转换为TIFF文件的方法。涉及的关键技术包括GDAL库、pyhdf库以及地理坐标系统的处理。代码详细解释了如何读取HDF文件中的经纬度和数据，计算分辨率，设置地理变换和投影，并最终写入TIFF文件。

HDF文件转tif文件，python代码

这个问题网上找了好多，基本用处不大，我专门找了地理学的同学给做的，对于计算机的我来说确实是个难点，话不多说

hdf中所用到的三个数据，我用的是co2浓度和经纬度，如图
在这里插入图片描述
我有两个代码，以下是代码1，

import os, sys
from osgeo import gdal,ogr,osr
import numpy as np
from pyhdf.SD import SD

os.environ['PROJ_LIB'] = r'D:\python\pythonData\venv\Lib\site-packages\osgeo\data\proj'
def readhdf(filename):

    file = SD(filename)   
    lon = (file.select('lon_grid')).get()
    lat = (file.select('lat_grid')).get()
    CO2Amount = (file.select('CO2Amount_grid')).get()
    lon_res = 1
    lat_res = -1
    return CO2Amount,lon_res,lat_res,lon,lat


def writetif(data,outname,geotransform):
    nl,ns = [data.shape[0],data.shape[1]] 
    bands = 1               
    driver = gdal.GetDriverByName("GTiff")
    out_tif = driver.Create(outname, ns, nl, bands, gdal.GDT_Flo

标签