uses internal extdata/error_ellipse.RDS params ordered by lc
to sample smaj, smin from log-normals & eor from a von Mises, then convert
ellipse params to covariance matrix and sample x,y errors from multivariate
normal. Returns both ellipse params and x,y errors