Arc Extract Breakline from DEM

Introduction

This page outlines the method used to take a line feature class, and extract elevations from DEM for input to TUFLOW breaklines. This is useful if you have a hydraulic control, for example a road embankment that you want to ensure is represented in the model. In this method the elevations will be extracted from the DEM. If detailed survey exists this should be used. This can also be used to process a thalweg / creek centre line.
TUFLOW breaklines (2d Z Shape), are a mixture of line and point feature classes. This allows points to be deleted, without deleting the vertices along the line. The process of converting a line feature class, to a point feature class and extracting the elevation. Is a multiple step process and we are looking to implement a python script to automate this. There are two ways of doing this:
Method 1

1. Convert line vertices to points
2. Extract elevations directly from DEM

This method requires that the line has been digitised at the high or low point that you are trying to represent. If the vertex on the line, is not at the crest of the ridge, the elevation may be slightly out.
Method 2 Method 2 creates a buffer around each points and can be used to extract the maximum, minimum or average elevation from the DEM within this buffer. This ensures that the control is represented, even if the line does not follow the levee crest or centreline perfectly. The steps involved are:

1. Convert line vertices to points
2. Create a buffer around each point (specified buffer distance)
3. Extract desired elevation statistic (max, min, mean) for each region
4. Assign elevation from buffer to point

The steps involved in each method are detailed below.

Method 1 (point inspection)

1. The first step is to create a point feature class from the line, to do this follow the method outlined in the convert lines to points tip.
2. To create a buffer (region) object around each of the points, firstly open a region feature class. Copy the 2d_zsh_empty_R template to a working directory and rename).
Make the layer