r - 为什么我的光栅堆栈显示所有单独的光栅图像而不是合并的光栅图像?
问题描述
我正在尝试绘制由一堆光栅图像组成的图像。使用可重现的数据:
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”徽标,并使用它们的平均值将“红色”、“绿色”和“蓝色”层堆叠在一起。
任何人都可以帮助实现这一目标吗?谢谢!
解决方案
获取 RasterStack(或 RasterBrick)的平均值
library(raster)
b <- brick(system.file("external/rlogo.grd", package = "raster"))
bm <- mean(b)
plot(bm)
推荐阅读
- python - numpy 2个数组加在一起导致相同的数组但移位
- java - 将代码转换为函数和链接问题
- .net - 如何从命令行获取 Visual Studio Intellisense 错误
- jakarta-ee - 无状态会话 bean 何时获得数据库连接?
- arrays - 将输入拆分为C中数组位置的问题
- c++ - 是否有可能让“命名构造函数”返回一个私有构造的、不可移动的、不可复制的 std::optional
? - algorithm - 是否有死锁检测算法的代码?
- python-3.x - 使用 Slack RTM API 收听消息时没有任何反应
- c++ - 在 C++ 中排列列
- java - 调用 WIN10 ChromeOption 时 Selenium Grid Hub/崩溃