如何仅为两个类创建包含分类栅格值的箱线图？

Question

我有一个分类栅格，包含12个名为“1”，“2”，“3”等的类，代表土地用途。我有第二个栅格，其中包含一些代表蒸发蒸腾速率的值。我正在尝试为2级和10级的蒸发蒸腾率创建一个箱线图。

我能够创建箱图，但它们包含所有类，但我只想获得类2和10。

代码下方有：r，分类栅格包含12个类并且蒸发，栅格包含蒸发蒸腾值

 boxplot(evapo, r)

任何帮助都会很棒！

Answer 1

以下是Paulo解决方案的变体

library(raster)
r <- raster(nc=10, nr=5)
r[] <- runif(ncell(r), min=10, max=20) * 2
s <- setValues(r, sample(c(1:4), replace = T, size=50))

rs <- stack(r, s)
names(rs) <- c('r', 's')

d <- as.data.frame(rs)
# all classes
boxplot(r~s, data= d)
# only class 2 and 4
boxplot(r~s, data=d[d$s %in% c(2,4), ])

如果由于栅格太大而无法制作d，您可以采用大样本代替近似结果

n <- 10000
d <- data.frame(sampleRegular(rs, n))

Answer 2

library(raster)
library(ggplot2)
library(dplyr)
r <- raster(nc=10, nr=5)
r[] <- runif(ncell(r), min=10, max=20) * 2
#plot(r)
s <- setValues(r,
               sample(c(1:4), replace = T, size=50)
               )
plot(s)
ct <- crosstab(r,s, useNA=TRUE, long=TRUE)
ct2 <- ct[ rep( seq(dim(ct)[1]), ct$Freq), ]

ggplot(filter(ct2, layer.2 %in% c(1,2)), aes(y=as.numeric(layer.1), x=layer.2)) +
  geom_boxplot()

Answer 3

如果无法看到您的数据集，我无法确定此解决方案是否适合您，但这应该可以帮到您。

boxplot(evapo[which(class == 2 | 10)], r[which(class == 2 | 10)])

或者，您可以从您的蒸发和r生成“新”栅格，其中只包含类别为2和10的信息，然后运行带有这些栅格的箱形图。