ArcGIS Reclassify gagal di Python Toolbox

Question

Saya mencoba membangun Kotak Alat Python untuk menganalisis tanah longsor bawah laut. Antara lain saya ingin menjalankan alat aspek dan kemudian mengklasifikasikan kembali aspek tersebut ke dalam kelompok dengan interval 45 derajat. Karenanya saya melakukan yang berikut:

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

Hal yang aneh adalah bahwa kode di atas berfungsi dengan baik selama saya menjalankannya dari Pycharm atau jendela ArcGIS Python, tetapi ketika saya mencoba memasukkannya ke dalam kotak alat python (.pyt) dan menjalankannya, kesalahan berikut muncul:

ERROR 000864: Reklasifikasi: Input tidak berada dalam domain yang ditentukan.

ERROR 000893: Nilai GP atau teks kosong dimasukkan.

Saya benar-benar tidak mengerti mengapa itu tidak berfungsi di kotak alat. Fakta bahwa itu berfungsi sebaliknya memberi tahu saya bahwa tidak ada masalah mendasar dalam kode saya, itulah sebabnya saya mengharapkannya menjadi beberapa hal latar belakang yang saya lewatkan. Saya benar-benar bingung di sini.

Saya menggunakan ArcGIS 10.6.1 di Windows 10

Answer 1

Oke, sebenarnya saya baru saja memperbaikinya sekarang. Untuk menjawab pertanyaan saya, saya perlu memberikan sedikit latar belakang tambahan:

Pada dasarnya, interval aspek tidak dikodekan dengan keras, itu disediakan melalui parameter. Namun, itu diakses melalui metode .ValueAsText, yang mengembalikan nilai sebagai string. Namun, setelah itu diperlakukan sebagai bilangan bulat, yang merupakan asal kesalahan. Pasti juga mengapa itu berfungsi dengan baik di PyCharm tetapi tidak sebagai .pyt. Lihat mwe di bawah ini:

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)