source: anuga_work/production/hobart_2009/Arc_Inundation_polygon.py @ 7355

Last change on this file since 7355 was 7355, checked in by kristy, 15 years ago
File size: 2.9 KB
Line 
1# ---------------------------------------------------------------------------
2# Inundation_line.py
3# Created on: Tue May 12 2009 01:50:45 PM
4#   (generated by ArcGIS/ModelBuilder)
5# Usage: Inundation_line <Bun_depth_M_E> <Bus_depth_M_E> <bunbury_1km> <busselton_1km>
6# ---------------------------------------------------------------------------
7
8# Import system modules
9import sys, string, os, arcgisscripting, glob
10
11# Create the Geoprocessor object
12gp = arcgisscripting.create()
13gp.overwriteoutput = 1
14
15# Set the necessary product code
16gp.SetProduct("ArcInfo")
17
18# Check out any necessary licenses
19gp.CheckOutExtension("spatial")
20
21# Load required toolboxes...
22gp.AddToolbox("C:/Program Files/ArcGIS/ArcToolbox/Toolboxes/Spatial Analyst Tools.tbx")
23gp.AddToolbox("C:/Program Files/ArcGIS/ArcToolbox/Toolboxes/Conversion Tools.tbx")
24gp.AddToolbox("C:/Program Files/ArcGIS/ArcToolbox/Toolboxes/Data Management Tools.tbx")
25gp.AddToolbox("C:/Program Files/ArcGIS/ArcToolbox/Toolboxes/Analysis Tools.tbx")
26
27# Script arguments...
28model = 'hobart'
29scenario_folder ="\\\\nas2\\gemd\\georisk_models\\inundation\\data\\tasmania\\" + model + "_tsunami_scenario_2009\\"
30output_folder = scenario_folder + "anuga\\outputs\\"
31
32time_dir1 = '20090714_162136_run_final_0_58260_None_kvanputt'
33time_dir2 = '20090714_162041_run_final_0.8_58292_None_kvanputt'
34time_dir3 = '20090714_161944_run_final_0.8_58280_None_kvanputt'
35time_dir4 = '20090714_161740_run_final_0.8_58260_None_kvanputt'
36
37time_dirs = [time_dir1, time_dir2, time_dir3, time_dir4]
38
39for time_dir in time_dirs:
40    print "time_dir = ",time_dir
41    raster = output_folder + time_dir + "\\raster.gdb\\"
42   
43    # Local variables...
44    locations = ['Hobart', 'NW', 'South']
45    # polygon clip
46    extent = scenario_folder + "ArcGIS\\data.gdb\\PriorityAreasHobart"
47    # name of rasters
48    for location in locations:
49        depth = raster + location +'_depth_M_E'
50        depth_raster = depth
51
52        # Local variables...
53        Output_raster = raster +"Temp_raster_"+ location
54        Output_polygon_features = raster +"Temp_polygon_"+ location
55        Output_Feature_Class = raster +"Temp_line_"+ location
56        Reclassification = "0 0.05 NODATA;0.05 4.5840001106262207 1"
57        Inundation = raster + "Inundation_" + location + "_polygon"
58        Inundation_line = raster + "Inundation_line"
59         
60        print 'Process: Reclassify'
61        gp.Reclassify_sa(depth, "Value", Reclassification, Output_raster, "DATA")
62       
63        print 'Process: Raster to Polygon'
64        gp.RasterToPolygon_conversion(Output_raster, Output_polygon_features, "SIMPLIFY", "VALUE")
65
66        print 'Process: Clip'
67        gp.Clip_analysis(Output_polygon_features, extent, Inundation, "")
68
69        print 'Process: Delete'
70        gp.Delete_management(Output_raster, "RasterDataset")
71        gp.Delete_management(Output_polygon_features, "FeatureClass")
72
Note: See TracBrowser for help on using the repository browser.