PostGIS'te kesişmeyen özellikleri seçin

Bu bana bu kadar basit bir soru gibi görünüyor (ve muhtemelen de öyle) ancak bana cevap veren bir örnek bulamıyorum. PostGIS'i kullanarak, sadece çokgenlerin dışına düşen noktaları seçmek istiyorum . Sonuçta bu, görebildiğim kadarıyla ST_Iersters'ın tersi.

Örnek: Bir taxlot katmanına ve bir adres noktası katmanına sahibim. ST_Iersters'ı kullanmam gerektiğini düşünüyorum, ama bunun tersini yapmasını nasıl söyleyebilirim? Aşağıdaki kodun önüne bir NOT ifadesi eklemeyi düşünmüştüm , ancak bu işe yaramadı.

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

“Kesişmeme” ile çalışmamasının nedeni, sadece geometrileri çiftler halinde karşılaştırmanızdır; ayrıklıkla aynı sorun olacaktır. Her ev noktası bir parselle kesişse bile bazı parselleri ayırır.

underdark'in önerisinde bu sorun yok. Ayrıca, endeksleri daha etkin kullanabilecek başka bir numara daha var:

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;

Fikir onlara st_intersects ile katılmak ve parsel kimliğinin bulunmadığı satırları almaktır.

Burada ihtiyaç duyulan endeksler uzamsal bir dizin ve parsellerdeki giddeki bir dizindir (parsel tablosundaki kimliğin de buna gideceğini varsayarsak).

ST_Disjoint arıyor olabilirsiniz

ST_Disjoint - Geometriler "uzamsal olarak kesişmiyorsa" TRUE döndürür - birlikte bir boşluk paylaşmazlarsa.

Özel bir fonksiyon yoksa:

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

Burada kullanıyoruz NOT EXISTSve 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 nasıl? - Geometriler "uzamsal olarak kesişmezse" - TRUE, birlikte bir boşluk paylaşmazlarsa TRUE döndürür.

Bazı durumlarda çok faydalıdır LATERAL JOIN, kullanımı çok hızlı olabilir.

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'ten önce basitçe NOT'u kullanmak hileyi yapmaz:

Bu, # 62 mahallesinde olmayan tüm adresleri alır:

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 sütunlarının sırasına dikkat edin - ilk önce çokgenler, ikinci, ST_Intersects kullanımından ters çevrilen puanlar.

Çabuk ve basit! Bir süredir bunun doğru şekilde nasıl yapıldığını merak ediyordum!

Bu en hızlı çözüm olmayabilir ... Ama genelde sadece diğer tablonun tüm özelliklerine katılarak hile yapıyorum.

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 tablo o kadar karmaşık değilse hoş ve çabuk .