Chọn các tính năng KHÔNG giao nhau trong PostGIS

Điều này đối với tôi giống như một câu hỏi đơn giản như vậy (và có lẽ là vậy) nhưng dường như tôi không thể tìm thấy một ví dụ cho tôi câu trả lời. Sử dụng PostGIS, tôi chỉ muốn chọn các điểm nằm ngoài đa giác. Cuối cùng, đây là nghịch đảo của ST_Intersects, theo như tôi có thể thấy.

Ví dụ: Tôi có một lớp taxlot và một lớp điểm địa chỉ. Tôi giả sử tôi nên sử dụng ST_Intersects, nhưng làm cách nào để bảo nó thực hiện lựa chọn ngược? Tôi nghĩ có thể thêm một câu lệnh KHÔNG vào trước đoạn mã bên dưới, nhưng điều đó không hiệu quả.

CREATE table t_intersect AS
SELECT 
  hp.gid, 
  hp.st_address, 
  hp.city, 
  hp.st_num,
  hp.the_geom
FROM 
  public.parcel as par,
  public.housepoints as hp
WHERE 
  ST_Intersects(hp.the_geom,par.the_geom);

Lý do nó không hoạt động với "không giao nhau" là vì bạn chỉ so sánh hình học theo cặp; sẽ có cùng một vấn đề với sự rời rạc. Mỗi housepoint sẽ phân biệt một số bưu kiện ngay cả khi nó giao nhau một bưu kiện.

đề nghị của underdark không có vấn đề đó. Ngoài ra còn có một mẹo khác có thể sẽ giúp sử dụng chỉ mục hiệu quả hơn:

CREATE TABLE t_intersect AS
SELECT 
  hp.gid, 
  hp.st_address, 
  hp.city, 
  hp.st_num,
  hp.the_geom
FROM 
  public.housepoints AS hp LEFT JOIN
  public.parcel AS par ON
  ST_Intersects(hp.the_geom,par.the_geom)
WHERE par.gid IS NULL;

Ý tưởng là tham gia cùng họ với st_intersects và nhận các hàng không có id bưu kiện.

Các chỉ mục cần thiết ở đây là một chỉ mục không gian và một chỉ số về gid tính theo bưu kiện (giả sử rằng id trong bảng bưu kiện cũng được gọi là gid).

Bạn có thể đang tìm kiếm ST_Disjoint

ST_Disjoint - Trả về ĐÚNG nếu Hình học không "giao nhau về mặt không gian" - nếu chúng không chia sẻ bất kỳ khoảng trống nào với nhau.

Trong trường hợp không có chức năng chuyên biệt:

CREATE table t_intersect AS
SELECT 
  hp.gid, 
  hp.st_address, 
  hp.city, 
  hp.st_num,
  hp.the_geom
FROM
  public.housepoints as hp
WHERE
  hp.gid NOT IN 
  (
    SELECT 
      h.gid
    FROM 
      public.parcel as p,
      public.housepoints as h
    WHERE 
      ST_Intersects(h.the_geom,p.the_geom)
  ) AS foo

Ở đây chúng tôi sử dụng NOT EXISTSvà CREATE TABLE AS SELECT(CTAS)

CREATE table t_intersect
AS
  SELECT 
    hp.gid,
    hp.st_address,
    hp.city, hp.st_num,
    hp.the_geom
  FROM public.housepoints AS hp
  WHERE NOT EXISTS (
    SELECT 1
    FROM public.parcel AS par 
    WHERE ST_Intersects(hp.the_geom,par.the_geom)
  );

Làm thế nào về ST_Disjoint? - Trả về ĐÚNG nếu Hình học không "giao nhau về mặt không gian" - nếu chúng không chia sẻ bất kỳ khoảng trống nào với nhau.

Trong một số trường hợp rất hữu ích, hãy sử dụng LATITH THAM GIA, nó có thể rất nhanh Nó trông giống như

SELECT * FROM houses h
LEFT JOIN LATERAL (
   SELECT True t FROM parcels p
   WHERE ST_Intersects(p.geom, h.geom)
   LIMIT 1
) p ON True
WHERE p.t IS NULL;

Đơn giản chỉ cần sử dụng KHÔNG trước khi ST_Intersects thực hiện thủ thuật:

Điều này nhận được tất cả các địa chỉ không nằm trong khu phố # 62:

select 
a.*
from denver.neighborhoods as n
join denver.addresses as a on not ST_Intersects(n.geom, a.geom)
where n.nbhd_id = '62'

Lưu ý thứ tự của các cột geom - đa giác đầu tiên, điểm thứ hai, được đảo ngược so với cách sử dụng thông thường của ST_Intersects.

Nhanh chóng và đơn giản! Đã tự hỏi làm thế nào để làm điều này một cách chính xác trong một thời gian!

Đây có thể không phải là giải pháp nhanh nhất ... Nhưng tôi thường chỉ gian lận bằng cách tham gia tất cả các tính năng của bảng khác.

Create table blah as
select
  d.*
from
  data_i_want d,
  (select st_union(geom) geom from not_in_here) n
where
  st_disjoint(d.geom,n.geom);

Đẹp và linh hoạt nếu bảng not_in_here không phức tạp.