source: anuga_work/production/new_south_wales/gosford/Arc_extract_by_mask.py @ 7252

# ---------------------------------------------------------------------------
# Arc_extract_by_mask.py
# Created on: Fri Jun 05 2009 12:09:58 
#   (generated by ArcGIS/ModelBuilder)
# ---------------------------------------------------------------------------

# Import system modules
import sys, string, os, arcgisscripting, glob, os.path

# Create the Geoprocessor object
gp = arcgisscripting.create()

# Check out any necessary licenses
gp.CheckOutExtension("spatial")

# Load required toolboxes...
gp.AddToolbox("C:/Program Files/ArcGIS/ArcToolbox/Toolboxes/Spatial Analyst Tools.tbx")
gp.overwriteoutput = 1

scenario_dir="\\\\nas2\\gemd\\georisk_models\\inundation\\data\\new_south_wales\\gosford_tsunami_scenario_2009\\"
output_dir = "anuga\\outputs\\"


variable = 'speed'
time_dir1 = '20090615_160527_run_final_1.0_51436_jgriffin'
time_dir2 = '20090615_160610_run_final_0.0_51436_jgriffin'
time_dir3 = '20090615_160900_run_final_0.0_58242_jgriffin'
time_dir4 = '20090615_160939_run_final_1.0_58242_jgriffin'
time_dir5 = '20090615_161039_run_final_0.0_58284_jgriffin'
time_dir6 = '20090615_161145_run_final_0.0_58349_jgriffin'
time_dir7 = '20090615_161318_run_final_0.0_58187_jgriffin'
time_dir8 = '20090615_161339_run_final_0.0_58113_jgriffin'
time_dir9 = '20090615_161420_run_final_0.0_58025_jgriffin'

time_dirs = [time_dir1, time_dir2, time_dir3, time_dir4, time_dir5, time_dir6, time_dir7, time_dir8, time_dir9]

for time_dir in time_dirs:

    print time_dir

    #get event number
    time_dir_split = time_dir.split('_')
    event_number = time_dir_split[5]
    print 'event_number', event_number

    # Local variables...
    folder =  scenario_dir + output_dir + time_dir + '\\'
    raster_gdb = folder + 'raster.gdb'

    gp.Workspace = raster_gdb

    print gp.Workspace
    
    #replication dictionary
    replicate = (('batemans_bay', ''),
                 ('_', ''),('Geordie', 'Geo'),('Sorrento', 'Sor'), ('max','M_'),
                 ('Fremantle', 'Fre'),('Rockingham', 'Roc'),('depth','_dep_'),
                 ('speed', '_spe_'), ('elevation', '_ele_'), ('stage','_stage'))

    generate_filename = []
    input_raster = raster_gdb + '\\All_spe_M_'
    #input_raster = glob.glob(raster_gbd + 'All_dep*')
    print 'input raster', input_raster

    output_raster = raster_gdb + '\speed_clip'
    print 'output_raster', output_raster
    
# Local variables...
#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"
#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"
#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_"
#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_"
##    aoi_land = "N:\\georisk_models\\inundation\\data\\new_south_wales\\gosford_tsunami_scenario_2009\\anuga\\polygons\\polygons.gdb\\aoi_umina_land_union"
    aoi_all = "N:\\georisk_models\\inundation\\data\\new_south_wales\\gosford_tsunami_scenario_2009\\anuga\\polygons\\polygons.gdb\\aoi_terrigal_umina_allarea_Union"

# Process: Extract by Mask...
    print 'extracting by mask'
    gp.ExtractByMask_sa(input_raster, aoi_all, output_raster)

# Process: Extract by Mask (2)...
#gp.ExtractByMask_sa(All_spe_M_, aoi_all, speed_58346_msl)