¿Existe una opción de python para "unir atributos por ubicación"?

Estoy tratando de realizar la función join attributes by locationque se encuentra en el menú QGIS Vector> Herramientas de administración de datos. Estoy buscando una opción de python de código abierto para esto. Sé que arcpy tiene una spatial joinfunción, pero estoy tratando de realizar esto fuera del entorno ESRI.

Es posible que desee echar un vistazo a Shapely y Fiona . Fiona es un contenedor de gdal para facilitar la importación y exportación de archivos espaciales. Shapely proporciona funcionalidad de geometría. Aquí hay un ejemplo muy simple para darle la idea. Une atributos de polígono a todos los puntos dentro de ese polígono.

Los datos de ejemplo que he usado son estos polígonos y estos puntos .

import fiona
from shapely.geometry import shape
from copy import deepcopy

with fiona.open("planning_neighborhoods.shp", "r") as n: 

    with fiona.open("Schools_Private_Pt.shp", "r") as s:

        # create a schema for the attributes
        outSchema =  deepcopy(s.schema)
        outSchema['properties'].update(n.schema['properties'])

        with fiona.open ("Schools_withNbhd.shp", "w", s.driver, outSchema, s.crs) as output:

            for school in s: 
                for neighborhood in n:
                    # check if point is in polygon and set attribute
                    if shape(school['geometry']).within(shape(neighborhood['geometry'])):  
                        school['properties']['neighborho'] = neighborhood['properties']['neighborho'] 
                    # write out
                        output.write({                                 
                            'properties': school['properties'], 
                            'geometry': school['geometry']
                        })

Aunque todavía es un poco difícil, especialmente cuando se trata de documentación y ejemplos, pero el futuro de las geopandas parece brillante. Básicamente combina el poder de los marcos de datos de pandas con las capacidades geoespaciales de Shapely.

la función que busca se llama sjoin

Asegúrese de que su máquina / instancia tenga suficiente memoria para realizar la operación

import geopandas as gpd
import pandas as pd
import os

gdfLeft = gpd.read_file(os.path.join(PATH,INPUT_FILE_NAME_1))
gdfRight = gpd.read_file(os.path.join(PATH,INPUT_FILE_NAME_2))

gdfJoined = gpd.sjoin(gdfLeft, gdfRight, how="left", op='intersects')