Подмножество ERA-промежуточная серия с шейп-файлом в R - PullRequest
Новая
       12

Подмножество ERA-промежуточная серия с шейп-файлом в R

0 голосов
/ 20 февраля 2019

Я учусь использовать ERA-interim в своей работе.Мне нужно проанализировать двухлетнюю серию для метеорологической станции, которая является точкой в ​​формате шейп-файла, а также визуализировать данные для бассейна.До сих пор я прочитал файл 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 МБ

...