Skip to contents

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().

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.

Details

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

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)
#> Simple feature collection with 1 feature and 3 fields
#> Geometry type: POINT
#> Dimension:     XY
#> Bounding box:  xmin: 48.34555 ymin: 48.89264 xmax: 48.34555 ymax: 48.89264
#> CRS:           NA
#>         sx       sy    theta                  geometry
#> 1 41.77136 43.83588 119.7011 POINT (48.34555 48.89264)