Виберіть функції, які НЕ перетинаються в PostGIS

Мені це здається таким простим питанням (і це, мабуть, є), але я не можу знайти приклад, який дає мені відповідь. Використовуючи PostGIS, я просто хочу виділити точки, які випадають поза полігонів. Зрештою, це обернення ST_Intersects, наскільки я це бачу.

Приклад: у мене є шар таксолота та рівень адресної точки. Я припускаю, що я повинен використовувати ST_Intersects, але як мені це зробити, щоб зробити зворотний вибір? Я подумав, що, можливо, додамо оператор NOT перед кодом нижче, але це не спрацювало.

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);

Причина, з якою вона не працює з "не перетинається", полягає в тому, що ви порівнюєте лише геометрії в парах; буде однакова проблема і з нерозбірливістю. Кожна домашня точка буде роз'єднана деякими посилками, навіть якщо вона перетинає одну посилку.

Пропозиція underdark не має такої проблеми. Також є ще одна хитрість, яка, ймовірно, дозволить ефективніше використовувати індекси:

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;

Ідея полягає в тому, щоб з'єднати їх із st_intersects та отримати рядки там, де ідентифікатор посилки відсутній.

Тут необхідні індекси - це просторовий індекс та індекс на gid у посилках (якщо припустити, що id у таблиці посилок також називається gid).

Можливо, ви шукаєте ST_Disjoint

ST_Disjoint - повертає ІСТИНА, якщо геометрії не "просторово перетинаються" - якщо вони не ділять простору разом.

У разі відсутності спеціалізованої функції:

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

Тут ми використовуємо NOT EXISTSі 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)
  );

Як щодо ST_Disjoint? - Повертає ІСТИНА, якщо Геометрії не "просторово перетинаються" - якщо вони не ділять простору разом.

У деяких випадках дуже корисно використовувати LATERAL JOIN, це може бути дуже швидко. Це повинно виглядати так

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;

Просто використовуючи НЕ, перш ніж ST_Intersects виконає трюк:

Це отримує всі адреси, які не знаходяться в околицях №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'

Зверніть увагу на порядок стовпчиків Geom - по-перше, багатокутників, пункти друге, яке відмінено від звичайного використання ST_Intersects.

Швидко і просто! Замислювались, як це правильно зробити на деякий час!

Це може бути не найшвидшим рішенням ... Але я зазвичай просто обманюю, приєднуючись до всіх особливостей іншої таблиці.

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);

Приємно і спритно, якщо таблиця not_in_here не така складна.