Shapely: parallel_offset 및 buffer가 윤곽선을 올바르게 닫지 않습니다.

예상되는 행동과 실제 행동.

parallel_offset 또는 buffer 를 사용할 때 결과 윤곽선이 올바르게 닫히지 않는 경우가 있습니다. 윤곽이 올바르게 닫히지 않은 지점은 항상 윤곽의 첫 번째 / 마지막 지점 인 것처럼 보입니다.

다음은 처음에 StackExchange에 게시 한 문제입니다.

문제를 재현하는 단계.

다음은 parallel_offset 사용할 때 닫히지 않는 윤곽선을 생성하는 최소한의 작업 예제입니다. buffer 하면 mitred 스타일에도 불구하고 둥근 모서리가 생성됩니다.
흥미롭게도 더 큰 데이터 세트에는 성공적으로 닫히는 다른 매우 유사한 윤곽이 있습니다.

import matplotlib.pyplot as plt
from shapely.geometry.polygon import LinearRing

def plot_line(ax, ob, color):
    x, y = ob.xy
    ax.plot(x, y, color=color, alpha=0.7, linewidth=3, 
            solid_capstyle='round', zorder=2)

polygon = [[-29.675, -30.675],
           [-28.4094, -29.4094],
           [-28.325, -29.325],
           [-28.325, -29.764],
           [-28.325, -29.7933],
           [-28.4587, -29.8274],
           [-28.4676, -29.8297],
           [-28.5956, -29.8814],
           [-28.6041, -29.8848],
           [-28.724, -29.953],
           [-28.732, -29.9576],
           [-28.8417, -30.0413],
           [-28.849, -30.0469],
           [-28.9466, -30.1445],
           [-28.9531, -30.151],
           [-29.0368, -30.2607],
           [-29.0424, -30.268],
           [-29.1106, -30.3879],
           [-29.1152, -30.3959],
           [-29.1669, -30.5239],
           [-29.1703, -30.5324],
           [-29.2044, -30.6661],
           [-29.2067, -30.675],
           [-29.6457, -30.675],
           [-29.675, -30.675]]

poly_line = LinearRing(polygon)
poly_line_offset = poly_line.parallel_offset(0.05, side="left", resolution=16, 
                                             join_style=2, mitre_limit=1)
# Alternative:
# poly_line_offset = poly_line.buffer(0.05, resolution=16, join_style=2, mitre_limit=1).exterior

fig = plt.figure()
ax = fig.add_subplot(111)
plot_line(ax, poly_line, "blue")
plot_line(ax, poly_line_offset, "green")
plt.show()

운영 체제

Windows 7, Python 3.6.2

매끈한 버전과 출처

Shapely 1.6.3 (PyPI에서 설치)