layer_us_county <-
dplyr::inner_join(
get_level_data_geo(gd_us, "county"),
get_level_layer(gd_us, "county"),
by = c("geoid", "statefp", "name", "type")
) |>
sf::st_as_sf()
us_county_point <-
coordinates_to_geometry(layer_us_county,
lon_lat = c("intptlon", "intptlat"))
Run the code above in your browser using DataLab