我试图在除经纬度以外的任何投影中用 geom_sf() 绘制多边形。
我正在使用 manual pages for geom_sf() 中的示例
导入数据集:
nc <- sf::st_read(system.file("shape/nc.shp", package = "sf"), quiet = TRUE)
从 latlong 转换为 epsg:3857nc_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 对象,其中包含除经纬度以外的任何投影,我们在Stack Overflow上找到一个类似的问题:https://stackoverflow.com/questions/44887635/