GIS/ArcPy: Landnutzungsdaten aus einem Raster extrahieren und Shigfile Polygonen zuweisen

Question

Ich habe zwei Karten für die gleiche Fläche (1) eine Raster Landnutzungskarte und (2) ein Shapefile mit Tausenden von Teilwassereinzugsgebieten. Ich versuche, den Landnutzungstyp aus dem Raster (Karte 1) jedem Teilwassereinzugsgebiet (Karte 2) basierend auf der Mehrheitsregel zuzuordnen. Ich habe versucht, räumliche Join, aber das Ergebnis scheint falsch zu sein. Was ist der beste Weg, um es entweder in ArcMap oder durch Arcpy zu tun?

Answer 1

Ich würde empfehlen, Zonal Statistics as Table (Spatial Analyst) zu verwenden, um diese Aufgabe zu erfüllen. Hier ist der allgemeine Arbeitsablauf:

1. Führen Sie Zonal Statistics as Table mithilfe der Statistik "Majority".
2. Treten Sie der Tabelle mit dem Wendepunkt Klasse Feature Join Field (Data Management)

---

import arcpy, os 
from arcpy.sa import * 
arcpy.CheckOutExtension("Spatial") 

# Your watershed feature class 
watersheds = r'C:\path\to\your\geodatabase.gdb\watersheds' 

# Your land cover raster 
raster = r'C:\path\to\your\landcover_raster.tif' 

# The workspace where the output table will go 
zone_table = r'C:\path\to\your\geodatabase.gdb' 

# Perform the zonal statistics and output a table 
arcpy.sa.ZonalStatisticsAsTable (watersheds, 'watershed_id', raster, zone_table, 'DATA', 'MAJORITY') 

# Join the table to the watershed feature class by the OBJECTID field (for feature class) 
arcpy.JoinField_management(watersheds, 'watershed_id', zone_table, 'OBJECTID', 'MAJORITY')