rendered paste body# ---------------------------------------------------------------------------# model.py# Created on: Wed Apr 05 2006 06:38:35 PM# (generated by ArcGIS/ModelBuilder)# ---------------------------------------------------------------------------# Import system modulesimport sys, string, os, win32com.client# Create the Geoprocessor objectgp = win32com.client.Dispatch("esriGeoprocessing.GpDispatch.1")# Check out any necessary licensesgp.CheckOutExtension("spatial")gp.CheckOutExtension("3D")# Load required toolboxes...gp.AddToolbox("E:/Program Files/ArcGIS/ArcToolbox/Toolboxes/Conversion Tools.tbx")gp.AddToolbox("E:/Program Files/ArcGIS/ArcToolbox/Toolboxes/Spatial Analyst Tools.tbx")gp.AddToolbox("E:/Program Files/ArcGIS/ArcToolbox/Toolboxes/Analysis Tools.tbx")gp.AddToolbox("E:/Program Files/ArcGIS/ArcToolbox/Toolboxes/3D Analyst Tools.tbx")# Local variables...eskers = "eskers"dist_eskers = "M:\\Geog255\\255final\\dist_eskers"Output_direction_raster = ""lakes = "lakes"lake_raster = "M:\\Geog255\\255final\\lake_raster"dem = "dem"slope_dem = "M:\\Geog255\\255final\\slope_dem"Reclass_slop1 = "M:\\Geog255\\255final\\Reclass_slop1"Reclass_lake1 = "M:\\Geog255\\255final\\Reclass_lake1"surfgeol = "surfgeol"Feature_surf1 = "M:\\Geog255\\termasn\\feature_surf1"Reclass_Feat1 = "M:\\Geog255\\termasn\\Reclass_Feat1"eskers__2_ = "M:\\Geog255\\255final\\eskers"Output_raster__4_ = "M:\\Geog255\\255final\\Plus_Reclass3"Reclass_dist1 = "M:\\Geog255\\255final\\Reclass_dist1"Output_raster__6_ = "M:\\Geog255\\255final\\Plus_Reclass1"Output_raster__7_ = "M:\\Geog255\\255final\\Plus_Plus_Re1"economic = "M:\\Geog255\\255final\\economic"streams = "streams"wolfdens = "wolfdens"archsites = "archsites"archsites_Buffer_shp = "M:\\Geog255\\termasn\\archsites_Buffer.shp"wolfdens_Buffer_shp = "M:\\Geog255\\termasn\\wolfdens_Buffer.shp"str_buff_shp = "M:\\Geog255\\255final\\str_buff.shp"reclass_arch2 = "M:\\Geog255\\255final\\reclass_arch2"Reclass_wolf1 = "M:\\Geog255\\255final\\Reclass_wolf1"Reclass_Feat1__2_ = "M:\\Geog255\\255final\\Reclass_Feat1"arch_rast = "M:\\Geog255\\255final\\arch_rast"wolfden_rast = "M:\\Geog255\\255final\\wolfden_rast"Feature_stre1 = "M:\\Geog255\\255final\\feature_str_1"Output_raster__15_ = "M:\\Geog255\\255final\\Plus_Reclass4"enviro = "M:\\Geog255\\255final\\enviro"EEcombo = "M:\\Geog255\\255final\\EEcombo"Reclass_eske1 = "M:\\Geog255\\255final\\Reclass_eske1"costdist = "M:\\Geog255\\255final\\costdist"blacklink = "M:\\Geog255\\255final\\blacklink"lupin = "lupin"expsite = "expsite"costpath = "M:\\Geog255\\255final\\costpath"# Process: Euclidean Distance...gp.EucDistance_sa(eskers, dist_eskers, "", "100", Output_direction_raster)# Process: Buffer (3)...gp.Buffer_analysis(streams, str_buff_shp, "500 Meters", "FULL", "ROUND", "NONE", "")# Process: Feature to Raster (6)...gp.FeatureToRaster_conversion(str_buff_shp, "LENGTH", Feature_stre1, "100")# Process: Reclassify (7)...gp.Reclassify_sa(Feature_stre1, "VALUE", "15.231550216674805 500 10;500 5050.0869140625 0", Reclass_Feat1__2_, "DATA")# Process: Buffer...gp.Buffer_analysis(archsites, archsites_Buffer_shp, "500 Meters", "FULL", "ROUND", "NONE", "")# Process: Feature to Raster (4)...gp.FeatureToRaster_conversion(archsites_Buffer_shp, "CONTENTS", arch_rast, "100")# Process: Reclassify (5)...gp.Reclassify_sa(arch_rast, "Value", "1 10", reclass_arch2, "DATA")# Process: Buffer (2)...gp.Buffer_analysis(wolfdens, wolfdens_Buffer_shp, "1000 Meters", "FULL", "ROUND", "NONE", "")# Process: Feature to Raster (5)...gp.FeatureToRaster_conversion(wolfdens_Buffer_shp, "DEN", wolfden_rast, "100")# Process: Reclassify (6)...gp.Reclassify_sa(wolfden_rast, "Value", "1 10", Reclass_wolf1, "DATA")# Process: Plus (5)...gp.Plus_sa(reclass_arch2, Reclass_wolf1, Output_raster__15_)# Process: Plus (6)...gp.Plus_sa(Reclass_Feat1__2_, Output_raster__15_, enviro)# Process: Slope...gp.Slope_sa(dem, slope_dem, "DEGREE", "1")# Process: Reclassify...gp.Reclassify_sa(slope_dem, "Value", "0 2 0;2 4 3;4 6 6;6 8 8;8 19.671970367431641 10", Reclass_slop1, "DATA")# Process: Feature to Raster (3)...gp.FeatureToRaster_conversion(eskers, "CODE", eskers__2_, "100")# Process: Reclassify (8)...gp.Reclassify_sa(eskers__2_, "Value", "1 10", Reclass_eske1, "DATA")# Process: Reclassify (4)...gp.Reclassify_sa(dist_eskers, "Value", "0 1000 0;1000 2000 3;2000 24011.167025654682 6", Reclass_dist1, "DATA")# Process: Plus...gp.Plus_sa(Reclass_eske1, Reclass_dist1, Output_raster__4_)# Process: Feature to Raster...gp.FeatureToRaster_conversion(lakes, "AREA", lake_raster, "100")# Process: Reclassify (2)...gp.Reclassify_sa(lake_raster, "VALUE", "4557 978594496 1", Reclass_lake1, "DATA")# Process: Feature to Raster (2)...gp.FeatureToRaster_conversion(surfgeol, "CODE", Feature_surf1, "100")# Process: Reclassify (3)...gp.Reclassify_3d(Feature_surf1, "Value", "1 10;2 3;3 10;4 1;5 10;6 10;7 10;8 10;9 10;10 2;11 4;12 10;13 10", Reclass_Feat1, "DATA")# Process: Plus (2)...gp.Plus_sa(Reclass_lake1, Reclass_Feat1, Output_raster__6_)# Process: Plus (3)...gp.Plus_sa(Output_raster__4_, Output_raster__6_, Output_raster__7_)# Process: Plus (4)...gp.Plus_sa(Reclass_slop1, Output_raster__7_, economic)# Process: Plus (7)...gp.Plus_sa(enviro, economic, EEcombo)# Process: Cost Distance...gp.CostDistance_sa(lupin, EEcombo, costdist, "", blacklink)# Process: Cost Path...gp.CostPath_sa(expsite, costdist, blacklink, costpath, "BEST_SINGLE", "Id")