This repository contains code for extracting terraces automatically from DEMs using slope and channel relief thresholds. The methodology is outlined in Clubb et al. (2017).
The user provides latitude and longitude coordinates of two points along a channel of interest. The code then extracts a swath profile along the channel between these two points following Hergarten et al. (2014). The terraces are extracted from the swath profile based on thresholds of local gradient and elevation compared to the baseline channel. These thresholds are calculated statistically from the DEM using quantile-quantile plots, and do not need to be set by the user for the landscape in question.
For more details on running the code please see the LSDTopoTools documentation