Changeset 6273 for anuga_work/production/patong/asc2raster_WGS84z47N.py

Timestamp:

Feb 4, 2009, 2:15:57 PM (15 years ago)

Author:

kristy

Message:

updating

File:

• anuga_work/production/patong/asc2raster_WGS84z47N.py (modified) (5 diffs)

anuga_work/production/patong/asc2raster_WGS84z47N.py

@@ -26 +26 @@
 
 
-time_dir1 = '20081126_141801_run_final_0.8_kvanputt'
+time_dir1 = '20090130_135044_run_final_0.8_buildings_onielsen'
+time_dir2 = '20090123_135355_run_final_0.8_buildings_onielsen'
 
-time_dirs = [time_dir1]
+time_dirs = [time_dir1] #, time_dir2]
 
 for time_dir in time_dirs:
@@ -35 +36 @@
     folder = scenario_dir + output_dir + time_dir + '\\'
     raster_gbd = folder + 'raster.gdb'
-   # land = scenario_dir + "map_work\\Patong.gdb\\poly_all_cut"
+    land = scenario_dir + "map_work\\patong.gdb\\initial_conditions_Dissolve"
     
-    print 'Process: Create File GDB'
-    gp.CreateFileGDB_management(folder, "raster")
+##    print 'Process: Create File GDB'
+##    gp.CreateFileGDB_management(folder, "raster")
 
     gp.Workspace = raster_gbd
@@ -49 +50 @@
 
     generate_filename = []
-    input_ascii = glob.glob(folder + '*.asc')
+    input_ascii = glob.glob(folder + '*speed.asc')
 
     for infile in input_ascii:
@@ -61 +62 @@
         generate_filename.append(output_DEM)
         print 'Output DEM ',output_DEM
-        output_extract = output_DEM + 'E'
-        print 'Output Extract ',output_extract
         
         print 'Process: ASCII to Raster'
@@ -75 +74 @@
                                                    ",PARAMETER['Latitude_Of_Origin',0.0],UNIT['Meter',1.0]]")
 
-
-
+##        output_extract = output_DEM + 'E'
+##        print 'Output Extract ',output_extract
 ##        print 'Process: Extract by Mask'
 ##        gp.ExtractByMask_sa(output_DEM, land, output_extract)