我在大伦敦地区有一个形状文件。我使用readShapePoly
包中的maptools
函数将其作为SpatialPolygonDataFrame
加载到R中。
我想绘制这些多边形..我已经通过使用R中的plot
函数的基础完成了这些工作。
输出如下图所示:
现在,我尝试使用ggplot2
绘制相同的形状文件,但对我不起作用。
我在图表中看到一些奇怪的线条,如下图所示:
我使用的代码是:
london.wards <- readShapePoly("~/TD/london_wards2013/london_wards2013.shp"
, proj4string=CRS(projString))
wards.count <- nrow(london.wards@data)
# assign id for each lsoa
london.wards@data$id <- 1:wards.count
wards.fort <- fortify(london.wards, region='id')
ggplot(wards.fort, aes(long, lat)) + geom_polygon(colour='black', fill='white')
其中projString是描述用于输入形状文件的投影的投影字符串。
最佳答案
您需要添加其他美观的group
。假设多边形ID为ID
,则synatx如下所示:
ggplot(wards.fort, aes(x = long, y = lat, group = ID)) +
geom_polygon(colour='black', fill='white')
关于r - 使用ggplot绘制SpatialPolygonsDataFrame,我们在Stack Overflow上找到一个类似的问题:https://stackoverflow.com/questions/18174703/