Changeset 3177 for inundation/pmesh/mesh.py

Timestamp:

Jun 19, 2006, 10:30:57 AM (18 years ago)

Author:

duncan

Message:

added function to pmesh.mesh for adding vertices by passing in a geospatial instance.

File:

• inundation/pmesh/mesh.py (modified) (4 diffs)

inundation/pmesh/mesh.py

@@ -13 +13 @@
    Geoscience Australia
 """
-import load_mesh.loadASCII
-import alpha_shape.alpha_shape
 
 import sys
@@ -29 +27 @@
 from coordinate_transforms.geo_reference import Geo_reference,DEFAULT_ZONE
 from utilities.polygon import point_in_polygon 
+import load_mesh.loadASCII
+import alpha_shape.alpha_shape
+from geospatial_data.geospatial_data import Geospatial_data, ensure_geospatial
 
 SET_COLOUR='red'
@@ -609 +610 @@
 
     
-    
+    def add_vertices(self, data):
+        data = ensure_geospatial(data)
+
+        # get points relative to the mesh geo_ref
+        points = data.get_data_points(geo_reference=self.geo_reference)
+        #FIXME finish....
+    
+        for point in points:
+            v=Vertex(point[0], point[1])
+            self.userVertices.append(v)
+            
     def add_hole_from_polygon(self, polygon, tags=None, geo_reference=None):
         """
@@ -853 +864 @@
         
     def getUserVertices(self):
+        """
+        Note: The x,y values will be relative to the mesh geo_reference
+        This returns a list of vertex objects
+        """
         return self.userVertices
+
+    def get_user_vertices(self, absolute=True):
+        """
+        This returns a list of (x,y) values
+        (maybe it should return a geospatical object?
+        It makes it a bit confusing though.)
+        """
+        pointlist=[]
+        for vertex in self.userVertices:
+            pointlist.append([vertex.x,vertex.y])
+        spat = Geospatial_data(pointlist, geo_reference=self.geo_reference)
+        return spat.get_data_points(absolute=absolute)
     
     def getUserSegments(self):