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

## 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.

## 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)
```