ArcGIS Reclassify nie działa w Python Toolbox

Question

Próbuję zbudować Python Toolbox do analizy osuwisk podmorskich. Między innymi chcę uruchomić narzędzie aspektu, a następnie przeklasyfikować aspekty na grupy z interwałem 45 stopni. Dlatego robię co następuje:

import arcpy
import arcpy.sa
arcpy.env.workspace = r'C:\FIP_Tools20191118\GIS_files\SBP_Toolbox.gdb'

arcpy.CheckOutExtension('Spatial')
        
aspect_ras = arcpy.sa.Aspect('Smooth_Bathymetry', 'PLANAR', 'METER')
aspect_ras.save('aspect_raster')

remap = RemapRange([[0, 45, 45], [45, 90, 90], [90, 135, 135], [135, 180, 180], [180, 225, 225], [225, 270, 270], [270, 315, 315], [315, 360, 360]])
reclass_ras = arcpy.sa.Reclassify('aspect_raster', 'Value', remap, 'NODATA')
arcpy.CopyRaster_management(reclass_ras, 'reclass_Ras')

Dziwne jest to, że powyższy kod działa dobrze, o ile uruchamiam go z Pycharm lub okna ArcGIS Python, ale kiedy próbuję umieścić go w przyborniku Pythona (.pyt) i uruchomić, generuje następujące błędy:

BŁĄD 000864: Zmiana klasyfikacji: dane wejściowe nie należą do zdefiniowanej domeny.

BŁĄD 000893: Wprowadzono pustą wartość lub tekst GP.

Naprawdę nie rozumiem, dlaczego nie działa w skrzynce narzędziowej. Fakt, że działa inaczej, mówi mi, że nie ma podstawowego problemu w moim kodzie, dlatego spodziewam się, że będzie to coś w tle, którego mi brakuje. Jestem tu naprawdę zagubiony.

Używam ArcGIS 10.6.1 w systemie Windows 10

Answer 1

OK, właściwie to teraz to naprawiłem. Aby odpowiedzieć na moje pytanie, muszę podać dodatkowe informacje:

Zasadniczo interwał proporcji nie jest zakodowany na stałe, jest dostarczany za pośrednictwem parametru. Jednak uzyskano do niego dostęp za pomocą metody .ValueAsText, która zwraca wartość jako ciąg. Jednak później jest traktowany jako liczba całkowita, z której pochodzą błędy. Na pewno też, dlaczego działało dobrze w PyCharm, ale nie jako .pyt. Zobacz mwe poniżej:

import arcpy
from arcpy.sa import *
arcpy.CheckOutExtension('Spatial')

class Toolbox(object):
    def __init__(self):
        """Define the toolbox (the name of the toolbox is the name of the
        .pyt file)."""
        self.label = "ShearBandPropagation"
        self.alias = "Shear Band Propagation"

    # List of tool classes associated with this toolbox
    self.tools = [CreateGeometry]
def getParameterInfo(self):
    """Define parameter definitions"""
    param0 = arcpy.Parameter(
        displayName='Workspace',
        name='workspace',
        datatype='DEWorkspace',
        parameterType='Required',
        direction='Input')

    param1 = arcpy.Parameter(
        displayName='Smooth Bathymetry',
        name='smooth_bathymetry',
        datatype='Raster Layer',
        parameterType='Required',
        direction='Input')

    param2 = arcpy.Parameter(
        displayName='Aspect Interval',
        name='aspect_interval',
        datatype='GPLong',
        parameterType='Required',
        direction='Input')

    param2.value = 45
    param2.filter.type = "Range"
    param2.filter.list = [1, 360]

    param3 = arcpy.Parameter(
        displayName='Output',
        name='output_name',
        datatype='Raster Layer',
        parameterType='Required',
        direction='Output')


    params = [param0, param1, param2, param3]
    return params

def isLicensed(self):
    """Set whether tool is licensed to execute."""
    return True

def updateParameters(self, parameters):
    """Modify the values and properties of parameters before internal
    validation is performed.  This method is called whenever a parameter
    has been changed."""
    return

def updateMessages(self, parameters):
    """Modify the messages created by internal validation for each tool
    parameter.  This method is called after internal validation."""
    return

def execute(self, parameters, messages):
    """The source code of the tool."""

    workspace = parameters[0].valueAsText
    smooth_bathymetry = parameters[1].valueAsText
    aspect_interval = parameters[2].valueAsText
    output = parameters[3].ValueAsText

aspect_ras = arcpy.sa.Aspect(smooth_bathymetry, 'PLANAR', 'METER')
aspect_ras.save('aspect_raster')

reclassified_list = [] # create values for remap range based on input parameter
# here is where it went wrong: the correct thing to do would have been: aspect_range = int(aspect_interval)
aspect_range = aspect_interval
while aspect_range <= 360:
    if aspect_range % 1 != 0:
        reclassified_list.append([aspect_range - aspect_interval, aspect_range, int(aspect_range*10)])
        aspect_range += aspect_interval
    else:
        reclassified_list.append([aspect_range - aspect_interval, aspect_range, int(aspect_range)])
        aspect_range += aspect_interval

reclass_ras = arcpy.sa.Reclassify('aspect_raster', 'Value', RemapRange(reclassified_list), 'NODATA')
arcpy.CopyRaster_management(reclass_ras, output)