Mean center calculates the geographic average center. One can specify the groups to calculate individual centers for groups and weights for each individual point. It is analagous to the ArcGIS Pro Mean Center tool.
If st_is_longlat(x)
, mean center is calculated assuming a
spherical Earth. Projected data is calculated assuming a
"flat" Earth.
Examples
df <- data.frame(
lon = c(20, 50, 30, 80, 10),
lat = c(25, 70, 30, 50, 30),
grp = c("a", "b", "a", "b", "a"),
wt = c(1, 5, 1, 3, 2)
)
x <- sf::st_as_sf(df, coords = c("lon", "lat"), crs = 4326)
mean_center(x, group = "grp", weight = "wt")
#> Simple feature collection with 2 features and 1 field
#> Geometry type: POINT
#> Dimension: XY
#> Bounding box: xmin: 17.51482 ymin: 29.00224 xmax: 65.92087 ymax: 63.32603
#> Geodetic CRS: WGS 84
#> group geometry
#> 1 a POINT (17.51482 29.00224)
#> 2 b POINT (65.92087 63.32603)