Changeset 6273
- Timestamp:
- Feb 4, 2009, 2:15:57 PM (16 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
anuga_work/production/patong/asc2raster_WGS84z47N.py
r6066 r6273 26 26 27 27 28 time_dir1 = '20081126_141801_run_final_0.8_kvanputt' 28 time_dir1 = '20090130_135044_run_final_0.8_buildings_onielsen' 29 time_dir2 = '20090123_135355_run_final_0.8_buildings_onielsen' 29 30 30 time_dirs = [time_dir1] 31 time_dirs = [time_dir1] #, time_dir2] 31 32 32 33 for time_dir in time_dirs: … … 35 36 folder = scenario_dir + output_dir + time_dir + '\\' 36 37 raster_gbd = folder + 'raster.gdb' 37 # land = scenario_dir + "map_work\\Patong.gdb\\poly_all_cut"38 land = scenario_dir + "map_work\\patong.gdb\\initial_conditions_Dissolve" 38 39 39 print 'Process: Create File GDB'40 gp.CreateFileGDB_management(folder, "raster")40 ## print 'Process: Create File GDB' 41 ## gp.CreateFileGDB_management(folder, "raster") 41 42 42 43 gp.Workspace = raster_gbd … … 49 50 50 51 generate_filename = [] 51 input_ascii = glob.glob(folder + '* .asc')52 input_ascii = glob.glob(folder + '*speed.asc') 52 53 53 54 for infile in input_ascii: … … 61 62 generate_filename.append(output_DEM) 62 63 print 'Output DEM ',output_DEM 63 output_extract = output_DEM + 'E'64 print 'Output Extract ',output_extract65 64 66 65 print 'Process: ASCII to Raster' … … 75 74 ",PARAMETER['Latitude_Of_Origin',0.0],UNIT['Meter',1.0]]") 76 75 77 78 76 ## output_extract = output_DEM + 'E' 77 ## print 'Output Extract ',output_extract 79 78 ## print 'Process: Extract by Mask' 80 79 ## gp.ExtractByMask_sa(output_DEM, land, output_extract)
Note: See TracChangeset
for help on using the changeset viewer.