子数据集
MODIS图像(中分辨率成像光谱仪)是一个分层数据格式(HDF)文件。如果数据包含子数据集,可用 GetSubDatasets() 函数得到他们的子数据集列表,再利用这个列表信息打开所需要的子数据集。

import os import numpy as np from osgeo import gdal data_dir = r'D:\modis' # 从2024文件中获取子数据集的名称和描述 # 并打印,输出NDVI(归一化植被指数) os.chdir(os.path.join(data_dir, '2024')) ds = gdal.Open('MOD0.A201.hdf') subdatasets = ds.GetSubDatasets() # GetSubDatasets函数返回一个元组列表 print('Number of subdatasets: {}'.format(len(subdatasets))) for sd in subdatasets: print('Name: {0}\nDescription:{1}\n'.format(*sd)) # 打开2024文件中的第一个子数据集:[0][0] # 第五个数据集为: [4][0] ndvi_ds = gdal.Open(subdatasets[0][0]) # 通过打印尺寸来确保正常工作 print('Dataset dimensions: {} {}'.format(ndvi_ds.RasterXSize, ndvi_ds.RasterYSize)) # 读取数据之前先获取波段 ndvi_band = ndvi_ds.GetRasterBand(1) print('Band dimensions: {} {}'.format(ndvi_band.XSize, ndvi_band.YSize))
,Best Wish 不负年华
【推荐】编程新体验,更懂你的AI,立即体验豆包MarsCode编程助手
【推荐】凌霞软件回馈社区,博客园 & 1Panel & Halo 联合会员上线
【推荐】抖音旗下AI助手豆包,你的智能百科全书,全免费不限次数
【推荐】博客园社区专享云产品让利特惠,阿里云新客6.5折上折
【推荐】轻量又高性能的 SSH 工具 IShell:AI 加持,快人一步
· CSnakes vs Python.NET:高效嵌入与灵活互通的跨语言方案对比
· DeepSeek “源神”启动!「GitHub 热点速览」
· 我与微信审核的“相爱相杀”看个人小程序副业
· Plotly.NET 一个为 .NET 打造的强大开源交互式图表库
· 上周热点回顾(2.17-2.23)