Có một tùy chọn python để tham gia các thuộc tính theo khu vực khác không?

Tôi đang cố gắng thực hiện chức năng join attributes by locationnhư được tìm thấy trên menu QGIS Vector> Công cụ quản lý dữ liệu. Tôi đang tìm kiếm một tùy chọn python mã nguồn mở cho việc này. Tôi biết arcpy có một spatial joinchức năng nhưng tôi đang cố gắng thực hiện điều này bên ngoài môi trường ESRI.

Bạn có thể muốn xem Shapely và Fiona . Fiona là một trình bao bọc cho gdal để làm cho việc nhập và xuất tệp không gian dễ dàng. Shapely cung cấp chức năng hình học. Dưới đây là một ví dụ rất đơn giản để cung cấp cho bạn ý tưởng. Nó tham gia các thuộc tính đa giác cho tất cả các điểm trong đa giác đó.

Dữ liệu ví dụ tôi đã sử dụng là những đa giác và những điểm này .

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

Mặc dù vẫn còn một chút khó khăn xung quanh các cạnh, đặc biệt là khi nói đến tài liệu và ví dụ, nhưng tương lai của geopandas có vẻ tươi sáng. Về cơ bản, nó kết hợp sức mạnh của các cơ sở dữ liệu gấu trúc với khả năng không gian địa lý của hình dạng.

chức năng bạn tìm kiếm được gọi là sjoin

Đảm bảo máy / thể hiện của bạn có đủ bộ nhớ để thực hiện thao tác

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