問題タブ [grib]

R を使用して気象データを操作しています。私は2つのステップで進みます：

1. ncar コマンド言語のコマンド ライン関数 ncl_convert2nc を使用して、grib を netcdf に変換します。
2. R のパッケージ ncdf を使用して、netcdf データをインポートします。

まだ問題が 1 つあります。

2- 一部の特定の grib ファイルでは、ncar ツールによる変換が機能しません。R で grib ファイルを読み取る他の方法やトリック (netcdf への転写以外) はありますか?

Dirk が回答した問題: 1- R 内で多数のファイルを自動的に処理したいと考えています。R 内で ncl_convert2nc を呼び出すことはできますか? （以下のDirk Eddelbuettelによる回答）

ここで気象学や地球科学の仕事をしている人はいますか？GRIB2形式のデータファイルを操作しようとしています。利用可能なライブラリ、特にUnidata GRIBJavaDecoderがあります。そして、大きな線形配列でデータを抽出する方法を正確に理解できますが、私がやりたいのは、緯度/経度で値にアクセスすることです。それを行う方法についての簡単な例を見つけることができません。提案？

TIA！

R 関数を使用して、ECMWF ERA-40 Web サイトから波高を含む GRIB ファイルwavedata.gribを読み込もうとしています。これまでの私のソースコードは次のとおりです。

wavedata.gribファイルには、1957 年 9 月 1 日から 2002 年 8 月 31 日までの期間の SWH 値がグリッド化されています。各バンドは緯度/経度のペアを参照し、各日の 00 時に一連の 16346 SWH 値を持ちます (1 バンド = 特定の緯度/経度での 16346 値)。

myframe のサイズは 16438 x 91 です。91 = 7 行 x 13 列であることに注意してください。そして、16438という数字はバンド数とほぼ同じです。追加の 2 つの行/バンドは経度と緯度の値で、他のすべての列は 16436 バンドに対応する波高でなければなりません。

問題は、緯度/経度 = 43.75,331.25 で SWH (波の高さ) を抽出したいのですが、同じ緯度/経度で Degrib ユーティリティを使用してファイルを読み取った値と一致しません。

また、必要な正しい値 (2.147、2.159、1.931、...) は myframe の列 3 ではなく列 4 にあります。 331.25 (長い)。どうしてこれなの？myframe[i,j] の値が実際にどの緯度/経度に対応しているか、またはプロセス中にデータ インポート エラーがあるかどうかを知りたいです。Degrib にはエラーがないと仮定しています。

グリッドポイント間の値を抽出したい場合、マトリックス内の値を簡単に補間する R ルーチンはありますか? より一般的には、次のような波の高さを抽出する効果的な R 関数を作成する際に助けが必要です。

助けてください。

GRIBデータを処理するための .NET ライブラリを知っている人はいますか?

GRIBファイルを処理するphpライブラリがあるかどうかは誰にもわかりませんか? Environment Canada で入手できる GRIB (または GRIB 2) ファイルを使用して、自分の Web サイトに世界の天気を表示したいと考えています。私は開発にphpを使用しているため、これを探しています。ありがとう....

grib2ファイルとして配布されている米国国立気象局の表面分析をインポートする必要があるアプリケーションがあります。それらをPostGIS2.0ラスターに取り込み、計算とモデリングを行い、データとモデルの結果をGeoServerに表示したいと思います。

grib2はGDALでサポートされている形式であるため、提供されているraster2pgsqlユーティリティはgrib2をPostGIS互換のSQLに直接スラップできるはずであり、そこに到達すると、GeoServerはそれを処理できるはずです。しかし、私は明白な解決策がない問題に直面しています-とにかく、私には明白ではありません！Raster2pgsqlは、明らかにエラーなしで実行され、SQLを生成し、SQLを実行すると、ラスターに非常によく似たものが作成されます。しかし、GeoServerはそれを表示できません。特に境界は奇妙に見え（0,0 -1、-1）、「プレビューレイヤー」はNullPointerExceptionをスローするだけです。

誰かがすでにこの道を進んでいますか？データのSRIDがどうあるべきかわからないのと同じくらい基本的な問題があります（おそらく4326？）。誰かが私の問題をデバッグすることは期待していませんが、誰かがすでにこのツールチェーンまたはその一部を機能させている場合は、既知の良いものをプラグインして、何を発見できるかを確認できます。

TIA、

rw

更新：Mikeによると、これがファイルの1つからの座標系のものです。「gdalinfo」からの出力で他の749バンドを削除しました。ファイル名が異なることに注意してください。元のファイルで「gdalinfo」を実行すると、何か問題があることがわかりました。gdalinfoはそれを読み取ることができませんでした。ここに新しい（35MB！）ファイルがあります。

Gdalinfo出力：

「高解像度初期条件」の気候予測データを 1 日分ダウンロードしました。拡張子は .tar.gz だったので、ローカル ディレクトリに展開すると、添付の画像のようなファイルが得られます。拡張子のないファイルはGRIBデータだと思います（最初の単語が「GRIB」であるため）。そのため、大きなファイル (グリッド内の温度や圧力などの気候データを含む GRIB および NetCDF 形式) からデータベースにデータを取得したいのですが、それらはバイナリです。これらのファイルからデータを取得する簡単な方法を教えてください。彼らのウェブサイトでのデータセットの取り扱いに関する情報を得ることができません。

これらのファイルを .csv に変換できればよいのですが、GRIB ファイルを変換するプログラムが見つかりません。

全世界の情報 (一部のパラメーター) を含む grib2 ファイルがあり、ユーザー(クライアント ソフトウェアからサーバー)によって指定された緯度と経度に基づいて wgrib2 を使用してデータを抽出したい場合。次のコマンドを試しましたが、完全な grib2 ファイルのみを取得しています:

wgrib2.exe input.grb -アウトボックスの定義を解除 10:90 -10:10 -grib output.grb

どこが間違っているのか教えてください。ありがとう。

私はwpfアプリケーションを開発しています。次のコードを使用してdegrib.exeを実行しています

degrib.exeは、gribファイルからcsvファイルを生成します。これがdegribhttp://www.nws.noaa.gov/mdl/degrib/index.phpの説明です。上記の関数fileNameには、gribファイルのパスがあります。上記の関数は、gribファイルがプログラムファイル以外のフォルダーに配置されている場合、gribファイルからcsvファイルを作成します。gribファイルをProgramFilesまたはProgramFilesの他のフォルダーに配置した場合、同じ関数はcsvファイルを生成しません。同じ関数がAppDataフォルダーでも機能しません。なぜこれが起こっているのですか？上記の問題の解決策を教えてください。

このgribファイルを読み込もうとしていますが、エラーが発生しました。この問題を解決する方法はありますか