shapely descartes와 matplotlib로 단절된 실체를 그려야합니다

Question

다른 목적을 위해 매끈하게 만든 분리 된 원의 목록을 그려야합니다.

나는 http://toblerity.org/shapely/manual.html#cascading-unions 쇼 ( code 참조)의 예로서 정확하게하려고 노력했지만 그 원이 겹치는 경우에만 작동 전반적인 것은 (내 경우에는하지 않은) 연결되어 있습니다. 당신은 일이 데카르트 (하여 다각형 없었던

polygons = [Point(i, 0).buffer(0.7) for i in (0,4)] 

함께 나누기

및 AssertionError와 라인

polygons = [Point(i, 0).buffer(0.7) for i in range(5)] 

를 대체하여 볼 수 있듯이 또는 하나는 데카르트의 주장을 주석 처리하면 assert vertices.ndim == 2 실패하기 matplotlib로 테스트로)

matplotlib.path.Path에 대한 문서를 보면 MOVETO를 사용하여이 목표를 달성 할 수 있지만 매끄럽게 지원하지는 않습니다. 이 올바른지? 어떤 해결 방법이 있습니까?

Answer 1

다음 코드는 작동 :

from shapely.ops import cascaded_union 
from shapely.geometry import Point 
import random 
from matplotlib.patches import Polygon 
import pylab as pl 
import numpy as np 

circles = [Point(random.random(), random.random()).buffer(random.random() * 0.1) 
      for i in range(100)] 

polygons = cascaded_union(circles) 

fig, ax = pl.subplots(figsize=(8, 8)) 

for polygon in polygons: 
    mpl_poly = Polygon(np.array(polygon.exterior), facecolor="g", lw=0, alpha=0.4) 
    ax.add_patch(mpl_poly) 

ax.relim() 
ax.autoscale() 

출력 :