Breeze's questions

Supondo que eu precise calcular a variância dos dados raster para a mesma longitude, usar a sappfunção é o caminho certo a seguir? Mas eu falhei.

Agora posso usar o método mais complexo de conversão para data.frame e agrupamento para calcular a variância:

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))

Existe uma função mais direta em terra? Por favor, me diga se você sabe.

Há também cenários mais complexos em que isso precisa ser feito em várias camadas de dados raster ao mesmo tempo . Como isso deve ser feito?

Ao usar as.data.frame para manipular múltiplas camadas, para a camada mesclada toda a nomenclatura é lyr.1. Existe uma maneira melhor de resolver esse problema de nomenclatura?

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

É possível trazer os nomes originais para as camadas (no exemplo: rme rn)?

Ao reproduzir uma imagem usando ggplot2, havia um requisito interessante. Adicione um segundo título à figura, defina o local para a direita (seria legal fazer alterações na cor e em outros recursos). Parece algo assim:

Usando a annotatefunção, essa abordagem não é muito conveniente devido às limitações de panel.

Após a busca, o seguinte fornece uma maneira de fazer isso: Adicionar anotação ggplot fora do painel? Ou dois títulos? Gostaria de saber se agora há uma maneira mais simples e direta de fazer isso.

library(tidyverse)

ggplot(mtcars, aes(wt, mpg)) + 
  geom_point()+
  ggtitle('left')+
  annotate("text", label="right", 
           x=Inf, y=Inf, vjust=1, hjust=1)

Por exemplo, quero usar as listras do ggpattern para preencher áreas com valores maiores que 50. Mas o parâmetro breaks (provavelmente meu problema) não atinge o resultado que desejo.

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'))

ou

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)

resultado: