Hot answers tagged spatial-analyst
This usually works: from arcpy import env env.overwriteOutput = True env.workspace = outFolder raster_obj=FocalStatistics(raster, neighborhood, 'MAXIMUM', "DATA") raster_obj.save("outZ") del raster_obj ExtractValuesToPoints(in_point_features, outFolder+os.sep+"outZ", out_point_features)
Just use the ArcGIS raster calculator with a CON statement. CON(("h2005" - "h2015") >= 30, 1, 0) This will result in a binary raster where  represents differences of >= 30m and  no change at this threshold. And yes, set your analysis environment for extent and snap raster.
There are number of tools from GRASS and SAGA which allows for cost analysis which you can access from the Processing Toolbox:
In case someone is interesting to accomplish this task, without needing it to be in ArcGIS for Desktop. In R, there is the automap package, function autoKrige. The new_data argument, which indicates where the interpolation will happen, accepts points (besides grid and polygons). The manual says: new_data: a sp object containing the prediction ...
Only top voted, non community-wiki answers of a minimum length are eligible