Focal Stats & Map algebra

CiaranRobb — Tue, 20 Nov 2012 11:58:27 GMT

Hi everyone,

I'm a Python novice trying to write some simple scripts to streamline operations for my research. I'm trying to perform some focal
stats to then use as variables in a map algebra expression to map topographic pits and peaks.

Here is the malfunctioning script - Can anyone tell me where I'm going wrong? The script works fine for the focal stats or map algebra in isolation but not when put together

# Name: SDVI
# Description: Calaculates neighbourhood stats
# Requirements: Spatial Analyst Extension
# Author: Ciaran Robb

# Import system modules
import arcpy
from arcpy import env
from arcpy.sa import *

# Set environment settings
env.workspace = "I:\Ciaran\CAMBRIDGE PHD\Batch"

# Set local variables (inRaster = elevation)
inRaster = "fil_tif1"
neighborhood = NbrCircle(3, "CELL")

# Check out the ArcGIS Spatial Analyst extension license
arcpy.CheckOutExtension("Spatial")

# Execute FocalStatistics
outFocalStatistics1 = FocalStatistics(inRaster, neighborhood, "MEAN",
"")
outFocalStatistics2 = FocalStatistics(inRaster, neighborhood, "RANGE",
"")

# map algebra process (mean elevation - elevation))/ elevation range)
stdv3 = (Raster("focalst_fil_1") - Raster("fil_tif1")) / Raster("focalst_fil_2")

topic Focal Stats &amp; Map algebra in ArcGIS Spatial Analyst Questions

Focal Stats & Map algebra

topic Focal Stats & Map algebra in ArcGIS Spatial Analyst Questions