我想使用geom_sf()函数绘制一个多边形,但不是在经纬度下的投影中。
我正在使用geom_sf()函数手册中找到的示例导入数据集。
最后使用ggplot2绘制图表,并定义绘图的crs。
我正在使用geom_sf()函数手册中找到的示例导入数据集。
nc <- sf::st_read(system.file("shape/nc.shp", package = "sf"), quiet = TRUE)
从经纬度转换为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)的地图,其纬度和经度轴。我希望将轴转换为米,因此需要ggplot
绘制投影多边形。我做错了什么?