powered by
From an sf object containing points, calculate the standard deviational ellipse.
std_dev_ellipse(geometry)
An sf object with three columns
sx: major axis radius in CRS units,
sx
sy: minor axis radius in CRS units,
sy
theta: degree rotation of the ellipse.
theta
sf object's geometry is the center mean point.
an sfc object. If a polygon, uses sf::st_point_on_surface().
sf::st_point_on_surface()
The bulk of this function is derived from the archived CRAN package aspace version 3.2.0.
#' # 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)
Run the code above in your browser using DataLab