Existe uma opção python para "juntar atributos por local"?

Estou tentando executar a função join attributes by locationconforme encontrado no menu QGIS Vector> Data Management Tools. Estou procurando uma opção python de código aberto para isso. Eu sei que o arcpy tem uma spatial joinfunção, mas estou tentando executar isso fora do ambiente ESRI.

Você pode dar uma olhada em Shapely e Fiona . O Fiona é um invólucro para o gdal para facilitar a importação e exportação de arquivos espaciais. Shapely fornece funcionalidade de geometria. Aqui está um exemplo muito simples para lhe dar a ideia. Ele une atributos de polígono a todos os pontos dentro desse polígono.

Os dados de exemplo que usei são esses polígonos e esses pontos .

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']
                        })

Embora ainda seja um pouco difícil, especialmente quando se trata de documentação e exemplos, o futuro das geopandas parece brilhante. Basicamente, combina o poder dos quadros de dados do panda com os recursos geoespaciais de bem torneados.

a função que você procura é chamada de junção

Verifique se a sua máquina / instância tem memória suficiente para executar a operação

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')