Guest User

Untitled

a guest
Nov 21st, 2017
95
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.56 KB | None | 0 0
  1. from shapely import wkt
  2. from shapely.ops import linemerge, unary_union, polygonize
  3.  
  4. POLY = "POLYGON ((34.67491149902344 31.59900710035676, 34.85000610351562 31.59900710035676, 34.85000610351562 31.73867905688433, 34.67491149902344 31.73867905688433, 34.67491149902344 31.59900710035676))"
  5. LINE = "LINESTRING (34.64401245117188 31.63292168314889, 34.80812072753906 31.75911546882192)"
  6.  
  7. poly = wkt.loads(POLY)
  8. line = wkt.loads(LINE)
  9.  
  10. merged = linemerge([poly.boundary, line])
  11. borders = unary_union(merged)
  12. polygons = polygonize(borders)
  13. for p in polygons:
  14. print(p)
Add Comment
Please, Sign In to add comment