Changeset 7011 for anuga_work/production/australia_ph2/carnarvon
- Timestamp:
- May 11, 2009, 3:43:32 PM (16 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
anuga_work/production/australia_ph2/carnarvon/Arc_asc2raster_GDA94z50.py
r6949 r7011 25 25 output_dir="anuga\\outputs\\" 26 26 27 time_dir1 = '20090429_170742_run_trial_0_27285_84_Tb__kvanputt' 28 time_dir2 = '' 27 time_dir1 = '20090504_114635_run_final_0_68754_84_Tb__mhingee' 28 time_dir2 = '20090504_114715_run_final_0_70881_84_Tb__mhingee' 29 time_dir3 = '20090505_100250_run_final_0_27285_84_Tb__mhingee' 30 time_dir4 = '20090505_173023_run_final_0_27285_84_Tb_internal_mhingee' 29 31 30 time_dirs = [time_dir1 ] #, time_dir2] #4, time_dir5] #, time_dir3, time_dir4, time_dir5, time_dir6]32 time_dirs = [time_dir1, time_dir2, time_dir3, time_dir4] 31 33 32 34 for time_dir in time_dirs: … … 35 37 folder = scenario_dir + output_dir + time_dir +'\\' 36 38 raster_gbd = folder + 'raster.gdb' 37 land = scenario_dir + "map_work\\Busselton.gdb\\Internal_polygons\\initial_conditions_extend_Di" 38 ocean = scenario_dir + "map_work\\Busselton.gdb\\input_boundaries\\Ocean" 39 40 41 print 'Process: Create File GDB' 42 gp.CreateFileGDB_management(folder, "raster") 39 40 ## print 'Process: Create File GDB' 41 ## gp.CreateFileGDB_management(folder, "raster") 43 42 44 43 gp.Workspace = raster_gbd … … 73 72 ",PROJECTION['Transverse_Mercator'],PARAMETER['False_Easting',500000.0],PARAMETER['False_Northing',10000000.0],PARAMETER['Central_Meridian',113.5],PARAMETER['Scale_Factor',0.9996]" 74 73 ",PARAMETER['Latitude_Of_Origin',0.0],UNIT['Meter',1.0]]") 75 ## output_extract = output_DEM + '_E'76 ## print 'Output Extract ',output_extract77 ## print 'Process: Extract by Mask'78 ## gp.ExtractByMask_sa(output_DEM, land, output_extract)79 74 80
Note: See TracChangeset
for help on using the changeset viewer.