SQL-Abfrage für eine vollständige Geojson-Funktion von PostGIS?

Ich möchte ein Geojson-Feature mit Eigenschaften von PostGIS erhalten. Ich habe ein Beispiel für eine Feature-Sammlung gefunden, kann sie jedoch nicht nur für ein Feature verwenden.

SELECT row_to_json(fc)
 FROM ( SELECT 'FeatureCollection' As type, array_to_json(array_agg(f)) As features
 FROM (SELECT 'Feature' As type
    , ST_AsGeoJSON(lg.geog)::json As geometry
    , row_to_json(lp) As properties
   FROM locations As lg 
         INNER JOIN (SELECT loc_id, loc_name FROM locations) As lp 
       ON lg.loc_id = lp.loc_id  ) As f )  As fc;

Bisher habe ich versucht, die Feature-Collection-Abfrage des Beispiels zu ändern. aber die Ausgabe ist ungültig.

Dies kann mit json_build_objectPostgreSQL 9.4+ ein bisschen einfacher gemacht werden, mit dem Sie einen JSON-Code erstellen können, indem Sie alternative Schlüssel- / Wertargumente angeben. Beispielsweise:

SELECT json_build_object(
    'type',       'Feature',
    'id',         gid,
    'geometry',   ST_AsGeoJSON(geom)::json,
    'properties', json_build_object(
        'feat_type', feat_type,
        'feat_area', ST_Area(geom)::geography
     )
 )
 FROM input_table;

Noch besser wird es in PostgreSQL 9.5+, wo einige neue Operatoren für den jsonbDatentyp ( docs ) hinzugefügt werden . Dies macht es einfach, ein "Eigenschaften" -Objekt einzurichten, das alles außer der ID und der Geometrie enthält .

SELECT jsonb_build_object(
    'type',       'Feature',
    'id',         gid,
    'geometry',   ST_AsGeoJSON(geom)::jsonb,
    'properties', to_jsonb(row) - 'gid' - 'geom'
) FROM (SELECT * FROM input_table) row;

Möchten Sie eine FeatureCollection erstellen? Einfach alles einpacken mit jsonb_agg:

SELECT jsonb_build_object(
    'type',     'FeatureCollection',
    'features', jsonb_agg(features.feature)
)
FROM (
  SELECT jsonb_build_object(
    'type',       'Feature',
    'id',         gid,
    'geometry',   ST_AsGeoJSON(geom)::jsonb,
    'properties', to_jsonb(inputs) - 'gid' - 'geom'
  ) AS feature
  FROM (SELECT * FROM input_table) inputs) features;

Diese Antwort könnte mit der PostgreSQL-Version vor 9.4 verwendet werden. Verwenden Sie die Antwort von dbaston für PostgreSQL 9.4+

Die Abfrage lautet wie folgt: (Wo 'GEOM'ist das Geometriefeld, iddas Feld, das in den JSON-Eigenschaften enthalten sein soll, shapefile_featureder Tabellenname und 489445die ID des gewünschten Features?)

SELECT row_to_json(f) As feature \
     FROM (SELECT 'Feature' As type \
     , ST_AsGeoJSON('GEOM')::json As geometry \
     , row_to_json((SELECT l FROM (SELECT id AS feat_id) As l)) As properties \
     FROM shapefile_feature As l WHERE l.id = 489445) As f;

Ausgabe:

{
   "geometry":{
      "type":"MultiPolygon",
      "coordinates":[
         [
            [
               [
                  -309443.24253826,
                  388111.579584133
               ],
               [
                  -134666.391073443,
                  239616.414560895
               ],
               [
                  -308616.222736376,
                  238788.813082666
               ],
               [
                  -309443.24253826,
                  388111.579584133
               ]
            ]
         ]
      ]
   },
   "type":"Feature",
   "properties":{
      "feat_id":489445
   }
}

Nur eine kleine Korrektur der Antwort von dbaston (ich würde kommentieren, aber ich habe keine Punkte). Sie müssen die Ausgabe von ST_AsGeoJSON als json (das ::jsonDing) umwandeln:

SELECT json_build_object(
  'type',       'Feature',
  'id',         gid,
  'geometry',   ST_AsGeoJSON(geom)::json,
  'properties', json_build_object(
    'feat_type', feat_type,
    'feat_area', ST_Area(geom)::geography
  )
)
FROM input_table;

Andernfalls ist das Geometriemitglied eine Zeichenfolge. Das ist kein gültiger GeoJSON

Die Antwort von @ dbaston wurde kürzlich von @ John Powell aka Barça geändert, und es entstehen auf meiner Seite ungültige Geojsons. In der geänderten Form gibt die Aggregation nach Features jedes Feature zurück, das in einem ungültigen json-Objekt verschachtelt ist.

Ich habe nicht den Ruf, die Antwort direkt zu kommentieren, aber das endgültige jsonb_agg sollte in der Spalte "feature" und nicht in der Unterabfrage "features" stehen. Wenn Sie den Spaltennamen aggregieren (oder "features.feature", wenn Sie ihn übersichtlicher finden), wird jedes Element nach der Aggregation direkt in das Array "features" eingefügt. Dies ist der richtige Weg.

Das Folgende, das der Antwort von @ dbaston bis vor einigen Wochen ziemlich ähnlich ist (plus der Korrektur von @ Jonh Powell für die Benennung von Unterabfragen), funktioniert also:

SELECT jsonb_build_object(
  'type',     'FeatureCollection',
  'features', jsonb_agg(feature)
)
FROM (
  SELECT jsonb_build_object(
    'type',       'Feature',
    'id',         gid,
    'geometry',   ST_AsGeoJSON(geom)::jsonb,
    'properties', to_jsonb(inputs) - 'gid' - 'geom'
  ) AS feature
  FROM (
    SELECT * FROM input_table
  ) inputs
) features;