Создание многомерного NetCDF в R
Я пытаюсь создать многомерный файл NetCDF с помощью пакета R Ncdf. Я работаю с климатическими ежедневными наблюдениями для набора из 1500 точек, количество наблюдений составляет ~ 18250 для каждой точки. Проблема в том, что структура файла NetCDF Create.ncdf) занимает 4 ГБ, и каждая точка увеличивает размер файла более чем на 3 ГБ Put.var.ncdf)
Это код, который я использую:
<code># Make a few dimensions we can use dimX <- dim.def.ncdf( "Long", "degrees", Longvector ) dimY <- dim.def.ncdf( "LAT", "degrees", Latvector ) dimT <- dim.def.ncdf( "Time", "days", 1:18250, unlim=FALSE ) # Make varables of various dimensionality, for illustration purposes mv <- -9999 # missing value to use var1d <- var.def.ncdf( "var1d", "units", dimX, mv,prec="double" ) var2d <- var.def.ncdf( "var2d", "units", list(dimX,dimY), mv,prec="double" ) var3d <- var.def.ncdf( "var3d", "units", list(dimX,dimY,dimT), mv,prec="double" ) # Create the test file nc <- create.ncdf( "writevals.nc", list(var1d,var2d,var3d) ) # !!Creates a nc file with + 4 Gb # Adding the complete time series for one point (the first point in the list of the dataset) put.var.ncdf( nc, var3d,dataset[[1]], start=c(Longvector[1],Latvector[1],1), count=c(1,1,-1)) </code>
Longvector и Latvector - векторы, взятые из матрицы с Long и Lat для каждой точки. Набор данных представляет собой формат списка, и для каждой точки у меня есть список числовых значений.
<code>dataset[[1]]=c(0,0,0,9.7,0,7.5,3.6,2.9,0,0.5,....) </code>
Я что-то пропустил или мне стоит попробовать другие пакеты ??