小编Mou*_*tot的帖子

我正在为大伦敦地区提供一个形状文件.我使用包中的readShapePoly函数将maptools其作为一个加载到R中SpatialPolygonDataFrame.

我想绘制那些多边形..我已经通过使用plotR中的函数基本已经完成了.输出看起来如下图所示:

现在,我正在尝试绘制相同的形状文件,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是描述用于输入形状文件的投影的投影字符串.