파이썬에서 여러 다각형의 교차로를 효율적으로 얻기

여러 다각형의 교차점을 얻고 싶습니다. 파이썬 shapely패키지를 사용하면 intersection함수를 사용하여 두 다각형의 교차점을 찾을 수 있습니다 . 여러 다각형의 교집합을 얻는 비슷한 효율적인 기능이 있습니까?

다음은 내가 의미하는 바를 이해하는 코드 스 니펫입니다.

from shapely.geometry import Point

coord1 = ( 0,0 )
point1 = Point(coord1)
circle1 = point1.buffer(1)

coord2 = ( 1,1 )
point2 = Point(coord2)
circle2 = point2.buffer(1)

coord3 = ( 1,0 )
point3 = Point(coord3)
circle3 = point3.buffer(1) 

두 서클의 교차점은에서 찾을 수 있습니다 circle1.intersection(circle2). 로 세 서클의 교차점을 모두 찾을 수 있습니다 circle1.intersection(circle2).intersection(circle3). 그러나이 방법은 점점 더 많은 코드가 필요하기 때문에 많은 수의 다각형에 적용 할 수 없습니다. 임의의 수의 다각형을 가져 와서 그 교차점을 반환하는 함수를 원합니다.

하나의 가능한 접근법은 다각형 쌍, 그 교차점 및 마지막으로 계단식 유니온을 통해 모든 교차점의 조합을 고려하는 것입니다 ( 여기서 제안 됨 ).

from shapely.geometry import Point
from shapely.ops import cascaded_union
from itertools import combinations

circles = [
    Point(0,0).buffer(1),
    Point(1,0).buffer(1),
    Point(1,1).buffer(1),
]

intersection = cascaded_union(
    [a.intersection(b) for a, b in combinations(circles, 2)]
)
print intersection

보다 효율적인 접근법은 많은 기하학을 처리하기 위해 Rtree 와 같은 공간 인덱스를 사용해야합니다 (세 개의 원이 아닌 경우).

from shapely.geometry import Point
from shapely.ops import cascaded_union
from rtree import index

circles = [
    Point(0,0).buffer(1),
    Point(1,0).buffer(1),
    Point(1,1).buffer(1),
]
intersections = []
idx = index.Index()

for pos, circle in enumerate(circles):
    idx.insert(pos, circle.bounds)

for circle in circles:
    merged_circles = cascaded_union([circles[pos] for pos in idx.intersection(circle.bounds) if circles[pos] != circle])
    intersections.append(circle.intersection(merged_circles))

intersection = cascaded_union(intersections)
print intersection

반복 또는 재귀를 사용하지 않는 이유는 무엇입니까? 같은 :

from shapely.geometry import Point

def intersection(circle1, circle2):
    return circle1.intersection(circle2)

coord1 = ( 0,0 )
point1 = Point(coord1)
circle1 = point1.buffer(1)

coord2 = ( 1,1 )
point2 = Point(coord2)    
circle2 = point2.buffer(1)


coord3 = ( 1,0 )
point3 = Point(coord3)
circle3 = point3.buffer(1)
circles = [circle1, circle2, circle3]
intersectionResult = None

for j, circle  in enumerate(circles[:-1]):

    #first loop is 0 & 1
    if j == 0:
        circleA = circle
        circleB = circles[j+1]
     #use the result if the intersection
    else:
        circleA = intersectionResult
        circleB = circles[j+1]
    intersectionResult = intersection(circleA, circleB)

result= intersectionResult

이 코드에 샷을 제공하십시오. 그것의 개념은 매우 간단하고 나는 당신이 찾고있는 것을 얻습니다.

from shapely.geometry import Point
from itertools import combinations
dic ={}
dic['coord1']=Point(0,0).buffer(1)
dic['coord2']=Point(1,1).buffer(1)
dic['coord3']=Point(1,0).buffer(1)
inter = {k[0]+v[0]:k[1].intersection(v[1]) for k,v in combinations(dic.items(),2)}
print inter

그리고 출력을 shapefile로 저장하려면 fiona를 사용하십시오.

from shapely.geometry import Point,mapping
import fiona
from itertools import combinations
schema = {'geometry': 'Polygon', 'properties': {'Place': 'str'}}
dic ={}
dic['coord1']=Point(0,0).buffer(1)
dic['coord2']=Point(1,1).buffer(1)
dic['coord3']=Point(1,0).buffer(1)
inter = {k[0]+v[0]:k[1].intersection(v[1]) for k,v in combinations(dic.items(),2)}
print inter
with fiona.open(r'C:\path\abid', "w", "ESRI Shapefile", schema) as output:
    for x,y in inter.items():
        output.write({'properties':{'Place':x},'geometry':mapping(y)})

이 출력-