r - R 中的比例和颜色光栅图例

Question

如何缩放光栅图例，使图例标签等距但值不规则？

## Generate Random Raster
rast <- raster( matrix(
    sample(c(0,10),9, replace=T) + runif(9),
    nrow=3, ncol=3) )
rast[4] <- NA

## legend parameters
zlim  <- c(0,15)
zlabs <- c(0,.1,.25,.5,1, 2.5, 5,10, 15)
zat <- seq(length(zlabs))

## Colors
greycols <- grey(seq(1, 0, length=length(zlabs)+1))

## Plot
raster::plot( rast,
    zlim=zlim,
    breaks=zlabs,
    lab.breaks=zat,
    col=greycols
   # axis.args=list(
   #     at=c(0,zat),
   #     labels=c("NA", zlabs))
)

Answer 1

如果我正确理解你想要什么，一个解决方案可能是为标签存储光栅的旧值，然后根据你对绘图的需要更改实际值(以便色标和实际值是匹配的，而标签是你喜欢的)。

为了清楚起见，我只是使用了另一个示例和其他颜色。在您的情况下，数据转换只是它们的 rank ，因为您需要等间隔:

library(raster)
library(scales)
library(ggplot2)

r = raster(matrix(c(0.75, 0.1, 1.05, 9.8, 9, 1), nrow=2, byrow=T))
zlim = c(0,10)
mycols = c("blue", "red", "green", "magenta", "cyan", "black")
# Convert to DF for ggplot
df = data.frame(rasterToPoints(r))
# Transform the data (and rescale to occupy the full zlim space)
df$layer = rescale(rank(r[]), zlim) 
# Plot the "wrong" data, but label them according to the original ones
ggplot() + geom_raster(data=df, aes(y=y, x=x, fill=layer)) +
  scale_fill_gradientn(colours = mycols,
                       breaks = df$layer, 
                       labels = r[], 
                       limits = zlim)

我建议使用 ggplot2，因为它可以很好地控制色标(参见 the online manual)，而且可以在这里找到许多问题和答案。此外，在 scale_fill_gradientn 中，values 参数在某些情况下可能很有趣 - 在您的情况下，我不太确定它是否会满足您的需求，但值得一提的是存在此功能。