Topo to Raster by File (3D Analyst)—ArcMap

Available with Spatial Analyst license.

Available with 3D Analyst license.

Summary
Usage
Syntax
Code sample
Environments
Licensing information

Summary

Interpolates a hydrologically correct raster surface from point, line, and polygon data using parameters specified in a file.

Learn more about how Topo to Raster works

Usage

The parameter file is structured with the input datasets listed first, followed by the various parameter settings, then the output options.

The input data identifies the input datasets and, where applicable, fields. There are nine types of input: Contours, Points, Sinks, Streams, Lakes, Boundaries, Cliffs, Exclusion, and Coastal polygons. You can use as many inputs as you want, within reason. The order in which the inputs are entered does not have any bearing on the outcome. <Path> indicates a path to a dataset, <Item> indicates a field name, and <#> indicates a value to be entered.

The following table lists all of the parameters, the definition of each, and their syntax.


Parameter	Definition	Syntax
Input datasets:
Contours	Contour line dataset with item containing height values.	Contour <Path> <Item>
Points	Point dataset with item containing height values.	Point <Path> <Item>
Sinks	Point dataset containing sink locations. If the dataset has elevation values for the sinks, specify that field name as the <Item>. If only the locations of the sinks are to be used, use NONE for <Item>.	Sink <Path> <Item>
Streams	Stream line dataset. Height values are not necessary.	Stream <Path>
Lakes	Lake polygon dataset. Height values are not necessary.	Lake <Path>
Boundary	Boundary polygon dataset. Height values are not necessary.	Boundary <Path>
Cliff	Line dataset of the cliffs. There is no Field option for Cliff.	Cliff <Path>
Exclusion	Exclusion polygon dataset of the areas in which the input data should be ignored. There is no Field option for Exclusion.	Exclusion <Path>
Coast	Coast polygon dataset containing the outline of a coastal area. There is no Field option for Coast.	Boundary <Path>
Parameter settings:
Enforce	Controls whether drainage enforcement is applied.	ENFORCE <ON \| OFF \| ON_WITH_SINK>
Datatype	Primary type of input data.	DATATYPE <CONTOUR \| SPOT>
Iterations	The maximum number of iterations the algorithm performs.	ITERATIONS <#>
Roughness penalty	The measure of surface roughness.	ROUGHNESS_PENALTY <#>
Profile curvature roughness penalty	The profile curvature roughness penalty is a locally adaptive penalty that can be use to partly replace total curvature.	PROFILE_PENALTY <#>
Discretisation error factor	The amount to adjust the data smoothing of the input data into a raster.	DISCRETE_ERROR_FACTOR <#>
Vertical standard error	The amount of random error in the z-values of the input data.	VERTICAL_STANDARD_ERROR <#>
Tolerances	The first reflects the accuracy of elevation data in relation to surface drainage, and the other prevents drainage clearance through unrealistically high barriers.	TOLERANCES <#> <#>
Z-Limits	Lower and upper height limits.	ZLIMITS <#> <#>
Extent	Minimum x, minimum y, maximum x, and maximum y coordinate limits.	EXTENT <#> <#> <#> <#>
Cell size	The resolution of the final output raster.	CELL_SIZE <#>
Margin	Distance in cells to interpolate beyond the specified output extent and boundary.	MARGIN <#>
Outputs:
Output stream features	The output line feature class of stream polyline features and ridge line features.	OUT_STREAM
Output sink features	The output point feature class of the remaining sink point features.	OUT_SINK
Output diagnostics file	The location and name of the diagnostics file.	OUT_DIAGNOSTICS <Path>
Output residual point features	The output point feature class of all the large elevation residuals as scaled by the local discretisation error.	OUT_RESIDUALS
Output stream and cliff point features	The output point feature class of locations where possible stream and cliff errors occur.	OUT_STREAM_CLIFF_ERRORS
Output contour error point features	The output point feature class of possible errors pertaining to the input contour data.	OUT_CONTOUR_ERRORS

Do not specify paths for the optional output feature datasets in the parameter file. Use the Output stream polyline features and Output remaining sink point features in the tool dialog box to identify these outputs.

The contents of an example parameter file are:

     Contour D:\data\contours2\arc HEIGHT
     Point D:\data\points2\point SPOTS
     Sink D:\data\sinks_200.shp 
     Stream D:\data\streams\arc 
     Lake D:\data\lakes\polygon 
     Boundary D:\data\clipcov\polygon
     Cliff D:\data\cliffs.shp 
     ENFORCE ON
     DATATYPE CONTOUR
     ITERATIONS 40
     ROUGHNESS_PENALTY 0.0
     PROFILE_PENALTY 0.5
     DISCRETE_ERROR_FACTOR 1.0
     VERTICAL_STANDARD_ERROR 0.0
     TOLERANCES 2.5 100.0
     ZLIMITS -2000.0 13000.0
     EXTENT -810480.625 8321785.0 810480.625 10140379.0
     CELL_SIZE 1800.00000000000
     MARGIN 20
     OUT_DIAGNOSTICS D:\data\ttr_diag.txt

Syntax

arcpy.ddd.TopoToRasterByFile(in_parameter_file, out_surface_raster, {out_stream_features}, {out_sink_features}, {out_residual_feature}, {out_stream_cliff_error_feature}, {out_contour_error_feature})

Parameter	Explanation	Data Type
in_parameter_file	The input ASCII text file containing the inputs and parameters to use for the interpolation. The file is typically created from a previous run of Topo to Raster with the optional output parameter file specified. In order to test the outcome of changing the parameters, it is easier to make edits to this file and rerun the interpolation than to correctly issue the Topo to Raster tool each time.	File
out_surface_raster	The output interpolated surface raster. It is always a floating-point raster.	Raster Dataset
out_stream_features (Optional)	Output feature class of stream polyline features. The polyline features are coded as follows: 1. Input stream line not over cliff. 2. Input stream line over cliff (waterfall). 3. Drainage enforcement clearing a spurious sink. 4. Stream line determined from contour corner. 5. Ridge line determined from contour corner. 6. Code not used. 7. Data stream line side conditions. 8. Code not used. 9. Line indicating large elevation data clearance.	Feature Class
out_sink_features (Optional)	Output feature class of remaining sink point features.	Feature Class
out_residual_feature (Optional)	The output point feature class of all the large elevation residuals as scaled by the local discretisation error. All the scaled residuals larger than 10 should be inspected for possible errors in input elevation and stream data. Large-scaled residuals indicate conflicts between input elevation data and streamline data. These may also be associated with poor automatic drainage enforcements. These conflicts can be remedied by providing additional streamline and/or point elevation data after first checking and correcting errors in existing input data. Large unscaled residuals usually indicate input elevation errors.	Feature Class
out_stream_cliff_error_feature (Optional)	The output point feature class of locations where possible stream and cliff errors occur. The locations where the streams have closed loops, distributaries, and streams over cliffs can be identified from the point feature class. Cliffs with neighboring cells that are inconsistent with the high and low sides of the cliff are also indicated. This can be a good indicator of cliffs with incorrect direction. Points are coded as follows: 1. True circuit in data streamline network. 2. Circuit in stream network as encoded on the out raster. 3. Circuit in stream network via connecting lakes. 4. Distributaries point. 5. Stream over a cliff (waterfall). 6. Points indicating multiple stream outflows from lakes. 7. Code not used. 8. Points beside cliffs with heights inconsistent with cliff direction. 9. Code not used. 10. Circular distributary removed. 11. Distributary with no inflowing stream. 12. Rasterized distributary in output cell different to where the data stream line distributary occurs. 13. Error processing side conditions—an indicator of very complex streamline data.	Feature Class
out_contour_error_feature (Optional)	The output point feature class of possible errors pertaining to the input contour data. Contours with bias in height exceeding five times the standard deviation of the contour values as represented on the output raster are reported to this feature class. Contours that join other contours with a different elevation are flagged in this feature class by the code 1; this is a sure sign of a contour label error.	Feature Class

Code sample

TopoToRasterByFile example 1 (Python window)

This example creates a hydrologically correct TIFF surface raster from a parameter file defining the input point, line, and polygon data.

import arcpy
from arcpy import env  
env.workspace = "C:/data"
arcpy.TopoToRasterByFile_3d("topotorasbyfile.txt", "c:/output/ttrbf_out.tif",
                            "c:/ouput/outstreams.shp", "#", "c:/ouput/outresid.shp")

TopoToRasterByFile example 2 (stand-alone script)

This example creates a hydrologically correct Grid surface raster from a parameter file defining the input point, line, and polygon data.

# Name: TopoToRasterByFile_3d_Ex_02.py
# Description: Interpolates a hydrologically correct surface from 
#    point, line, and polygon data using parameters specified in a file.
# Requirements: 3D Analyst Extension

# Import system modules
import arcpy
from arcpy import env

# Set environment settings
env.workspace = "C:/data"

# Set local variables
inParameterFile = "topotorasterbyfile.txt"
outRaster =  "C:/output/ttrbyfout02"
outResidFC = "C:/ouput/outresid.shp"
# Check out the ArcGIS 3D Analyst extension license
arcpy.CheckOutExtension("3D")

# Execute TopoToRasterByFile
arcpy.TopoToRasterByFile_3d(inParameterFile, outRaster, "#", "#", outResidFC)

Licensing information

Basic: Requires 3D Analyst or Spatial Analyst
Standard: Requires 3D Analyst or Spatial Analyst
Advanced: Requires 3D Analyst or Spatial Analyst