Используйте R для извлечения временных рядов из данных netcdf

Вопрос новичка, связанный с R. Как извлечь данные временных рядов для определенного местоположения с помощью R из файла netdcf. Так, например, следующий снимок показывает, что временной ряд для местоположения (1,2) равен 13,28,43.

Образец netcdf

Заранее спасибо.


person umbersar    schedule 16.09.2013    source источник
comment
Пожалуйста, начните с воспроизводимого примера: stackoverflow.com/questions/5963269/   -  person zero323    schedule 16.09.2013
comment
Очевидно, что эти данные находятся в Excel (или аналогичном) — как вы планировали импортировать их в R и в каком виде? Это будет иметь огромное значение для потенциальных ответов. Скриншоты данных также практически бесполезны для дальнейшего использования.   -  person thelatemail    schedule 16.09.2013
comment
Удалось ли вам вообще получить данные netcdf в R? Взгляните на пакет ncdf.   -  person Scott Ritchie    schedule 16.09.2013


Ответы (2)


Это может сделать это, где "my.variable" — это имя интересующей вас переменной:

library(survival)
library(RNetCDF)
library(ncdf)
library(date)

setwd("c:/users/mmiller21/Netcdf")

my.data <- open.nc("my.netCDF.file.nc");

my.time <- var.get.nc(my.data, "time")

n.latitudes  <- length(var.get.nc(my.data, "lat"))
n.longitudes <- length(var.get.nc(my.data, "lon"))

n.dates <- trunc(length(my.time))
n.dates

my.object <- var.get.nc(my.data, "my.variable")

my.array  <- array(my.object, dim = c(n.latitudes, n.longitudes, n.dates))
my.array[,,1:5]
my.vector <- my.array[1, 2, 1:n.dates]  # first latitude, second longitude
head(my.vector)

baseDate <- att.get.nc(my.data, "NC_GLOBAL", "base_date")
bdInt    <- as.integer(baseDate[1])

year     <- date.mdy(seq(mdy.date(1, 1, bdInt), by = 1,
                     length.out = length(my.vector)))$year 

head(year)
person Mark Miller    schedule 16.09.2013

your_data <- read.csv("")

#Subsetting your data
location12 <- subset(your_data, latitude == 1 & column2_value == 2)

Data_location12 <- table(location12)

timeseries12 <- ts(Data_location12)

Это должно работать.

person Hari Narayan Singh    schedule 25.02.2017