r - 使用 geom_sf() 绘制 sf 对象，其中包含除经纬度以外的任何投影

Question

我试图在除经纬度以外的任何投影中用 geom_sf() 绘制多边形。
我正在使用 manual pages for geom_sf() 中的示例
导入数据集:

nc <- sf::st_read(system.file("shape/nc.shp", package = "sf"), quiet = TRUE)

从 latlong 转换为 epsg:3857

nc_3857 <- sf::st_transform(nc, "+init=epsg:3857")

最后用 ggplot2 定义绘图的 crs 绘图:

ggplot() +
 geom_sf(data = nc_3857, colour = "red", fill = NA) +
 coord_sf(crs=st_crs(3857))

我一直在 wgs84(即 epsg:4326)中得到一张带有经纬度轴的 map 。我想要以米为单位的轴，所以我需要 ggplot 来绘制投影多边形。我究竟做错了什么？

Answer 1

另请参阅 https://github.com/tidyverse/ggplot2/issues/2200 并尝试

ggplot() +  geom_sf(data = nc_3857, colour = "red", fill = NA) +
   coord_sf(datum=st_crs(3857))

这使