问题描述
我正在尝试从从cordex 气候模型下载的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 的数据。
解决方法
该数组由 1826 个矩阵组成,每个矩阵有 106 行和 103 列。
data_v1[,1826]
应该提取您想要的矩阵。请参阅Grouping functions (tapply,by,aggregate) and the *apply family,简要概述一些处理数组的有用函数。