Đếm các đa giác chồng chéo từ một Shapefile bằng ArcGIS Desktop?

Tôi có một shapefile với 16.400 đa giác. Mỗi đa giác cho thấy sự mở rộng của một loài chim cho toàn thế giới.

Bây giờ tôi phải đếm các đa giác chồng chéo. Tôi đã thử nó với công đoàn và giải thể (đếm công đoàn), nhưng công đoàn không hoạt động cho rất nhiều đa giác.

Sau đó, tôi đã cố gắng cắt các lục địa, nhưng điều này cũng không hoạt động vì số lượng đa giác rất lớn. Hơn nữa tôi đã thử phương pháp này , cũng không thành công.

Vì vậy, tôi đang hỏi bạn rằng có cách nào để đếm các đa giác chồng chéo nếu 16400 đa giác nằm trong một shapefile không?

Tôi đang làm việc với 10.0 và có thể làm việc với 10.2. Một giải pháp ArcPy cũng tuyệt vời.

Hiện tại tôi đang suy nghĩ về việc tạo một fishnet và lặp qua các hàng của shp với 16400 đa giác và viết 1 vào trường giá trị của một ô fishnet nếu đa giác nằm trong ô này và hơn là lấy hàng tiếp theo (đa giác) và nếu đây cũng là trong tế bào fishnet đếm +1.

Nhưng tôi không biết liệu đây có phải là một giải pháp tốt hay không và làm thế nào để nhận ra nó. Hoặc tôi phải học R để sử dụng phương pháp này .

Kết quả: Nó phải là một hình dạng trong đó bạn có các đa giác mới trong số các đa giác chồng chéo và một trường trong đó các phần trùng lặp được tính.

Vì vậy, cuối cùng nên có một shapefile nơi bạn có thể thấy có bao nhiêu loài chim được tìm thấy ở cùng một nơi.

Tôi sẽ khuyên bạn nên sử dụng công cụ Đa giác chồng chéo tùy chỉnh .

Sự miêu tả:

Mẫu này chứa một hộp công cụ với một công cụ, Đếm các đa giác chồng chéo. Đưa ra một lớp đối tượng hoặc lớp chứa đa giác chồng lấp, xuất ra một lớp đối tượng mới với các lớp phủ bị loại bỏ và trường Join_Count chứa số lượng đa giác chồng lấp.

Sử dụng mã thông báo hình học arcpy, bạn có thể thử một cái gì đó như thế này:

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

Đối với dữ liệu mẫu ở trên, mã sẽ trả về số lượng trùng lặp sau:

Mã như sẽ chỉ trả về số lượng cho đa giác có ít nhất một chồng lấp.

Một phương pháp rất đơn giản là:

1. Liên minh các shapefile với chính nó;
2. Chuyển đổi đầu ra nhiều phần thành một phần ;
3. Sử dụng công cụ nối không gian để đếm các lớp phủ (sử dụng tùy chọn đối sánh IS_IDENTICAL_TO);
4. Tượng trưng bằng cách sử dụng các join_countlĩnh vực.

Tôi đoán bạn đã thử phương pháp này: Đếm và quét các đa giác chồng chéo trong ArcGIS Desktop?

16.400 đa giác không phải là nhiều. Tuy nhiên, một giải pháp tiềm năng là chỉ cần thực hiện Tham gia không gian thông thường. Trong ArcMap toolbox, > Analysis Tools -> Overlap -> Spatial Join.

Đặt cả hai targetvà joincác tính năng cho cùng một tập dữ liệu và chỉ định một đầu ra. Để lại phần còn lại của cài đặt.

Sau một vài phút, bạn sẽ lấy lại một shapefile có chứa cột "đếm số tham gia". Trừ 1 từ điều này (vì rõ ràng mỗi tính năng nên giao nhau) và đó phải là số lượng "chồng chéo" (thực sự giao nhau) cho mỗi đa giác.

Tôi chỉ thực hiện nó trên

Tôi đã tải xuống và dùng thử công cụ "Đếm đa giác chồng chéo". Nó có thể hoạt động, nhưng phải mất một thời gian rất dài (có thể vì kích thước tệp, nhưng FC đầu vào của tôi chỉ có <5.000 bản ghi).

Trong khi chờ công cụ đó chạy, tôi đã mở một cửa sổ ArcMap khác và chỉ mất vài bước để có được thứ tôi muốn. 1) Tham gia không gian - sử dụng cùng một lớp tính năng như Tính năng mục tiêu và Tham gia và chọn tùy chọn "Tham gia một đến nhiều". 2) Hòa tan - sử dụng đầu ra từ bước cuối cùng. Sử dụng "TARGET_FID" làm trường hòa tan và để thống kê, bạn có thể SUM trường "Join_Count" hoặc COUNT trường "THAM GIA". 3) Trong tệp đầu ra từ bước 2, sử dụng máy tính trường trừ 1 từ trường thống kê ("SUM_Join_Count" hoặc "COUNT_JOIN_FID") - vì mỗi tính năng tự giao nhau.

Tôi đề nghị sử dụng phương pháp này qua công cụ "Đếm đa giác chồng chéo". Tôi đã bắt đầu chạy công cụ COP ~ 5 phút trước khi bắt đầu phương thức Tham gia-> Hòa tan này và nó cho tôi kết quả có đủ thời gian để viết nó lên trước khi công cụ "Đếm chồng chéo đa giác" thậm chí đã kết thúc.

Hi vọng điêu nay co ich!