compute_curvature: Compute curvature over a raster

In VincentGodard/gtbox: A Geomorphology Toolbox With Various Functions For The Analysis of Topography and Surface Processes

Compute curvature over a raster

Description

Fits quadratic surfaces at each pixel

Usage

compute_curvature(rast, win, gisbase, central = TRUE, planar = FALSE)

Arguments

rast	a 'SpatRaster' representing the raster (usually a Digital Elevation Model)
win	an odd integer (>=3) specifying the size of the computation window (in pixels)
gisbase	The directory path to GRASS binaries and libraries, containing bin and lib sub-directories among others
central	Constrain fitted surface through central window cell (default TRUE)
planar	Compute planar curvature instead (default FALSE)

Details

The surface is parametrized as : z = ax^2 + by^2 +cxy +dx + ey + f

Curvature is calculated as C = 2a + 2b

the computation relies on GRASS GIS 'r.param.scale' function

Value

a curvature 'SpatRaster' (in 1/m)

VincentGodard/gtbox documentation built on Sept. 4, 2022, 3:46 a.m.