4

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出力：

Driver: GRIB/GRIdded Binary (.grb)
Files: ruc2.t00z.bgrb13anl.grib2
Size is 451, 337
Coordinate System is:
PROJCS["unnamed",
    GEOGCS["Coordinate System imported from GRIB file",
        DATUM["unknown",
            SPHEROID["Sphere",6371229,0]],
        PRIMEM["Greenwich",0],
        UNIT["degree",0.0174532925199433]],
    PROJECTION["Lambert_Conformal_Conic_2SP"],
    PARAMETER["standard_parallel_1",25],
    PARAMETER["standard_parallel_2",25],
    PARAMETER["latitude_of_origin",0],
    PARAMETER["central_meridian",265],
    PARAMETER["false_easting",0],
    PARAMETER["false_northing",0]]
Origin = (-3332155.288903323933482,6830293.833488883450627)
Pixel Size = (13545.000000000000000,-13545.000000000000000)
Corner Coordinates:
Upper Left  (-3332155.289, 6830293.833) (139d51'22.04"W, 54d10'20.71"N)
Lower Left  (-3332155.289, 2265628.833) (126d 6'34.06"W, 16d 9'49.48"N)
Upper Right ( 2776639.711, 6830293.833) ( 57d12'21.76"W, 55d27'10.73"N)
Lower Right ( 2776639.711, 2265628.833) ( 68d56'16.73"W, 17d11'55.33"N)
Center      ( -277757.789, 4547961.333) ( 98d 8'30.73"W, 39d54'5.40"N)
Band 1 Block=451x1 Type=Float64, ColorInterp=Undefined
  Description = 1[-] HYBL="Hybrid level"
  Metadata:
    GRIB_UNIT=[Pa]
    GRIB_COMMENT=Pressure [Pa]
    GRIB_ELEMENT=PRES
[Etc., Etc., for all 750 bands]
4

2 に答える 2

1

少なくともこのスレッドに参加する人たちにとって、これがお役に立てば幸いです。

GeoServerはPostGISからラスターデータをロードできますが、デフォルトのPostGIS「インポート」モジュールはベクターデータでのみ使用できるため、これらの奇数の境界（-1 -1 0 0）を取得することに注意してください。

ジオサーバーのインストールにImageMosaicJDBCプラグインを追加する必要があります。ここの手順に従ってください！

http://docs.geoserver.org/latest/en/user/tutorials/imagemosaic-jdbc/imagemosaic-jdbc_tutorial.html

于 2012-10-19T11:46:48.920 に答える
0

ここで私の問題に対する優れた答えを得ました。別の答えとしてそれを入れます。

彼は、gdalwarpを使用してGRIB2ファイルを既知のSRIDにプルすることを推奨しました。

gdalwarp -t_srs EPSG:4326 original_file.grib2 4326_file.grib2

次に、raster2pgsqlは問題なく機能します。

raster2pgsql -M -a 4326_file.grib2 some_sql.sql
于 2013-09-27T22:09:29.830 に答える