我有一个简单的 shapefile,我想用一般的 plot() 来绘制(我注意到 ggplot 在绘制地图时非常慢)。
我可以用代码正确地绘制形状
library(maptools)
map_shp <- readShapePoly(map_filepath)
map <- fortify(map_shp)
plot(map)
但是我如何定义线条的颜色和宽度呢?
如果您不打算对所有内容都使用ggplot2,则可以回避fortify()
并使用sp包的完善工具来绘制类对象SpatialPolygons*
:
library(rgdal) # For readOGR(), generally preferable to readShapePoly() IMHO
library(spdep) # For the "columbus.shp" file that ships with the package
map <- readOGR(dsn = system.file("etc/shapes", package="spdep"),
layer = "columbus")
plot(map, border="red", lwd=3)