Question

我将 NetCDF 文件读取为光栅。

library(raster)
r <- brick("./DATA/consecutive_wet_days_index_50_.nc",
varname="consecutive_wet_days_index_per_time_period")

现在，如何从栅格 varname 中提取指定的 r？我稍后需要 varname 来命名图。另外，如果在阅读时不是物种，如何提取varname？

此处为示例 NetCDF 文件：

Answer 1

其中一个可能？

gsub("(.*)\\\\","", r@file@name, perl=TRUE)
# [1] "consecutive_wet_days_index_50_.nc"

r@file@datanotation
# [1] "FLT4S"

您可以使用 str(r) 检查结构。

或者使用 capture.output hack。

x <- capture.output(r)
trimws(gsub("(.*)\\:","", x[9], perl=TRUE))
# [1] "consecutive_wet_days_index_per_time_period"