I'd expect something like this to work ```r library(ggplot2) library(sf) library(albersusa) usa_sf <- st_as_sf(usa_composite("laea")) centroids <- st_centroid(usa_sf) ggplot() + geom_sf(data = usa_sf) + stat_sf(data = centroids, aes(geometry = geometry, label = name), geom = "text") ```