ArcGIS Desktop kullanarak bir Shapefile'dan çakışan çokgenleri sayıyor musunuz?

16.400 çokgen içeren bir şekil dosyam var. Her çokgen, bir kuş türünün tüm dünya için genişlemesini gösterir.

Şimdi örtüşen çokgenleri saymalıyım. Sendika ile denedim ve dağıttım (birliği sayın), ancak sendika çok fazla çokgen için çalışmıyor.

Sonra kıtaları kırpmaya çalıştım, ancak çok sayıda çokgen nedeniyle bu da çalışmıyor. Üstelik bu yöntemi de denedim , başarılı olamadı.

Bu nedenle, size 16400 çokgen bir şekil dosyasındaysa, çakışan çokgenleri saymanın bir yolu olduğunu soruyorum.

10.0 ile çalışıyorum ve 10.2 ile çalışabiliyorum. Bir ArcPy çözümü de harika.

Şu anda bir ağ oluşturmayı düşünüyorum ve 16400 çokgenleri ile shp satırları üzerinde yineleme yapıyorum ve çokgen bu hücrede ise ve bir sonraki satırı (çokgen) almak için bir fishnet hücresinin değer alanına 1 yazıyorum ve bu da fishnet hücre sayısında +1 varsa.

Ama bunun iyi bir çözüm olup olmadığını ve nasıl gerçekleştirileceğini bilmiyorum. Ya da bu yaklaşımı kullanabilmek için R öğrenmem gerekiyor .

Sonuç: Üst üste binenlerden yeni çokgenlere sahip olduğunuz bir şekil ve örtüşmelerin sayıldığı bir alan olmalıdır.

Sonuçta, aynı yerde kaç kuş türünün bulunduğunu görebileceğiniz bir şekil dosyası olmalıdır.

Sayılan Çakışan Özellikler (Analiz) aracını kullanmanızı öneririm .

Giriş özelliklerinden düzlemsel örtüşen özellikler üretir. Çakışan özelliklerin sayısı çıktı özelliklerine yazılır.

Arcpy geometri belirteçlerini kullanarak, şöyle bir şey deneyebilirsiniz:

import os
import arcpy

arcpy.env.workspace = r"" #path to workspace
arcpy.env.overwriteOutput = 1

polygon_fc = r"" #path to polygon fc

base = [row for row in arcpy.da.SearchCursor(polygon_fc,["OID@","SHAPE@"])]
compare = base

overlaps_stats = {}

for b in base:
    for c in compare:
        if b[1].overlaps(c[1]):
            #print "{0} overlaps {1}".format(b[0],c[0])
            if overlap_stats.has_key(b[0]):
                overlap_stats[b[0]].append(c[0])
            else:
                overlap_stats[b[0]] = [c[0]]

for key,value in overlap_stats.iteritems():
    print "Polygon {0}:  Overlaps: {1}.".format(key,len(value))

Yukarıdaki örnek veriler için, kod aşağıdaki çakışma sayılarını döndürür:

Olduğu gibi kod, yalnızca en az bir çakışması olan çokgenler için sayıları döndürür.

Çok basit bir yöntem:

1. Şekil dosyasını kendisiyle birleştirin;
2. Çok parçalı çıktıyı tek parçaya dönüştürün ;
3. Çakışmaları saymak için uzamsal birleştirme aracını kullanın (ARE_IDENTICAL_TO eşleme seçeneğini kullanın);
4. join_countAlanı kullanarak sembolize edin .

Sanırım şu yöntemi denediniz: ArcGIS Desktop'ta çokgen çakışmalarını sayma ve rasterleştirme?

16.400 çokgen o kadar da değil. Ancak, potansiyel çözümlerden biri, düzenli bir Mekansal Katılma yapmaktır. İçinde ArcMap toolbox, > Analysis Tools -> Overlap -> Spatial Join.

Hem targetve joinözelliklerini aynı veri kümesine ayarlayın ve bir çıktı belirtin. Ayarların geri kalanını bırakın.

Birkaç dakika sonra "birleştirme sayısı" sütunu içeren bir şekil dosyasını geri almanız gerekir. Bundan 1 çıkarın (açıkçası her özelliğin kendisi kesişmesi gerekir) ve bu her çokgen için "örtüşmeler" (aslında kesişir) sayısı olmalıdır.

Az önce yaptım

"Çakışan Çokgenleri Say" aracını indirip denedim. Çalışabilir, ancak çok uzun zaman alır (muhtemelen dosya boyutu nedeniyle, ancak FC girişimin yalnızca <5.000 kaydı vardı).

Bu aracın çalışmasını beklerken başka bir ArcMap penceresi açtım ve istediğim şeyi elde etmek için sadece birkaç hızlı adım attı. 1) Uzamsal Birleştirme - Hedef ve Birleştirme Özellikleri ile aynı özellik sınıfını kullanarak ve "Bire Çok Katıl" seçeneğini seçin. 2) Çözün - son adımdaki çıktıyı kullanarak. Çözme alanı olarak "TARGET_FID" kullanın ve istatistikler için "Join_Count" alanını TOPLA veya "JOIN_FID" alanını COUNT yapabilirsiniz. 3) 2. adımdaki çıktı dosyasında, her özellik birbiriyle kesiştiği için alan hesap makinesini istatistik alanından ("SUM_Join_Count" veya "COUNT_JOIN_FID") 1 çıkarın.

Bu yöntemi "Çakışan Çokgen Say" aracı üzerinden kullanmanızı öneririm. Bu Join-> Dissolve yöntemine başlamadan önce COP aracını ~ 5 dakika çalıştırmaya başladım ve bana "Üst üste binen çokgen sayımı" aracı bile bitmeden bunu yazmak için yeterli zaman verdi.

Bu yardımcı olur umarım!