blocks|key|2374344|text|如果不能看到您的数据集，我不能确定此解决方案是否适用于您，但这应该会对您有所帮助。|type|unstyled|depth|inlineStyleRanges|entityRanges|data|2374345|boxplot(evapo[which(class+==+2+%7C+10)],+r[which(class+==+2+%7C+10)])|code-block|syntax|javascript|2374346|或者，你也可以从你的raster和r生成'new‘栅格，它只包含类2和10的信息，然后用它们运行一个盒子图。|2374347|entityMap^0|0|0|0^^$0|@$1|2|3|4|5|6|7|K|8|@]|9|@]|A|$]]|$1|B|3|C|5|D|7|L|8|@]|9|@]|A|$E|F]]|$1|G|3|H|5|6|7|M|8|@]|9|@]|A|$]]|$1|I|3|-4|5|6|7|N|8|@]|9|@]|A|$]]]|J|$]]

Without being able to see your dataset I can not be certain this solution will work for you, but this should help you out.

<pre><code>boxplot(evapo[which(class == 2 | 10)], r[which(class == 2 | 10)])
</code></pre>

Alternatively you can generates 'new' raster from your evapo and r, which only contain information where the classes are 2 and 10, then run a boxplot with those.

blocks|key|2374410|text|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+%25in%25+c(1,2)),+aes(y=as.numeric(layer.1),+x=layer.2))+%2B
++geom_boxplot()|type|code-block|depth|inlineStyleRanges|entityRanges|data|syntax|javascript|2374411|​|unstyled|2374412|📷|atomic|offset|length|2374413|2374414|entityMap|0|IMAGE|mutability|IMMUTABLE|imageUrl|https://ask.qcloudimg.com/http-save/yehe-900000/e6c9ec62f9f69c6268dedc4213ec280a.png|imageAlt^0|0|0|0|1|0|0|0^^$0|@$1|2|3|4|5|6|7|V|8|@]|9|@]|A|$B|C]]|$1|D|3|E|5|F|7|W|8|@]|9|@]|A|$]]|$1|G|3|H|5|I|7|X|8|@]|9|@$J|Y|K|Z|1|10]]|A|$]]|$1|L|3|E|5|F|7|11|8|@]|9|@]|A|$]]|$1|M|3|-4|5|F|7|12|8|@]|9|@]|A|$]]]|N|$O|$5|P|Q|R|A|$S|T|U|-4]]]]

<pre><code>library(raster)
library(ggplot2)
library(dplyr)
r &lt;- raster(nc=10, nr=5)
r[] &lt;- runif(ncell(r), min=10, max=20) * 2
#plot(r)
s &lt;- setValues(r,
 sample(c(1:4), replace = T, size=50)
 )
plot(s)
ct &lt;- crosstab(r,s, useNA=TRUE, long=TRUE)
ct2 &lt;- 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()
</code></pre>

<a href="https://i.stack.imgur.com/pNKGo.png" rel="nofollow noreferrer"><img src="https://i.stack.imgur.com/pNKGo.png" alt="enter image description here"></a>

blocks|key|2585558|text|这是保罗解的一个变种|type|unstyled|depth|inlineStyleRanges|entityRanges|data|2585559|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+%25in%25+c(2,4),+])|code-block|syntax|javascript|2585560|如果由于栅格太大而无法生成d，则可以采用较大的采样来获得近似结果|offset|length|style|CODE|2585561|n+<-+10000
d+<-+data.frame(sampleRegular(rs,+n))|2585562|entityMap^0|0|0|D|1|0|0^^$0|@$1|2|3|4|5|6|7|Q|8|@]|9|@]|A|$]]|$1|B|3|C|5|D|7|R|8|@]|9|@]|A|$E|F]]|$1|G|3|H|5|6|7|S|8|@$I|T|J|U|K|L]]|9|@]|A|$]]|$1|M|3|N|5|D|7|V|8|@]|9|@]|A|$E|F]]|$1|O|3|-4|5|6|7|W|8|@]|9|@]|A|$]]]|P|$]]

Here is a variation on Paulo's solution

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

rs &lt;- stack(r, s)
names(rs) &lt;- c('r', 's')

d &lt;- 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), ])
</code></pre>

If you cannot make <code>d</code> because the rasters are too large, you could take a large sample instead for an approximate result

<pre><code>n &lt;- 10000
d &lt;- data.frame(sampleRegular(rs, n))
</code></pre>

I have a classified raster containing 12 classes named "1", "2", "3", etc. representing the land use. I have a second raster containing some values representing the evapotranspiration rate.
I am trying to create a boxplot for the evapotranspiration rate for the classes 2 and 10. 

I was able to create boxplots but they contain all the classes but I only want to get the classes 2 and 10.

Below the code with:
r, the classified raster containing 12 classes and evapo, the raster containing the evapotranspiration values

<pre><code> boxplot(evapo, r)
</code></pre>

Any help would be great!

How to create a boxplot containing the values of a classified raster for only two classes?

翻译质量差，导致语言生硬或混乱。

没有提供实际的解决方法或示例。

解答不清晰，无法理解或解决问题。

页面排版不美观，阅读体验差。

文章

问答

视频

学习中心

腾讯云实验室

直播

竞赛

腾讯云代码分析专区

腾讯iOA零信任安全管理系统专区

腾讯云架构师技术同盟交流圈

腾讯云数据库专区

腾讯云顾问专区

腾讯云原生专区

腾讯混元专区

腾讯云TCE专区

腾讯云Lighthouse专区

腾讯云HAI专区

腾讯云Edgeone专区

腾讯云存储专区

腾讯云智能专区

腾讯轻联专区 

腾讯云开发专区

TAPD专区

腾讯轻量云游戏服专区

腾讯云最具价值专家

腾讯云架构师技术同盟

腾讯云创作之星

腾讯云开发者先锋

腾讯云代码助手

云原生构建

TAPD 敏捷项目管理

Cloud Studio

SDK中心

API中心

命令行工具

涵盖代码开发、场景应用、自动测试全流程，助你从零构建专属AI助手

一站式MCP教程库，解锁AI应用新玩法

我有一个分类栅格，其中包含12个类别，分别命名为"1“、"2”、"3“等，代表土地使用情况。我有第二个栅格，其中包含一些表示蒸散速率的值。我正在尝试为2类和10类的蒸散速率创建一个boxplot。我能够创建盒子图，但是它们包含所有的类，但我只想得到类2和类10。下面的代码包含：r (包含12个类的分类栅格)和(包含蒸散...

问如何创建一个只包含两个类的分类栅格值的箱线图？
EN

社区

活动

圈层

关于

腾讯云开发者

热门产品

热门推荐

更多推荐

问如何创建一个只包含两个类的分类栅格值的箱线图？EN

社区

活动

圈层

关于

腾讯云开发者

热门产品

热门推荐

更多推荐

问如何创建一个只包含两个类的分类栅格值的箱线图？
EN