""" Spatial Data Modeller for ESRI* ArcGIS 9.2 Copyright 2007 Gary L Raines, Reno, NV, USA: production and certification Don L Sawatzky, Spokane, WA, USA: Python software development """ # Create the Geoprocessor object # Import system modules import sys, string, os #GP = win32com.client.Dispatch("esriGeoprocessing.GpDispatch.1") import arcgisscripting GP = arcgisscripting.create() try: Evidence_Raster = GP.GetParameterAsText(0) #desc = GP.Describe(Evidence_Raster) #Evidence_Raster = desc.CatalogPath #GP.AddMessage(Evidence_Raster) MidPoint = GP.GetParameterAsText(1) #GP.AddMessage(Mean) Spread = GP.GetParameterAsText(2) #GP.AddMessage(Spread) Output_Raster = GP.GetParameterAsText(3) if Output_Raster == "": Output_Raster = Evidence_Raster + "_Small" GP.SetParameterAsText(3,Output_Raster) #GP.AddMessage(Output_Raster) # Check out ArcGIS Spatial Analyst extension license GP.CheckOutExtension("Spatial") # Process GP.OverwriteOutput = 1 #1 / (1 + Pow(Float(Evidence_Raster) / MidPoint,Spread)) InExpression = "1.0 / (1.0 + Pow(Float(" + Evidence_Raster + ") / " + MidPoint + ", " + Spread + "))" #GP.AddMessage(InExpression) # Process: MapAlgebraStatement GP.SingleOutputMapAlgebra_sa(InExpression,Output_Raster,"#") except: GP.AddError("Error in Execution" + GP.GetMessages()) GP.AddMessage("Exit Small Fuzzy.")