Topo to raster / Cliff function problem

[fsbusschers]

Spoiler

Dear members,

I'm using topo to raster to create a base of a geological surface in combination with faults. I thought this would be doable using the Cliff option in topo to raster.

But:

* If I run topo to raster using a sinlge faultline (see result_polyline2_single_fault.pdf) it works.

* If I run topo to raster using a shape with multiple faultlines (see result_polyline3_multple_fault.pdf) the faults are not taken into account (it just runs the ''normal topo to raster function).

The shapes are the same. Orientations are ok. Changing id's does not help, dissolving the multiple faults into one id also does not.

Is it possible anyway to run multiple cliffs at once?

I'm running topo to raster via pyhton 2.7 and arcmap 1.6 (see below).

Can anyone help me out?

Cheers, Freek, Geological Survey of the Netherlands  !

p.s.: input datasets are gven in the attachment (file with the single fault, file with the multiple faults and the input borehole data from which I use the field "bottom")

CODE:

#---------------------------------------------------------------------------------------------------
# ArcGIS Geoprocessing instellen
#---------------------------------------------------------------------------------------------------

import sys, string, os, time, arcpy
from arcpy.sa import *
import ReadIniFile

#---------------------------------------------------------------------------------------------------
# ArcGIS overwrite aanzetten
#---------------------------------------------------------------------------------------------------

arcpy.env.overwriteOutput = True

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

#-------------------------------------------------------------------------------------
# Ini-file inlezen
#-------------------------------------------------------------------------------------
print "Ini-file lezen ..."
IniFile = ReadIniFile.LoadConfig("verrasteren.ini")

# Gebiedsbegrenzing
xmin = int(IniFile.get("extent.xmin"))
xmax = int(IniFile.get("extent.xmax"))
ymin = int(IniFile.get("extent.ymin"))
ymax = int(IniFile.get("extent.ymax"))

#---------------------------------------------------------------------------------------------------
# Environment settings
#---------------------------------------------------------------------------------------------------
## General settings
print "Environment settings ..."
arcpy.Workspace = IniFile.get("workspace.directory")
arcpy.ScratchWorkspace = IniFile.get("scratchworkspace.directory")

# Extent goed instellen : left, bottom, right, top (!!!)
extent_string = str(xmin) + " " + str(ymin) + " " + str(xmax) + " " + str(ymax)
arcpy.env.extent = extent_string

## Raster Analysis settings
arcpy.env.cellSize = int(IniFile.get("cellsize.cellsize"))

print "-----------------------------------------------------/n"
print "Settings"
print "Workspace : " + arcpy.Workspace
print "ScratchWorkspace : " + arcpy.ScratchWorkspace
print "Cellsize : " + arcpy.env.cellSize
print "-----------------------------------------------------/n"

# Output
out_raster = "c:/Temp/be_base11"

# Input coverages
pointCov1 = "c:/Temp/Be_doorboord.shp "
attributeCov1 = "bottom "
typeCov1 = "PointElevation "

# breuken!!!!
be_breuk_utm = "c:/Temp/polyline2.shp "

outTTR = TopoToRaster([TopoPointElevation([[pointCov1, attributeCov1]]),
TopoCliff([be_breuk_utm])],
"250",
"682000 5685000 715000 5710000",
"20",
"#",
"#",
"NO_ENFORCE",
"SPOT",
"40",
"",
"1",
"0",
"#",
"#",
"#",
"#",
"#",
"#")
outTTR.save(out_raster)

[DanPatterson]

Topo to Raster (3D Analyst)—ArcGIS Pro | Documentation

out_stream_cliff_error_feature 

out_contour_error_feature

Did you try to run it with those optional parameters?

[fsbusschers]

Dan, I (now 🙂 ) did (see below). About in the middle it says: "WINDOWED CLIFF SIZE 7 EXCEEDS MAXIMUM 35". You have any idea what that means? Cheers Freek!

ANUDEM VERSION 5.3.0 (32BIT) 19/06/2011
COPYRIGHT AUSTRALIAN NATIONAL UNIVERSITY

=== TOPOtoRASTER DIAGNOSTICS ===

--- PARAMETER PROCESSING ---

DRAINAGE OPTION (0,1 OR 2): 0
0 NO DRAINAGE ENFORCEMENT
1 DRAINAGE ENFORCEMENT
2 DRAINAGE ENFORCEMENT WITH SINK DIAGNOSTICS

CONTOUR DATA OPTION (0 OR 1): 0
0 - ELEVATION DATA CONSISTS MAINLY OF SPOT HEIGHTS
1 - ELEVATION DATA CONSISTS MAINLY OF CONTOURS

NON-NEGATIVE RMS FACTOR: 1.00000
VERTICAL STANDARD ERROR: 0.00000
CLIFF ERROR: 0.20000
STREAM ERROR: 0.20000

POTENTIAL (0.0-1.0), PROFILE CURV (0.0-0.9):
0.0 - CONTOUR DATA 0.0 - DEFAULT
0.5 - POINT DATA 0.5 - WEAK PROF CURV
1.0 - NOT RECOMMENDED 0.8 - STRONG PROF CURV
0.500 0.000

ELEVATION TOLERANCE (HALF DATA CONTOUR INTERVAL):
0.000 200.000

MAXIMUM NO. OF ITERATIONS (NORMALLY 20): 40

ELEVATION LIMITS: -34.14000 36.84000

X (LONGITUDE) LIMITS: 682125.00000000 714875.00000000

Y (LATITUDE) LIMITS: 5685125.0000000 5709875.0000000

RASTER CELL SIZE: 250.00000000000
NUMBER OF ROWS IN OUTPUT RASTER = 100
NUMBER OF COLS IN OUTPUT RASTER = 132

RASTER MARGIN (NON-NEGATIVE): 5000.0000000000

NO. OF ROWS IN WORKING GRID = 140
NO. OF COLS IN WORKING GRID = 172

NUMBER OF DATA FILES: 2
TYPE 1 = POINT DATA FILE
TYPE 2 = SINK POINT FILE
TYPE 3 = STREAM LINE FILE
TYPE 4 = COAST POLYGON FILE
TYPE 5 = CONTOUR FILE
TYPE 6 = LAKE BOUNDARY FILE
TYPE 7 = CLIFF FILE
TYPE 8 = DATA MASK FILE

DATA FILE TYPES
---------------
1= POINT
2= CONTOUR
3= SINK
4= STREAM
5= COAST
6= LAKE
7= CLIFF
8= DATAMASK

INPUT FEATURE DATA 1:
NAME: C:\TEMP\BE_DOORBOORD
TYPE: POINT
ELEVATION ITEM NAME: bottom

INPUT FEATURE DATA 2:
NAME: C:\TEMP\POLYLINE3
TYPE: CLIFF

OUTPUT ELEVATION RASTER NAME: C:\TEMP\TOPOTOR_BE_DOOR1.TIF

INPUT ELEVATION POINT DATA 1: C:\TEMP\BE_DOORBOORD
NUMBER OF ELEVATION DATA POINTS = 342
NUMBER OF POINTS WITHIN LIMITS = 342
NUMBER OF POINTS ACCEPTED = 342

INPUT CLIFF LINE DATA 2: C:\TEMP\POLYLINE3

WINDOWED CLIFF SIZE 7 EXCEEDS MAXIMUM 35
NUMBER OF CONTOUR LINE POINTS = 342
NUMBER OF POINTS WITHIN LIMITS = 342
NUMBER OF POINTS ACCEPTED = 342

NUMBER OF RASTER RESOLUTIONS = 5

RESOLUTION 1
SPACING = 0.4000E+04 NCOLS = 11 NROWS = 9
RSETIW K = 1 DFAC = 6.84 EPS = 0.0000

NUMBER OF DISTINCT DATA POINTS = 50
RHO = 0.195E-01 RMS1 = 0.100E+01 RMS2 = 0.248E+01 SLOPE = 0.00 %
NUMBER OF DATA POINT SINKS = 2
NUMBER OF NON-DATA POINT SINKS = 0
TOTAL NUMBER OF SINKS = 2

RESOLUTION 2
SPACING = 0.2000E+04 NCOLS = 22 NROWS = 18
RSETIW K = 2 DFAC = 2.31 EPS = 0.0000

NUMBER OF DISTINCT DATA POINTS = 148
RHO = 0.159E-01 RMS1 = 0.100E+01 RMS2 = 0.232E+01 SLOPE = 0.00 %
NUMBER OF DATA POINT SINKS = 8
NUMBER OF NON-DATA POINT SINKS = 2
TOTAL NUMBER OF SINKS = 10

RESOLUTION 3
SPACING = 0.1000E+04 NCOLS = 43 NROWS = 35
RSETIW K = 3 DFAC = 1.28 EPS = 0.0000

NUMBER OF DISTINCT DATA POINTS = 268
RHO = 0.228E-01 RMS1 = 0.100E+01 RMS2 = 0.214E+01 SLOPE = 0.00 %
NUMBER OF DATA POINT SINKS = 29
NUMBER OF NON-DATA POINT SINKS = 4
TOTAL NUMBER OF SINKS = 33

RESOLUTION 4
SPACING = 0.5000E+03 NCOLS = 86 NROWS = 70
RSETIW K = 4 DFAC = 1.08 EPS = 0.0000

NUMBER OF DISTINCT DATA POINTS = 316
RHO = 0.520E-01 RMS1 = 0.100E+01 RMS2 = 0.138E+01 SLOPE = 0.00 %
NUMBER OF DATA POINT SINKS = 43
NUMBER OF NON-DATA POINT SINKS = 8
TOTAL NUMBER OF SINKS = 51

RESOLUTION 5
SPACING = 0.2500E+03 NCOLS = 172 NROWS = 140
RSETIW K = 5 DFAC = 1.03 EPS = 0.0000

NUMBER OF DISTINCT DATA POINTS = 331
RHO = 0.118E+00 RMS1 = 0.100E+01 RMS2 = 0.102E+01 SLOPE = 0.00 %
NUMBER OF DATA POINT SINKS = 66
NUMBER OF NON-DATA POINT SINKS = 11
TOTAL NUMBER OF SINKS = 77

NUMBER OF SPURIOUS DATA POINT SINKS = 64
NUMBER OF SPURIOUS NON-DATA POINT SINKS = 9
TOTAL NUMBER OF REMAINING SINKS = 73

-----------------------------------------
INPUT DATA SUMMARY
-----------------------------------------
ELEVATION DATA POINTS 331
SINK DATA POINTS 0
STREAM LINES 0
DISTRIBUTARIES 0
CLIFFS ON STREAMS 0
COAST POLYGONS 0
ELEVATION CONTOURS 0
MISSING ELEVATION CONTOURS 0
LAKES 0
CLIFF LINES 1
DATA MASK POLYGONS 0

-----------------------------------------
SUMMARY STATISTICS AND DATA DIAGNOSTICS
-----------------------------------------
GRID SPACING 250.000000
HORIZONTAL ERROR FACTOR 1.000
VERTICAL STANDARD ERROR 0.000
POTENTIAL ROUGHNESS PENALTY 0.500
PROFILE CURVATURE PENALTY 0.000
REMAINING SPURIOUS SINKS 73
DRAINAGE ENFORCEMENTS 0
LARGE RESIDUALS 9
FLAT STREAM SEGMENTS 0
CIRCULAR STREAM ERRORS 0
CLIFF HEIGHT ERRORS 6
COAST POLYGON DATA ERRORS 0
CONTOUR DATA ERRORS 0
LAKE POLYGON ERRORS 0
LAKE HEIGHT REVERSES 0
MULTIPLE DOWNSTREAM LAKES 0
DATA MASK POLYGON ERRORS 0
-----------------------------------------

=== END OF DIAGNOSTICS ===

[DanPatterson]

Can't find anything in the help documentation, but since you said you had 7 cliffs but there are CLIFF HEIGHT ERRORS 6 my guess would be ...

[JayantaPoddar]

Could you test if Spline with Barriers gives you the desired result?

[fsbusschers]

Dear Jayanta, spline works like a rocket..... a very exact interpolator by the way. Not sure how it deviates from topo to raster but this seems very good... still I can't compare it with topo to raster... Cheers Freek