我正在尝试创建一个拥有1981 - 2014年月度数据的.nc文件。但我从.dat文件中获取数据,每个文件代表一年中的一个月。
从本教程http://geog.uoregon.edu/GeogR/topics/netcdf-to-raster.html我正在使用此代码:
library(ncdf4)
nlon <- 360
nlat<- 180
dlon <- 360.0/nlon
dlat <- 180.0/nlat
lon <- seq(-180.0+(dlon), +180.0,by=dlon)
lat <- seq(-90.0+(dlat),+90.0,by=dlat)
#This is how I add one data set
DD <- as.matrix(read.table("~/1981_01_NAM.dat", header=FALSE))
datos1 <- t(matrix(t(DD),ncol=360,nrow=180, byrow=TRUE))
londim <- ncdim_def("lon", "degrees_east", as.double(lon))
latdim <- ncdim_def("lat", "degrees_north", as.double(lat))
varname ="humedad"
units = "z-scores"
dlname <- "test-variable -- original"
fillvalue <- 1e20
tmp.def <- ncvar_def(varname, units, list(londim, latdim), fillvalue, dlname, prec = "single")
ncfname <- "test-humedad.nc"
ncout <- nc_create(ncfname, list(tmp.def), force_v4 = TRUE)
ncvar_put(ncout, tmp.def, datos1)
ncatt_put(ncout, "lon","axis","X")
ncatt_put(ncout, "lat","axis","Y")
title <- "Prueba Humedad"
ncatt_put(ncout, 0, "title", title)
nc_close(ncout)
那么如何将每个.dat文件与.nc文件中的一个月对齐以便我可以获得完整的时间序列?