Mais

Criação de boxplots de dois rasters usando o pacote Raster em R, usando apenas um subconjunto de classes de um

Criação de boxplots de dois rasters usando o pacote Raster em R, usando apenas um subconjunto de classes de um


Em R usando o pacote raster, você pode criar boxplots a partir de dois rasters (x e y), assim:

boxplot (x, y)

Onde:

x = Raster * objeto y = Se x for um objeto RasterLayer, y pode ser um RasterLayer adicional para agrupar os valores de x por 'zona'

O que funciona perfeitamente quando quero agrupar os valores de x por todas as classes em meu raster temático y.

Mas se eu quiser apenas agrupar os valores em x por um subconjunto das classes temáticas em y, há uma maneira de fazer isso? Eu poderia substituir quaisquer valores em y fora dos valores que estou interessado pelos valores NODATA, mas esperava um método mais limpo e simples.


Descobri que o operador de correspondência de valor% em% funciona perfeitamente, para substituir valores raster que eu não quero (por exemplo, 4, 7, 10) por NA, e então posso chamar o método boxplot novamente.

y [y% em% c (4,7,10)] <- NA

Assista o vídeo: Wykres pudełkowy z wąsami