使用 ggplot 从 SpatialPolygonsDataFrame 创建分区统计图

我有一个来自 UScensus2010 包的 SpatialPolygonsDataFrame 。我正在尝试创建一个分区统计图。当我这样做时，这将按预期工作：

data(colorado.county10)
choropleth(colorado.county10,
           "P0010001",
           color = list(fun = "rainbow", 
                        attr = list(4)),
           main="2010 US Counties",
           type="plot",
           border="black")

但由于未找到“P0010001”而失败

data(colorado.county10)
ggplot(colorado.county10, aes(long, lat, group = group)) +
  geom_polygon(aes(fill = P0010001), colour = alpha("white", 1/2), size = 0.2) +
  scale_fill_brewer(pal = "PuRd")

当我试图弄清楚这一点时，我注意到 colorado.county10$P0010001 返回一个数字数组，但 colorado.county10["P0010001"] 返回一个 SpatialPolygonsDataFrame。

对正在发生的事情有什么见解吗？

如果你想使用ggplot你需要强制SpatialPolygonsDataFrame to a data.frame.

ggplot2提供了一些fortify将创建正确格式的数据的方法。

目前fortify.SpatialPolygonsDataFrame方法不保留data组件，它确实提供了一个列id其中包含来自data.frame内data原来的插槽SpatialPolygonsDataFrame.

注意data.frames是存储此信息的低效方法（每个多边形的每个顶点 1 行）。

因此，以下内容可以工作，但速度很慢并且可能会导致内存问题

c10 <- fortify(colorado.county10)

c10d <- cbind(c10, colorado.county10@data[c10$id,])

ggplot(c10d, aes(long, lat, group = group)) +
   geom_polygon(aes(fill = factor(P0010001)), colour = alpha("white", 1/2), size = 0.2) +
   scale_fill_brewer(pal = "PuRd") 

Using base绘图功能会更快，并且不会占用太多资源。