This function creates a three-band raster stack from an input digital terrain
model (DTM) of bare earth surface elevations using torch. This implementation is
faster than the terra-based implementation, especially when using GPU-based computation.
The first band is a topographic position index (TPI) calculated using a moving
window with a 50 m circular radius. The second band is the square root of slope
calculated in degrees. The third band is a TPI calculated using an annulus moving
window with an inner radius of 2 and outer radius of 5 meters. The TPI values are
clamped to a range of -10 to 10 then linearly rescaled from 0 and 1. The square
root of slope is clamped to a ange of 0 to 10 then linearly rescaled from 0 to 1.
Values are provided in floating point.
The stack is described in the following publication and was originally proposed by
William Odom of the United States Geological Survey (USGS):
Maxwell, A.E., W.E. Odom, C.M. Shobe, D.H. Doctor, M.S. Bester, and T. Ore,
2023. Exploring the influence of input feature space on CNN-based geomorphic
feature extraction from digital terrain data, Earth and Space Science,
10: e2023EA002845. https://doi.org/10.1029/2023EA002845.