我试图在除经纬度以外的任何投影中用 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 来绘制投影多边形。我究竟做错了什么?

最佳答案

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

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

这使

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

关于r - 使用 geom_sf() 绘制 sf 对象,其中包含除经纬度以外的任何投影,我们在Stack Overflow上找到一个类似的问题:https://stackoverflow.com/questions/44887635/

10-12 23:21