Changeset 2376 for inundation

Timestamp:

Feb 10, 2006, 4:37:46 PM (19 years ago)

Author:

ole

Message:

Geo

File:

• inundation/pmesh/create_mesh.py (modified) (5 diffs)

inundation/pmesh/create_mesh.py

@@ -13 +13 @@
 
         if refzone is None:
-            zone = refzone
+            refzone = zone
         else:
             assert zone == refzone
@@ -68 +68 @@
                              resolution,
                              filename = None,
-                             interior_regions = None):
+                             interior_regions = None,
+                             UTM = False,
+                             refzone = None):
     """Create mesh from bounding polygon, tags for all segments and resolution.
 
@@ -99 +101 @@
     from coordinate_transforms.redfearn import redfearn
     from utilities.polygon import populate_polygon
-    
-    m = Mesh()
+    from utilities.numerical_tools import ensure_numeric
+    from coordinate_transforms.geo_reference import Geo_reference
 
 
+    bounding_polygon = ensure_numeric(bounding_polygon)
+    #print 'refzone', refzone
+
     #Convert to UTM
-    bounding_polygon, refzone = convert_points_from_latlon_to_utm(bounding_polygon)    
+    if not UTM:
+        bounding_polygon, refzone = convert_points_from_latlon_to_utm(bounding_polygon)
+    else:
+       assert refzone is not None, 'Refzone must be specified, got %s' %refzone  
+
+
+    #FIXME: Must give georeference here.
+    #Why can't Mesh work out the ll corners?
+
+    xllcorner = min(bounding_polygon[:,0])
+    yllcorner = min(bounding_polygon[:,1])    
+    
+    geo = Geo_reference(xllcorner = xllcorner,
+                        yllcorner = yllcorner,
+                        zone = refzone)
+
+    print 'geo reference derived from bounding polygon', geo
+    m = Mesh(geo_reference=geo)
+    
     
     
@@ -116 +139 @@
         excluded_polygons = []        
         for P, res in interior_regions:
-            polygon, _ = convert_points_from_latlon_to_utm(P, refzone)            
+            if not UTM:            
+                polygon, _ = convert_points_from_latlon_to_utm(P, refzone)
+            else:
+                polygon = P
             excluded_polygons.append( polygon )
     else:
@@ -133 +159 @@
     if interior_regions is not None:    
         for P, res in interior_regions:
-            polygon, _ = convert_points_from_latlon_to_utm(P, refzone)                        
+            if not UTM:                        
+                polygon, _ = convert_points_from_latlon_to_utm(P, refzone)
+            else:
+                polygon = P
+                
             region_dict = create_region(polygon, None, refzone)