利用gdalUtils库从MODIS LAI/FPAR hdf格式中提取sds

Question

我正在尝试使用gdalUTILs库提取MODIS LAI/FPAR图像的所有图层，我遇到了这个错误，有人能帮我解决这个问题吗我正在尝试使用

获取SDS名称列表

sds <- get_subdatasets('MOD15A2H.A2016289.h28v08.006.2016302012314.hdf')

隔离第一个SDS的名称

name <- sds[1]`enter code here`
filename <- 'name/of/output/file.tif'
gdal_translate(sds[1], dst_dataset = filename)

将创建的Geotiff加载到R中

r <- raster(filename)

错误是

> sds <- get_subdatasets('MOD15A2H.A2016289.h28v08.006.2016302012314.hdf')

文件[1]中出错:下标超出界限此外:警告消息:正在运行的命令'"C:\Program Files (X86) "MOD15A2H.A2016289.h28v08.006.2016302012314.hdf"‘“split1的状态为1

Answer 1

检查hdf文件大小。如果它是1或2kb，那么它的html错误文件。我在MOD15A2H产品上也遇到了同样的问题。