使用R从netcdf数据中提取时间序列

时间:2013-09-16 01:04:53

标签: r netcdf

与R相关的新手问题如何使用netdcf文件中的R提取特定位置的时间序列数据。例如,以下快照显示位置(1,2)的时间序列为13,28,43。

Sample netcdf

提前致谢。

2 个答案:

答案 0 :(得分:1)

这可能会这样做,其中"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)

答案 1 :(得分:0)

{{1}}

这应该有效。