ST_Intersection nasıl kullanılır?

İşte yapmaya çalıştığım şeylerle ilgili kısa bir özet: Postgres'te 3 tablo var, 'a' ve 'b', her birinin bir Çokgen sütunu var ve 'c' bir Nokta sütununa sahip. Burada yapmaya çalıştığım şey 'a', 'b' ve 'c' arasındaki geometri kesişimlerini elde etmek ve bu geometrileri bir OpenLayers vektör katmanında görüntülemek.

Zaten OpenLayers bir String herhangi bir geometri görüntülemek nasıl biliyorum, ama PostGIS 'ST_Intersection işlevi ile sorun yaşıyorum, bunu yapıyorum:

SELECT ST_Intersection(a.geom, b.geom) as inter from a, b;

Burada a.geom ve b.geom her ikisi de geometri sütunlarıdır ve bu hata iletisini alıyorum:

NOTICE:  TopologyException: found non-noded intersection between 515172 2.14408e+06, 497067 2.13373e+06 and 501321 2.13546e+06, 471202 2.14843e+06 500621 2.13576e+06 
ERROR:  GEOS Intersection() threw an error!

Ayrıca böyle geometri ST_AsText kullanarak metin olarak ortaya çıkan geometri ifade etmeye çalıştım:

SELECT ST_AsText(ST_Intersection(a.geom, b.geom)) as inter from a, b;

ama bana şu hata mesajını gönderdi:

HINT: No function matches the given name and argument types. You might need to add explicit type casts.

Ne yaptığımı bilmiyorum, çokgenlerin WKT'sini OpenLayers'da görüntülemesini istiyorum, bir WKT'den bir geometriyi nasıl görüntüleyeceğim:

                    var in_options = {
                        'internalProjection': new OpenLayers.Projection("EPSG:4326"),
                        'externalProjection': new OpenLayers.Projection("EPSG:4326")
                    }; 

                    var fea= new OpenLayers.Format.WKT(in_options).read(data); //data is the string with the WKT
                    vectorLayer.addFeatures([fea]); //this piece of code works great
                    map.zoomToExtent(bounds);

GÜNCELLEME: Sonrakini denedim:

SELECT ST_Intersection(a.geom, b.geom) as intersect_ab FROM a INNER JOIN b ON 
ST_Intersection(a,b) WHERE ST_Overlaps(a.geom, b.geom) 
AND ST_isvalid(a.geom)='t' AND ST_isvalid(b.geom)='t';

ama bir sonraki hata mesajını alıyorum:

ERROR: Function st_intersection(a,b) does not exist.
HINT: No function matches the given name and argument types. You might need to add explicit type casts.

Yalnızca geçerli çokgenlerin değerlendirildiğini doğrulamak için geçersiz ekledim, ancak hatanın ST_Intersection'da (a, b) olduğunu söylüyor, hem a, b hem de c aynı SRID'ye sahip olduğunu söylüyor, bu yüzden gerçekten kafam karıştı, üzgünüm çok soruyorum, ama PostGIS ile oldukça yeniyim, umarım seni fazla rahatsız etmiyorum. Teşekkürler.

Benim tahminim kavşak NULL döndürürse başarısız olurdu. Bu nedenle, WKT'yi oluşturmaya çalışmadan önce gerçekten bir kavşak olup olmadığını kontrol etmek için where cümlesi eklemelisiniz.

İpucu

ERROR: Function st_intersection(a,b) does not exist.
HINT: No function matches the given name and argument types. You might need to add explicit type casts.

Hata mesajının dediği gibi st_intersection'ı bu şekilde kullanamazsınız. Diğer cevapları özetlerken şöyle bir şey kullanmalısınız:

SELECT ST_Intersection(a.geom, b.geom) as intersect_ab 
FROM a INNER JOIN b ON ST_Intersects(a.geom,b.geom)
WHERE ST_isvalid(a.geom)='t' AND ST_isvalid(b.geom)='t';

AFAIK, st_overlaps ve st_intersects öğelerini oldukça benzer oldukları cümlede kullanmanın bir anlamı yoktur .

Farklı çokgen katmanları arasında test yapıyorum ve katmanın birinde en azından geçersiz bir geometri varsa başarısız oldu. Çokgenlerinizin geçerliliğini ST_isvalid (the_geom) kullanarak kontrol ettiniz mi? Anahtar olabilir.

Bunun gibi bir şey deneyin:

SELECT  ST_Intersection(a.geom, b.geom) As intersect_ab
    FROM a INNER JOIN b ON ST_Intersection(a,b)
    WHERE ST_Overlaps(a.geom, b.geom)
    ;

Kaynak

Geçersiz geometrileri hariç tutmaya çalıştım ama işe yaramadı, bu yüzden sonunda her geçersiz geometriyi silmek ve bunu kullanmak zorunda kaldım:

SELECT ST_AsText(ST_Intersection(a.geom, b.geom)) as intersect_ab FROM a,b 
WHERE ST_Overlaps(a.geom, b.geom) AND ST_isvalid(a.geom)='t' AND ST_isvalid(b.geom)='t';

Gördüğünüz gibi ST_Intersection (a, b) kısmını atladım ve bu harika çalıştı, üzgünüm çünkü herhangi bir geçersiz geometri seçimimden hariç tutmanın bir yolunu bulamadım. Buradan dışarı.

Bu sorunu bir kez yaşadım.

<pre>NOTICE:  TopologyException: found non-noded intersection between xxx, xxxx and xxx, xxx  ERROR:  GEOS Intersection() threw an error!</pre>

Bu hatayı bu yöntemi kullanarak çözmeyi başardım.
- QGIS kullanın
- Veritabanınızdan vektör katmanı ekleyin
- Hata mesajından noktayı alın ve QGIS içinde arayın
   "QuickWKT" (eklenti) bulmak için kullanılabilir
- O zaman sorun linestring göreceksiniz
- Düzenleme modunu açın
- Yeşil düğümü göstermek için "düğüm aracı" nı seçin (düğüm sorunu)
- Düğümü çakışan düğümden uzağa taşı
- Değişiklikleri kaydet