geodiv (version 1.0.5)

fitplane: Calculate a Least Squares Polynomial Surface

Description

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

Usage

fitplane(x, order)

Value

A matrix with values predicted from the polynomial fit.

Arguments

x

A raster or matrix.

order

Numeric. Indicates the polynomial order to be fit.

Examples

library(raster)

# import raster image
data(orforest)

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

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

# plot the fit
plot(x)