Как поместить файл netCDF в поле в R? - PullRequest
Новая
       76

Как поместить файл netCDF в поле в R?

0 голосов
/ 03 марта 2020

У меня есть файл netCDF, и я хочу поместить его в блок, в котором у меня есть координаты блока. Я хочу использовать индексы файла netCDF для подмножества файла. Я использую пакет ncdf4. Раньше в старой версии R была функция n c .get.var.subset.by.axes, которая выполняла бы эту работу. Но теперь, в более новой версии, это было удалено. Вот что я пробовал: 

library(ncdf4)

ncfile <- nc_open(myfile)

LonIdx <- which( ncfile$dim$lon$vals >= x1 & ncfile$dim$lon$vals <= x2)
LatIdx <- which( ncfile$dim$lat$vals >= y1 & ncfile$dim$lat$vals <= y2)

data <- ncvar_get( ncfile, "var1", subset(ncfile, list(X = LonIdx,Y = LatIdx)))

Но похоже, что это не работает. Можете ли вы придумать что-то, что работает так? Мы использовали это в более старой версии R: 

data <- nc.get.var.subset.by.axes(ncfile, "var1", axis.indices = list(X = LonIdx,Y = LatIdx))
...