การค้นหารูปทรงเรขาคณิตที่ใกล้เคียงที่สุดใน PostGIS

ฉันได้ดูฟังก์ชั่น PostGIS "API" และฉันสังเกตว่าส่วนใหญ่ใช้องค์ประกอบสองอย่างเพื่อเปรียบเทียบ ตัวอย่างเช่นฟังก์ชัน ST_Distance ใช้องค์ประกอบทางเรขาคณิต / ภูมิศาสตร์สององค์ประกอบเพื่อค้นหาระยะทาง

ไม่มีฟังก์ชั่นที่จะทำสิ่งใดเช่น: "ให้รูปทรงเรขาคณิต G ให้รูปทรงเรขาคณิตที่ใกล้ที่สุดกับ GClosest ในตาราง T ที่ 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 ได้เขียนฟังก์ชั่น Neighbor 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;

... เลือกโหนดที่ใกล้ที่สุดในตาราง nw_node สำหรับทุกรายการใน my_point_table

นอกจากนี้ยังมีฟังก์ชั่นทั่วไปมากขึ้นบนเว็บไซต์ GIS บอสตัน