Question

我有一个shapefile，特别是一个名为cdtract的SpatialPolygonsDataFrame。包含一个名为varcount的区域的变量。 varcount是＆＃34; NA＆＃34;或者1.当我做的时候

plot(cdtract)

...我看到地图已打印出来，但我看不到大片是＆＃39; 1＆＃39;在varcount中标记的任何不同于那些＆NA; NA。＆＃39;理想情况下，我希望在varcount中有一系列值，并查看以不同颜色反映的值。我考虑过使用ggplot，但根据this post here来说，这是非常耗费资源的，使用plot（）代替ggplot（）来创建等值区可能会更好。但我真的不知道如何去做。我正在使用来自州/国家的不同shapefile，所以我不确定choroplethr包是否可行。

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

Answer 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")

从SpatialPolygonsDataFrame-ggplot vs ggplot2 vs plot创建一个等值线？

1 个答案: