# 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_Rasters = GP.GetParameterAsText(0) #GP.AddMessage(str(Evidence_Rasters)) Output_Raster = GP.GetParameterAsText(1) if Output_Raster == "": outp_name = "" #GP.AddMessage(str(string.split(Evidence_Rasters,";"))) for inp_name in Evidence_Rasters: outp_name = outp_name + inp_name[0:3] Output_Raster = outp_name + "_Or" GP.SetParameterAsText(1,Output_Raster) GP.AddMessage("Output_Raster: "+ Output_Raster) # Check out ArcGIS Spatial Analyst extension license GP.CheckOutExtension("Spatial") # Process GP.OverwriteOutput = 1 #1.0 / (1.0 + Pow(Float(Evidence_Raster) / MidPoint, -Spread)) InExpression = "max(" + string.replace(Evidence_Rasters,';',',') + ")" #GP.AddMessage(InExpression) # Process: MapAlgebraStatement GP.SingleOutputMapAlgebra_sa(InExpression,Output_Raster,"#") except: GP.AddError("Error in Execution" + GP.GetMessages()) GP.AddMessage("Exit Fuzzy Or.")