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

问题描述

我正在尝试从从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]) 

我得到了这些数据:

enter image description here

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

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

https://easyupload.io/ch55ic

解决方法

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

data_v1[,1826]

应该提取您想要的矩阵。请参阅Grouping functions (tapply,by,aggregate) and the *apply family,简要概述一些处理数组的有用函数。