我正试图从哥白尼中提取一个气候变量。
devtools::install_github("r-spatial/stars")
library(stars)
library(ncdf4)
library(RNetCDF)
pp <- read_ncdf("~/climate_data/pp_ens_mean_0.1deg_reg_v25.0e.nc", proxy = TRUE)
> print(pp)
netcdf source stars proxy object from:
[1] "[...]/pp_ens_mean_0.1deg_reg_v25.0e.nc"
Available nc variables:
pp
dimension(s):
from to offset delta refsys point values x/y
longitude 1 705 -25.0001 0.1 WGS 84 NA NULL [x]
latitude 1 465 24.9999 0.1 WGS 84 NA NULL [y]
time 1 26298 1950-01-01 UTC 1 days POSIXct NA NULL
这就是我作为代理人得到的。
我想把经度在10°E到25°E之间,纬度在55°N到70°N之间。
我试过了,但不起作用,我真的不明白为什么。
pp <- read_ncdf(pp,
var = pp,
ncsub = cbind(start = c(10, 55, 1),
count = c(25, 70, 26298)))
有什么想法吗?非常感谢你的宝贵帮助。