Выберите объекты, которые НЕ пересекаются в 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);

Причина, по которой он не работает с «не пересекаются», состоит в том, что вы сравниваете геометрию только парами; будет такая же проблема с дизъюнктом. Каждая точка дома будет разъединять некоторые посылки, даже если она пересекает одну посылку.

Предложение Подземья не имеет этой проблемы. Существует также еще одна хитрость, которая, вероятно, позволит более эффективно использовать индексы:

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 - возвращает TRUE, если геометрия не «пространственно пересекается» - если они не разделяют пространство вместе.

В случае, если нет специализированной функции:

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? - Возвращает TRUE, если геометрия не "пространственно пересекается" - если они не разделяют пространство вместе.

В некоторых случаях очень полезно использовать 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;

Простое использование NOT до того, как 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 не так сложна.