¿Cómo es que arcpy.Clip_analysis() es tan rápido frente a implementaciones que no son de arcpy?

Question

Estoy utilizando una versión de la solución 2 de la respuesta seleccionada aquí para construir un nuevo shapefile de todos los puntos que caen dentro de un objeto poligonal dado (es decir, un clip).

Cuando se ejecuta contra 27K puntos, tarda unos 9 minutos. En arcpy tarda menos de 10 segundos.

No puedo encontrar ninguna documentación sobre el Clip_analysis() (¿es propietario?) así que tengo curiosidad por saber cómo puede ejecutar las comprobaciones de 27K tan rápido?

¿Es multiproceso? ¿Está llamando a algún tipo de map() en la iterable de puntos?

Answer 1

Tengo entendido que los algoritmos de intersección de polígonos utilizados por la herramienta Clip están escritos a un nivel mucho más bajo que el de Python.

La herramienta Clip está disponible dentro del paquete del sitio ArcPy a través de una envoltura de Python.

Answer 2

Falta de un índice espacial. Estás evaluando cada punto para la intersección si es probable que sea un candidato o no.

Si quieres hacerlo en Python, debes usar algo como Rtree para eliminar rápidamente los puntos que podrían intersecar su polígono de recorte. Para un índice espacial puro de Python puedes probar PyQuadTree .