0

名前付きのポリゴン形状ファイルCROPS.shpと次のキーがあります-

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

rgdalパッケージを使用して形状ファイルを という名前の変数に読み込みますmy_crops。現在、、、 、、のmy_crops@data6 つの異なるフィールドがあります。これらの各フィールドには、1 から 4 までの値があります。Y1Y2Y3Y4Y5Y6

たとえば、 の最初の数行は次の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

私は何をしたいですか？

最初のプロットは field に基づいてY1おり、最後のプロットは field に基づいていますY6。プロット (ポリゴン形状) の色は、各シェープファイル プロットの値 1 ～ 4 に基づいている必要があります。

プロットには、ページの下部に単一の凡例が付随する必要があります。

4

1 に答える 1

0

このようなことを意味しますか？

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)

ここに画像の説明を入力

于 2015-06-18T19:31:30.973 に答える