4

我偶然发现了fastshp库，根据描述（以及我的快速粗略测试），与其他三种方法相比，它确实提供了读取大型 shapefile 的时间。

我正在使用read.shp函数从 maptools 包中加载示例数据集：

library("maptools")

setwd(system.file("shapes", package="maptools"))

shp <- read.shp("columbus.shp", format="polygon")

我选择了“多边形”格式，因为根据文档：

这通常是绘图的首选格式。

我的问题是如何使用 ggplot2 包绘制这些多边形？

4

1 回答 1

4

由于read.shp在fastshp包中以列表列表的形式返回多边形数据，因此只需将其减少为绘图所需的单个数据框ggplot2。

library(fastshp)
library(ggplot2)

setwd(system.file("shapes", package="maptools"))

shp <- read.shp("columbus.shp", format="polygon")
shp.list <- sapply(shp, FUN = function(x) do.call(cbind, x[c("id","x","y")]))
shp.df <- as.data.frame(do.call(rbind, shp.list))
shp.gg <- ggplot(shp.df, aes(x = x, y=y, group = id))+geom_polygon()

编辑：基于@otsaw 关于多边形孔的评论，以下解决方案需要更多步骤，但确保最后绘制孔。它利用 shp.df$hole 是合乎逻辑的，并且带有 hole==TRUE 的多边形将最后绘制。

shp.list <- sapply(shp, FUN = function(x) Polygon(cbind(lon = x$x, lat = x$y)))
shp.poly <- Polygons(shp.list, "area")
shp.df <- fortify(shp.poly, region = "area")
shp.gg <- ggplot(shp.df, aes(x = long, y=lat, group = piece, order = hole))+geom_polygon()
于 2012-04-25T00:27:55.473 回答