Підрахунок перекриваючих полігонів з одного Shapefile за допомогою ArcGIS Desktop?

У мене формуляр з 16 400 полігонами. Кожен багатокутник показує розширення виду птахів на весь світ.

Тепер мені належить порахувати перекриваються багатокутники. Я спробував це з об'єднанням і розпустити (порахуйте союз), але союз не працює для такої кількості багатокутників.

Тоді я спробував обрізати континенти, але це теж не працює через величезну кількість полігонів. Більше того, я спробував цей метод , також без успіху.

Тому я запитую вас, чи існує спосіб підрахунку полігонів, що перекриваються, якщо 16400 багатокутників знаходяться в одному форматі?

Я працюю з 10.0 і можу працювати з 10.2. Рішення ArcPy також чудове.

На даний момент я замислююся над тим, щоб створити риб’ячу мережу і проітеравати над рядами SHP з полігонами 16400 і записати 1 у поле значення осередку ажурної мережі, якщо полігон знаходиться в цій комірці, а потім візьмемо наступний рядок (багатокутник) і якщо це також є в осередку ажурної мережі +1.

Але я не знаю, чи це хороше рішення, і як його реалізувати. Або я повинен навчитися R використовувати цей підхід.

Результат: Це має бути форма, де у вас є нові багатокутники поза перекриваються, і поле, де рахуються перекриття.

Отже, зрештою має бути файл форми, де можна побачити, скільки видів птахів зустрічається на одному місці.

Я рекомендую скористатися інструментом " Перерахування функцій (Аналіз)" .

Створює заплановані функції перекриття із вхідних функцій. Кількість функцій, що перекриваються, записується на вихідні функції.

Використовуючи маркери геометрії arcpy, ви можете спробувати щось подібне:

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

Для наведених вище зразків даних код поверне наступні підрахунки перекриття:

Код, як є, повертає рахунки лише для полігонів, що мають принаймні одне перекриття.

Дуже простий метод:

1. Союз шейпфайл з самим собою;
2. Перетворення багаточастинного виводу в одну частину ;
3. Використовуйте інструмент просторового з'єднання для підрахунку перекриттів (використовуйте параметр відповідності ARE_IDENTICAL_TO);
4. Символізуйте за допомогою join_countполя.

Напевно, ви спробували цей метод: підрахунок та растрування перекриттів полігонів у ArcGIS Desktop?

16 400 багатокутників не так вже й багато. Однак одне потенційне рішення - просто зробити звичайний просторовий приєднання. В ArcMap toolbox, > Analysis Tools -> Overlap -> Spatial Join.

Встановіть targetі joinтой самий набір даних та функції та вкажіть вихід. Залиште інші налаштування.

Через кілька моментів вам слід повернути файл форми, який містить стовпчик "кількість приєднань". Відніміть 1 з цього (як очевидно, кожна ознака повинна перетинатися), і це має бути кількість "перекриттів" (фактично перетинається) для кожного багатокутника.

Я щойно виконував це

Я завантажив і спробував інструмент "Полігони, що перекриваються". Це може спрацювати, але це займає жахливо тривалий час (можливо, тому що розмір файлу, але мій вхідний FC лише мав <5000 записів).

Поки я чекав запуску цього інструменту, я відкрив ще одне вікно ArcMap, і для того, щоб отримати те, що потрібно, було потрібно лише кілька швидких кроків. 1) Просторове приєднання - використовуючи той самий клас функцій, що і цільові та приєднатися, та вибравши опцію "Приєднатись до багатьох". 2) Dissolve - використовуючи вихід з останнього кроку. Використовуйте "TARGET_FID" як поле для розчинення, а для статистики ви можете ПІДПИСАТИ поле "Join_Count" або ЗРОБИТИ поле "JOIN_FID". 3) У вихідному файлі з кроку 2 за допомогою калькулятора поля віднімайте 1 з поля статистики ("SUM_Join_Count" або "COUNT_JOIN_FID") - оскільки кожна функція перетинається сама.

Я пропоную використати цей метод через інструмент "Полігон, що перекривається". Я почав запускати інструмент COP за 5 хвилин до того, як запустити цей метод Join-> Dissolve, і це дало мені результат з достатньою кількістю часу, щоб написати це до того, як інструмент "Count Overlopping Polygon" ще закінчився.

Сподіваюся, це допомагає!