Changeset 906

Timestamp:

10/13/07 10:41:40

Author:

seang

Message:

Provide geo interface from multipoints, multilines, and multipolygons (#125)

Files:

• PCL/trunk/PCL-Core/cartography/geometry/geometry.py (modified) (3 diffs)

PCL/trunk/PCL-Core/cartography/geometry/geometry.py

@@ -500 +500 @@
-            x = createCollection(GEOS_MULTIPOINT, points)
-            self._copyFromGeom(x)
-
+
+
+
+
+
+
+
+
+
+
+
+
-
-class MultiLineString(GeometryCollection):
@@ -526 +537 @@
-            self._copyFromGeom(x)
-
+    def _provideGeoInterface(self):
+        allcoords = []
+        for i in range(self.numGeometries()):
+            line = self.geometryN(i)
+            coords = []
+            for j in range(line.numPoints()):
+                p = line.pointN(j)
+                coords.append((p.x, p.y, p.z))
+            allcoords.append(tuple(coords))
+        return {
+            'type': 'MultiLineString',
+            'coordinates': tuple(allcoords) 
+            }
+    __geo_interface__ = property(_provideGeoInterface)
+
-
-class MultiPolygon(GeometryCollection):
@@ -550 +576 @@
-            x = createCollection(GEOS_MULTIPOLYGON, polygons)
-            self._copyFromGeom(x)
+
+    def _provideGeoInterface(self):
+        allcoords = []
+        for i in range(self.numGeometries()):
+            g = self.geometryN(i)
+            geomcoords = []
+            exring = g.exteriorRing()
+            ringcoords = []
+            for k in range(exring.numPoints()):
+                p = exring.pointN(k)
+                ringcoords.append((p.x, p.y, p.z))
+            geomcoords = [tuple(ringcoords)]
+            for j in range(g.numInteriorRings()):
+                inring = g.interiorRingN(j)
+                ringcoords = []
+                for k in range(inring.numPoints()):
+                    p = inring.pointN(k)
+                    ringcoords.append((p.x, p.y, p.z))
+                geomcoords.append(tuple(ringcoords))
+            allcoords.append(tuple(geomcoords))
+        return {
+            'type': 'MultiPolygon',
+            'coordinates': tuple(allcoords) 
+            }
+    __geo_interface__ = property(_provideGeoInterface)
-
-