GRIB2到PostGIS的光栅 - 任何人得到这个工作？(Grib2 to PostGIS ras

我有我需要进口美国国家气象局分析表面，分布作为应用GRIB2文件。我想拉那些到了PostGIS 2.0栅格，做一些计算和建模，并显示在GeoServer的数据和模型结果。

由于GRIB2是GDAL支持的格式，提供的raster2pgsql工具应该能够正确啜一GRIB2到PostGIS的兼容SQL，一旦它的存在，的GeoServer应该能够处理它。不过，我运行到其没有明显的解决方案的问题 - 不是很明显对我来说，无论如何！ Raster2pgsql运行，显然没有错误，产生SQL和运行SQL创建看起来非常像一个光栅。但是GeoServer的不能显示它 - 的界限，特别是出来找怪（0,0，-1）和“预览层”只是抛出一个NullPointerException异常。

有没有人在这条路上了吗？我有问题，基本不知道SRID应该是什么样的数据（4326，也许？）。我不希望任何人来调试我的问题，但对我来说，如果有人已经得到了这个工具链的工作，或它的一部分，我可以插入已知好的东西，看看我能发现。

TIA，

更新时间：每迈克，这里是从文件中的一个坐标系统的东西; 我省略掉其他749个乐队从“gdalinfo”的输出。需要注意的是，文件名是不同的 - 我发现通过运行在我的原文件中“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]

Answer 1:

我希望这可以帮助，至少那些正在添加此主题。

记住，GeoServer的，同时能够从PostGIS的加载栅格数据，默认PostGIS的“进口”模块仅适用于矢量数据，这就是为什么你获得这些奇怪的界限（-1 -1 0 0）。

你必须ImageMosaicJDBC插件添加到您的GeoServer安装，请按照步骤在这里！

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

Answer 2:

有一个很好的回答了我的问题在这里。把它作为一个单独的答案。

他建议使用gdalwarp到GRIB2文件拉入一个已知的SRID，从而：

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

然后，raster2pgsql工作得很好，如

raster2pgsql -M -a 4326_file.grib2 some_sql.sql