Использование R для открытия файлов Grib

Question

Я использую R для работы с метеорологическими данными.Я продолжаю в два этапа:

1. преобразовать grib в netcdf, используя функцию командной строки ncl_convert2nc из языка команд ncar
2. , используйте пакет ncdf в R для импорта данных netcdf.

У меня все еще есть одна проблема:

2 - Для некоторых конкретных файлов grib преобразование с помощью инструмента ncar не работает.Есть ли другие способы или хитрости (кроме транскрипции в netcdf) для чтения файлов grib в R ?

Проблема, на которую ответил Дирк: 1- Я хотел бы сделать автоматическую обработку многих файлов в R. Могу ли я вызвать ncl_convert2nc в R?(отвечает Дирк Эддельбюттель ниже)

Answer 1

По вопросу 1 ответ «Да» - см. help(system) и вариант internal=TRUE, если вы хотите получить результаты.

Answer 2

rgdal также может это сделать, но он менее гибок и требует большей тщательности и детализации, чем ncdf или RNetCDF - и зависит от вашего GDAL / rgdal, включая драйвер GRIB.

Answer 3

ncl_convert2nc кажется лучшим решением.Однако, если структура данных немного сложнее, я использую GrADS для преобразования файла GRIB в ASCII (например, .csv), а затем можно создать файл NetCDF, используя пакет ncdf4, выделенный для R. GrADS также предоставляет поддержкузапись GRIB в NetCDF, но есть ограничение только на 1 переменную.

Answer 4

Другое быстрое и простое решение для командной строки - использовать cdo для преобразования:

cdo -f nc copy file.grb file.nc

Если вы хотите вывести файл netcdf4, укажите "-f nc4".

Answer 5

Другое решение состоит в том, чтобы использовать программное обеспечение wgrib / wgrib2 (http://www.cpc.ncep.noaa.gov/products/wesley/wgrib2/) и выгрузить файл GRIB-1 / GRIB-2 непосредственно в формат CSV, например:

/path/to/your/wgrib2 input_file.grb -csv output_file.csv

. Тогда оно можетчитать прямо в R ...