首页 > 解决方案 > 为什么我的光栅堆栈显示所有单独的光栅图像而不是合并的光栅图像?

问题描述

我正在尝试绘制由一堆光栅图像组成的图像。使用可重现的数据:

library(raster)
b <- brick(system.file("external/rlogo.grd", package = "raster"))
b
#class      : RasterBrick 
#dimensions : 77, 101, 7777, 3  (nrow, ncol, ncell, nlayers)
#resolution : 1, 1  (x, y)
#extent     : 0, 101, 0, 77  (xmin, xmax, ymin, ymax)
#crs        : +proj=merc +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0 
#source     : /Library/Frameworks/R.framework/Versions/4.0/Resources/library/raster/external/rlogo.grd 
#names      : red, green, blue 
#min values :   0,     0,    0 
#max values : 255,   255,  255 

library(rasterVis)
levelplot(b)

这会并排生成 3 个栅格层的图像

我想要做的是绘制它,以便图像仅显示一次“R”徽标,并使用它们的平均值将“红色”、“绿色”和“蓝色”层堆叠在一起。

任何人都可以帮助实现这一目标吗?谢谢!

标签: rstackraster

解决方案


获取 RasterStack(或 RasterBrick)的平均值

library(raster)
b <- brick(system.file("external/rlogo.grd", package = "raster"))

bm <- mean(b)
plot(bm)

推荐阅读