Bir şekil dosyasındaki çokgen diğerinin çokgen içindeki yüzdesi

Ben bir acemi, bu açık / özür dilerim zaten sorulmuş ve cevap ama hiçbir şey bulamadım.

İki şekil dosyaları var: 1. İngiltere'de 500 küçük bölgeleri olan bir LSOA sınırı olarak bilinen bir ilçe için bir idari sınır katmanı 2. bir sel bölgesi.

İdeal olarak, küçük LSOA bölgelerinden hangisinin sel bölgesi içinde ≥% 50 olduğunu bulmak ve 500 LSOA bölgesinin her biri için evet / hayır veya 1/0 ile sonuçlanmak istiyorum.

Ama bunu nasıl yapacağımı bilmiyorum. İki şekil dosyasına katılabileceğimi düşündüm, ancak aralarında ortak bir özellik yok. Daha sonra, sel bölgesinde hangi LSOA'nın bulunduğunu gösteren ve bana göre gösterilen Özelliğe Göre Katıl Özelliği işlevini kullanabileceğimi düşündüm, ancak neredeyse hepsi bu (bkz. Resim 2).

Bu bir SQL sorunu olduğunu düşünüyorum ama bilmiyorum. QGIS'te yeniyim ve hiç PostgreSQL kullanmadım.

Herhangi bir yardım büyük mutluluk duyacağız. Güzel insanların bana yardım etmek için ihtiyaç duydukları bilgileri sağlayabilirim.

Bu, QGIS'de bulunan coğrafi işlem araçlarını kullanan nispeten basit bir görevdir.

1. LSOA bölgelerinizin alanını hesaplayın.

   • LSOA katmanı öznitelik tablosunu açın.
   • Düzenleme modunu etkinleştirin.
   • Alan hesap makinesini açın.
   • "$ Area" ifadesi ile "Ondalık sayı (gerçek)" türünde yeni bir alan oluşturun.
   • Düzenleme modunu devre dışı bırakın (düzenlemeleri kaydedin).

2. Taşkın bölgesi katmanını tek bir çok parçalı özelliğe birleştirin.

   • Vector > Geometry Tools > Singleparts to Multipart.
   • Benzersiz Kimlik alanı için "--- Tümünü birleştir ---" seçeneğini belirleyin.

3. LSOA bölge katmanını çok parçalı sel bölgesi katmanıyla kesiştirin.

   • Vector > Geoprocessing Tools > Intersect.
   • Giriş katmanı LSOA bölgeleridir, kesişen katman taşkın bölgeleridir.

4. Ortaya çıkan katman, sel alanları katmanıyla örtüşen LSOA bölgelerinin (LSOA bölgeleri katmanındaki özniteliklerle) parçaları olacaktır. Taşkın bölgesi içindeki her LSOA bölgesinin oranını hesaplamak için:

   • Kesişen özelliklerin alanını hesaplayın (1. adımdaki gibi), ardından
   • Orijinal (toplam) alanı kesişen alana bölerek başka bir alan ekleyin. Sonuç, 0 ile 1 arasında bir ondalık sayıdır. Yüzde vermek için 100 ile çarpın.

5. Her iki katmanın paylaştığı benzersiz kimliği kullanarak orijinal LSOA katmanını kesişen katmana birleştirin.

6. Birleştirilen katmanı yeni bir şekil dosyası olarak dışa aktarın.

7. Çoğaltılan öznitelikleri silin.

Et voilà!

Adım # 2 olmadan, her LSOA özelliği için her farklı taşkın bölgesi özelliği için ayrı bir özellik oluşturulacaktır. Her LSOA bölgesi için yalnızca toplam kapsama alanı ile ilgileniyorsanız, bu muhtemelen istediğiniz bir şey değildir. Fluvial / tidal / pluvial sel arasında ayrım yapmak istiyorsanız (ve sel bölgesi verileri bunu destekliyorsa), singlepart s'yi "IDPE" alanını Unique ID alanı olarak belirten multipart'a dönüştürebilirsiniz.

Uzamsal ve bazı uzamsal SQL işlevlerini kullanabilirsiniz.

Select t1.geometry, t1.ID, area(t1.geometry), area(t2.geometry) ...... (anything you need to have in the table results)

(area(intersection(t1.geometry,t2.geometry))) as "Commun_AREA"

, ("Commun_AREA"*100/(area(t1.geometry))) as "Percent_AREA"

From lsoa as t1, flood_zone as t2

Where Intersects( t1.geometry,t2.geometry ) = 1

Bu, verilen cevaplardan çok daha kolay yapılabilecek bir şey gibi görünüyor. Şahsen basit bir python betiği kullanırdım:

floodName = "the layer name here"
boundryName = "the layer name here"
fieldName = "the name of the field to contain the output 1/0"
minCoverage = 0.5 # the minimum amount of area covered to write 1
updateMap = [] # this will store values to be written    

# get layers
floodLayer = QgsMapLayerRegistry.instance().mapLayersByName(floodName)[0]
boundryLayer = QgsMapLayerRegistry.instance().mapLayersByName(boundryName)[0]
fieldIndex = boundryLayer.dataProvider().fieldNameIndex(fieldName)    

# iterate through boundries
for b in boundryLayer.getFeatures():
    # get only flood features that intersect with this feature's bounding box
    # this will make the script go way faster than it would otherwise
    request = QgsFeatureRequest().setFilterRect(b.geometry().boundingBox())
    floodGeom = geometry()
    floodFeat = QgsFeature()
    iter = floodLayer.getFeatures(request)
    iter.nextFeature(feat)
    while iter.nextFeature(feat):
        floodGeom = floodGeom.combine(feat.geometry())
    intersectGeom = b.geometry().intersection(feat.geometry())
    if intersectGeom.area() > minCoverage * b.geometry().area():
        updateMap[b.id()] = {fieldIndex : 1}
    else:
        updateMap[b.id()] = {fieldIndex : 0}

boundryLayer.dataProvider().changeAttributeValues(updateMap)

bu, her sınır katmanının sınırlayıcı kutusuyla kesişen taşkın poligonlarını değerlendirir, böylece çalıştırılması oldukça hızlı olmalıdır, daha sonra mevcut katmandaki yalnızca bir alanı günceller (tamamen yeni bir katman oluşturma ve eski değerleri kopyalama karmaşık bir işlem yerine) sonra silme)

Adım 3'te "Kesişim" yöntemini kullanarak Snorfalorpagus'un talimatlarını takip eden KJ ile aynı problemi yaşadım Hesaplamak biraz zaman aldı ve neyle kaldığım boştu.

QGIS içinde Intersect yerine "Clip" yöntemini kullanmak dışında aynı adımları izlemeye çalıştım - yani, örneğinizde kalan alanın sel bölgesi tarafından kaplanmayan kısımları olurdu. Bu bir nedenle işe yaramış gibi görünüyordu ve bir önceki adımdan "Alan" alan hesaplamasını ve her bir poligonun geri kalan kısımlarında yeni bir "Alan" hesaplamasını, her alanın% OLMADIĞINI anlamak için kullanabildim diğer Çokgen katmanı ile kaplanmıştır.

Bu teknik olarak istediğin tersi. Ama oradan ne elde etmek 1'den her değerin çıkarılmasından sadece bir meselesi olduğunu taşkın sınırının içindedir.