如何使用 R 中的 ncdf4 包正确读取气候研究数据（CRU）？

【问题标题】：How to read climate research data (CRU) correctly using ncdf4 package in R?如何使用 R 中的 ncdf4 包正确读取气候研究数据（CRU）？
【发布时间】：2017-10-18 00:38:14
【问题描述】：

CRU netCDF 数据具有月平均值。我的文件有 1901-1910 年期间的数据。当我读取变量 (tmp) 时，它只显示了 10 年 120 个月中的 12 个级别。

library(ncdf4)
nc<- nc_open(file.choose())
lat<- ncvar_get(nc, "latitude")
lon<- ncvar_get(nc, "longitude")
tm<- ncvar_get(nc, "tmp")

dim(tm)
[1] 720 360  12

数据文件的链接在这里：- https://drive.google.com/open?id=0B9ACBRSUXCjnWWxZclR3alNIdEE

【问题讨论】：

我认为我的回答解决了您的问题 - 请您勾选绿色箭头以标记为已接受的答案，以便可以关闭此问题？

标签： r netcdf

【解决方案1】：

如果你查看文件的标题

ncdump -h cru_tmp_clim_1901-1910.nc

您会看到该文件只有 12 个时间片。

netcdf cru_tmp_clim_1901-1910 {
dimensions:
time = 12 ;
latitude = 360 ;
longitude = 720 ;
bounds = 2 ;

如果你使用 cdo：

cdo showdate cru_tmp_clim_1901-1910.nc 
1905-01-16  1905-02-15  1905-03-16  1905-04-16  1905-05-16  1905-06-16  1905-07-16  1905-08-16  1905-09-16  1905-10-16  1905-11-16  1905-12-16

您会看到数据仅适用于 1905 年。

我认为有人向您传递了一个已被删除的数据文件。我建议从源重新下载原始 CRU 数据。请记住，如果您有疑问，最好使用 ncdump 来示例文件头和元数据以查看其中的内容。

【讨论】：