如何从 R 中的 netcdf 文件中提取降水数据？

Question

我正在尝试从从 cordex 气候 model 下载的 netcdf 文件中提取降水数据。 我应用了以下代码：

library(ncdf4)    
v1 <- nc_open('pr_EUR-44_ICHEC-EC-EARTH_rcp45_r12i1p1_SMHI-RCA4_v1_day_20960101-21001231.nc')
data_v1<-ncvar_get(v1,attributes(v1$var)$names[3]) 

我得到了这个数据：

在这里，我需要获取行号为 1826 的数据。

我在这里上传了 20 mb 大小的 nc 文件：

https://easyupload.io/ch55ic

Answer 1

该数组由 1826 个矩阵组成，每个矩阵有 106 行和 103 列。

data_v1[,,1826]

应该提取你所追求的矩阵。 有关处理 arrays 的一些有用函数的简要概述，请参阅分组函数（tapply、by、aggregate）和 *apply 系列。