页面中的多个 spplots 和单个索引答案

【问题标题】：Multiple spplots in a page and with a single index页面中的多个 spplots 和单个索引
【发布时间】：2015-08-24 02:37:38
【问题描述】：

我有一个名为 CROPS.shp 的多边形形状文件和以下键 -

1 - Rice
2 - Wheat 
3 - Corn 
4 - Millets

我使用rgdal 包将形状文件读入名为my_crops 的变量中。现在my_crops@data 有六个不同的字段，即Y1、Y2、Y3、Y4、Y5 和Y6。这些字段中的每一个都有从 1 到 4 的值。

例如，crops@data 的前几行会是这样的

Y1   Y2   Y3   Y4   Y5    Y6 
1    2    4    1    1    2
3    4    1    1    1    2
4    2    2    2    1    3

我想做什么？

我想要一个包含 6 个图（2 行 3 列）的页面，其中第一个图将基于字段 Y1，最后一个图将基于字段 Y6。在每个 shapefile 图中，绘图（多边形形状）的颜色应基于值 1 到 4。

图应在页面底部附有一个图例。

【问题讨论】：

标签： r plot rgdal sp

【解决方案1】：

你的意思是这样的吗？

library(maptools)
nc.sids <- readShapePoly(system.file("etc/shapes/sids.shp", package="spdep")[1],
  ID="FIPSNO", proj4string=CRS("+proj=longlat +ellps=clrk66"))

set.seed(1)
nc.sids$a = factor(sample(letters[1:4], 100, replace=T))
nc.sids$b = factor(sample(letters[1:4], 100, replace=T))
nc.sids$c = factor(sample(letters[1:4], 100, replace=T))
nc.sids$d = factor(sample(letters[1:4], 100, replace=T))
nc.sids$e = factor(sample(letters[1:4], 100, replace=T))
nc.sids$f = factor(sample(letters[1:4], 100, replace=T))
spplot(nc.sids[c("a", "b", "c", "d", "e", "f")], 
    colorkey = list(space = "bottom"), layout = c(3, 2),
    as.table = TRUE)

【讨论】：