Pregunta:
Estoy usando ArcMap 10. Tengo dos conjuntos de datos, uno de los cuales consta de muchos polígonos contiguos pequeños y otro que tiene muchos polígonos no contiguos pequeños. El primer conjunto de datos son bloques censales y el segundo es un esquema de polígono de áreas urbanas generado a partir de los conjuntos de datos ráster que proporciona NLCD. Lo que quiero hacer es recortar cada bloque censal de modo que cada polígono de bloque censal solo contenga el área urbana especificada por el polígono NLCD. Ejecutar el clip con los bloques del censo como características de entrada me da esto. Sin embargo, existen restricciones adicionales: si el área de la intersección entre un bloque censal y el polígono NLCD está por debajo de un cierto umbral (por ejemplo, 1000 pies cuadrados), entonces el bloque censal completo debe mantenerse en la salida.
En este ejemplo, el polígono verde son los datos del NLCD y los contornos rojos son los bloques censales. Puedes ver el clip resultante.
El problema surge cuando no hay intersección entre el bloque censal y el polígono NLCD; no se guarda nada en la salida de ese bloque censal.
Necesito hacer esto para una gran cantidad de datos, probablemente en arcpy. ¿Qué estrategia sugerirías?
Respuesta:
Me acercaría a esto utilizando primero la herramienta Unión .
El resultado tendrá atributos que le permitirán distinguir tres conjuntos de polígonos, es decir, aquellos en los que las dos clases de entidades de entrada se superponen y aquellos en los que existe una u otra pero no se superponen con la otra.
Desde esta base, puede comenzar a aplicar sus otras reglas, tal vez usando una combinación de Seleccionar , Disolver y Agregar / Combinar .
Una vez que tenga su procedimiento, debería poder convertirlo en un modelo o secuencia de comandos de Python.