AnsweredAssumed Answered

Error in executing tool (Focal Statistics)

Question asked by krasovecn on Feb 24, 2018
Latest reply on Feb 25, 2018 by Dan_Patterson

Hi, I'm really new at using python and I'm trying to create a script for terrain ruggedness index but I'm using 4 shapefiles. I want to convert them to rasters and then merge them and do the TRI on this one raster. My script keeps failing because of the RuntimeError: Object: Error in executing tool. I would be really grateful if anybody helps me.

 

I got this:

Traceback (most recent call last):
File "F:\GIS_mag\py_nmg\test.py", line 60, in <module>
arcpy.gp.FocalStatistics_sa(dmr, v3x3max, "Rectangle 3 3 CELL", "MAXIMUM", "DATA")
File "c:\program files (x86)\arcgis\desktop10.5\arcpy\arcpy\geoprocessing\_base.py", line 510, in <lambda>
return lambda *args: val(*gp_fixargs(args, True))
RuntimeError: Object: Error in executing tool

 

This is my script:

# Import arcpy module
import arcpy

# Script arguments
dmr1 = arcpy.GetParameterAsText(0)
if dmr1 == '#' or not dmr1:
    dmr1 = "dmr1" # provide a default value if unspecified

dmr2 = arcpy.GetParameterAsText(1)
if dmr2 == '#' or not dmr2:
    dmr2 = "dmr2" # provide a default value if unspecified

dmr3 = arcpy.GetParameterAsText(2)
if dmr3 == '#' or not dmr3:
    dmr3 = "dmr3" # provide a default value if unspecified

dmr4 = arcpy.GetParameterAsText(3)
if dmr4 == '#' or not dmr4:
    dmr4 = "dmr4" # provide a default value if unspecified

TRI = arcpy.GetParameterAsText(4)
if TRI == '#' or not TRI:
    TRI = "F:\\GIS_mag\\py_nmg\\tri" # provide a default value if unspecified

# Local variables:
dmr1_ras = "F:\\GIS_mag\\py_nmg\\dmr1_ras"
dmr2_ras = "F:\\GIS_mag\\py_nmg\\dmr2_ras"
dmr3_ras = "F:\\GIS_mag\\py_nmg\\dmr3_ras"
dmr4_ras = "F:\\GIS_mag\\py_nmg\\dmr4_ras"
dmr = dmr1_ras, dmr2_ras, dmr3_ras, dmr4_ras
py_nmg = "F:\\GIS_mag\\py_nmg"
v3x3max = "F:\\GIS_mag\\py_nmg\\3x3max"
v3x3min = "F:\\GIS_mag\\py_nmg\\3x3min"

# Process: Point to Raster
arcpy.PointToRaster_conversion(dmr1, "Field3", dmr1_ras, "MOST_FREQUENT", "NONE", "1")

# Process: Point to Raster (2)
arcpy.PointToRaster_conversion(dmr2, "Field3", dmr2_ras, "MOST_FREQUENT", "NONE", "1")

# Process: Point to Raster (3)
arcpy.PointToRaster_conversion(dmr3, "Field3", dmr3_ras, "MOST_FREQUENT", "NONE", "1")

# Process: Point to Raster (4)
arcpy.PointToRaster_conversion(dmr4, "Field3", dmr4_ras, "MOST_FREQUENT", "NONE", "1")

# Process: Mosaic To New Raster
arcpy.MosaicToNewRaster_management("F:\\GIS_mag\\py_nmg\\dmr1_ras;F:\\GIS_mag\\py_nmg\\dmr2_ras;F:\\GIS_mag\\py_nmg\\dmr3_ras;F:\\GIS_mag\\py_nmg\\dmr4_ras", py_nmg, "dmr", "", "8_BIT_UNSIGNED", "", "1", "LAST", "FIRST")

# Process: Focal Statistics (2)
arcpy.gp.FocalStatistics_sa(dmr, v3x3max, "Rectangle 3 3 CELL", "MAXIMUM", "DATA")

# Process: Focal Statistics
arcpy.gp.FocalStatistics_sa(dmr, v3x3min, "Rectangle 3 3 CELL", "MINIMUM", "DATA")

# Process: Raster Calculator
arcpy.gp.RasterCalculator_sa("SquareRoot(Abs(Square(\"3x3max\") - Square(\"3x3min\")))", TRI)

Outcomes