michelstos

Script for "Topo to Raster"

Discussion created by michelstos on Nov 17, 2011
I have a problem in the script to generate the DEM to Raster using the topo to Raster.

I wrote the script using the "ArcGIS Help" website. I do not know what is wrong, but using the Debug gives "ERROR 000582: Error occurred during execution."

Does anyone know how to help me?

Thanks....
...i'm using ArcGIS 9.3...

# TopoToRaster_IGC.py
# Description: Interpolar os dados do IGC da Bacia do Cavalheiro
# Interpolate a series of point features onto a rectangular raster using TopoToRaster
# Requirements: Ponto Cotado, Curvas de Nível e Hidrografia...
# Author: ESRI, modificado por Michel Metran da Silva
# Date: February 19, 2011

# Import system modules
import sys, string, os, arcgisscripting

# Create the Geoprocessor object
gp = arcgisscripting.create()

try:
    # Output Workspace
    out_workspace = "E:/GIS_Mestrado/Bacia_Cavalheiro/DEM/IGC/Topo_Raster_Testes/Arquivos_Teste/"

    # Versão da Rodada
    versao = "01"

    # Output Features...    
    #out_raster = out_workspace + "t" + versao + "_dem.tif"
    #out_shp_stream = out_workspace + "t" + versao + "_stream.shp"
    #out_shp_sink = out_workspace + "t" + versao + "_sink.shp"
    #dia_txt = out_workspace + "t" + versao + "_diag.txt"
    #out_txt = out_workspace + "t" + versao + "_parameter.txt"

    out_raster = "E:/GIS_Mestrado/Bacia_Cavalheiro/DEM/IGC/Topo_Raster_Testes/Arquivos_Teste/dem.tif"
    out_shp_stream = "E:/GIS_Mestrado/Bacia_Cavalheiro/DEM/IGC/Topo_Raster_Testes/Arquivos_Teste/stream.shp"
    out_shp_sink = "E:/GIS_Mestrado/Bacia_Cavalheiro/DEM/IGC/Topo_Raster_Testes/Arquivos_Teste/sink.shp"
    dia_txt = "E:/GIS_Mestrado/Bacia_Cavalheiro/DEM/IGC/Topo_Raster_Testes/Arquivos_Teste/diag.TXT"
    out_txt = "E:/GIS_Mestrado/Bacia_Cavalheiro/DEM/IGC/Topo_Raster_Testes/Arquivos_Teste/parameter.TXT"

    # Input Data...
    ptos_cotados = "E:/GIS_Mestrado/Bacia_Cavalheiro/Geodata_Cavalheiro/Geodata_Cavalheiro.mdb/IGC/Ptos_Cotados "
    ptos_cotados_atributos = "Elevacao "
    ptos_cotados_typeof = "PointElevation "

    represas = "E:/GIS_Mestrado/Bacia_Cavalheiro/Geodata_Cavalheiro/Geodata_Cavalheiro.mdb/IGC/Represa "
    represas_atributos = "# "
    represas_typeof = "Lake "

    hidrografia = "E:/GIS_Mestrado/Bacia_Cavalheiro/Geodata_Cavalheiro/Geodata_Cavalheiro.mdb/IGC/Hidrografia "
    hidrografia_atributos = "# "
    hidrografia_typeof = "Stream "

    curvas_nivel = "E:/GIS_Mestrado/Bacia_Cavalheiro/Geodata_Cavalheiro/Geodata_Cavalheiro.mdb/IGC/Curvas_Nivel "
    curvas_nivel_atributos = "Elevacao "
    curvas_nivel_typeof = "Contour "

    limite = "E:/GIS_Mestrado/Bacia_Cavalheiro/Geodata_Cavalheiro/Geodata_Cavalheiro.mdb/Bacia_Cavalheiro/MB_250_Contorno "
    limite_atributos = "# "
    limite_typeof = "Boundary "

    # Check out ArcGIS Spatial Analyst extension license
    gp.CheckOutExtension("Spatial")

    # Process: TopoToRaster
    gp.TopoToRaster_sa (ptos_cotados + ptos_cotados_atributos + ptos_cotados_typeof + ";" +
        represas + represas_atributos + represas_typeof + ";" +
        hidrografia + hidrografia_atributos + hidrografia_typeof + ";" +
        curvas_nivel + curvas_nivel_atributos + curvas_nivel_typeof + ";" +
        limite + limite_atributos + limite_typeof, out_raster,
        "5",            #cellsize
        "222035.12296793022 7550213.11435197950 226272.14127189049 7555020.42812331390",
                        #x-min x-max y-min y-max
        "20",           #grid margin
        "",             #Smallest z value to be used in interpolation (optional)
        "",             #Largest z value to be used in interpolation (optional)
        "ENFORCE",      #Drainage option
        "CONTOUR",      #Contour data option
        "40",           #Maximum number of iterations (optional)
        "",             #Roughness penalty (optional)
        "1",            #Discretisation error factor (optional)
        "0",            #Vertical standard error (optional)
        "",             #Tolerance 1 (optional)
        "",             #Tolerance 2 (optional)
        out_shp_stream, #Output stream polyline features (optional)
        out_shp_sink,   #Output remaining sink point features (optional)
        dia_txt,        #Output diagnostic file (optional)
        out_txt)        #Output parameter file (optional)

except:
    # If an error occurred while running a tool, then print the messages
    print "Error in script"
    print gp.GetMessages()

Outcomes