benchmark-geom_sf-points.R

library(sf)
library(ggplot2)

do_bench <- function(n) {
  d <- data.frame(
    x = runif(n) * 360,
    y = runif(n) * 180
  )
  
  d_sf <- st_multipoint(as.matrix(d)) %>%
    st_sfc() %>% 
    st_cast("POINT")
  
  
  result <- microbenchmark::microbenchmark(
    plot(d),
    print(ggplot(d) + geom_point(aes(x, y))),
    plot(d_sf),
    print(ggplot(d_sf) + geom_sf()),
    times = 5
  )
  
  # as data.frame
  result <- summary(result)
  
  result$n <- n
  result
}

x <- purrr::map_dfr(c(10, 50, 100, 500, 1000), do_bench)
x

ggplot(x, aes(n, median, colour = expr)) +
  geom_linerange(aes(ymin = min, ymax = max)) +
  geom_point() +
  geom_line() +
  theme_minimal() +
  theme(legend.position = "top") +
  ggtitle("Why geom_sf() is so slow...?")

Result: