source: anuga_work/production/bunbury_storm_surge_2009/Arc_extract_by_mask.py @ 7647

Last change on this file since 7647 was 7613, checked in by fountain, 14 years ago

Bunbury storm surge modelling

File size: 3.8 KB
Line 
1# ---------------------------------------------------------------------------
2# Arc_extract_by_mask.py
3# Created on: Fri Jun 05 2009 12:09:58
4#   (generated by ArcGIS/ModelBuilder)
5# ---------------------------------------------------------------------------
6
7# Import system modules
8import sys, string, os, arcgisscripting, glob, os.path
9
10# Create the Geoprocessor object
11gp = arcgisscripting.create()
12
13# Check out any necessary licenses
14gp.CheckOutExtension("spatial")
15
16# Load required toolboxes...
17gp.AddToolbox("C:/Program Files/ArcGIS/ArcToolbox/Toolboxes/Spatial Analyst Tools.tbx")
18gp.overwriteoutput = 1
19
20scenario_dir="\\\\nas2\\gemd\\georisk_models\\inundation\\data\\queensland\\gold_coast_tsunami_scenario_2009\\"
21output_dir = "anuga\\outputs\\"
22
23
24variable = 'speed'
25time_dir1 = '20090508_150215_run_final_0_51469_lfountai'
26time_dir2 = '20090511_165539_run_final_0_50863_lfountai'
27time_dir3 = '20090511_161526_run_final_1.1_51469_kvanputt'
28time_dir4 = '20090518_154710_run_final_0_50994_lfountai'
29time_dir5 = '20090519_160510_run_final_1.1_50994_lfountai'
30time_dir6 = '20090521_220101_run_final_1.1_50863_kvanputt'
31time_dir7 = '20090522_164526_run_final_0_51392_lfountai'
32time_dir8 = '20090522_164640_run_final_1.1_51392_lfountai'
33time_dir9 = '20090522_164948_run_final_0_51423_lfountai'
34time_dir10 = '20090522_165600_run_final_1.1_51423_lfountai'
35
36time_dirs = [time_dir1, time_dir2, time_dir3, time_dir4, time_dir5,
37             time_dir6, time_dir7, time_dir8, time_dir9, time_dir10]
38
39
40for time_dir in time_dirs:
41
42    print time_dir
43
44    #get event number
45    time_dir_split = time_dir.split('_')
46    event_number = time_dir_split[5]
47    print 'event_number', event_number
48
49    # Local variables...
50    folder =  scenario_dir + output_dir + time_dir + '\\'
51    raster_gdb = folder + 'raster.gdb'
52
53    gp.Workspace = raster_gdb
54
55    print gp.Workspace
56   
57    #replication dictionary
58    replicate = (('batemans_bay', ''),
59                 ('_', ''),('Geordie', 'Geo'),('Sorrento', 'Sor'), ('max','M_'),
60                 ('Fremantle', 'Fre'),('Rockingham', 'Roc'),('depth','_dep_'),
61                 ('speed', '_spe_'), ('elevation', '_ele_'), ('stage','_stage'))
62
63    generate_filename = []
64    input_raster = raster_gdb + '\\All_spe_M_'
65    #input_raster = glob.glob(raster_gbd + 'All_dep*')
66    print 'input raster', input_raster
67
68    output_raster = raster_gdb + '\speed_'+event_number+'_clip'
69    print 'output_raster', output_raster
70   
71# Local variables...
72#depth_58346_msl = "N:\\georisk_models\\inundation\\data\\new_south_wales\\batemans_bay_tsunami_scenario_2009\\anuga\\outputs\\20090528_153025_run_final_0.0_58346_jgriffin\\raster.gdb\\depth_58346_msl"
73#speed_58346_msl = "N:\\georisk_models\\inundation\\data\\new_south_wales\\batemans_bay_tsunami_scenario_2009\\anuga\\outputs\\20090528_153025_run_final_0.0_58346_jgriffin\\raster.gdb\\speed_58346_msl"
74#All_dep_M_ = "N:\\georisk_models\\inundation\\data\\new_south_wales\\batemans_bay_tsunami_scenario_2009\\anuga\\outputs\\20090528_153025_run_final_0.0_58346_jgriffin\\raster.gdb\\All_dep_M_"
75#All_spe_M_ = "N:\\georisk_models\\inundation\\data\\new_south_wales\\batemans_bay_tsunami_scenario_2009\\anuga\\outputs\\20090528_153025_run_final_0.0_58346_jgriffin\\raster.gdb\\All_spe_M_"
76    aoi_land = "N:\\georisk_models\\inundation\\data\\queensland\\gold_coast_tsunami_scenario_2009\\anuga\\polygons\\polygons.gdb\\areas_of_interest"
77#aoi_terrigal_umina_allarea_Union = "N:\\georisk_models\\inundation\\data\\new_south_wales\\batemans_bay_tsunami_scenario_2009\\anuga\\polygons\\polygons.gdb\\aoi_terrigal_umina_allarea_Union"
78
79# Process: Extract by Mask...
80    print 'extracting by mask'
81    gp.ExtractByMask_sa(input_raster, aoi_land, output_raster)
82
83# Process: Extract by Mask (2)...
84#gp.ExtractByMask_sa(All_spe_M_, aoi_all, speed_58346_msl)
85
Note: See TracBrowser for help on using the repository browser.