MODIS 数据格式为hdf,然而,gdal 并不能直接读取hdf数据,而是要通过hdf所包含的子数据集合来读取。
datasets = gdal.Open(r"MCD12C1.A2003001.006.2018053185458.hdf")
print(datasets.GetSubDatasets())
subdataset_one = datasets.GetSubDatasets()[0][0] # 第一个子数据集合
raster = gdal.Open(subdataset_one)
raster_arrlike = raster.ReadAsArray() # 读入arr类型