0

我正在尝试在 R 中制作一个包含许多不同栅格的多面板图。我正在使用 Windows 机器。

栅格将按以下格式显示：

第一行应包含一个栅格（居中），带有可见的“主”标题和 ylabel。

第 2 行到第 6 行包含两个相邻的栅格（即两列）。仅在第 2 行需要“主要”标题，因为标题适用于其余行。

第 7 行应包含分类图例。类似于此处提供的答案中描述的内容：

https://gis.stackexchange.com/questions/73143/legend-of-a-raster-map-with-categorical-data

间距很重要。我想最小化第 2 行和第 8 行之间的距离，同时在第 1 行和第 2 行之间保留一些额外的空间以允许标题。每个栅格应具有相同的大小。

请在下面找到我正在处理的可复制示例。我不确定为什么会产生“主要”标题的彩色背景，并且是在我尝试使用“top.padding”和“bottom.padding”命令减少空白之后发生的。

library(raster)

library(rasterVis)

f <- system.file("external/test.grd", package="raster")

r <- raster(f)

rmat = as.matrix(r) 

pc = c(0,0.6,0.8,0.9,0.92,0.95,0.97,0.98,0.99,1) 

digits_n = c(0,1,c(rep(0,length(pc)-2)))

color_breaks = NULL

for (i in 1:length(pc)){color_breaks = c(color_breaks, signif(quantile(c(rmat[which(!is.na(rmat))]),pc[i]), digits = digits_n[i]))}

color_breaks = unique(color_breaks)

color_type = c(colorRampPalette(c("snow1","snow3","seagreen","orange","firebrick"))(length(color_breaks))) 

ss = stack(r)

p0 = levelplot(ss, col.regions=color_type, colorkey=FALSE,margin=FALSE, at=color_breaks,names.attr=c(2000), xlab=NULL, ylab=c("A0"), scales=list(draw=FALSE), par.settings=list(layout.heights=list(bottom.padding=-2)))

ss = stack(r,r)

p1 = levelplot(ss,col.regions=color_type, colorkey=FALSE, margin=FALSE, at=color_breaks, names.attr=c(2050,2100), xlab=NULL, ylab=c("A1"), scales=list(draw=FALSE), layout=c(2,1), par.settings=list(layout.heights=list(top.padding=-2,bottom.padding=-2)))

ss = stack(r,r)

p2 = levelplot(ss,col.regions=color_type, colorkey=FALSE, margin=FALSE, at=color_breaks, names.attr=c("",""), xlab=NULL, ylab=c("A2"), scales=list(draw=FALSE), layout=c(2,1), par.settings=list(layout.heights=list(top.padding=-2,bottom.padding=-2))) 

ss = stack(r,r)

p3 = levelplot(ss,col.regions=color_type, colorkey=FALSE, margin=FALSE, at=color_breaks, names.attr=c("",""), xlab=NULL, ylab=c("A3"), scales=list(draw=FALSE), layout=c(2,1), par.settings=list(layout.heights=list(top.padding=-2,bottom.padding=-2))) 

ss = stack(r,r)

p4 = levelplot(ss,col.regions=color_type, colorkey=FALSE, margin=FALSE, at=color_breaks, names.attr=c("",""), xlab=NULL, ylab=c("A4"), scales=list(draw=FALSE), layout=c(2,1), par.settings=list(layout.heights=list(top.padding=-2,bottom.padding=-2)))

ss = stack(r,r)

p5 = levelplot(ss,col.regions=color_type, colorkey=FALSE, margin=FALSE, at=color_breaks, names.attr=c("",""), xlab=NULL, ylab=c("A5"), scales=list(draw=FALSE), layout=c(2,1), par.settings=list(layout.heights=list(top.padding=-2,bottom.padding=-2))) 

windows()

print(p0, split= c(1,1,1,7), between = c(3,3))

print(p1, split= c(1,2,1,7), newpage=FALSE)

print(p2, split= c(1,3,1,7), newpage=FALSE)

print(p3, split= c(1,4,1,7), newpage=FALSE)

print(p4, split= c(1,5,1,7), newpage=FALSE)

print(p5, split= c(1,6,1,7), newpage=FALSE)

然后，我想在第 7 行的图例中添加类似于以下内容的内容：

legend("center", 
    ncol = length(color_breaks)/2,
    legend=c(color_breaks),
    fill=color_type,cex=0.6, 
    title=expression(bold("Magnitude [ UNITS ]")), box.col = FALSE, xjust = 0, title.adj = 0
    )

我会发布我获得的图像，但我想我需要声誉积分。

为什么标题中的背景颜色出现了？如何减小第 2 行和第 6 行之间的间距，同时增加第 1 行和第 2 行之间的间距？现在如何在第 7 行添加分类图例？

谢谢你。

4

1 回答 1

1

我将尝试简要回答您的问题：

  1. 条带区域的背景颜色是lattice. 将其rasterTheme设置为白色，但您必须使用这里描述par.settings的而不是。col.regions
  2. 您应该尝试 的panel.width和panel.height参数print.trellis。
  3. levelplot能够绘制带有适当图例的分类栅格，但您必须相应地定义栅格。在这里你会找到一个例子。
于 2015-01-24T15:55:18.780 回答