在PostGIS中找到最接近的几何

我一直在研究PostGIS的“ API”功能，并且我注意到它们中的大多数需要两个元素进行比较。例如，ST_Distance函数采用两个几何/地理元素来查找距离。

没有函数可以执行以下操作：“给定几何G，在表T中给我最接近的几何GClosest，其中G.id <> GClosest.id”

我意识到我可以编写一个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表中最近的节点。

波士顿GIS站点上还有一个更通用的功能。