Changeset 7546 for anuga_work/production/wa/mandurah_2009/build_elevation.py

Timestamp:

Oct 9, 2009, 5:17:27 PM (15 years ago)

Author:

monica

Message:

mandurah tsunami file setup - input elevation and polygons

File:

• anuga_work/production/wa/mandurah_2009/build_elevation.py (modified) (3 diffs)

anuga_work/production/wa/mandurah_2009/build_elevation.py

@@ -64 +64 @@
         geospatial_data[filename] = G_grid.clip(project.bounding_polygon)
 
-# Create Geospatial data from TXT files
+#Create Geospatial data from TXT files
 if not project.point_filenames == []:
     for filename in project.point_filenames:
@@ -79 +79 @@
 print 'Add geospatial objects' # except', project.offshore_name5
 G = None
+G1 = None
+G2 = None
 for key in geospatial_data:
-    #if key == project.point_filenames[4] or key == project.ascii_filenames[1]:
-    #    # Skip these files for now
-    #    continue
+   if key == 'm_peel_1m':
+        G1 = geospatial_data[key].clip(project.clip_box)
+
+   if key == 'ph10mex_clipped':
+        G2 = geospatial_data[key].clip_outside(project.clip_box)
+
+    G += geospatial_data[key]
+G = G + G1 + G2
     
-    G += geospatial_data[key]
-        
-#print 'Clip combined geospatial data'
-##G_clip = G.clip_outside(project.poly_aoi1)
-##G_all = G_clip + geospatial_data[project.point_filenames[4]] 
-#G_clipped = G_all.clip(project.poly_all)
-#G_clip = G.clip(project.bounding_polygon)
-
-
 print 'Export combined DEM file'
 G.export_points_file(project.combined_elevation + '.pts')
@@ -98 +96 @@
 # Use for comparision in ARC
 G.export_points_file(project.combined_elevation + '.txt')
-