PostGIS'de en yakın Geometri'yi bulma

PostGIS "API" fonksiyonlarını inceledim ve çoğunun karşılaştırmak için iki öğe aldığını fark ettim. Örneğin, ST_Distance işlevi mesafeyi bulmak için iki geometri / coğrafya öğesi alır.

"G geometrisi verildiğinde bana Tablo T'de en yakın geometri GClosest'i verin, burada G.id <> GClosest.id" ifadesini kullanın.

Ben tablo üzerinde yineleme ve her öğe üzerinde ST_Distance aramak için bir PL / PgSQL işlevi yazabilirim farkında, ama daha iyi, daha verimli bir çözüm var umuyorum.

Sorunuz ayrıca, tüm kaydı ve referans geometriye olan mesafeyi döndüren aşağıdaki gibi tek (karmaşık da olsa) bir sorgu ile yanıtlanabilir. Birden fazla kayıt minimum mesafeyle eşleşirse bunların hepsinin döndürüldüğünü lütfen unutmayın .

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;

Bu sorgu adresler tablosunda test ettik ve çalışıyor. Yukarıdaki sorguda ben id = 3 ile en yakın noktayı arıyorum.

George MacKerron, oldukça kullanışlı bulduğum basit bir En Yakın Komşu işlevi yazdı . Bu işlev, verilen bir özelliğe en yakın komşunun kimliğini döndürür:

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;

Kullanım örneği:

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

... my_point_table içindeki her giriş için nw_node tablosundaki en yakın düğümü seçer.

Boston GIS sitesinde daha genel bir işlev de var .