PostGIS ile “uzun ve dar” çokgenlerin belirlenmesi

Şehir mahallelerini, geniş alanları temsil eden çokgenler setim var. Aralarındaki büyük örtüşen alanları tanımlamak istiyorum.

Ancak bir sorun var: bazen bu çokgenler çevreleri boyunca örtüşür (çünkü az bir hassasiyetle çizilirler). Bu, umursamadığım uzun ve dar çakışmalar oluşturacaktır.

Ancak diğer zamanlarda, güçlü çokgenlerle büyük örtüşmeler olacak, yani bir mahallenin çokgeninin diğeriyle örtüştüğü geniş alanlar olacak. Sadece bunları seçmek istiyorum.

Yalnızca çakışmaların aşağıdaki resmine bakın. Sadece sol alt köşedeki mavi çokgeni seçmek istediğimi düşünün.

Alanlara bakabilirdim, ama bazen dar olanları o kadar uzun olurlar ki mavi poligon kadar büyük alanlara sahip olurlar. Alan / çevre oranı yapmaya çalıştım, ama bu da karışık sonuçlar verdi.

Kullanmayı bile denedim ST_MinimumClearance, ancak bazen geniş alanların dar bir kısmı veya iki çok yakın köşesi olacak.

Başka yaklaşımlarla ilgili fikirleriniz var mı?

Sonunda benim için en iyi olan şey, aşağıdaki @Cyril ve @FGreg tarafından önerildiği gibi negatif bir tampon kullanmaktı.

Gibi bir şey kullandım:

ST_Area(ST_Buffer(geom, -10)) as neg_buffer_area

Benim durumumda, birimler metre idi, bu yüzden 10 m negatif tampon.

Dar çokgenler için bu alan sıfır döndürdü (ayrıca geometri boş olurdu). Sonra bu sütunu dar çokgenleri filtrelemek için kullandım.

Negatif bir tampon oluşturmaya çalışırdım, eğer ince çokgenler yer, o zaman iyi, çokgen yemiyorsa, o zaman benim ... :-)

daha önce doğrusal çokgenlerin genişliğinin 2 / 3'ünü ayarlayan bu komut dosyasını çalıştırın ...

create table name_table as
SELECT ST_Buffer(
(ST_Dump(
(ST_Union(
ST_Buffer(
(geom),-0.0001))))).geom,
0.0001)) as geom from source_table

İŞLETİM SİSTEMİ :-)...

Alan / çevre yerine, alanı çevrenin karesine (veya tersine) bölmek daha iyidir.

Buna "şekil dizini" de denir. Çevrenin alana bölünen karesi minimum 4 * Pi () değerine sahiptir (en kompakt 2D geometrisi olan bir disk durumunda), böylece kolay bir şekilde 4 * Pi () ile normalleştirilebilir yorumlama (1'e yakın normalleştirilmiş değerler, çok kompakt nesneleriniz ve kareleriniz yaklaşık 1,27 değerine sahip olduğunuz anlamına gelir).

DÜZENLEME: Alandaki bir eşik, kompakt olabilen çok küçük eserleri kaldırmak için yararlı olacaktır. O zaman şekil endeksi daha iyi kontrast gösterecektir. DÜZENLEME: Bu yanıta ek olarak, ST_Snap kullanımı sorunu oluşmadan çözmenize yardımcı olabilir.

Bir seçenek, çokgen alanının uçları kullanılarak çizilebilecek en uzun çizgiye oranını kullanmak olacaktır. Uzun dar çokgenlerin belirlenmesi.

select * from polygons where ST_Length(ST_LongestLine(geom, geom)) < ST_Area(geom) * 4

Bu şerit çokgenler için oldukça iyi çalışır. İhtiyaçlarınıza ve projeksiyonunuza göre oranın (alanı çoğalttığınız şey) ayarlayabilirsiniz.

Bunun kullanım durumunuzla eşleşebileceği anlaşılıyor : Seçili çokgenleri ortadan kaldırın

Giriş katmanının seçilen çokgenlerini, ortak sınırlarını silerek belirli bitişik çokgenlerle birleştirir. Bitişik çokgen, en büyük veya en küçük alana sahip olan veya çıkarılacak çokgen ile en büyük ortak sınırı paylaşan poligon olabilir.

Eliminat, normal olarak şerit çokgenlerden, yani girdilerin sınırlarının benzer olduğu ancak aynı olmadığı çokgen kesişme işlemlerinin bir sonucu olan küçük çokgenlerden kurtulmak için kullanılır.

"En Büyük Ortak Sınır" seçeneğini denemek istediğiniz anlaşılıyor.

Bu bana PostGIS topolojisi uzantısı için mükemmel bir kullanım örneği gibi görünüyor . Topolojinin tolerans parametresi, köşelerin diğer mevcut çokgenlere yapışmasına, kaynak verilerin düşük hassasiyetiyle başa çıkmasına ve temizlenmesine ne kadar izin vereceğinizi belirleyecektir.

Kısacası, strateji:

1. Topoloji uzantısını etkinleştirin

CREATE EXTENSION postgis_topology;

2. Yeni bir boş topoloji oluşturun

SELECT topology.CreateTopology('neighborhoods_topo', 4326, 1e-7);

Üçüncü parametre, CRS birimlerinde toleranstır; akıllıca seçin. İdeal olarak, birimin metre olduğu bir CRS istiyorsunuz. CRS ünitesi metre değilse, WGS 84 aka 4326'da olduğu gibi, ST_Transformçokgenlerinizi yeniden projekte etmek için kullanın .

3. Çokgenler tablosuna bir TopoGeometry sütunu ekleyin

SELECT topology.AddTopoGeometryColumn('neighborhoods_topo', 'public', 'neighborhoods', 'topogeom', 'POLYGON');

Bu yeni bir sonuç verir layer_id. Kaydet, daha sonra gerekecek. 1Başlangıçtan itibaren sıfırdan katman olacak ve her yeni çağrıda artırılacaktır .

4. Tüm çokgenleri topolojiye ekleyin

UPDATE public.neighborhoods
SET topogeom = topology.toTopoGeom(geom, 'neighborhoods_topo', 1, 1e-7);

Bu, büyük bir veri kümesi için birkaç saat sürebilir, sabırlı olun. 1daha önce döndürülen katman_kididir.

5. Çeşitli mahallelerde görünen yüzleri bulun

Topolojiden 2 veya daha fazla topogeometride bulunan tüm yüzleri bulun. Sorguyu bir alıştırma olarak bırakacağım. En kolayı muhtemelen GetTopoGeomElementsişleve göre, sonra yüz kimliğine göre gruplandırın ve 2 veya daha fazla sayıya sahip olanlara bakın. Alternatif olarak, topogeom sütunundan temizlenmiş geometri ile yeni bir tablo oluşturabilir, sadece standart geometriye atabilir topogeom::geometryve şu anda sahip olduğunuz şeyi tekrarlayabilirsiniz, ancak şimdi şerit üst üste gelmeden temiz bir veri kümesiyle.