compute_curvature.RdFits quadratic surfaces at each pixel
compute_curvature(rast, win, gisbase, central = TRUE, planar = FALSE)
| 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) |
a curvature `SpatRaster` (in 1/m)
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