std_dev_ellipse {sfdep}

Calculation Standard Deviational Ellipse

Description

From an sf object containing points, calculate the standard deviational ellipse.

Usage

std_dev_ellipse(geometry)

Arguments

geometry

an sfc object. If a polygon, uses sf::st_point_on_surface().

Details

The bulk of this function is derived from the archived CRAN package aspace version 3.2.0.

Value

An sf object with three columns

• sx: major axis radius in CRS units,

• sy: minor axis radius in CRS units,

• theta: degree rotation of the ellipse.

sf object's geometry is the center mean point.

Examples

#' # Make a grid to sample from
grd <- sf::st_make_grid(n = c(1, 1), cellsize = c(100, 100), offset = c(0,0))

# sample 100 points
pnts <- sf::st_sample(grd, 100)
std_dev_ellipse(pnts)