在线条与多边形相交的每一点上拆分线条

Question

我有一个横跨各种多边形的线串，存储为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))

我期望有三条线，一条存在于正方形多边形内，然后两条分别存在于多边形外。

Answer 1

看起来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))