I am trying to compute a moving window maximum on a raster, such that the output is the max within a 5 mile radius. This should result in an output grid where the maximum is equal to the maximum of the input grid. I am using the FocalStatistics function in spatial analyst,
Here is the code:
import arcpy
from arcpy import env
from arcpy.sa import *
env.extent = "MAXOF"
arcpy.CheckOutExtension("Spatial")
localMax = FocalStatistics(InRaster, NbrCircle(5, "MAP"), "MAXIMUM")
Pretty straight forward. But if I check that stats: arcpy.GetRasterProperties(localMax, “MAXIMUM”) The value returned is 3.4e+38 (out to lunch), whereas:
arcpy.GetRasterProperties_management(inRaster, “MAXIMUM”)
returns 5.73
I have used this function without problems before – I am out of ideas. Anyone?
edit: inRaster is projected, and the linear unit is Miles (US)