Sélectionnez les fonctionnalités qui ne se croisent pas dans PostGIS

Cela me semble être une question si simple (et c'est probablement le cas), mais je n'arrive pas à trouver un exemple qui me donne la réponse. À l’aide de PostGIS, je souhaite simplement sélectionner des points situés hors des polygones. En fin de compte, c'est l'inverse de ST_Intersects, autant que je puisse le voir.

Exemple: j'ai une couche taxlot et une couche de points d'adresses. Je suppose que je devrais utiliser ST_Intersects, mais comment puis-je lui demander de faire la sélection inverse? Je pensais peut-être ajouter une instruction NOT devant le code ci-dessous, mais cela n'a pas fonctionné.

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

La raison pour laquelle cela ne fonctionne pas avec "pas d'intersection" est que vous ne comparez les géométries que par paires; il y aura le même problème avec disjoint. Chaque point de séparation divisera certaines parcelles même s'il en coupe une.

La suggestion de Underdark n'a pas ce problème. Il existe également une autre astuce qui permettra probablement une utilisation plus efficace des index:

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;

L'idée est de les joindre avec st_intersects et d'obtenir les lignes où l'ID de parcelle n'est pas présent.

Les index nécessaires ici sont un index spatial et un index sur gid dans les parcelles (en supposant que l'id de la table parcels est également appelé gid).

Vous cherchez peut-être ST_Disjoint

ST_Disjoint - Retourne VRAI si les géométries ne "se croisent pas spatialement" - si elles ne partagent aucun espace ensemble.

S'il n'y a pas de fonction spécialisée:

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

Ici nous utilisons NOT EXISTSet 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)
  );

Que diriez-vous de ST_Disjoint? - Retourne VRAI si les géométries ne "se croisent pas spatialement" - si elles ne partagent aucun espace.

Dans certains cas, il est très utile d’utiliser LATERAL JOIN, cela peut être très rapide.

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;

Il suffit de ne pas utiliser NOT avant ST_Intersects:

Cela donne toutes les adresses qui ne sont pas dans le voisinage n ° 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'

Notez l'ordre des colonnes geom - les polygones en premier, les points en second, qui sont inversés par rapport à l'usage habituel de ST_Intersects.

Rapide et simple! Je me demandais comment faire cela correctement depuis un moment!

Ce n'est peut-être pas la solution la plus rapide ... Mais d'habitude, je ne fais que tricher en joignant toutes les fonctionnalités de l'autre table.

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

Sympa et accrocheur si la table not_in_here n'est pas si complexe.