Existe-t-il une option python pour «joindre les attributs par emplacement»?

J'essaie d'exécuter la fonction join attributes by locationtelle que trouvée dans le menu QGIS Vector> Outils de gestion des données. Je cherche une option python open source pour cela. Je sais qu'arcpy a une spatial joinfonction, mais j'essaie de réaliser cela en dehors de l'environnement ESRI.

Vous voudrez peut-être jeter un œil à Shapely et Fiona . Fiona est un wrapper pour gdal pour faciliter l'importation et l'exportation de fichiers spatiaux. Shapely fournit des fonctionnalités de géométrie. Voici un exemple très simple pour vous donner l'idée. Il joint les attributs du polygone à tous les points de ce polygone.

Les données d'exemple que j'ai utilisées sont ces polygones et ces points .

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

Bien que toujours un peu rude sur les bords, surtout en ce qui concerne la documentation et les exemples, mais l'avenir des géopandas s'annonce brillant. Il combine essentiellement la puissance des trames de données pandas avec des capacités géospatiales de galbe.

la fonction que vous recherchez s'appelle sjoin

Assurez-vous que votre machine / instance dispose de suffisamment de mémoire pour effectuer l'opération

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