Используя R, чтобы открыть файлы GRIB
-
27-09-2019 - |
Вопрос
Я использую R для работы с метеорологическими данными. Я продолжаю два шага:
- Конвертировать GRIB на NetCDF с помощью функции командной строки NCL_CONVERT2NC из языка команд NCAR
- Используйте пакет NCDF в R, чтобы импортировать данные NetCDF.
У меня все еще есть одна проблема:
2- Для некоторых конкретных файлов GRIB преобразование с инструментом NCAR не работает. Есть ли другие способы или трюк (кроме транскрипции в NetCDF) для Прочитайте файлы GRIB в R ?
Проблема, отвечающая Дирком: 1- Я хотел бы сделать автоматическую обработку многих файлов в R. Могу ли я позвонить NCL_Convert2nc внутри R? (Ответил Дирк Эдделбуетель ниже)
Решение
Что касается вопроса 1, ответ «да» - см. help(system)
и то internal=TRUE
Вариант Если вы хотите захватить результаты.
Другие советы
RGDal также может это сделать, но менее гибко и требует большего ухода и детализации, чем NCDF или RNETCDF - и зависит от вашего GDAL / RGDAL, включая водитель GRIB.
ncl_convert2nc, кажется, лучшее решение. Однако, если структура данных немного более сложнее, я использую GRADS для преобразования файла GRIB в ASCII (например, .CSV), а затем можно создать файл NetCDF с использованием пакета NCDF4, предназначенный для R. GALDS, также обеспечивает поддержку Написание GRIB на NetCDF, но есть ограничение только 1 переменной.
Другое решение - использовать программное обеспечение 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 ...
Еще одним быстрым и простым решением командной строки является использование CDO для преобразования:
cdo -f nc copy file.grb file.nc
Если вы хотите вывести файл NetCDF4, который вы указываете «-F NC4».