Changeset 6063 for anuga_work/production/busselton/asc2raster_GDA94z50.py
- Timestamp:
- Dec 11, 2008, 4:10:29 PM (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
anuga_work/production/busselton/asc2raster_GDA94z50.py
r6020 r6063 25 25 output_dir="anuga\\outputs\\" 26 26 27 ##time_dir1 = '20081209_155610_run_final_0_27255_250m_none_lfountai' 28 ##time_dir2 = '20081209_160607_run_final_0_68693_250m_none_lfountai' 27 29 28 time_dir1 = '20081009_184721_run_final_0.6_68693_alpha0.1_kvanputt' 30 31 ##time_dir1 = '20081009_184721_run_final_0.6_68693_alpha0.1_kvanputt' 29 32 time_dir2 = '20081031_120515_run_final_0_27255_alpha0.1_kvanputt' 30 time_dir3 = '20081031_120753_run_final_0.6_27283_alpha0.1_kvanputt'31 time_dir4 = '20081031_120642_run_final_0_27283_alpha0.1_kvanputt'33 ##time_dir3 = '20081031_120753_run_final_0.6_27283_alpha0.1_kvanputt' 34 ##time_dir4 = '20081031_120642_run_final_0_27283_alpha0.1_kvanputt' 32 35 time_dir5 = '20081031_120316_run_final_0_68693_alpha0.1_kvanputt' 33 time_dir6 = '20081111_113832_run_final_0.6_27255_alpha0.1_kvanputt'36 ##time_dir6 = '20081111_113832_run_final_0.6_27255_alpha0.1_kvanputt' 34 37 35 time_dirs = [time_dir 6] #1, time_dir2, time_dir3, time_dir4, time_dir5]38 time_dirs = [time_dir2, time_dir5]#, time_dir3, time_dir4, time_dir5, time_dir6] 36 39 37 40 for time_dir in time_dirs: … … 41 44 raster_gbd = folder + 'raster.gdb' 42 45 land = scenario_dir + "map_work\\Busselton.gdb\\Internal_polygons\\initial_conditions_extend" 46 ocean = scenario_dir + "map_work\\Busselton.gdb\\input_boundaries\\Ocean" 43 47 44 # Process: Create File GDB... 45 gp.CreateFileGDB_management(folder, "raster") 48 49 ## print 'Process: Create File GDB' 50 ## gp.CreateFileGDB_management(folder, "raster") 46 51 47 52 gp.Workspace = raster_gbd 48 53 49 #print gp.Workspace54 print time_dir 50 55 51 56 #replication dictionary 52 replicate = (('busselton', ''),('_', ''),('max',' M_'),57 replicate = (('busselton', ''),('_', ''),('max','_M'), 53 58 ('Busselton', 'Bus'),('Bunbury', 'Bun'), 54 59 ('depth','_dep_'),('speed', '_spe_'), 55 ('elevation', '_ele_'), ('stage','_sta _'))60 ('elevation', '_ele_'), ('stage','_stage')) 56 61 57 62 generate_filename = [] 58 input_ascii = glob.glob(folder + '* max.asc')63 input_ascii = glob.glob(folder + '*stage_max.asc') 59 64 60 65 for infile in input_ascii: … … 69 74 70 75 print 'Output DEM ',output_DEM 71 output_extract = output_DEM + 'E'72 print 'Output Extract ',output_extract73 76 74 77 print 'Process: ASCII to Raster' … … 79 82 ",PROJECTION['Transverse_Mercator'],PARAMETER['False_Easting',500000.0],PARAMETER['False_Northing',10000000.0],PARAMETER['Central_Meridian',117.0],PARAMETER['Scale_Factor',0.9996]" 80 83 ",PARAMETER['Latitude_Of_Origin',0.0],UNIT['Meter',1.0]]") 81 82 print 'Process: Extract by Mask' 83 gp.ExtractByMask_sa(output_DEM, land, output_extract) 84 ## output_extract = output_DEM + 'E' 85 ## print 'Output Extract ',output_extract 86 ## print 'Process: Extract by Mask' 87 ## gp.ExtractByMask_sa(output_DEM, ocean, output_extract) 84 88 85 89
Note: See TracChangeset
for help on using the changeset viewer.