Trovare la geometria più vicina in PostGIS

Ho esaminato le "API" PostGIS delle funzioni e noto che la maggior parte di esse richiede due elementi per il confronto. Ad esempio, la funzione ST_Distance accetta due elementi geometria / geografia per trovare la distanza.

Non esiste una funzione per fare qualcosa del tipo: "Data una geometria G, dammi la geometria più vicina GClosest nella tabella T dove G.id <> GClosest.id"

Mi rendo conto di poter scrivere una funzione PL / PgSQL per iterare sul tavolo e chiamare ST_Distance su ciascun elemento, ma spero che ci sia una soluzione migliore, più efficiente.

Alla tua domanda si può anche rispondere con una singola query (anche se complessa) come la seguente che restituisce l'intero record e la distanza alla geometria di riferimento. Se più di un record corrisponde alla distanza minima, vengono tutti restituiti.

SELECT 
  i.*,
  md.min_distance
FROM
  address AS i, 
  (SELECT 
     ga.address_geom,
     min( ST_Distance(
            ga.address_geom,
            gb.address_geom)
        ) AS min_distance
   FROM
     address AS ga,
     address AS gb 
   WHERE 
     ga.id <> gb.id 
   AND 
     ga.id = 3
   GROUP BY 
     ga.address_geom
  ) AS md 
WHERE 
  ST_Distance( i.address_geom, md.address_geom) = md.min_distance;

Ho testato questa query sulla tabella degli indirizzi e funziona. Nella query sopra sto cercando il punto più vicino a quello con id = 3.

George MacKerron ha scritto una semplice funzione di vicinato più vicino che ho trovato abbastanza utile. Questa funzione restituisce l'ID del vicino più vicino a una determinata funzione:

create or replace function 
  nn(nearTo                   geometry
   , initialDistance          real
   , distanceMultiplier       real 
   , maxPower                 integer
   , nearThings               text
   , nearThingsIdField        text
   , nearThingsGeometryField  text)
 returns integer as $$
declare 
  sql     text;
  result  integer;
begin
  sql := ' select ' || quote_ident(nearThingsIdField) 
      || ' from '   || quote_ident(nearThings)
      || ' where st_dwithin($1, ' 
      ||   quote_ident(nearThingsGeometryField) || ', $2 * ($3 ^ $4))'
      || ' order by st_distance($1, ' || quote_ident(nearThingsGeometryField) || ')'
      || ' limit 1';
  for i in 0..maxPower loop
     execute sql into result using nearTo              -- $1
                                , initialDistance     -- $2
                                , distanceMultiplier  -- $3
                                , i;                  -- $4
    if result is not null then return result; end if;
  end loop;
  return null;
end
$$ language 'plpgsql' stable;

Esempio di utilizzo:

SELECT id, nn(pt_geom,0.00001,2,100,'nw_node','node_id','node_geom') FROM my_point_table;

... seleziona il nodo più vicino nella tabella nw_node per ogni voce in my_point_table.

Esiste anche una funzione più generica sul sito GIS di Boston .