假设我需要计算同一经度的栅格数据的方差，使用该sapp函数是正确的方法吗？但我失败了。
我现在可以使用更复杂的方法来转换为 data.frame 并分组来计算方差：
library(terra)
library(tidyverse)
r <- rast(nrows = 5, ncols = 10, vals = c(1:50))
plot(r)
r |>
terra::as.data.frame(xy= TRUE) |>
group_by(x) |>
summarise(lonsd = sd(lyr.1))
有没有更直接的函数terra？如果你知道的话请告诉我。
还有更复杂的场景，需要同时对多层栅格数据进行此操作，该怎么做？
当使用 as.data.frame 处理多个图层时，对于合并的图层，所有的命名都是lyr.1。有没有更好的方法来解决这个命名问题？
library(terra)
m <- matrix(1:25, nrow=5, ncol=5)
rm <- rast(m)
n <- matrix(rep(5,time=25),nrow = 5,ncol = 5)
rn <- rast(n)
dt <- c(rm,rn)
data <- terra::as.data.frame(dt,xy=TRUE)
head(data)
> head(data)
x y lyr.1 lyr.1
1 0.5 4.5 1 5
2 1.5 4.5 6 5
3 2.5 4.5 11 5
4 3.5 4.5 16 5
5 4.5 4.5 21 5
6 0.5 3.5 2 5
是否可以将原始名称带入图层中（例如：rm和rn）？
在使用 ggplot2 重现图像时，有一个有趣的要求。为图添加第二个标题，将位置设置为右侧（最好对颜色和其他特征进行更改）。它看起来像这样：
使用annotate函数，由于的限制，这种方法不太方便panel。
经过搜索，下面提供了一种方法： 在面板外添加 ggplot 注释？或者两个标题？ 我想知道现在有没有更简单、更直接的方法来实现这一点。
library(tidyverse)
ggplot(mtcars, aes(wt, mpg)) +
geom_point()+
ggtitle('left')+
annotate("text", label="right",
x=Inf, y=Inf, vjust=1, hjust=1)
例如，我想使用 ggpattern 的条纹来填充值大于 50 的区域。但是 breaks 参数（可能是我的问题）无法达到我想要的结果。
library(tidyverse)
d <- expand.grid(x = 1:10 ,y=1:10)
d$value <- 1:100
# The target area is concentrated in the upper right corner
ggplot(data = d,aes(x = x,y = y,fill = value))+
geom_tile()
ggplot(d)+
stat_contour_filled(aes(x = x, y = y, z = value),
breaks = c(-Inf, 50, Inf),
geom = ggpattern::GeomPolygonPattern,
pattern_angle = 60,
fill = 'transparent', color = NA,
pattern = 'stripe',
pattern_fill = 'black',
pattern_density = 0.001,
pattern_spacing = 0.03,
pattern_size = 0.3) +
scale_pattern_manual(values = c(NA,'stripe'))
或者
ggplot(d)+
stat_contour_filled(aes(x = x, y = y, z = value>50),
geom = ggpattern::GeomPolygonPattern,
pattern_angle = 60,
fill = 'transparent', color = NA,
pattern = 'stripe',
pattern_fill = 'black',
pattern_density = 0.001,
pattern_spacing = 0.03,
pattern_size = 0.3)
结果：
