r - 使用 ggplot 绘制 SpatialPolygonDataFrame

Question

我有大伦敦地区的形状文件。我使用包中的readShapePoly函数将maptools其作为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 是描述用于输入形状文件的投影的投影字符串。

Answer 1

您需要添加额外的美感，group. 假设多边形 id 被调用ID，则 synatx 将如下所示：

ggplot(wards.fort, aes(x = long, y = lat, group = ID)) + 
   geom_polygon(colour='black', fill='white')

Answer 2

或者，过渡到通过几何结构很好地集成的sf包是很好的。ggplot2geom_sf

library(sf)
library(ggplot2)

# Download the London shapefile.
# Links at Statistical GIS Boundary Files for London:
# https://data.london.gov.uk/dataset/statistical-gis-boundary-files-london
dataset_url <- "https://data.london.gov.uk/download/statistical-gis-boundary-files-london/b381c92b-9120-45c6-b97e-0f7adc567dd2/London-wards-2014.zip"
download.file(dataset_url, destfile = "London-wards-2014.zip")
unzip("London-wards-2014.zip", exdir = "London-wards-2014")

# Read the shapefile
polys <- st_read("./London-wards-2014/London-wards-2014 (1)/London-wards-2014_ESRI/London_Ward.shp")
#> Reading layer `London_Ward' from data source `~\London-wards-2014\London-wards-2014 (1)\London-wards-2014_ESRI\London_Ward.shp' using driver `ESRI Shapefile'
#> Simple feature collection with 654 features and 7 fields
#> geometry type:  POLYGON
#> dimension:      XY
#> bbox:           xmin: 503568.2 ymin: 155850.8 xmax: 561957.5 ymax: 200933.9
#> epsg (SRID):    NA
#> proj4string:    +proj=tmerc +lat_0=49 +lon_0=-2 +k=0.999601272 +x_0=400000 +y_0=-100000 +datum=OSGB36 +units=m +no_defs

# Fast plot/map
ggplot(polys) +
  geom_sf()

^{由reprex 包（v0.2.1）于 2019-05-20 创建}