Estoy tratando de implementar una herramienta de script en ArcGIS 10 que completa algunas ecuaciones con conjuntos de datos ráster.
Parece que siempre obtengo el mismo error sin importar lo que intente:
: ERROR 000732: Entrada Ráster: El conjunto de datos 5 no existe o no es compatible
No se pudo ejecutar (TWI).
TWI es el nombre de la herramienta de script.
Mi código está debajo:
import arcpy
from arcpy import env
import os
from arcpy.sa import *
arcpy.CheckOutExtension("spatial")
DEM = arcpy.GetParameterAsText(0)
cellsize = arcpy.GetParameterAsText(1)
FlowAcc = arcpy.GetParameterAsText(2)
K_Sat = arcpy.GetParameterAsText(3)
OutRaster = arcpy.GetParameterAsText(4)
OutSpace = arcpy.GetParameterAsText(5)
slope = (Slope(Raster(DEM)) * 1.570796 ) / 90
tan_slp = Con(slope > 0, Tan(slope), 0.001)
FlowAcc_scaled = (Raster(FlowAcc) + 1 ) * cellsize
TWI = Ln(FlowAcc_scaled / (Raster(K_Sat) * tan_slp))
arcpy.env.workspace = OutSpace
TWI.save(OutRaster)
Tengo la sensación de que el error está relacionado con la configuración correcta de los parámetros. Es algo con lo que parezco tener dificultades para entender, pero según entiendo, "DEM", "FlowAcc" y "K_Sat" son capas ráster, "cellsize" es un entero largo (también probé con doble), "OutSpace" es un espacio de trabajo y "OutRaster" es una cadena.
El error parece estar en la forma en que se guarda/nombra el ráster.
EDICIÓN: Por favor, vea una captura de pantalla de la configuración de parámetros a continuación (tenga en cuenta que he probado el espacio de trabajo de salida y el ráster de salida tanto como entradas y salidas y obtengo el mismo error)
EDICIÓN 2 - Respuesta
Ok, el problema era que el segundo parámetro 'cellsize' necesita ser llamado con la función 'arcpy.GetParameter' en lugar de la función 'arcpy.GetParameterAsText' ya que es un valor de entrada en lugar de una cadena.