¿Intercambio de categorías en un shapefile con ArcGIS for Desktop?

Question

Tengo un shapefile compuesto por polígonos de cicatrices de incendios identificables como categorías únicas basadas en el año del incendio.

Me interesa medir la superficie de las cicatrices de incendios superpuestas de diferentes años.

¿Es posible realizar una intersección de categorías dentro de un shapefile o tendré que dividir primero cada categoría en un shapefile separado?

Answer 1

La respuesta corta es que sí, pero la forma de hacerlo varía según el software que utilices (QGIS, GDAL & OGR, SAGA, etc). Voy a suponer que está utilizando ArcGIS 10.x:

Si se carga el shapefile dos veces en la misma sesión (es decir, dos capas, ambas leyendo de la misma fuente de datos shapefile), y se establecen diferentes Consultas de Definición en ellas para diferentes años de incendio, se puede entonces realizar la Intersección en estas dos capas. Se puede hacer todo esto a través de la interfaz gráfica de usuario, por supuesto, o en Python podría ser algo como lo siguiente:

import arcpy

# Two layers from the same shapefile
arcpy.MakeFeatureLayer_management(r"C:\aDir\myshapefile.shp","test_lyr1")
arcpy.MakeFeatureLayer_management(r"C:\aDir\myshapefile.shp","test_lyr2")
lyr1 = arcpy.mapping.Layer("test_lyr1")
lyr2 = arcpy.mapping.Layer("test_lyr2")

# Name the layers
lyr1.name = "lyr1"
lyr2.name = "lyr2"

# Set definition queries for each layer
lyr1.definitionQuery = '"FIREYEAR" = ' + "'2010'"
lyr2.definitionQuery = '"FIREYEAR" = ' + "'2011'"

# Perform the intersection on the named layers
# You can intersect among more than two layers at a time if 
# you've got the full ArcGIS license.
arcpy.Intersect_analysis(["lyr1", "lyr2"], r"c:\aDir\Intersection.shp")

Answer 2

También puede utilizar simplemente el Intersect con un conjunto de datos como entrada, como se explica en Cómo funciona Intersect :

Intersect puede funcionar con una sola entrada. En este caso, en lugar de descubrir las intersecciones entre las características de las diferentes clases de características o capas, descubrirá las intersecciones entre características dentro de la entrada única. Esto puede ser útil para descubrir superposición de polígonos e intersecciones de líneas (como puntos o líneas).