Разделите линию в каждой точке, в которой она пересекает многоугольник.

У меня есть строка, которая охватывает различные многоугольники, хранящиеся как GeoJsons. Я хочу разбить линию на отдельные части в пределах каждой области многоугольника. Однако мне пока не удалось этого добиться. Это воспроизводимый пример того, что у меня есть:

from shapely.geometry import Point, LineString, Polygon
from shapely.ops import split
from matplotlib import pyplot as plt

poly = {
    'type': "Feature",
    'geometry': {
        "type": "Polygon",
        "coordinates": ([(2,2),(2,4),(4,4),(4,2)]),
    },
    'properties': {
        'id': 'A'
    }
}

line = {
    'type': "Feature",
    'geometry': {
        "type": "Linestring",
        "coordinates": ([(3,3),(5,1)]),
    },
    'properties': {
        'id': 'A'
    }
}

poly = Polygon(poly['geometry']['coordinates'])
line = LineString(line['geometry']['coordinates'])

x,y = poly.exterior.xy
x2,y2 = line.coords

plt.plot(x, y, x2, y2)
plt.show()

Этот код создает следующий квадратный многоугольник с пересекающей его линией: введите описание изображения здесь

Затем я попытался разделить линию через многоугольник, вот так:

new_lines = split(line, poly)

Но я получаю следующий результат, который не кажется правильным:

GEOMETRYCOLLECTION (LINESTRING (3 3, 4 2), LINESTRING (4 2, 5 1))

Я ожидал, что будут три линии, одна будет внутри квадратного многоугольника, а две - отдельно за пределами многоугольника.


python geometry shapely
person Thirst for Knowledge    schedule 03.02.2019    source источник

Ответы (1)


arrow_upward
2
arrow_downward

Похоже, что split работает так, как вы ожидаете, но координаты линии нанесены неверно - следует использовать .xy, чтобы получить линию xs и ys, как с многоугольником:

from shapely.geometry import Point, LineString, Polygon
from shapely.ops import split
from matplotlib import pyplot as plt

poly = Polygon([(2,2),(2,4),(4,4),(4,2)])
original_line = LineString([(3,3),(5,1)])
line = LineString([(3,1), (3,5)])
x, y = poly.exterior.xy
x2, y2 = line.xy
x3, y3 = original_line.xy

plt.plot(x, y, x2, y2)
plt.show()

plt.plot(x, y, x3, y3)
plt.show()

изображение с исходными координатами

изображение с ожидаемыми координатами

>>> print(split(line, poly))
GEOMETRYCOLLECTION (LINESTRING (3 1, 3 2), LINESTRING (3 2, 3 4), LINESTRING (3 4, 3 5))
person tlrss    schedule 04.02.2019