如何在光栅堆栈中绘制超过 16 个图形？答案

【问题标题】：how to plot more than 16 graphics in a raster stack?如何在光栅堆栈中绘制超过 16 个图形？
【发布时间】：2019-07-22 16:59:12
【问题描述】：

我有一个 32 层的光栅堆栈，我想用一些层制作一个图，但我可以在图中获得的最大层数是 16 (4 x 4)。

covs<-stack(NDVI,BTA,biovars)
plot(covs)

我需要在同一个图表中再插入 3 层有人可以帮帮我吗？

【问题讨论】：

标签： r plot stack raster

【解决方案1】：

将maxnl参数设置为32（见raster中plot的help）：

plot(covs, maxnl=32)

这是一个例子：

library(raster)
fn <- system.file("external/test.grd", package="raster")
r <- raster(fn)
stk <- stack(r, fn, r, fn, r, fn, r, fn, r, fn, r, fn, r, fn, r, fn, r, fn, r, fn, 
             r, fn, r, fn, r, fn, r, fn, r, fn, r, fn, r, fn, r, fn, r, fn, r, fn) 
plot(stk, nc=6, maxnl=32)

【讨论】：

没有解决问题，我还有16个图形。用levelplot()试试，但是图形的比例不一样。
@franciscocorvalan 我在我的电脑上用 32 层测试了这个解决方案，它可以工作。你能分享你的 covs 对象（或 NDVI、BTA、biovars）吗？

【解决方案2】：

Marco Sadri 的解决方案应该可行，你能show(covs) 吗？但您也可以“手动”执行此操作

par(mfrow=c(6,6))
for (i in 1:nlayers(stk)) {
   plot(stk[[i]])
}

【讨论】：