如何从光栅中提取变量名?
How to extract varname from raster?
我将 NetCDF 文件读取为栅格。
library(raster)
r <- brick("./DATA/consecutive_wet_days_index_50_.nc",
varname="consecutive_wet_days_index_per_time_period")
现在,如何从光栅r
中提取指定的varname
?稍后我需要 varname 来命名地块。另外,如果不是物种,如何在阅读时提取 varname?
示例 NetCDF 文件位于此处:
https://mega.nz/file/HE8X3SbB#06DxcP2_jJkTItyealz9UACKrkQ31_4_Xmkx_r-fB28
也许是其中之一?
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"
我将 NetCDF 文件读取为栅格。
library(raster)
r <- brick("./DATA/consecutive_wet_days_index_50_.nc",
varname="consecutive_wet_days_index_per_time_period")
现在,如何从光栅r
中提取指定的varname
?稍后我需要 varname 来命名地块。另外,如果不是物种,如何在阅读时提取 varname?
示例 NetCDF 文件位于此处:
https://mega.nz/file/HE8X3SbB#06DxcP2_jJkTItyealz9UACKrkQ31_4_Xmkx_r-fB28
也许是其中之一?
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"