library(sf)
library(gdalUtils)
library(rgdal)
library(tidyverse)
library(raster)
file_A1 <- "F:/mcd19A1/h19v08/MCD19A1.A2000058.h19v08.006.2018013034346.hdf"
sds = get_subdatasets(file_A1)
tif <- "F:/tif.tif"
tmp <- gdal_translate(sds[1], dst_dataset = tif)
mcd19a1数据集在sds[1]中有不不止一个波段,但最后得到的是两个相同的band
相似问题