Хранить GeoJSON FeatureCollection в PostgreSQL с помощью PostGIS?

21

Я новичок в GeoJSON. У меня есть коллекция объектов GeoJSON, как показано, и я хотел бы сохранить ее в таблице postgres (testtable). Моя таблица postgres имеет серийный идентификатор и столбец геометрии.

{
    "type": "FeatureCollection",
    "features": [
        {
            "type": "Feature",
            "geometry": {
                "type": "Point",
                "coordinates": [
                    2565453.1826721914,
                    -3835048.659760314
                ]
            }
        },
        {
            "type": "Feature",
            "geometry": {
                "type": "LineString",
                "coordinates": [
                    [
                        2727584.7219710173,
                        -3713449.1942418693
                    ],
                    [
                        2732476.691781269,
                        -3992291.473426192
                    ]
                ]
            }
        },
        {
            "type": "Feature",
            "geometry": {
                "type": "Polygon",
                "coordinates": [
                    [
                        [
                            2442627.9025405287,
                            -3705499.954308534
                        ],
                        [
                            2425506.008204649,
                            -3886502.837287831
                        ],
                        [
                            2425506.008204649,
                            -3886502.837287831
                        ],
                        [
                            2555143.2081763083,
                            -3910962.686339088
                        ],
                        [
                            2442627.9025405287,
                            -3705499.954308534
                        ]
                    ]
                ]
            }
        }
    ]
}

Я хотел бы вставить данные GeoJSON в таблицу тестов.

Как мне это сделать?

Я использую Postgres версии 9.3.5 с Postgis версии 2.1.3


Я был направлен на ранее заданные вопросы, которые отвечают, как хранить одну особенность, например, точку или многоугольник. Мой вопрос спрашивает, как сохранить несколько объектов в файле GeoJSON. Под несколькими объектами я подразумеваю сочетание точек, линий и типов объектов полигонов в одном файле.

сойка
источник
Возможный дубликат Как вставить полигон GeoJSON в таблицу PostGIS?
Рикардо Оливейра
Привет Рикардо, я видел этот вопрос, но он не решает мою проблему. Я хочу сохранить список функций, а не только один тип функций. Пожалуйста, взгляните на мою коллекцию возможностей GeoJSON в моем вопросе.
Джей
@Jay Так что теперь или ваш вопрос: «Как я могу разбить коллекцию геоджонов на отдельные объекты, или вам нужно добавить больше информации (может быть, хранить информацию о том, что эти геометрии принадлежат какой-то коллекции?)
Jakub Kania
1
Спасибо @Джон за ваш ответ. Поскольку я новичок в ГИС и GeoJSON, я хотел несколько указателей с моей проблемой под рукой. Предпосылки вопроса: пользователь рисует объекты на карте, а я фиксирую коллекцию нарисованных объектов. Я хотел бы сохранить эту коллекцию в БД с уникальным идентификатором. Позже сохраненные данные могут быть получены для предоставленного идентификатора. Тестируемый в postgres имеет 2 столбца. столбец gid, являющийся последовательным типом, для хранения идентификатора, и столбец geom, имеющий геометрический тип.
Джей
1
@Jay Да, вы можете хранить JSON, но тогда это не будет геометрия, поэтому вы не сможете легко запросить ближайшего соседа и т. Д.
Якуб Кания

Ответы:

26

Предполагая, что у вас есть по крайней мере PostgreSQL версии 9.3, вы можете использовать несколько функций и операторов JSON для извлечения соответствующих частей спецификации GeoJSON, требуемой ST_GeomFromGeoJSON для создания геометрий.

Попробуйте следующее, где вы можете заменить JSON в верхней части:

WITH data AS (SELECT '{ "type": "FeatureCollection",
    "features": [
      { "type": "Feature",
        "geometry": {"type": "Point", "coordinates": [102.0, 0.5]},
        "properties": {"prop0": "value0"}
        },
      { "type": "Feature",
        "geometry": {
          "type": "LineString",
          "coordinates": [
            [102.0, 0.0], [103.0, 1.0], [104.0, 0.0], [105.0, 1.0]
            ]
          },
        "properties": {
          "prop0": "value0",
          "prop1": 0.0
          }
        },
      { "type": "Feature",
         "geometry": {
           "type": "Polygon",
           "coordinates": [
             [ [100.0, 0.0], [101.0, 0.0], [101.0, 1.0],
               [100.0, 1.0], [100.0, 0.0] ]
             ]
         },
         "properties": {
           "prop0": "value0",
           "prop1": {"this": "that"}
           }
         }
       ]
     }'::json AS fc)

SELECT
  row_number() OVER () AS gid,
  ST_AsText(ST_GeomFromGeoJSON(feat->>'geometry')) AS geom,
  feat->'properties' AS properties
FROM (
  SELECT json_array_elements(fc->'features') AS feat
  FROM data
) AS f;

Находит три геометрии. geomСтолбец имеет геометрический объект, и gidэто номер функции. В ST_AsTextфункции показывает WKT эквивалент каждой геометрии. Я также включил propertiesатрибуты or, которые могут быть определены для каждой геометрии, как показано в спецификации.

 gid |                   geom                   |              properties
-----+------------------------------------------+--------------------------------------
   1 | POINT(102 0.5)                           | {"prop0": "value0"}
   2 | LINESTRING(102 0,103 1,104 0,105 1)      | {                                   +
     |                                          |           "prop0": "value0",        +
     |                                          |           "prop1": 0.0              +
     |                                          |           }
   3 | POLYGON((100 0,101 0,101 1,100 1,100 0)) | {                                   +
     |                                          |            "prop0": "value0",       +
     |                                          |            "prop1": {"this": "that"}+
     |                                          |            }
(3 rows)

Вы должны назначить SRID для геометрии, используя ST_SetSRID.

Или, если вам просто нужна одна гетерогенная GEOMETRYCOLLECTION, вы можете сделать ее такой компактной:

SELECT ST_AsText(ST_Collect(ST_GeomFromGeoJSON(feat->>'geometry')))
FROM (
  SELECT json_array_elements('{ ... put JSON here ... }'::json->'features') AS feat
) AS f;

GEOMETRYCOLLECTION(POINT(2565453.18267219 -3835048.65976031),LINESTRING(2727584.72197102 -3713449.19424187,2732476.69178127 -3992291.47342619),POLYGON((2442627.90254053 -3705499.95430853,2425506.00820465 -3886502.83728783,2555143.20817631 -3910962.68633909,2442627.90254053 -3705499.95430853)))

См. Также Создание коллекций объектов GeoJSON с помощью функций JSON и PostGIS из журнала Postgres OnLine, который делает обратное.

Майк Т
источник