Shapely: 폴리 라벨 버그

https://github.com/mapbox/polylabel 자체에서 제안하는 좌표가 다른지 확인해 주시겠습니까?

TypeScript의 원래 알고리즘:

import polylabel from "polylabel";

let exterior = [[0,0],[10,0],[10,10],[0,10],[0,0]];
let hole = [[2,2],[6,2],[6,6],[2,6],[2,2]];
let polygon = [exterior, hole]

console.log(polylabel(polygon, 0.01));

PS> .\node_modules\.bin\tsc test.ts
PS> node .\test.js
[ 7.65625, 7.65625 ]

매끈한:

>>> from shapely.ops import polylabel
>>> from shapely.geometry import Polygon
>>> exterior = [[0,0],[10,0],[10,10],[0,10],[0,0]]
>>> hole = [[2,2],[6,2],[6,6],[2,6],[2,2]]
>>> p = Polygon(exterior, [hole]) 
>>> polylabel(p, 0.01).coords[:]
[(6.25, 3.75)]

당신이 맞습니다 @benediktbrandt. Shapely 구현은 폴리곤 외부까지의 거리만 고려하므로 구멍을 올바르게 처리하지 않습니다. 이것은 원래 알고리즘의 동작과 일치하도록 수정되어야 하는 버그라고 생각합니다.

https://github.com/Toblerity/Shapely/blob/cfa66d418fe3afe2fcb48438fcc9156142abaeb8/shapely/algorithms/polylabel.py#L44 -L53

폴리곤 외부를 라인스트링으로 변환하는 것은 생성된 모든 셀에 대해 수행되기 때문에 여기에서도 약간 비효율적입니다.