st_coordinates()
sf
inla()
library(mapview) mapview(nc, col.regions = "navajowhite3", label = "NAME") + mapview(nc_points, col.regions = "red", cex = 5)