我正在学习在我的工作中使用ERA过渡。我需要分析一个气象站的两年序列，这是shapefile格式的一个点，并且还需要可视化流域的数据。 到目前为止，我已经阅读了netcdf文件并声明了变量
rm(list=ls())
graphics.off()
library(raster)
library(ncdf4)
library(maps)
library(maptools)
library(rgdal)
library(tidyverse)
library(reticulate)
library(ggthemes)
library(viridis)
setwd("C:/Users/TPPC/Desktop")
rutas <- list.files(path=getwd())
ncin<-nc_open("interim_2t_20152016.nc")
lat=ncvar_get(ncin,'latitude')
lon=ncvar_get(ncin,'longitude')
t<- ncvar_get(ncin, "time")
tunits<-ncatt_get(ncin,'time')
tunits$units
tustr<- strsplit(tunits$units, " ")
timestamp = as.POSIXct(t*3600,tz='GMT',origin=tustr[[1]][3])
data<-data_frame(name=attributes(ncin$var)$names) %>%
bind_cols(map_df(.$name,ncatt_get,nc=ncin)) %>%
mutate(values=map(name,ncvar_get,nc=ncin))
nc_close(ncin)
data
df<-expand.grid(lon=lon-360,lat=lat,timestamp=timestamp,name=data$name) %>%
mutate(coord=factor(paste(lon,lat,'/')))
*但是在最后一行中，我得到下一个错误 错误：无法分配大小为645.2 Mb的向量