Add randomPoints function to the geom module

Author: jericks

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
+from geom import prepare, simplify, densify, transform, buffer, delaunay, voronoi, randomPoints
 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

@@ -14,6 +14,7 @@
 from org.locationtech.jts.operation.buffer import BufferOp, BufferParameters
 from org.geotools.geometry.jts import JTS
 from org.geotools.referencing.operation.transform import AffineTransform2D
+from org.locationtech.jts.shape.random import RandomPointsBuilder
 from geoscript.geom.bounds import Bounds
 
 _factory = GeometryFactory()
@@ -153,6 +154,19 @@ def buffer(g, distance, singleSided=False):
 
   return BufferOp.bufferOp(g, distance, bp) 
 
+def randomPoints(g, number):
+  """
+  Create a number of random Points within the given Geometry
+
+  *g* The constraining Geometry
+
+  *number* The number of Points to create
+  """
+  builder = RandomPointsBuilder(GeometryFactory())
+  builder.setExtent(g)
+  builder.numPoints = number
+  return builder.getGeometry()
+
 def _bounds(g):
   return Bounds(env=g.getEnvelopeInternal())
 

tests/test_geom.py

@@ -227,4 +227,8 @@ def testSubLine(self):
     subLine = line.subLine(0.33, 0.67)
     assert "LINESTRING (1156010.153864557 649246.3016361536, 1175115.6870342216 648021.5879714314)" == str(subLine)
 
-
+  def testRandomPoints(self):
+    polygon = geom.buffer(geom.Point(1153461.34, 649950.30), 100)
+    pts = geom.randomPoints(polygon, 10)
+    print(str(pts))
+    assert pts.getNumGeometries() == 10