如何解决HDF格式变量提取
我有一个包含hdf文件的文件夹,试图从变量列表中提取变量(NO2的列数)。提取变量后,堆叠包含该特定变量的所有文件,并取所提取变量的平均值。 Data link
P.S:我可以使用以下代码从单个文件中提取文件,但是,不能从文件列表中提取文件。
a <- raster("E:2020/jan.hdf",var="Data Fields/ColumnAmountNO2Trop",ncdf=TRUE)
library(raster)
library(ncdf4)
library(rgdal)
a <- raster("E:2020/jan.hdf",ncdf=TRUE)
dput(a)
new("RasterLayer",file = new(".RasterFile",name = "E:2020/jan.hdf",datanotation = "FLT4S",byteorder = "little",nodatavalue = -1.26765060022823e+30,NAchanged = FALSE,nbands = 1L,bandorder = "BIL",offset = 0L,toptobottom = TRUE,blockrows = 0L,blockcols = 0L,driver = "netcdf",open = FALSE),data = new(".SingleLayerData",values = logical(0),offset = 0,gain = 1,inmemory = FALSE,fromdisk = TRUE,isfactor = FALSE,attributes = list(),haveminmax = FALSE,min = Inf,max = -Inf,band = 1L,unit = "",names = "Data.Fields.ColumnAmountNO2Trop"),legend = new(".RasterLegend",type = character(0),color = logical(0),names = logical(0),colortable = logical(0)),title = character(0),extent = new("Extent",xmin = 0.5,xmax = 1440.5,ymin = -0.000695410292072429,ymax = 1441.00069541029),rotated = FALSE,rotation = new(".Rotation",geotrans = numeric(0),transfun = function ()
NULL),ncols = 1440L,nrows = 720L,crs = new("CRS",projargs = NA_character_),history = list(),z = list(
NULL))
r <- list.files("E:/2020",pattern = "*.hdf",full.names = TRUE)
当我想要堆叠变量时,出现以下错误消息
ext_r <- raster::stack(r,varname = "Data Fields/ColumnAmountNO2Trop")
rgdal :: GDALinfo中的错误(文件名,无声=无声,returnRAT = RAT,returnCategoryNames = RAT): 找不到对象“ RATlist” h(simpleError(msg,call))中的错误: 选择函数'stack'的方法时评估参数'x'的错误:选择函数'stack'的方法时评估参数'x'的错误:无法从此文件创建RasterLayer对象。 另外:警告消息: 在dim(x)中:数据集中没有条带
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。