Гриб2 для растра PostGIS у кого-нибудь заставить это работать?

У меня есть приложение, для которого мне нужно импортировать анализы поверхности Национальной метеорологической службы США, которые распространяются как grib2 файлы. Я хочу перенести их в растры PostGIS 2.0, выполнить некоторые вычисления и моделирование, а также отобразить данные и результаты моделирования в GeoServer.

Поскольку grib2 является форматом, поддерживаемым GDAL, поставляемая утилита raster2pgsql должна иметь возможность захлебнуть grib2 прямо в PostGIS-совместимый SQL, и как только он появится, GeoServer сможет с этим справиться. Однако я сталкиваюсь с проблемами, у которых нет очевидных решений - во всяком случае, не очевидных для меня! Raster2pgsql запускается, по-видимому, без ошибок, создавая SQL, а запуск SQL создает то, что очень похоже на растр. Но GeoServer не может его отобразить - границы, в частности, выглядят странно (0,0 -1, -1), а «слой предварительного просмотра» просто выдает исключение NullPointerException.

Кто-нибудь уже прошел по этому пути? У меня есть такие простые проблемы, как незнание, каким должен быть SRID для данных (возможно, 4326?). Я не ожидаю, что кто-то отладит мои проблемы за меня, но если у кого-то уже есть эта цепочка инструментов или ее часть, я могу подключить заведомо хорошие вещи и посмотреть, что я могу обнаружить.

TIA,

rw

Обновлено: Пер Майк, вот информация о системе координат из одного из файлов; Я исключил остальные 749 полос в выводе "gdalinfo". Обратите внимание, что имя файла другое - я обнаружил, запустив "gdalinfo" в моем исходном файле, что с ним что-то не так, gdalinfo не может его прочитать. Новый (35 МБ!) Файл здесь.

Вывод 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]

postgresql postgis gis geoserver grib
person Rick Wayne    schedule 17.05.2012    source источник
comment
не могли бы вы запустить gdalinfo your.grb на растре и обновить свой вопрос? Вы можете перенести это на gis.SE   -  person Mike T    schedule 18.05.2012
comment
Майк, спасибо за ответ! Я действительно поместил это на gis.se, стесняюсь сказать. Но я клянусь, что буду обрабатывать ОБЕИ темы, чтобы никто не ответил на вопрос, на который уже был дан ответ! Итак, мой следующий шаг - пролистать списки EPSG в поисках чего-то подходящего, чтобы я мог ввести SRID для raster2pgsql? Не знаете, как перепроецировать файл GRIB2! TIA ^ 2!   -  person Rick Wayne    schedule 19.05.2012

Ответы (2)


arrow_upward
1
arrow_downward

Я надеюсь, что это поможет, по крайней мере, тем, кто зашел в эту ветку.

Имейте в виду, что GeoServer, хотя и способен загружать растровые данные из PostGIS, модуль «импорта» PostGIS по умолчанию доступен ТОЛЬКО для векторных данных, поэтому вы получаете эти нечетные границы (-1 -1 0 0).

Вам нужно будет добавить плагин ImageMosaicJDBC к вашей установке геосервера, выполните следующие действия!

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

person Mikel G. Gainza    schedule 19.10.2012

arrow_upward
0
arrow_downward

Получил отличный ответ на мою проблему здесь . Выкладываю как отдельный ответ.

Он рекомендовал использовать gdalwarp для извлечения файла GRIB2 в известный SRID, таким образом:

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

Тогда raster2pgsql работает нормально, например.

raster2pgsql -M -a 4326_file.grib2 some_sql.sql
person Rick Wayne    schedule 27.09.2013