layer_us_state <- get_level_layer(gd_us, "state")
state <-
geolevel(name = "state",
layer = layer_us_state,
key = "statefp",
snake_case = TRUE) |>
complete_point_geometry()
Run the code above in your browser using DataLab