# FireIncidentRiskAnalysis.py # Created on: 2022-03-12 12:32:56.00000 # (generated by ArcGIS/ModelBuilder) # Usage: FireIncidentRiskAnalysis # Description: # --------------------------------------------------------------------------- # Import arcpy module import arcpy from arcpy import env from arcpy.sa import * arcpy.env.overwriteOutput = True arcpy.env.workspace = arcpy.env.scratchGDB StationProximity = arcpy.GetParameter(0) if StationProximity == '#' or not StationProximity: StationProximity = 25 Landuse = arcpy.GetParameter(1) if Landuse == '#' or not Landuse: Landuse = 25 HazardousSiteProximity = arcpy.GetParameter(2) if HazardousSiteProximity == '#' or not HazardousSiteProximity: HazardousSiteProximity = 25 HydrantProximity = arcpy.GetParameter(3) if HydrantProximity == '#' or not HydrantProximity: HydrantProximity = 25 # Script arguments #Fire_Incident_Risk_Analysis = arcpy.GetParameterAsText(0) #if Fire_Incident_Risk_Analysis == '#' or not Fire_Incident_Risk_Analysis: Fire_Incident_Risk_Analysis = "C:\\CapstoneProject\\Project_Toolbox\\FireIncidentRiskAssessment\\FireIncidentRiskAssessment.gdb\\AnalysisOutput" # provide a default value if unspecified # Local variables: Proximity_to_Fire_Station_Raster = "ParcelsNearStationRaster" Landuse_Raster = "LanduseReclass" Hazardous_Material_Site_Buffer_Raster = "HazardousMaterialSites" Fire_Hydrant_Buffer_Raster = "FireHydrantDistanceBuffers" remapstation = RemapValue([[1,4],[2,3],[3,2],[4,1],["NODATA","NODATA"]]) remaplanduse = RemapValue([[1,4],[2,3],[3,2],[4,1],["NODATA","NODATA"]]) remaphazardous = RemapValue([[1,4],[2,3],[3,2],[4,1],["NODATA","NODATA"]]) remaphydrant = RemapValue([[1,5],[2,4],[3,3],[4,1],["NODATA","NODATA"]]) myWOTable = WOTable([[Proximity_to_Fire_Station_Raster, StationProximity, "VALUE", remapstation], [Landuse_Raster, Landuse, "VALUE", remaplanduse], [Hazardous_Material_Site_Buffer_Raster, HazardousSiteProximity, "VALUE", remaphazardous], [Fire_Hydrant_Buffer_Raster, HydrantProximity, "VALUE", remaphydrant]], [1, 5, 1]) # Execute WeightedOverlay WeightedOverlayOutput = WeightedOverlay(myWOTable) arcpy.CopyRaster_management(WeightedOverlayOutput, Fire_Incident_Risk_Analysis) #WeightedOverlayOutput.save(Fire_Incident_Risk_Analysis)