Похоже, у вас две проблемы.
Первая связана с ожиданием той же структуры, что и файл netCDF в R, что является проблемой. само по себе, потому что, когда вы переводите многомерную структуру массива netCDF в двухмерный фрейм данных. Формат NetCDF нуждается в некотором изменении формы в R для того, чтобы манипулировать им, как в python (см .: http://geog.uoregon.edu/bartlein/courses/geog490/week04-netCDF.html).
Второй заключается в том, что вы используете значения вместо индексов при подборе данных.
umed <- (uwnd[ , , 10, 421] + uwnd[ , , 10, 422] + uwnd[ , , 10, 423])/3
Решение , которое я вижу для этого, начинается с создания индексов измерений, которые вы хотите поднастроить. В этом примере я устанавливаю уровень уверенности 10 миллибар и все, что идет между долготой 230 и 300 и широтой 25 и 40.
nc <- nc_open("uwnd.mon.mean.nc")
LonIdx <- which( nc$dim$lon$vals > 230 & nc$dim$lon$vals <300 )
## [1] 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113
## 114 115 116 117 118 119 120
LatIdx <- which( nc$dim$lat$vals >25 & nc$dim$lat$vals < 40)
## [1] 22 23 24 25 26
LevIdx <- which( nc$dim$level$vals==10)
## [1] 17
Тогда вам нужно будет применять индексы для каждого измерения, кроме времени, которое я бы Предположим, вы не хотите подмножество. Подстановка lon и latitude важна, потому что R сохраняет все в памяти, поэтому оставление всего их диапазона потребует значительного объема оперативной памяти.
lat <- ncvar_get(nc,"lat")[LatIdx]
lon <- ncvar_get(nc,"lon")[LonIdx]
lev <- ncvar_get(nc,"level")[LevIdx]
time <- ncvar_get(nc,"time")
После этого вы можете получить переменную, которую искали uwnd Ежемесячный U-ветер на уровнях давления и окончание sh чтение файла netCDF с помощью nc_close(nc).
uwnd <- ncvar_get(nc,"uwnd")[LonIdx,LatIdx,LevIdx,]
nc_close(nc)
В конце вы можете расширить сетку всеми четырьмя измерения: долгота, широта, уровень уверенности и время.
uwndf <- data.frame(as.matrix(cbind(expand.grid(lon,lat,lev,time))),c(uwnd))
names(uwndf) <- c("lon","lat","level","time","U-wind")
Свяжите его с кадром данных с помощью переменной U-wind и преобразуйте переменную времени netcdf в объект времени R.
uwndf$time_final<-convertDateNcdf2R(uwndf$time, units = "hours", origin =
as.POSIXct("1800-01-01", tz = "UTC"),time.format="%Y-%m-%d %Z %H:%M:%S")
В конце у вас будет дата-кадр, который вы ищете в период с января 1979 года по март 2020 года.
max(uwndf$time_final)
## [1] "2020-03-01 UTC"
min(uwndf$time_final)
## [1] "1979-01-01 UTC"
head(uwndf)
## lon lat level time U-wind time_final
## 1 232.5 37.5 10 1569072 3.289998 1979-01-01
## 2 235.0 37.5 10 1569072 5.209998 1979-01-01
## 3 237.5 37.5 10 1569072 7.409998 1979-01-01
## 4 240.0 37.5 10 1569072 9.749998 1979-01-01
## 5 242.5 37.5 10 1569072 12.009998 1979-01-01
## 6 245.0 37.5 10 1569072 14.089998 1979-01-01
Надеюсь, это полезно! Приветствия!
Примечание. Для преобразования временной переменной netcdf в объект времени R убедитесь, что у вас установлена библиотека ncdf.tools.