Seleccionar por ubicación en Python

Question

Tengo una capa de polígonos y una capa de puntos (con coordenadas para cada punto). La capa de polígonos se creó con el plugin MMQGIS (polígonos hexagonales). Quiero seleccionar todos los puntos dentro de cada polígono.

Antes trabajaba con ArcGIS pero quiero cambiar a código abierto. Además, quiero abordarlo a través de un script de Python ya que quiero analizar los datos extraídos. Tengo instalado QGIS 2.8. Estoy un poco abrumado por todas las librerías GIS que hay para Python.

¿Puedo abordar ese tipo de selección con una biblioteca SIG que ya esté instalada con QGIS? Si es así, ¿podría alguien indicarme la dirección correcta?

Answer 1

Suponiendo que quiera ejecutar su script dentro de QGIS (desde un archivo de script o la consola de Python), puede utilizar lo siguiente:

import processing

processing.runalg("qgis:selectbylocation", INPUT, INTERSECT, METHOD, OUTPUT)

Aquí está la descripción de ayuda proporcionada por la consola de Python que define cada parámetro:

processing.alghelp("qgis:selectbylocation")
ALGORITHM: Select by location
    INPUT <ParameterVector>
    INTERSECT <ParameterVector>
    METHOD <ParameterSelection>
    OUTPUT <OutputVector>

METHOD(Modify current selection by)
    0 - creating new selection
    1 - adding to current selection
    2 - removing from current selection

Para crear un script completamente desde cero que no llame a las funciones de QGIS, los siguientes enlaces pueden ayudar a desarrollar una funcionalidad similar:

• Desde el plugin Processing

• Desde el plugin fTools

---

EDITAR :

Hay un par de posts que podrían ser de alguna utilidad para escribir scripts fuera de QGIS para Mac:

• ¿Cómo ejecutar un simple script de python para QGIS desde fuera (por ejemplo, Sublime Text)?

• Fijación de la ruta de Python para los scripts de QGIS en Mac OSX

Answer 2

La respuesta de Joseph es correcta, aunque es necesario incluir otro parámetro: "PREDICADO".

es decir:

processing.run("qgis:selectbylocation", {
    "INPUT":lyr_input,\
    "PREDICATE":0,\
    "INTERSECT":lyr_intersect,\
    "METHOD":0,\
    "OUTPUT":path}
)

Desde processing.algorithmHelp("qgis:selecybylocation") :

PREDICATE: Where the features (geometric predicate)

    Parameter type: QgsProcessingParameterEnum

    Available values:
        - 0: intersect
        - 1: contain
        - 2: disjoint
        - 3: equal
        - 4: touch
        - 5: overlap
        - 6: are within
        - 7: cross

    Accepted data types:
        - int
        - str: as string representation of int, e.g. '1'
        - QgsProperty