Shapely: Ошибка с полилиблической этикеткой

Созданный на 8 янв. 2020 · 3Комментарии · Источник: Toblerity/Shapely

Ожидаемое поведение и фактическое поведение.

Отверстия не учитываются правильно в реализации алгоритма полилиблевой метки shapely.

Шаги по воспроизведению проблемы.

from shapely.geometry import Polygon
from shapely.algorithms.polylabel import polylabel
poly1=Polygon(shell=[(0,0),(10,0),(10,10),(0,10),(0,0)][::-1], holes=[[(2,2),(6,2),(6,6),(2,6),(2,2)]])
import geopandas
gdf2 = geopandas.GeoDataFrame(geometry=[poly1])
gdf2.plot()
label_location = polylabel(poly1, tolerance=0.05) 
print(label_location.x)
print(label_location.y)

Результатом этого сценария является:

6.25
3.75

Но правильная позиция - что-то близкое к 7,65625, 7,65625. В этом можно убедиться, посмотрев на график этого многоугольника.

Возможно, я делаю что-то не так в определении многоугольника, но я пробовал много разных способов создания экземпляра этого многоугольника и всегда получаю один и тот же результат.

Операционная система

Linux

Элегантная версия и происхождение

Протестировано с shapely версии 1.6 и 1.7a3

bug

Источник

benediktbrandt

Самый полезный комментарий

Вы правы @benediktbrandt. Реализация Shapely неправильно обрабатывает отверстия, поскольку учитывает только расстояние до внешней части многоугольника. Я думаю, что это ошибка, и ее следует исправить, чтобы она соответствовала поведению исходного алгоритма.

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

Преобразование внешнего полигона в линейную линию здесь тоже немного неэффективно, так как оно будет выполняться для каждой созданной ячейки.

snorfalorpagus 21 февр. 2020

❤1 👍1

Все 3 Комментарий

Не могли бы вы убедиться, что координаты, предлагаемые самим https://github.com/mapbox/polylabel , отличаются?

kannes 6 февр. 2020

Оригинальный алгоритм в 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)]

snorfalorpagus 21 февр. 2020

👍2

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

snorfalorpagus 21 февр. 2020

❤1 👍1

Была ли эта страница полезной?

0 / 5 - 0 рейтинги