Add sierpinski carpet method to geometry

Author: jericks

examples/sierpinskiCarpet.py

@@ -0,0 +1,6 @@
+from geoscript import geom
+from geoscript import viewer
+
+bounds = geom.Bounds(-180,-90,180,90)
+sierpinskiCarpet = geom.createSierpinskiCarpet(bounds, 50)
+viewer.draw(sierpinskiCarpet)

geoscript/geom/__init__.py

@@ -11,7 +11,7 @@
 from compoundring import CompoundRing
 from bounds import Bounds
 from geom import Geometry
-from geom import prepare, simplify, densify, transform, buffer, delaunay, voronoi, randomPoints
+from geom import prepare, simplify, densify, transform, buffer, delaunay, voronoi, randomPoints, createSierpinskiCarpet
 from io.wkt import readWKT, writeWKT, fromWKT
 from io.wkb import readWKB, writeWKB, fromWKB, toWKB
 from io.json import writeJSON, readJSON

geoscript/geom/geom.py

@@ -15,6 +15,7 @@
 from org.geotools.geometry.jts import JTS
 from org.geotools.referencing.operation.transform import AffineTransform2D
 from org.locationtech.jts.shape.random import RandomPointsBuilder
+from org.locationtech.jts.shape.fractal import SierpinskiCarpetBuilder
 from geoscript.geom.bounds import Bounds
 
 _factory = GeometryFactory()
@@ -167,6 +168,19 @@ def randomPoints(g, number):
   builder.numPoints = number
   return builder.getGeometry()
 
+def createSierpinskiCarpet(b, numberOfPoints):
+  """
+  Create a Sierpinski Carpet
+
+  *g* The Bounds
+
+  *number* The number of Points
+  """
+  builder = SierpinskiCarpetBuilder(GeometryFactory())
+  builder.extent = b
+  builder.numPoints = numberOfPoints
+  return builder.getGeometry()
+
 def _bounds(g):
   return Bounds(env=g.getEnvelopeInternal())
 

tests/test_geom.py

@@ -231,3 +231,8 @@ def testRandomPoints(self):
     polygon = geom.buffer(geom.Point(1153461.34, 649950.30), 100)
     pts = geom.randomPoints(polygon, 10)
     assert pts.getNumGeometries() == 10
+
+  def testCreateSierpinskiCarpet(self):
+    bounds = geom.Bounds(-180,-90,180,90)
+    sierpinskiCarpet = geom.createSierpinskiCarpet(bounds, 50)
+    assert sierpinskiCarpet.getNumGeometries() > 0