如何将三角形添加到 R 中的 ggplot2 颜色条以指示超出范围的值？

我想用ggplot2其中一些填充值被剪裁，即高于或低于色阶限制的值显示为最小/最大颜色。我可以使用以下组合让它像这样工作limit and oob（超出范围）：

library(ggplot2)
library(scales)
ggplot() + ... + scale_fill_viridis(na.value="white", limit=c(0, 10), oob=squish)

但颜色栏中没有任何信息表明存在超出限制的值。 如何在 ggplot 中重现这个 matplotlib 示例：https://stackoverflow.com/a/32072348

具体来说，如何获取颜色条末尾的三角形？

据我所知，没有一个包可以在 ggplot2 中实现颜色条的三角形末端（但如果有，请告诉我！）。但是，我们可以实现我们自己的。我们需要一个自定义指南的构造函数以及绘制它的方法。大部分内容已经在guide_colourbar()以及他们的类的方法，所以我们需要做的就是标记我们自己的类并扩展guide_gengrob方法。下面的代码应该适用于垂直方向的颜色条。你需要了解一些关于grid包装和gtable包跟随。

library(ggplot2)
library(gtable)
library(grid)

my_triangle_colourbar <- function(...) {
  guide <- guide_colourbar(...)
  class(guide) <- c("my_triangle_colourbar", class(guide))
  guide
}

guide_gengrob.my_triangle_colourbar <- function(...) {
  # First draw normal colourbar
  guide <- NextMethod()
  # Extract bar / colours
  is_bar <- grep("^bar$", guide$layout$name)
  bar <- guide$grobs[[is_bar]]
  extremes <- c(bar$raster[1], bar$raster[length(bar$raster)])
  # Extract size
  width  <- guide$widths[guide$layout$l[is_bar]]
  height <- guide$heights[guide$layout$t[is_bar]]
  short  <- min(convertUnit(width, "cm",  valueOnly = TRUE),
                convertUnit(height, "cm", valueOnly = TRUE))
  # Make space for triangles
  guide <- gtable_add_rows(guide, unit(short, "cm"),
                           guide$layout$t[is_bar] - 1)
  guide <- gtable_add_rows(guide, unit(short, "cm"),
                           guide$layout$t[is_bar])
  
  # Draw triangles
  top <- polygonGrob(
    x = unit(c(0, 0.5, 1), "npc"),
    y = unit(c(0, 1, 0), "npc"),
    gp = gpar(fill = extremes[1], col = NA)
  )
  bottom <- polygonGrob(
    x = unit(c(0, 0.5, 1), "npc"),
    y = unit(c(1, 0, 1), "npc"),
    gp = gpar(fill = extremes[2], col = NA)
  )
  # Add triangles to guide
  guide <- gtable_add_grob(
    guide, top, 
    t = guide$layout$t[is_bar] - 1,
    l = guide$layout$l[is_bar]
  )
  guide <- gtable_add_grob(
    guide, bottom,
    t = guide$layout$t[is_bar] + 1,
    l = guide$layout$l[is_bar]
  )
  
  return(guide)
}

然后您可以使用您的自定义指南作为guide争论的规模。

g <- ggplot(mtcars, aes(mpg, wt)) +
  geom_point(aes(colour = drat))

g + scale_colour_viridis_c(
    limits = c(3, 4), oob = scales::oob_squish,
    guide = my_triangle_colourbar()
  )

确实没有一种自然的方法可以对越界值进行不同的着色，但您可以在极值附近制作非常小的切片，并使用不同的颜色。

g + scale_colour_gradientn(
    colours = c("red", scales::viridis_pal()(255), "hotpink"),
    limits = c(3, 4), oob = scales::oob_squish,
    guide = my_triangle_colourbar()
  )

^{Created on 2021-07-19 by the reprex package (v1.0.0)}