Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- from shapely.geometry import LineString
- line1 = LineString([(0, 0), (2, 2),(3,1)])
- line2 = LineString([(2, 0), (2, 1),(1,2)])
- print line1.intersection(line2)
- POINT (1.5 1.5)
- for line in line1.union(line2):
- print line
- LINESTRING (0 0, 1.5 1.5)
- LINESTRING (1.5 1.5, 2 2, 3 1)
- LINESTRING (2 0, 2 1, 1.5 1.5)
- LINESTRING (1.5 1.5, 1 2)
- # or
- from shapely.ops import unary_union
- for line in unary_union([line1,line2]):
- print line
- LINESTRING (0 0, 1.5 1.5)
- LINESTRING (1.5 1.5, 2 2, 3 1)
- LINESTRING (2 0, 2 1, 1.5 1.5)
- LINESTRING (1.5 1.5, 1 2)
- from shapely.ops import linemerge
- new_line1 = linemerge([line.union(line2)[0],line.union(line2)[1]])
- new_line2 = linemerge([line.union(line2)[2],line.union(line2)[3]])
- print new_line1
- LINESTRING (0 0, 1.5 1.5, 2 2, 3 1)
- print new_line2
- LINESTRING (2 0, 2 1, 1.5 1.5, 1 2)
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement