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:
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)