การนับรูปหลายเหลี่ยมที่ทับซ้อนกันจากหนึ่ง Shapefile โดยใช้ ArcGIS Desktop?

ฉันมีไฟล์รูปร่างที่มี 16,400 รูปหลายเหลี่ยม รูปหลายเหลี่ยมแต่ละอันแสดงให้เห็นถึงการขยายพันธุ์ของนกทั้งโลก

ตอนนี้ฉันต้องนับรูปหลายเหลี่ยมที่ทับซ้อนกัน ฉันลองกับสหภาพและยุบ (นับยูเนี่ยน) แต่สหภาพไม่ทำงานสำหรับรูปหลายเหลี่ยมจำนวนมาก

จากนั้นฉันพยายามที่จะคลิปทวีป แต่ก็ไม่ได้ผลเพราะรูปหลายเหลี่ยมจำนวนมาก ยิ่งกว่านั้นฉันลองวิธีนี้แต่ก็ไม่ประสบความสำเร็จ

ดังนั้นฉันขอให้คุณฉันมีวิธีการนับรูปหลายเหลี่ยมที่ทับซ้อนกันถ้ารูปหลายเหลี่ยม 16400 มีอยู่ในรูปร่างเดียวหรือไม่

ฉันทำงานกับ 10.0 และสามารถทำงานกับ 10.2 โซลูชันของ ArcPy นั้นยอดเยี่ยมเช่นกัน

ในขณะนี้ฉันกำลังคิดเกี่ยวกับการสร้างแหอวนและวนซ้ำแถวของ shp ด้วย 16400 polygons และเขียน 1 ไปยังเขตข้อมูลค่าของเซลล์ fishnet หากรูปหลายเหลี่ยมอยู่ในเซลล์นี้และใช้แถวถัดไป (รูปหลายเหลี่ยม) และ ถ้าสิ่งนี้อยู่ในเซลล์ fishnet ด้วย +1

แต่ฉันไม่รู้ว่านี่เป็นวิธีแก้ปัญหาที่ดีหรือไม่และจะรู้ได้อย่างไร หรือฉันต้องเรียนรู้ R เพื่อใช้วิธีนี้

ผลลัพธ์: ควรเป็นรูปร่างที่คุณมีรูปหลายเหลี่ยมใหม่จากส่วนที่ทับซ้อนกันและเขตข้อมูลที่มีการนับการทับซ้อน

ดังนั้นในที่สุดควรมีไฟล์รูปร่างที่คุณสามารถดูจำนวนนกที่พบในสถานที่เดียวกัน

ฉันจะแนะนำให้ใช้จำนวนที่ทับซ้อนกันคุณสมบัติ (วิเคราะห์)เครื่องมือ

สร้างคุณลักษณะการทับซ้อนของ planarized จากคุณสมบัติอินพุต จำนวนของคุณลักษณะที่ทับซ้อนกันถูกเขียนไปยังคุณลักษณะของเอาต์พุต

การใช้โทเค็นเรขาคณิตของ 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. ยูเนี่ยน shapefile กับตัวมันเอง;
2. แปลงการส่งออกหลายส่วนที่จะเป็นส่วนหนึ่งเดียว ;
3. ใช้เครื่องมือเข้าร่วมเชิงพื้นที่เพื่อนับการทับซ้อน (ใช้ตัวเลือกการจับคู่ ARE_IDENTICAL_TO)
4. สัญญลักษณ์ใช้join_countฟิลด์

ฉันเดาว่าคุณได้ลองวิธีนี้แล้ว: การนับและการซ้อนทับหลายเหลี่ยมแรสเตอร์ใน ArcGIS Desktop หรือไม่

รูปหลายเหลี่ยม 16,400 ไม่มากนัก อย่างไรก็ตามวิธีแก้ปัญหาที่อาจเกิดขึ้นอย่างหนึ่งคือเพียงเข้าร่วม Spatial แบบปกติ ArcMap toolbox, > Analysis Tools -> Overlap -> Spatial Joinใน

ตั้งค่าทั้งสองtargetและjoinคุณสมบัติเป็นชุดข้อมูลเดียวกันและระบุเอาท์พุท ออกจากการตั้งค่าที่เหลือ

หลังจากนั้นซักครู่คุณควรจะเรียกคืน shapefile ที่มีคอลัมน์ "เข้าร่วมนับ" ลบ 1 จากสิ่งนี้ (เห็นได้ชัดว่าแต่ละคุณลักษณะควรตัดกันเอง) และนั่นควรเป็นจำนวน "ทับซ้อน" (จริง ๆ ตัดกัน) สำหรับแต่ละรูปหลายเหลี่ยม

ฉันทำมันต่อไป

ฉันดาวน์โหลดและทดลองใช้เครื่องมือ "Count Overlapping Polygons" อาจใช้งานได้ แต่ใช้เวลานานมาก (อาจเป็นเพราะขนาดไฟล์ แต่อินพุต FC ของฉันมีบันทึกน้อยกว่า 5,000 รายการ)

ในขณะที่ฉันกำลังรอให้เครื่องมือทำงานฉันเปิดหน้าต่าง ArcMap อื่นและใช้เวลาเพียงไม่กี่ก้าวในการรับสิ่งที่ฉันต้องการ 1) Spatial Join - ใช้คลาสคุณลักษณะเดียวกับ Target และ Join Features และเลือกตัวเลือก "เข้าร่วมหนึ่งถึงหลายคน" 2) Dissolve - การใช้เอาต์พุตจากขั้นตอนสุดท้าย ใช้ฟิลด์ "TARGET_FID" เป็นฟิลด์ละลายและสำหรับสถิติคุณสามารถ SUM ฟิลด์ "Join_Count" หรือ COUNT ในฟิลด์ "JOIN_FID" 3) ในไฟล์เอาต์พุตจากขั้นตอนที่ 2 ให้ใช้เครื่องคิดเลขลบฟิลด์ 1 จากฟิลด์สถิติ ("SUM_Join_Count" หรือ "COUNT_JOIN_FID") - เนื่องจากแต่ละฟีเจอร์จะตัดกัน

ฉันแนะนำให้ใช้วิธีนี้กับเครื่องมือ "Count Overlapping Polygon" ฉันเริ่มใช้เครื่องมือ COP ~ 5 นาทีก่อนที่จะเริ่มวิธีการเข้าร่วม -> นี้และมันให้ผลลัพธ์ฉันด้วยเวลาเพียงพอที่จะเขียนสิ่งนี้ก่อนที่เครื่องมือ "Count Overlapping Polygon" ก็เสร็จสิ้นแล้ว

หวังว่านี่จะช่วยได้!