topic Re: Cant get ArcGIS 10.3.1 “topo to raster” to work in python script in Python Questions

Cant get ArcGIS 10.3.1 “topo to raster” to work in python script

GeoffreyRawling — Sun, 12 Dec 2021 02:53:21 GMT

I copied the sample python script from this webpage and modified for my own data and paths.

Topo to Raster—Help | ArcGIS for Desktop

# Name: TopoToRaster_Ex_02.py
# Description: Interpolates a hydrologically correct surface 
#    from point, line, and polygon data.
# Requirements: Spatial Analyst Extension

# testing topo to raster gcr

# Import system modules
import arcpy
from arcpy import env
from arcpy.sa import *

# Set environment settings
env.workspace = "D:/clovis_portales/cp_spatial/hydrology/water_level_analysis_v7_py"

# Set local variables
##inPointElevations = TopoPointElevation(['wells_in.shp', 'usgs_wl'])

inContours = TopoContour(['rwl_krg_cntr_smth.shp', 'CONTOUR'])

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

# Execute TopoToRaster
outTTR = TopoToRaster(inContours)

### Save the output 
##outTTR.save("D:/clovis_portales/cp_spatial/hydrology/water_level_analysis_v7_py/ttrout03")
‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍

I keep getting the following errors:

Traceback (most recent call last):
  File "D:\clovis_portales\cp_spatial\hydrology\water_level_analysis_v7_py\ord_krige_2010_2015_wls_v2_for_question_post.py", line 25, in <module>
    outTTR = TopoToRaster(inContours)
  File "C:\Program Files (x86)\ArcGIS\Desktop10.3\ArcPy\arcpy\sa\Functions.py", line 2777, in TopoToRaster
    out_contour_error_feature)
  File "C:\Program Files (x86)\ArcGIS\Desktop10.3\ArcPy\arcpy\sa\Utils.py", line 53, in swapper
    result = wrapper(*args, **kwargs)
  File "C:\Program Files (x86)\ArcGIS\Desktop10.3\ArcPy\arcpy\sa\Functions.py", line 2753, in Wrapper
    out_contour_error_feature)
  File "C:\Program Files (x86)\ArcGIS\Desktop10.3\ArcPy\arcpy\geoprocessing\_base.py", line 504, in <lambda>
    return lambda *args: val(*gp_fixargs(args, True))
RuntimeError: Object: Error in executing tool
‍‍‍‍‍‍‍‍‍‍‍‍

It doesnt work with point or contour data. Topo to raster works fine with same data and paths when using the GUI in arc toolbox

thanks in advance for any help

Geoff

Re: Cant get ArcGIS 10.3.1 “topo to raster” to work in python script

JoshuaBixby — Mon, 15 Aug 2016 17:55:31 GMT

https://community.esri.com/community/developers/gis-developers/python?sr=search&searchId=bace911a-a053-438f-868e-b672ad5e8969&searchIndex=0‌

Re: Cant get ArcGIS 10.3.1 “topo to raster” to work in python script

JoshuaBixby — Mon, 15 Aug 2016 17:55:53 GMT

/blogs/dan_patterson/2016/08/14/script-formatting?sr=search&searchId=0379bdd5-849e-48b8-983a-89b54c4d7003&searchIndex=0‌

Re: Cant get ArcGIS 10.3.1 “topo to raster” to work in python script

XanderBakker — Mon, 15 Aug 2016 18:35:06 GMT

Apart from the tip to use Code formatting to make your code more readable (as Joshua suggested), there are some things I see to start with:

You use Spatial Analyst functionality before checking out the extension. If this script is executed as standalone script, it will yield an error (although different from the one you have)
Your workspace is pretty long and has a lot of underscores. Some functionality of Spatial Analyst will not "like" that. You could try it with a shorter path
Most importantly, first try to do this manually (using the tools in ArcMap). The results window will provide access to Python Snippets that you may copy to examine the way the tools should be called to make this work.

See: Using the Results window—Help | ArcGIS for Desktop

Look at the option "Copy As Python Snippet"

Re: Cant get ArcGIS 10.3.1 “topo to raster” to work in python script

XanderBakker — Mon, 15 Aug 2016 18:36:15 GMT

And since this may well be data related I would be could to provide some data you are using to see if we can reproduce the error.

Re: Cant get ArcGIS 10.3.1 “topo to raster” to work in python script

XanderBakker — Mon, 15 Aug 2016 19:15:04 GMT

When you copy the Python Snippet you will get something like this:

arcpy.gp.TopoToRaster_sa("contours_v2 Elevation Contour", "C:/Users/Xander/Documents/ArcGIS/Default.gdb/DEM10m", "10", "824700 1173050 826900 1175000", "20", "", "", "ENFORCE", "CONTOUR", "20", "", "1", "0", "2,5", "100", "", "", "", "", "", "", "", "")

The first parameter is a string and consists of 3 items:

contours_v2 (the featureclass with the contourlines)
Elevation (the field with the height values
Contour (the way you want to use the lines to create the topo raster)

The second parameter is the output raster and the other parameters are optional.

In this case the tool was executed from ArcMap and the featureclass is present in my TOC. If you want to execute the tool from a standalone script, you will either need to have the featureclass available in a previously defined workspace or use a complete path to access the data.

Also note the for some strange reason the python snippet does not use arcpy.sa.TopoToRaster, but arcpy.gp.TopoToRaster_sa. Since the Interpolation tools are also available in 3D Analyst you also have "arcpy.TopoToRaster_3d" and "arcpy.ddd.TopoToRaster" which should all do the same thing.

Re: Cant get ArcGIS 10.3.1 “topo to raster” to work in python script

GeoffreyRawling — Sun, 12 Dec 2021 02:53:23 GMT

Thank you for replies. This bit of code modified from the help page on the 3d analyst works. I had tried all sorts of versions using the spatial analysts topo to raster tool, including exporting the code snippet from arctoolbox and they did not work, kept getting the same error. Gui interface in arc toolbox always worked fine.

Geoff

# V3  - modified by gcr

# Description: trying to recreate workflow of v7 analysis w py scripts
# ---------------------------------------------------------------------------

# import necessary modules etc.
import arceditor
import arcpy
from arcpy import env
from arcpy.sa import *
from arcpy.cartography import *


# set workspace
env.workspace = "D:\\clovis_portales\\cp_spatial\\hydrology\\water_level_analysis_v7_py"

# check out license
arcpy.CheckOutExtension("Spatial")
arcpy.CheckOutExtension("3d")

# set to overwrite outputs
arcpy.env.overwriteOutput = True

# Process: Topo to Raster using original data points and the smoothed kriged contours, 3d version, works

inPointElevations = "usgs_CR_wells_w_good_2010_2015_wls_matching_T_08_wells_33_rise_removed.shp usgs_wl Pointelevation"
inContours = "rwl_krg_cntr_smth.shp CONTOUR Contour"
inFeatures = (inPointElevations + ";" + inContours)
outRaster =  "wl_surf_smth_kgr_cntrs_and_points.tif"

# Execute TopoToRaster
arcpy.TopoToRaster_3d(inFeatures, outRaster)
‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍

Re: Cant get ArcGIS 10.3.1 “topo to raster” to work in python script

DanPatterson_Retired — Mon, 15 Aug 2016 20:20:46 GMT

I formatted the original code.

The last line, where you saved the raster, is commented out.... if you uncomment it, does it work? In a standalone script you need to save it... to a much shorter path, it is almost too long

Re: Cant get ArcGIS 10.3.1 “topo to raster” to work in python script

GeoffreyRawling — Mon, 15 Aug 2016 20:47:40 GMT

I uncommented the saving of the raster and moved the "check out extension" line up. Still does not work, same errors. Paths can be 128 chars for a grid but longer for other file types correct?

Re: Cant get ArcGIS 10.3.1 “topo to raster” to work in python script

DanPatterson_Retired — Mon, 15 Aug 2016 21:04:36 GMT

can't see anything obvious so it is time for print statements particularly the output filename... I know, you set the env, but stranger things have happened... try c:/temp/x.tif

The print statements will narrow down where it failed since your inputs 'look' good, but not confirmed yet.

Re: Cant get ArcGIS 10.3.1 “topo to raster” to work in python script

XanderBakker — Sun, 12 Dec 2021 02:53:26 GMT

If I look at the code of example 2 at Topo to Raster—Help | ArcGIS for Desktop , I notice that they are using "PointElevation" instead of "Pointelevation". This difference might be something that could create the error.

When using the TopoContour—Spatial Analyst module | ArcGIS for Desktop and TopoPointElevation—Spatial Analyst module | ArcGIS for Desktop classes and using the sample code from those pages. your code would be something like:

import arcpy
import os

# set workspace (still think the path is very long)
ws = "D:\\clovis_portales\\cp_spatial\\hydrology\\water_level_analysis_v7_py"
arcpy.env.workspace = ws

# check out license
arcpy.CheckOutExtension("Spatial")
arcpy.CheckOutExtension("3D")

# set to overwrite outputs
arcpy.env.overwriteOutput = True

# Topo inputs (assuming all input data is store in the ws)
inPointElevations = arcpy.sa.TopoPointElevation([["usgs_CR_wells_w_good_2010_2015_wls_matching_T_08_wells_33_rise_removed.shp","usgs_wl"]])
inContours = arcpy.sa.TopoContour([["rwl_krg_cntr_smth.shp","CONTOUR"]])

# Execute TopoToRaster
outTopoToRaster = arcpy.TopoToRaster_3d([inPointElevations, inContours])

# save the result
outRaster = "wl_surf_smth_kgr_cntrs_and_points.tif"
outTopoToRaster.save(os.path.join(ws, outTopoToRaster))‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍

I still believe you should shorten the names of your paths and input data, since this may also be a cause for errors.