从 SpatialPolygonsDataFrame- ggplot vs ggplot2 vs plot 制作一个等值线？答案

【问题标题】：making a choropleth from SpatialPolygonsDataFrame- ggplot vs ggplot2 vs plot?从 SpatialPolygonsDataFrame- ggplot vs ggplot2 vs plot 制作一个等值线？
【发布时间】：2015-06-19 15:39:28
【问题描述】：

我有一个 shapefile，特别是 SpatialPolygonsDataFrame，称为 cdtract。它包含每个区的一个变量，称为 varcount。 varcount 是“NA”或 1。当我这样做时

plot(cdtract)

...我看到打印出来的地图，但我没有看到 varcount 中标记为“1”的区域与标记为“NA”的区域有任何不同。理想情况下，我希望在 varcount 中有一系列值，并查看以不同颜色反映的值。我考虑过使用 ggplot，但根据this post here，这是非常耗费资源的，使用 plot() 而不是 ggplot() 来创建等值线可能会更好。但我真的不知道该怎么做。我正在使用来自州/国家/地区的不同 shapefile，所以不确定 choroplethr 包是要走的路。

有人可以解释如何有效地将 SpatialPolygonsDataFrame 转换为等值线吗？谢谢！

【问题讨论】：

你考虑过spplot和lattice吗？
我没有，我现在正在调查。
您能否提供有关如何使用 spplot 为变量 varcount @RomanLuštrik 创建这种等值线的任何指导？谢谢。

标签： r plot ggplot2

【解决方案1】：

可耻地取自Ari's post on gis.stackexchange.com。请参阅图片链接。还有更多例子here。

library(sp)

Srs1 = Polygons(list(Polygon(cbind(c(2,4,4,1,2),c(2,3,5,4,2)))), "s1")
Srs2 = Polygons(list(Polygon(cbind(c(5,4,2,5),c(2,3,2,2)))), "s2")

SpDF <- SpatialPolygonsDataFrame( SpatialPolygons(list(Srs1,Srs2)), 
                                  data.frame(z=1:2, row.names=c("s1","s2") ) )
spplot(SpDF, zcol="z")

【讨论】：

谢谢！做到了。我也不知道 gis.stackexchange.com。太棒了！