MODIS数据HDF批量转化为TIF

mmmmmyyyyyhhhhh

已于 2023-03-08 16:30:51 修改

阅读量1.2k

点赞数

文章标签： java javascript 开发语言 Powered by 金山文档

于 2023-03-08 11:23:57 首次发布

本文链接：https://blog.csdn.net/qq_42655653/article/details/129397397

版权

原始数据MCD12Q1，空间分辨率500m，转化为WGS84.

hdf ->tif

GDAL

from osgeo import gdal
import os
import glob
os.chdir("D:\MCD12Q1\\try")
file_list = glob.glob("*.hdf")
for i in file_list:
    name = (i.split('\\'))[-1]
    namee = name[:-4]
    datasets = gdal.Open(i)
    SubDatasets = datasets.GetSubDatasets()
    Metadata = datasets.GetMetadata()
    data = datasets.GetSubDatasets()[0][0]
    Raster_DATA = gdal.Open(data)
    DATA_Array = Raster_DATA.ReadAsArray()
    TifName = "D:\MCD12Q1\\try"+'\\'+namee+'.tif'
    geoData = gdal.Warp(TifName, Raster_DATA,
                    dstSRS = 'EPSG:4326', format = 'GTiff',
                    resampleAlg = gdal.GRA_Bilinear)
    del geoData

得到结果之后发现问题，每个tif文件的像元大小都不一致

0.002997 x 0.002997 Degrees

0.003661 x 0.003661 Degrees

0.102315 x 0.102344 Degrees

。。。。。。。。。。

不知道为什么，望解答