"위치별로 속성에 가입"하는 Python 옵션이 있습니까?

join attributes by locationQGIS 메뉴 Vector> Data Management Tools에 있는 기능을 수행하려고합니다 . 이를 위해 오픈 소스 파이썬 옵션을 찾고 있습니다. arcpy에 spatial join기능이 있지만 ESRI 환경 외부 에서이 작업을 수행하려고합니다.

Shapely 와 Fiona를 살펴볼 수 있습니다 . Fiona는 공간 파일 가져 오기 및 내보내기를 쉽게하기 위해 gdal의 래퍼입니다. 지오메트리 기능을 제공합니다. 다음은 아이디어를 제공하는 매우 간단한 예입니다. 다각형 속성을 해당 다각형 내의 모든 점에 결합합니다.

내가 사용한 예제 데이터는 이러한 다각형 및 이러한 점 .

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

특히 문서 및 예제와 관련하여 여전히 가장자리에서 약간 거칠지 만 geopandas의 미래는 밝게 보입니다. 기본적으로 팬더 데이터 프레임의 기능과 지형 공간 기능을 매끄럽게 결합합니다.

찾는 기능을 sjoin 이라고 합니다

기계 / 인스턴스에 작업을 수행하기에 충분한 메모리가 있는지 확인하십시오

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