Consulta SQL para ter um recurso geojson completo do PostGIS?

Gostaria de obter um recurso geojson com propriedades do PostGIS. Encontrei um exemplo para ter uma coleção de recursos, mas não posso fazê-lo funcionar apenas para um recurso.

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;

até agora, tentei modificar a consulta de coleção de recursos do exemplo. mas a saída não é válida.

Isso pode ser feito de maneira um pouco mais simples json_build_objectno PostgreSQL 9.4+, que permite criar um JSON fornecendo argumentos alternativos de chave / valor. Por exemplo:

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;

As coisas ficam ainda melhores no PostgreSQL 9.5+, onde alguns novos operadores são adicionados ao jsonbtipo de dados ( docs ). Isso facilita a configuração de um objeto "propriedades" que contém tudo, exceto o ID e a geometria .

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;

Deseja fazer uma FeatureCollection? Basta encerrar tudo com 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;

Esta resposta pode ser usada na versão PostgreSQL anterior à 9.4. Use a resposta do dbaston para o PostgreSQL 9.4+

A consulta é a seguinte: (onde 'GEOM'está o campo de geometria, ido campo a ser incluído nas propriedades json, shapefile_featureo nome da tabela e 489445o ID do recurso desejado)

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;

saída:

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

Apenas uma ligeira correção na resposta do dbaston (gostaria de comentar, mas não tenho pontos). Você precisa converter a saída de ST_AsGeoJSON como json (the ::jsonthingie):

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;

Caso contrário, o membro da geometria será uma sequência. GeoJSON inválido

A resposta do @ dbaston foi modificada ultimamente por @ John Powell, também conhecido como Barça, e produz geojsons inválidos para mim. Conforme modificada, a agregação nos recursos retorna cada recurso aninhado dentro de um objeto json, que é inválido.

Não tenho reputação de comentar diretamente a resposta, mas o jsonb_agg final deve estar na coluna "feature" e não na subconsulta "features". A agregação no nome da coluna (ou "features.feature", se você achar melhor) coloca todos os elementos diretamente na matriz "features" após a agregação, que é o caminho certo a seguir.

Portanto, o seguinte, que é bastante semelhante à resposta do @ dbaston, como era até algumas semanas atrás (mais a correção de @ Jonh Powell para nomear subconsultas) funciona:

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;