Usando R para abrir arquivos GRIB

Question

Estou usando R para trabalhar com dados meteorológicos. Eu prossegui em duas etapas:

1. Converta GRIB em NETCDF usando a função de linha de comando ncl_convert2nc da linguagem de comando ncar
2. Use o pacote NCDF no R para importar os dados do NetCDF.

Eu ainda tenho um problema:

2- Para alguns arquivos GRIB específicos, a conversão com a ferramenta NCAR não funciona. Existe outras maneiras ou truques (exceto a transcrição para netcdf) para Leia os arquivos GRIB em R ?

Problema respondido por Dirk: 1- Gostaria de fazer tratamento automático de muitos arquivos dentro de R. Posso ligar para NCL_CONVERT2NC dentro de r? (respondido por Dirk Eddelbuettel abaixo)

Answer 1

Em relação à pergunta 1, a resposta é 'sim' - veja help(system) e a internal=TRUE opção se você deseja capturar resultados.

Answer 2

O RGDAL também pode fazê -lo, mas é menos flexível e requer mais cuidados e detalhes que o NCDF ou RNETCDF - e depende do seu GDAL/RGDAL construído, incluindo o driver GRIB.

Answer 3

NCL_CONVERT2NC parece ser a melhor solução. No entanto, se a estrutura dos dados for um pouco mais complicada, eu uso graduados para converter o arquivo GRIB em ASCII (por exemplo .CSV) e, em seguida, é possível criar um arquivo NETCDF usando o pacote NCDF4 dedicado para R. grads também fornece suporte para recuperação Escrevendo GRIB para NetCDF, mas há limitação para apenas 1 variável.

Answer 4

Outra solução é usar o software WGRIB/WGRIB2 (http://www.cpc.ncep.noaa.gov/products/wesley/wgrib2/) e despeje seu arquivo GRIB-1/GRIB-2 diretamente no formato CSV, por exemplo:

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

Então pode ser lido diretamente em r ...

Answer 5

Outra solução rápida e fácil da linha de comando é usar o CDO para fazer a conversão:

cdo -f nc copy file.grb file.nc

Se você deseja produzir um arquivo netcdf4 que você especificar "-f nc4".