Changeset 2396

Timestamp:

Feb 14, 2006, 2:21:30 PM (19 years ago)

Author:

duncan

Message:

did a bit more work on changing the georef

Location:

inundation/pmesh

Files:

• mesh.py (modified) (1 diff)
• test_mesh.py (modified) (1 diff)

inundation/pmesh/mesh.py

@@ -801 +801 @@
         user vertex values.
 
+        Not to be used with the graphical interface
+        
         Not implemented.
         Don't implement now.  using the new georeferenced points class
-        will change this
+        will change this?
         """
         #to do
         # find the lower left hand corner
+        [xmin, ymin, xmax, ymax] = self.boxsize()
 
         # set all points to that lower left hand corner.
         # use change_geo_reference
-        pass
-
+        new_geo = Geo_reference(self.geo_reference.get_zone(), xmin, ymin)
+        self.change_geo_reference(new_geo)
+        
     def change_geo_reference(self, new_geo_reference):
         """
-        """
-        # chang egeoreference of;
-        #self.userVertices
-        #self.holes
-        #self.regions
-        pass
+        Change from one geo_reference to another.
+        Not to be used with the graphical interface
+        """
+        # FIXME
+        # change georeference of;
+        #self.userVertices = self.geo_reference.change_points_geo_ref( \
+        #self.userVertices)
+        #self.holes = self.geo_reference.change_points_geo_ref(self.holes)
+        #self.regions = self.geo_reference.change_points_geo_ref(self.regions)
+        # The above will not work.
+        # since userVertices (etc) is a list of point objects,
+        #not a list of lists.
+        # add a method to the points class to fix this up.
         
     def addVertsSegs(self, outlineDict):

inundation/pmesh/test_mesh.py

@@ -1894 +1894 @@
         m.export_mesh_file("aaat.msh")
         
+    def FIXMEtest_auto_set_geo_reference(self):
+        x=50000
+        y=1000
+        m=Mesh(geo_reference=Geo_reference(56,x,y))
+        polygon = [[0,0],[1,0],[1,1],[0,1]]
+        
+        m.add_region_from_polygon(polygon,
+                               {'tagin':[0,1],'bom':[2]},
+                                  max_triangle_area=10)
+        m.auto_set_geo_reference()
+        
        
 def list_comp(A,B):