""" 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) Minimum = GP.GetParameterAsText(1) #GP.AddMessage(Mean) Maximum = GP.GetParameterAsText(2) #GP.AddMessage(Spread) Output_Raster = GP.GetParameterAsText(3) if Output_Raster == "": Output_Raster = Evidence_Raster + "_Lin" GP.SetParameterAsText(3,Output_Raster) #GP.AddMessage(Output_Raster) # Check out ArcGIS Spatial Analyst extension license GP.CheckOutExtension("Spatial") # Process GP.OverwriteOutput = 1 #con((Float(Evidence_Raster) - Minimum ) / ( Maximum - Minimum) < 0, 0, con((Float(Evidence_raster) - Minimum ) / ( Maximum - Minimum) > 1,1,(Float(Evidence_Raster) - Minimum ) / ( Maximum - Minimum))) InExpression = "con((Float(" + Evidence_Raster + ") - " + Minimum + ") / (" + Maximum + " - " + Minimum + ") < 0, 0, con((Float(" + Evidence_Raster + ") - " + Minimum + " ) / ( " + Maximum + " - " + Minimum + ") > 1,1,(Float(" + Evidence_Raster + ") - " + Minimum + " ) / ( " + Maximum + " - " + Minimum + ")))" #GP.AddMessage(InExpression) # Process: MapAlgebraStatement GP.SingleOutputMapAlgebra_sa(InExpression,Output_Raster,"#") except: GP.AddError("Error in Execution" + GP.GetMessages()) GP.AddMessage("Exit Linear Fuzzy.")