Changeset 5943

Timestamp:

Nov 12, 2008, 11:20:44 AM (15 years ago)

Author:

kristy

Message:

Updated script, works like a charm!

File:

• anuga_work/production/perth/asc2raster_GDA94z50.py (modified) (4 diffs)

anuga_work/production/perth/asc2raster_GDA94z50.py

@@ -16 +16 @@
 gp.AddToolbox("C:/Program Files/ArcGIS/ArcToolbox/Toolboxes/Conversion Tools.tbx")
 gp.AddToolbox("C:/Program Files/ArcGIS/ArcToolbox/Toolboxes/Data Management Tools.tbx")
+gp.overwriteoutput = 1
 
 output_dir = "\\\\nas2\\gemd\\georisk_models\\inundation\\data\\western_australia\\"
+scenario_dir="perth_tsunami_scenario\\anuga\\outputs\\"
 
-##Perth
-scenario_dir="perth_tsunami_scenario\\anuga\\outputs\\"
+
 time_dir1 = '20081031_133353_run_final_0.6_68693_alpha0.1_kvanputt'
 time_dir2 = '20081031_133511_run_final_0_68693_alpha0.1_kvanputt'
@@ -34 +35 @@
 
     # Local variables...
-    folder = output_dir + scenario_dir + time_dir
-    raster_gbd = output_dir + scenario_dir + time_dir + 'raster.gdb'
+    folder = output_dir + scenario_dir + time_dir + '\\'
+    raster_gbd = folder + 'raster.gdb'
     
     # Process: Create File GDB...
-    gp.CreateFileGDB_management(time_dir, "raster")
+    gp.CreateFileGDB_management(folder, "raster")
 
     gp.Workspace = raster_gbd
 
+    #print gp.Workspace
+    
     #replication dictionary
-    replicate = (('perth', ''),('Geordie', 'Geo'),('Sorrento', 'Sor'),
-                 ('Fremantle', 'Fre'),('Rockingham', 'Rock'),
-                 ('time_39900_0', 'b'), ('time_79800_0', 'c'), ('_', ''))
+    replicate = (('perth', ''),('time_39900_0', 'b'), ('time_79800_0', 'c'),
+                 ('_', ''),('Geordie', 'Geo'),('Sorrento', 'Sor'),
+                 ('Fremantle', 'Fre'),('Rockingham', 'Roc'),('depth','_depth'),
+                 ('speed', '_speed'), ('elevation', '_elev_'), ('stage','_stage'))
 
     generate_filename = []
-    input_ascii = glob.glob(folder + '*.asc')
+    input_ascii = glob.glob(folder + '*max.asc')
 
     for infile in input_ascii:
@@ -60 +64 @@
         generate_filename.append(output_DEM)
 
+        #print 'input %s and output %s' %(infile,output_DEM)
+
         # Process: ASCII to Raster (4)...
         gp.ASCIIToRaster_conversion(infile, output_DEM, "FLOAT")
@@ -69 +75 @@
 
 
+