source: anuga_work/production/geraldton/geraldton_rerun/Inundation_polygon.py @ 7118

Last change on this file since 7118 was 7118, checked in by kristy, 15 years ago
File size: 2.8 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...
28time_dir1 = '20090324_134504_run_final_0.6_27255_Bt_kvanputt'
29time_dir2 = '20090324_140411_run_final_0.6_27283_Bt_kvanputt'
30time_dir3 = '20090324_140741_run_final_0.6_68693_Bt_kvanputt'
31time_dir4 = '20090327_155230_run_final_0_68693_Bt_kvanputt'
32time_dir5 = '20090327_160022_run_final_0_27283_Bt_kvanputt'
33time_dir6 = '20090327_160043_run_final_0_27255_Bt_kvanputt'
34
35time_dirs = [time_dir1, time_dir2, time_dir3, time_dir4, time_dir5, time_dir6]   
36
37community = 'geraldton'
38scenario = "\\\\nas2\\gemd\\georisk_models\\inundation\\data\\western_australia\\" + community +"_tsunami_scenario\\"
39output_folder = scenario + "anuga\\outputs\\"
40location = 'CBD'
41extent = scenario + "map_work\\Geraldton.gdb\Internal_Polygons\CBD_500m_clip"
42depth = 'CDB_depth_E'
43
44for time_dir in time_dirs:
45   
46    raster = output_folder + time_dir + "\\raster.gdb\\"
47
48    print 'time_dir', time_dir
49   
50    # Local variables...
51    Output_raster = raster +"Temp_raster_"+ location
52    Output_polygon_features = raster +"Temp_polygon_"+ location
53    Output_Feature_Class = raster +"Temp_line_"+ location
54    Reclassification = "0 0.05 NODATA;0.05 4.5840001106262207 1"
55    Inundation = raster + "Inundation_" + location + "_polygon"
56    depth_loc = raster + depth
57   
58    print 'Process: Reclassify'
59    gp.Reclassify_sa(depth_loc, "Value", Reclassification, Output_raster, "DATA")
60   
61    print 'Process: Raster to Polygon'
62    gp.RasterToPolygon_conversion(Output_raster, Output_polygon_features, "SIMPLIFY", "VALUE")
63
64    print 'Process: Clip'
65    gp.Clip_analysis(Output_polygon_features, extent, Inundation, "")
66
67    print 'Process: Delete'
68    gp.Delete_management(Output_raster, "RasterDataset")
69    gp.Delete_management(Output_polygon_features, "FeatureClass")
Note: See TracBrowser for help on using the repository browser.