Changeset 6499 for anuga_work

Timestamp:

Mar 12, 2009, 1:40:40 PM (16 years ago)

Author:

kristy

Message:

reading multiple internal polygons from the one csv file

Location:

anuga_work/production/hobart_2009

Files:

• project.py (modified) (3 diffs)
• setup_model.py (modified) (2 diffs)

anuga_work/production/hobart_2009/project.py

@@ -59 +59 @@
 # Used in build_elevation.py
 # Format for ascii grids, as produced in ArcGIS + a projection file
-ascii_grid_filenames = ['hob3_project',   # Topo
-                        'grid_250m_project'] # Busselton Topo
+ascii_grid_filenames = ['hob3_pro_extract', 'grid_250m_project', 'ldr_bbay',
+                        'ldr_bellerive', 'ldr_brunyislandneck',
+                        'ldr_calvertsbeach', 'ldr_carltonbeach',
+                        'ldr_connelly', 'ldr_cremorne', 'ldr_dodge1',
+                        'ldr_dodge2', 'ldr_kingston', 'ldr_lauderdale',
+                        'ldr_linersfarne', 'ldr_oppossumbay', 'ldr_primrose',
+                        'ldr_rosny', 'ldr_sandybay', 'ldr_selfspoint',
+                        'ldr_sevenmile', 'ldr_southarm1', 'ldr_southarm2',
+                        'ldr_sullivanscove', 'ldr_tranmere', 'tasdem_mask']
 
 # Format for point is x,y,elevation (with header)
@@ -75 +82 @@
 # Format for points easting,northing (no header)
 bounding_polygon_filename = 'bounding_polygon.csv'
-bounding_polygon_maxarea = 100000
+bounding_polygon_maxarea = 1000000
 
 # INTERIOR REGIONS -  for designing the mesh
 # Used in run_model.py
 # Format for points easting,northing (no header)                    
-interior_regions_data = []
+interior_regions_data = [['aos1.csv', 1500],
+                         ['aos2.csv', 1500],
+                         ['sw.csv', 30000]]
+PriorityArea_filename = 'PriorityAreas.csv'
 
 # LAND - used to set the initial stage/water to be offcoast only
@@ -237 +247 @@
 mux_input = join(event_folder, mux_input_filename)
 
+#Multiple polygons in one CSV file to make internal polygons
+PriorityAreas = join(polygons_folder, PriorityArea_filename)
+
+

anuga_work/production/hobart_2009/setup_model.py

@@ -11 +11 @@
 from os.path import join, exists
 from anuga.utilities.polygon import read_polygon, number_mesh_triangles
-
+from anuga.shallow_water.data_manager import csv2polygons
 import project
 
@@ -157 +157 @@
 # Create list of interior polygons with scaling factor
 project.interior_regions = []
+
+polygons, maxareas = csv2polygons(project.PriorityAreas)
+print 'Creating %d internal polygons' % len(polygons)
+#def create_polygon_function(polygons, geo_reference=None):
+project.interior_regions = []
+for i, key in enumerate(polygons):
+    if i%100==0: print i
+    poly = polygons[key]
+    maxarea = float(maxareas[key])
+    project.interior_regions.append([poly,
+                                     maxarea*project.scale_factor])
+
 for filename, maxarea in project.interior_regions_data:
     polygon = read_polygon(join(project.polygons_folder, filename))