Changeset 2391 for inundation/pmesh/mesh_interface.py

Timestamp:

Feb 13, 2006, 6:35:18 PM (19 years ago)

Author:

ole

Message:

Ideas

File:

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

inundation/pmesh/mesh_interface.py

@@ -27 +27 @@
 
     Boundary tags is a dictionary of symbolic tags. For every tag there
-    is a list of
-    indices referring to segments associated with that tag
+    is a list of indices referring to segments associated with that tag
 
     Resolution is the maximal area per triangle for the bounding polygon
@@ -51 +50 @@
     # check the segment indexes - throw an error if they are out of bounds
 
+    #IDEA (Ole): Let us work out what the mesh_geo_reference is automatically:
+    #If bounding_polygon is absolute it would look as follows:
+    #xllcorner = min(bounding_polygon[:,0])
+    #yllcorner = min(bounding_polygon[:,1])    
+    #
+    #mesh_geo_reference = Geo_reference(xllcorner = xllcorner,
+    #                                   yllcorner = yllcorner,
+    #                                   zone = "zone from bounding polygon")
+    #
+    #In addition I reckon the polygons could be of class Geospatial_data 
+
     m = Mesh(geo_reference=mesh_geo_reference)
 
     #Do bounding polygon
     m.add_region_from_polygon(bounding_polygon,
-                                tags=boundary_tags,
-                                geo_reference=poly_geo_reference)
+                              tags=boundary_tags,
+                              geo_reference=poly_geo_reference)
 
     #Find one point inside region automatically
@@ -62 +72 @@
         excluded_polygons = []        
         for polygon, res in interior_regions:
-            #polygon = convert_points_from_latlon_to_utm(P, refzone)            
             excluded_polygons.append( polygon )
     else:
@@ -98 +107 @@
             region = m.add_region(inner_point[0], inner_point[1])
             region.setMaxArea(res)
+            
     m.generate_mesh(minimum_triangle_angle=minimum_triangle_angle,
-                   maximum_triangle_area=maximum_triangle_area)
+                    maximum_triangle_area=maximum_triangle_area)
 
     if filename is None: