Changeset 7354 for anuga_work/production/hobart_2009
- Timestamp:
- Aug 12, 2009, 2:54:13 PM (16 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
anuga_work/production/hobart_2009/Arc_asc2raster_GDA94z55.py
r7268 r7354 26 26 output_dir = "anuga\\outputs\\" 27 27 28 ##time_dir1 = '20090601_133348_run_final_0.8_58280_None_kvanputt' 29 ##time_dir2 = '20090601_133503_run_final_0.8_58292_None_kvanputt' 30 ##time_dir3 = '20090505_150711_run_final_0_58280_None_kvanputt' 31 ##time_dir4 = '20090505_150517_run_final_0_58292_None_kvanputt' 32 ##time_dir5 = '20090505_151322_run_final_0.8_64477_None_kvanputt' 33 ##time_dir6 = '20090505_151447_run_final_0_64477_None_kvanputt' 34 ## 35 ##time_dirs = [time_dir4] #1, time_dir2] #, time_dir3, time_dir4, time_dir5, time_dir6] 28 time_dir1 = '20090714_162136_run_final_0_58260_None_kvanputt' 29 time_dir2 = '20090714_162041_run_final_0.8_58292_None_kvanputt' 30 time_dir3 = '20090714_161944_run_final_0.8_58280_None_kvanputt' 31 time_dir4 = '20090714_161740_run_final_0.8_58260_None_kvanputt' 36 32 37 time_dir1 = '20090619_165107_run_final_0_58260_None_kvanputt' 38 time_dir2 = '20090619_165136_run_final_0.8_58260_None_kvanputt' 33 time_dirs = [time_dir1, time_dir2, time_dir3, time_dir4] 39 34 40 time_dirs = [time_dir1, time_dir2]41 35 for time_dir in time_dirs: 42 36 # Local variables... … … 44 38 raster_gbd = folder + 'raster.gdb' 45 39 land_folder = scenario_dir + "ArcGIS\\data.gdb\\Land" 46 areas = ['NW', 'South', 'Hobart'] 40 areas = ['NW', 'South', 'Hobart'] 47 41 print time_dir 48 42 … … 55 49 #replication dictionary 56 50 replicate = (('hobart', ''), ('fitting_problem', ''), 57 ('_', ''), (' max','_M'), ('depth','_depth'),51 ('_', ''), ('Hobart','H'), ('max','_M'), ('depth','_depth'), 58 52 ('speed', '_speed'), ('elevation', '_ele'), ('stage','_stage')) 59 53 … … 70 64 sys.exit(10) 71 65 generate_filename.append(output_DEM) 66 output_DEM = raster_gbd + "\\" + output_DEM 72 67 print 'Output DEM ',output_DEM 73 68 … … 82 77 ",PARAMETER['Central_Meridian',147.0],PARAMETER['Scale_Factor',0.9996]" 83 78 ",PARAMETER['Latitude_Of_Origin',0.0],UNIT['Meter',1.0]]") 84 79 85 80 ## land = land_folder + '_' + area 86 81 ## print 'Process: Extract by Mask' … … 88 83 ## print 'Output Extract ',output_extract 89 84 ## gp.ExtractByMask_sa(output_DEM, land, output_extract) 90 ## 91 ## 85 ## print 'Process: Delete' 86 ## gp.Delete_management(output_DEM, "RasterDataset") 87
Note: See TracChangeset
for help on using the changeset viewer.