don't use the sf implementation: library(tidyverse) shape <- sf::st_sf(tibble(key = 1:2, GEOM = sf::st_sfc(sf::st_point(1:2)))) shape #> Simple feature collection with 2 features and 1 field #> geo...