SQL-запит, щоб мати повну функцію geojson від PostGIS?

Я хотів би отримати функцію geojson з властивостями від PostGIS. Я знайшов приклад наявності колекції функцій, але не можу змусити її працювати лише для функції.

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;

поки що я намагався змінити запит колекції функцій на прикладі. але вихід недійсний.

Це можна зробити трохи простіше за допомогою json_build_objectPostgreSQL 9.4+, що дозволяє створити JSON, надаючи змінні аргументи ключ / значення. Наприклад:

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;

Все стає ще краще в PostgreSQL 9.5+, де для jsonbтипу даних додаються нові оператори ( docs ). Це полегшує налаштування об’єкта "властивості", який містить усе, крім id та геометрії .

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;

Хочете зробити FeatureCollection? Просто оберніть все це 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;

Цю відповідь можна використати у версії PostgreSQL перед 9.4. Використовуйте відповідь dbaston для PostgreSQL 9.4+

Запит такий: (де 'GEOM'поле геометрії, idполе для включення у властивості json, shapefile_featureназва таблиці та 489445ідентифікатор потрібної функції)

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;

вихід:

{
   "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
   }
}

Лише невелика корекція відповіді dbaston (я б прокоментував, але у мене немає балів) Вам потрібно призначити вихід ST_AsGeoJSON як json ( ::jsonріч):

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;

Інакше член геометрії буде рядком. Це не дійсно GeoJSON

@ Dbaston в відповідь був змінений в останній час по @John Пауелл ака Барсою, і він виробляє неприпустимі geojsons на моєму кінці. Як змінено, агрегація за функціями повертає кожну функцію, вкладену в об'єкт json, який є недійсним.

У мене немає репутації коментувати відповідь безпосередньо, але остаточний jsonb_agg повинен знаходитись у стовпці "особливість", а не в підзапиті "особливості". Агрегуючи по імені стовпця (або "Features.feature", якщо ви вважаєте його акуратнішим), кожен елемент проставляє в масив "особливості" після агрегації, що є правильним шляхом.

Отже, наступне, що дуже схоже на відповідь @ dbaston, як це було до декількох тижнів тому назад (плюс виправлення @Jonh Powell для іменування підзаписів) працює:

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;