First commit

Author: thiago

rotate.py

@@ -0,0 +1,58 @@
+from matplotlib import pyplot
+from shapely.geometry import LineString
+from shapely import affinity
+
+#from figures import SIZE, BLUE, GRAY
+
+GRAY='#858a8c'
+BLUE='#00b4fc'
+def add_origin(ax, geom, origin):
+    x, y = xy = affinity.interpret_origin(geom, origin, 2)
+    ax.plot(x, y, 'o', color=GRAY, zorder=1)
+    ax.annotate(str(xy), xy=xy, ha='center',
+                textcoords='offset points', xytext=(0, 8))
+
+
+def plot_line(ax, ob, color):
+    x, y = ob.xy
+    ax.plot(x, y, color=color, alpha=0.7, linewidth=3,
+            solid_capstyle='round', zorder=2)
+
+fig = pyplot.figure(1, dpi=90)
+
+line = LineString([(1, 3), (1, 1), (4, 1)])
+
+xrange = [0, 5]
+yrange = [0, 4]
+
+# 1
+ax = fig.add_subplot(121)
+
+plot_line(ax, line, GRAY)
+plot_line(ax, affinity.rotate(line, 90, 'center'), BLUE)
+add_origin(ax, line, 'center')
+
+ax.set_title(u"90\N{DEGREE SIGN}, default origin (center)")
+
+ax.set_xlim(*xrange)
+ax.set_xticks(range(*xrange) + [xrange[-1]])
+ax.set_ylim(*yrange)
+ax.set_yticks(range(*yrange) + [yrange[-1]])
+ax.set_aspect(1)
+
+# 2
+ax = fig.add_subplot(122)
+
+plot_line(ax, line, GRAY)
+plot_line(ax, affinity.rotate(line, 90, 'centroid'), BLUE)
+add_origin(ax, line, 'centroid')
+
+ax.set_title(u"90\N{DEGREE SIGN}, origin='centroid'")
+
+ax.set_xlim(*xrange)
+ax.set_xticks(range(*xrange) + [xrange[-1]])
+ax.set_ylim(*yrange)
+ax.set_yticks(range(*yrange) + [yrange[-1]])
+ax.set_aspect(1)
+
+pyplot.show()

skew.py

@@ -0,0 +1,78 @@
+from matplotlib import pyplot
+from shapely.wkt import loads as load_wkt
+from shapely import affinity
+from descartes.patch import PolygonPatch
+
+#from figures import SIZE, BLUE, GRAY
+GRAY='#858a8c'
+BLUE='#00b4fc'
+
+def add_origin(ax, geom, origin):
+    x, y = xy = affinity.interpret_origin(geom, origin, 2)
+    ax.plot(x, y, 'o', color=GRAY, zorder=1)
+    ax.annotate(str(xy), xy=xy, ha='center',
+                textcoords='offset points', xytext=(0, 8))
+
+fig = pyplot.figure(1, dpi=90)
+
+# Geometry from JTS TestBuilder with fixed precision model of 100.0
+# Using CreateShape > FontGlyphSanSerif and A = triangle.wkt from scale.py
+R = load_wkt('''\
+POLYGON((2.218 2.204, 2.273 2.18, 2.328 2.144, 2.435 2.042, 2.541 1.895,
+  2.647 1.702, 3 1, 2.626 1, 2.298 1.659, 2.235 1.777, 2.173 1.873,
+  2.112 1.948, 2.051 2.001, 1.986 2.038, 1.91 2.064, 1.823 2.08, 1.726 2.085,
+  1.347 2.085, 1.347 1, 1 1, 1 3.567, 1.784 3.567, 1.99 3.556, 2.168 3.521,
+  2.319 3.464, 2.441 3.383, 2.492 3.334, 2.536 3.279, 2.604 3.152,
+  2.644 3.002, 2.658 2.828, 2.651 2.712, 2.63 2.606, 2.594 2.51, 2.545 2.425,
+  2.482 2.352, 2.407 2.29, 2.319 2.241, 2.218 2.204),
+ (1.347 3.282, 1.347 2.371, 1.784 2.371, 1.902 2.378, 2.004 2.4, 2.091 2.436,
+  2.163 2.487, 2.219 2.552, 2.259 2.63, 2.283 2.722, 2.291 2.828, 2.283 2.933,
+  2.259 3.025, 2.219 3.103, 2.163 3.167, 2.091 3.217, 2.004 3.253, 1.902 3.275,
+  1.784 3.282, 1.347 3.282))''')
+
+xrange = [0, 5]
+yrange = [0, 4]
+
+# 1
+ax = fig.add_subplot(121)
+
+patch1a = PolygonPatch(R, facecolor=GRAY, edgecolor=GRAY,
+                       alpha=0.5, zorder=1)
+skewR = affinity.skew(R, xs=20, origin=(1, 1))
+patch1b = PolygonPatch(skewR, facecolor=BLUE, edgecolor=BLUE,
+                       alpha=0.5, zorder=2)
+ax.add_patch(patch1a)
+ax.add_patch(patch1b)
+
+add_origin(ax, R, (1, 1))
+
+ax.set_title("a) xs=20, origin(1, 1)")
+
+ax.set_xlim(*xrange)
+ax.set_xticks(range(*xrange) + [xrange[-1]])
+ax.set_ylim(*yrange)
+ax.set_yticks(range(*yrange) + [yrange[-1]])
+ax.set_aspect(1)
+
+# 2
+ax = fig.add_subplot(122)
+
+patch2a = PolygonPatch(R, facecolor=GRAY, edgecolor=GRAY,
+                       alpha=0.5, zorder=1)
+skewR = affinity.skew(R, ys=30)
+patch2b = PolygonPatch(skewR, facecolor=BLUE, edgecolor=BLUE,
+                       alpha=0.5, zorder=2)
+ax.add_patch(patch2a)
+ax.add_patch(patch2b)
+
+add_origin(ax, R, 'center')
+
+ax.set_title("b) ys=30")
+
+ax.set_xlim(*xrange)
+ax.set_xticks(range(*xrange) + [xrange[-1]])
+ax.set_ylim(*yrange)
+ax.set_yticks(range(*yrange) + [yrange[-1]])
+ax.set_aspect(1)
+
+pyplot.show()

union.py

@@ -0,0 +1,62 @@
+from matplotlib import pyplot
+from shapely.geometry import Point
+from descartes import PolygonPatch
+
+#from figures import SIZE, BLUE, GRAY
+
+fig = pyplot.figure(1, dpi=90)
+
+a = Point(1, 1).buffer(1.5)
+b = Point(2, 1).buffer(1.5)
+GRAY='#858a8c'
+BLUE='#00b4fc'
+
+# 1
+ax = fig.add_subplot(121)
+
+patch1 = PolygonPatch(a, fc=GRAY, ec=GRAY, alpha=0.2, zorder=1)
+ax.add_patch(patch1)
+patch2 = PolygonPatch(b, fc=GRAY, ec=GRAY, alpha=0.2, zorder=1)
+ax.add_patch(patch2)
+c = a.union(b)
+patchc = PolygonPatch(c, fc=BLUE, ec=BLUE, alpha=0.5, zorder=2)
+ax.add_patch(patchc)
+
+ax.set_title('a.union(b)')
+
+xrange = [-1, 4]
+yrange = [-1, 3]
+ax.set_xlim(*xrange)
+ax.set_xticks(range(*xrange) + [xrange[-1]])
+ax.set_ylim(*yrange)
+ax.set_yticks(range(*yrange) + [yrange[-1]])
+ax.set_aspect(1)
+
+def plot_line(ax, ob, color=GRAY):
+    x, y = ob.xy
+    ax.plot(x, y, color, linewidth=3, solid_capstyle='round', zorder=1)
+
+#2
+ax = fig.add_subplot(122)
+
+plot_line(ax, a.exterior)
+plot_line(ax, b.exterior)
+
+u = a.exterior.union(b.exterior)
+if u.geom_type in ['LineString', 'LinearRing', 'Point']:
+    plot_line(ax, u, color=BLUE)
+elif u.geom_type is 'MultiLineString':
+    for p in u:
+        plot_line(ax, p, color=BLUE)
+
+ax.set_title('a.boundary.union(b.boundary)')
+
+xrange = [-1, 4]
+yrange = [-1, 3]
+ax.set_xlim(*xrange)
+ax.set_xticks(range(*xrange) + [xrange[-1]])
+ax.set_ylim(*yrange)
+ax.set_yticks(range(*yrange) + [yrange[-1]])
+ax.set_aspect(1)
+
+pyplot.show()