使用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')

split1 [[1]]出错：下标越界 另外：警告信息： 运行命令'“C：\ Program Files（x86）\ Boundless \ OpenGeo \ bin \ gdalinfo.exe”“MOD15A2H.A2016289.h28v08.006.2016302012314.hdf”'状态为1

Answer 1

检查hdf文件大小。如果它是1或2kb，则其html错误文件。我对MOD15A2H产品有同样的问题。