Shapely: разбить LineString в произвольной точке вдоль края - PullRequest
Новая
       41

Shapely: разбить LineString в произвольной точке вдоль края

0 голосов
/ 14 мая 2018

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

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

Вот что я сделал до сих пор: 

from shapely.ops import split
from shapely.geometry import Point, LineString

line = LineString([(0, 0), (5,8)])
point = Point(2,3)

# Find coordinate of closest point on line to point
d = line.project(point)
p = line.interpolate(d)
print(p)
# >>> POINT (1.910112359550562 3.056179775280899)

# Split the line at the point
result = split(line, p)
print(result)
# >>> GEOMETRYCOLLECTION (LINESTRING (0 0, 5 8))

Спасибо!

1 Ответ

0 голосов
/ 15 мая 2018

Как оказалось, ответ, который я искал, был изложен в документации как метод cut: 

def cut(line, distance):
    # Cuts a line in two at a distance from its starting point
    if distance <= 0.0 or distance >= line.length:
        return [LineString(line)]
    coords = list(line.coords)
    for i, p in enumerate(coords):
        pd = line.project(Point(p))
        if pd == distance:
            return [
                LineString(coords[:i+1]),
                LineString(coords[i:])]
        if pd > distance:
            cp = line.interpolate(distance)
            return [
                LineString(coords[:i] + [(cp.x, cp.y)]),
                LineString([(cp.x, cp.y)] + coords[i:])]

Теперь я могу использовать project ed расстояниевырезать LineString: 

...
d = line.project(point)
# print(d) 3.6039927920216237

cut(line, d)
# LINESTRING (0 0, 1.910112359550562 3.056179775280899)
# LINESTRING (1.910112359550562 3.056179775280899, 5 8)
...