我在大伦敦地区有一个形状文件。我使用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/

10-10 15:31