Guest User

Untitled

a guest
Jan 18th, 2019
84
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 1.42 KB | None | 0 0
  1. SELECT id, ST_Buffer ((geom), 0.0001,'endcap=round join=round') as geom
  2. FROM roads_split;
  3.  
  4. SELECT a.id,
  5. ST_LineMerge(
  6. ST_Difference(ST_ExteriorRing(a.geom), b.geom)
  7. ) geom
  8. FROM roads_union_byf a;
  9. LEFT JOIN roads_union_byf b ON (ST_Intersects(a.geom, b.geom)
  10. AND a.id < b.id);
  11.  
  12. create table roads_cut_poly_l as
  13. SELECT id, ST_MakePolygon(ST_AddPoint(geom, ST_StartPoint(geom))) geom
  14. FROM roads_cut_l;
  15.  
  16. create table roads_cut_r as
  17. SELECT a.id,
  18. ST_LineMerge(
  19. ST_Difference(ST_ExteriorRing(a.geom), b.geom)
  20. ) geom
  21. FROM roads_union_byf a
  22. LEFT JOIN roads_union_byf b ON (ST_Intersects(a.geom, b.geom)
  23. AND a.id > b.id);
  24.  
  25. create table roads_cut_poly_r as
  26. SELECT id, ST_MakePolygon(ST_AddPoint(geom, ST_StartPoint(geom))) geom
  27. FROM roads_cut_r;
  28.  
  29. create table roads_union_split as
  30. SELECT ST_Intersection(a.geom, b.geom) as geom
  31. FROM roads_cut_poly_r a, roads_cut_poly_l b
  32. WHERE ST_Intersects (a.geom, b.geom)
  33. AND a.id=b.id;
  34.  
  35. create table roadsegment_pt_split1 AS
  36. SELECT ST_LineInterpolatePoint (roadsegments_marge.ggeom, 0.1) as geom
  37. FROM roadsegments_marge
  38.  
  39. create table roadsegment_part1 as
  40. SELECT ST_MakeLine (geom) as geom
  41. FROM roadsegments_pt_split;
  42.  
  43. create table roadsegment_part2 as
  44. SELECT ST_MakeLine (geom) as geom
  45. FROM roadsegment_pt_split2;
  46.  
  47. create table roadsegment_buffer_part1 as
  48. SELECT ST_Buffer ((geom), 0.0001, 'endcap = flat join = round') as geom
  49. FROM roadsegment_part1
Add Comment
Please, Sign In to add comment