RNomad для обработки файлов grib2 в R занимает слишком много времени - PullRequest
Новая
       50

RNomad для обработки файлов grib2 в R занимает слишком много времени

0 голосов
/ 03 мая 2018

Мы пытаемся прочитать grib2, содержащий прогноз погоды с 13 переменными на одном уровне давления для 6 лат / лонг в течение одного дня. Файл grib2 составляет 62 МБ. Чтобы прочитать одну переменную почасовых данных за один день, требуется около часа, особенно ReadGrib.

Мы используем wgrib.exe, версия v2.0.3. Мы используем R версии 3.4 в Windows 10.

Если чтение займет много времени или мы делаем что-то не так. Пожалуйста, сообщите.

Код, который мы используем: 

file <- "./xxx.grib2"

### define the geographical coordinated to be extracted
lat_list <- c(36.8200, 43.2800, 52.7200)
lon_list <- c(85.8000, 87.3500, 95.8000)

var_list <- c("APCP", "DSWRF", "PRES")
lev_list = c("surface")

### read data
data <- ReadGrib(file, variable = var_list,  levels = lev_list)

for (i in 1: length(lat_list)) {
   cat("\n", "         ... extracting point",  i,  "/",  length(lat_list))
   data_extract <- BuildProfile(data, spatial.average = F,
                    lat = lat_list[i], lon = lon_list[i])
}

1 Ответ

0 голосов
/ 27 июня 2018

Можно попытаться определить параметр домена в пределах географического интереса: 

domain=c(85.8000,95.8000,52.7200,36.8200)

### read data
data <- ReadGrib(file, variable = var_list,  levels = lev_list,domain=domain)
...