如何解决从多个 hdf 文件中提取变量
我有这 3 个文件 here。试图堆叠它们并提取一个名为“OzoneTropColumn”的变量。我可以为单个文件提取数据,但无法为多个文件提取数据。
library(raster)
library(ncdf4)
list_col1 <- list.files("E:/TES",pattern = "*.hdf",full.names = TRUE)
ncin1 <- raster::stack(list_col1,varname = "Data Fields/OzoneTropColumn",ncdf=TRUE)
任何帮助将不胜感激。
谢谢
解决方法
我不知道这是否真的有效,因为我不是 raster
或导入 hdf
文件的专家。但是这段代码确实为我加载了 R 中的三个文件,但仍然给出了一些警告:-
library(raster)
library(ncdf4)
list_col1 <- as.list(list.files("E:/TES",pattern = "*.hdf",full.names = TRUE))
ncin1 <- raster::stack(list_col1,varname = "Data Fields/OzoneTropColumn",ncdf=TRUE)
[1] "vobjtovarid4: **** WARNING **** I was asked to get a varid for dimension named Data Fields/phony_dim_1 BUT this dimension HAS NO DIMVAR! Code will probably fail at this point"
[1] "vobjtovarid4: **** WARNING **** I was asked to get a varid for dimension named Data Fields/phony_dim_0 BUT this dimension HAS NO DIMVAR! Code will probably fail at this point"
[1] "vobjtovarid4: **** WARNING **** I was asked to get a varid for dimension named Data Fields/phony_dim_1 BUT this dimension HAS NO DIMVAR! Code will probably fail at this point"
[1] "vobjtovarid4: **** WARNING **** I was asked to get a varid for dimension named Data Fields/phony_dim_0 BUT this dimension HAS NO DIMVAR! Code will probably fail at this point"
[1] "vobjtovarid4: **** WARNING **** I was asked to get a varid for dimension named Data Fields/phony_dim_1 BUT this dimension HAS NO DIMVAR! Code will probably fail at this point"
[1] "vobjtovarid4: **** WARNING **** I was asked to get a varid for dimension named Data Fields/phony_dim_0 BUT this dimension HAS NO DIMVAR! Code will probably fail at this point"
[1] "vobjtovarid4: **** WARNING **** I was asked to get a varid for dimension named Data Fields/phony_dim_1 BUT this dimension HAS NO DIMVAR! Code will probably fail at this point"
[1] "vobjtovarid4: **** WARNING **** I was asked to get a varid for dimension named Data Fields/phony_dim_0 BUT this dimension HAS NO DIMVAR! Code will probably fail at this point"
但确实告诉我它已经读取了 3 个文件
ncin1[[1]]
class : RasterLayer
dimensions : 83,90,7470 (nrow,ncol,ncell)
resolution : 1,1 (x,y)
extent : 0.5,90.5,0.5,83.5 (xmin,xmax,ymin,ymax)
crs : NA
source : /home/shawn/Downloads/TES/TES-Aura_L3-O3-M2004m09_F01_12.hdf
names : Data.Fields.OzoneTropColumn.1
zvar : Data Fields/OzoneTropColumn
ncin1[[2]]
class : RasterLayer
dimensions : 83,ymax)
crs : NA
source : /home/shawn/Downloads/TES/TES-Aura_L3-O3-M2004m10_F01_12.hdf
names : Data.Fields.OzoneTropColumn.2
zvar : Data Fields/OzoneTropColumn
ncin1[[3]]
class : RasterLayer
dimensions : 83,ymax)
crs : NA
source : /home/shawn/Downloads/TES/TES-Aura_L3-O3-M2004m11_F01_12.hdf
names : Data.Fields.OzoneTropColumn.3
zvar : Data Fields/OzoneTropColumn
,
Robert 发布的脚本运行良好,唯一的区别是使用 as.list 和不使用 as.list]1]1。我在没有as.list的情况下使用之前尝试过。由于未知原因,提取的数据丢失了其地理参考 。所以必须在光栅堆栈后给出范围和投影。谢谢罗伯特。
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。