Changeset 6561
- Timestamp:
- Mar 20, 2009, 4:03:56 PM (14 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
anuga_work/production/busselton/Arc_asc2raster_GDA94z50.py
r6201 r6561 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' 29 30 ##time_dir1 = '20081002_111432_run_final_0_27283_250m_none_kvanputt' 31 ##time_dir2 = '20081217_115336_run_final_0_27283_250m_none_dp_kvanputt' 32 33 34 time_dir1 = '20081211_154006_run_final_0.6_27255_alpha0.1_kvanputt' 27 time_dir1 = '20090130_165532_run_final_0.6_27255_extend_dt_kvanputt' 35 28 ####time_dir2 = '20081211_162311_run_final_0_27255_alpha0.1_kvanputt' 36 29 ####time_dir3 = '20081211_162346_run_final_0_68693_alpha0.1_kvanputt' … … 46 39 folder = scenario_dir + output_dir + time_dir +'\\' 47 40 raster_gbd = folder + 'raster.gdb' 48 land = scenario_dir + "map_work\\Busselton.gdb\\Internal_polygons\\initial_conditions_extend _Di"41 land = scenario_dir + "map_work\\Busselton.gdb\\Internal_polygons\\initial_conditions_extend" 49 42 ocean = scenario_dir + "map_work\\Busselton.gdb\\input_boundaries\\Ocean" 50 43 51 44 52 #print 'Process: Create File GDB'53 #gp.CreateFileGDB_management(folder, "raster")45 print 'Process: Create File GDB' 46 gp.CreateFileGDB_management(folder, "raster") 54 47 55 48 gp.Workspace = raster_gbd … … 58 51 59 52 #replication dictionary 60 replicate = (('busselton', ''),('_', ''),('max','_M _5'),53 replicate = (('busselton', ''),('_', ''),('max','_M'), 61 54 ('Busselton', 'Bus'),('Bunbury', 'Bun'), 62 55 ('depth','_depth'),('speed', '_speed'), … … 64 57 65 58 generate_filename = [] 66 input_ascii = glob.glob(folder + '* Bunbury_depth_max.asc')59 input_ascii = glob.glob(folder + '*.asc') 67 60 68 61 for infile in input_ascii: … … 85 78 ",PROJECTION['Transverse_Mercator'],PARAMETER['False_Easting',500000.0],PARAMETER['False_Northing',10000000.0],PARAMETER['Central_Meridian',117.0],PARAMETER['Scale_Factor',0.9996]" 86 79 ",PARAMETER['Latitude_Of_Origin',0.0],UNIT['Meter',1.0]]") 87 output_extract = output_DEM + '_E'88 print 'Output Extract ',output_extract89 print 'Process: Extract by Mask'90 gp.ExtractByMask_sa(output_DEM, land, output_extract)80 ## output_extract = output_DEM + '_E' 81 ## print 'Output Extract ',output_extract 82 ## print 'Process: Extract by Mask' 83 ## gp.ExtractByMask_sa(output_DEM, land, output_extract) 91 84 92 85
Note: See TracChangeset
for help on using the changeset viewer.