10

我有一个简单的 shapefile,我想用一般的 plot() 来绘制(我注意到 ggplot 在绘制地图时非常慢)。

我可以用代码正确地绘制形状

library(maptools)    
map_shp <- readShapePoly(map_filepath)
map <- fortify(map_shp)
plot(map)

但是我如何定义线条的颜色和宽度呢?

4

1 回答 1

21

如果您不打算对所有内容都使用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)

在此处输入图像描述

于 2013-09-25T00:23:45.530 回答