fitplane {geodiv}

Calculate a Least Squares Polynomial Surface

Description

Fits a polynomial surface of order n to a raster or matrix.

Usage

fitplane(x, order)

Arguments

x

A raster or matrix.

order

Numeric. Indicates the polynomial order to be fit.

Value

A matrix with values predicted from the polynomial fit.

Examples

# import raster image
data(orforest)
orforest <- terra::unwrap(orforest)

# find the 2nd order least squares polynomial surface
polyfit <- fitplane(orforest, order = 2)

# create raster of polyfit
x <- terra::setValues(orforest, polyfit)

# plot the fit
terra::plot(x)