Apakah ada opsi python untuk "menggabungkan atribut berdasarkan lokasi"?

Saya mencoba menjalankan fungsi join attributes by locationseperti yang ditemukan pada menu QGIS, Vektor> Alat Manajemen Data. Saya mencari opsi python open source untuk ini. Saya tahu arcpy memiliki spatial joinfungsi tetapi saya mencoba melakukan ini di luar lingkungan ESRI.

Anda mungkin ingin melihat Shapely dan Fiona . Fiona adalah pembungkus untuk gdal untuk memudahkan impor dan ekspor file spasial. Shapely menyediakan fungsionalitas geometri. Ini adalah contoh yang sangat sederhana untuk memberi Anda ide. Ini menggabungkan atribut poligon ke semua titik dalam poligon itu.

Contoh data yang saya gunakan adalah poligon dan titik-titik ini .

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

Meskipun masih agak kasar, terutama dalam hal dokumentasi dan contoh, tetapi masa depan geopanda terlihat cerah. Ini pada dasarnya menggabungkan kekuatan panda dataframe dengan kemampuan geospasial dari rupawan.

fungsi yang Anda cari disebut sjoin

Pastikan mesin / instance Anda memiliki memori yang cukup untuk melakukan operasi

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