PostGIS에서 가장 가까운 지오메트리 찾기

함수의 PostGIS "API"를 살펴 보았는데 대부분 비교하기 위해 두 가지 요소가 필요하다는 것을 알았습니다. 예를 들어 ST_Distance 함수는 거리를 찾기 위해 두 개의 지오메트리 / 지리 요소를 사용합니다.

"지오메트리 G를 지정하고 G.id <> GClosest.id 인 표 T에서 가장 가까운 지오메트리 G를 제공하십시오."와 같은 기능은 없습니다.

PL / PgSQL 함수를 작성하여 테이블을 반복하고 각 요소에서 ST_Distance를 호출 할 수 있다는 것을 알고 있지만 더 좋고 효율적인 솔루션이 있기를 바랍니다.

질문은 전체 레코드와 참조 지오메트리까지의 거리를 반환하는 다음과 같은 단일 (복잡하지만) 쿼리로도 대답 할 수 있습니다. 둘 이상의 레코드가 최소 거리와 일치하면 모두 반환됩니다.

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;

주소 테이블 에서이 쿼리를 테스트했으며 작동합니다. 위의 쿼리에서 id = 3 인 가장 가까운 지점을 찾고 있습니다.

George MacKerron은 내가 찾은 간단한 Nearest Neighbor 함수 를 작성했습니다 . 이 함수는 주어진 지형지 물에 가장 가까운 이웃의 ID를 반환합니다.

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;

사용 예 :

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

... my_point_table의 모든 항목에 대해 nw_node 테이블에서 가장 가까운 노드를 선택합니다.

Boston GIS 사이트 에는 보다 일반적인 기능 도 있습니다 .