Tìm hình học gần nhất trong PostGIS

Tôi đã xem qua "API" các chức năng của PostGIS và tôi nhận thấy rằng hầu hết chúng đều có hai yếu tố để so sánh. Ví dụ: hàm ST_Distance lấy hai thành phần hình học / địa lý để tìm khoảng cách.

Không có chức năng để làm một cái gì đó như: "Đưa ra một hình học G, cho tôi hình học GClosest gần nhất trong Bảng T trong đó G.id <> GClosest.id"

Tôi nhận ra rằng tôi có thể viết một hàm PL / PGQuery để lặp lại bảng và gọi ST_Distance trên mỗi phần tử, nhưng tôi hy vọng có một giải pháp tốt hơn, hiệu quả hơn.

Câu hỏi của bạn cũng có thể được trả lời bằng một truy vấn (mặc dù phức tạp) như sau đây trả về toàn bộ bản ghi và khoảng cách đến hình học tham chiếu. Xin lưu ý rằng nếu có nhiều hơn một bản ghi khớp với khoảng cách tối thiểu thì tất cả chúng sẽ được trả về.

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;

Tôi đã kiểm tra truy vấn này trên bảng địa chỉ và nó hoạt động. Trong truy vấn trên tôi đang tìm kiếm điểm gần nhất với id = 3.

George MacKerron đã viết một hàm Neighbor Neighbor đơn giản mà tôi thấy khá hữu ích. Hàm này trả về ID của hàng xóm gần nhất với một tính năng nhất định:

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;

Ví dụ sử dụng:

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

... Chọn nút gần nhất trong bảng nw_node cho mỗi mục trong my_point_table.

Ngoài ra còn có một chức năng chung hơn trên trang web Boston GIS .