Access individual points in Shapely MultiPoint

Question

I am working with the Shapely library in Python. I find the intersection of two lines, the return value is given as a MultiPoint object.

How do I deconstruct the object to get the individual points in the intersection?

Here is the code:

from shapely.geometry import LineString, MultiLineString
a = LineString([(0, 1), (0, 2), (1, 1), (2, 0)])
b = LineString([(0, 0), (1, 1), (2, 1), (2, 0)])
x = a.intersection(b)

Output:

print(x) 
MULTIPOINT (1 1, 2 0)

So, in this case, I'd be looking for a way to extract the intersection points (1,1) and (2,0).

score 6 · Accepted Answer · answered Jun 27 '18 at 10:50

6

You can index the resulting MultiPoint:

>>> str(x)
'MULTIPOINT (1 1, 2 0)'
>>> print(len(x))
2
>>> print(x[0].x)
1.0
>>> print(x[0].y)
1.0

If you want a new list with the coordinates, you can use:

>>> [(p.x, p.y) for p in x]
[(1.0, 1.0), (2.0, 0.0)]

answered Jun 27 '18 at 10:50

jjmontes

24,679
4
39
51

What about if there are 3 intersection points, you obtain a MultiPoint and a LineString? – Chris B Jun 27 '18 at 17:56
5

I have installed `shapely==2.0.0` and `MultiPoint` is not iterable or subscriptable now. :/ – s.paszko Jan 12 '23 at 11:13
4

Try x.geoms[0] (see [here](https://shapely.readthedocs.io/en/latest/reference/shapely.MultiPoint.html) ) – Max Peglar-Willis Jan 13 '23 at 21:08

score 1 · Answer 2 · answered Jun 23 '23 at 13:34

Use .geoms:

from shapely.geometry import LineString
a = LineString([(0, 1), (0, 2), (1, 1), (2, 0)])
b = LineString([(0, 0), (1, 1), (2, 1), (2, 0)])

multipoint = a.intersection(b)
print(multipoint)
#MULTIPOINT (2 0, 1 1)
points = [p for p in multipoint.geoms]
print(points)
#[<POINT (2 0)>, <POINT (1 1)>]

Access individual points in Shapely MultiPoint

2 Answers2