在filled.contour 中创建一个围绕零的颜色渐变

Question

我是 R 新手，正在努力寻找我认为相对常见的问题的答案。我正在使用 fill.contour 创建一个变量的世界地图。例如：

z=matrix(rnorm(7008),nrow=96)
x=seq(-176.25,180, by=3.75)
y=seq(-90,90, by=2.5)
filled.contour(x,y,z, plot.axes={axis(1); axis(2); map(add=TRUE, interior=FALSE)} )

其中x & y 为经纬度，z 为数据矩阵。我已经花时间应用我自己的颜色和级别，但是我想要一个颜色渐变，白色分配为零。通过绿色和黄色，负数分级为深蓝色，正数分级为深红色。

我已尝试使用“plotrix”包中的 color.scale 函数

cellcol=matrix(rep("#000000",7008),nrow=96) # replicating the size of my matrix z
cellcol[z>0]=color.scale(z[z>0], c(0,1,1),c(1,1,0),0) # values above zero grading to red 
cellcol[z<0]=color.scale(z[z<0], 0, 0,c(0,1)) # values below zero grading to blue

但是现在我被卡住了。这是最好的方法吗？如果是这样，我如何将其输入到上面的filled.contour 代码中？我敢肯定它很简单，但不能让它工作。

提前感谢您的帮助。

Answer 1

我只是使用ggplot2 实现了你想要的。

您可以尝试以下方法（我以淹没火山数据为例）：

library(ggplot2)
library(reshape2)

## Just an example, I subtract the mean to have positive and negative values for z 
dd <- volcano-mean(volcano)

## Creates a data.frame with columns x, y, z
dd <- melt(dd)
names(dd) <- c('x','y','z')

## Does the contour plot
d <- ggplot(dd, aes(x,y,z=z))
d + geom_tile(aes(fill=z))  + scale_fill_gradient2(low="blue", high="red")

---

我写了一个小函数来做你想要实现的：

myFilled.contour <- function(x = seq(0, 1, length.out = nrow(z)),
                             y = seq(0, 1, length.out = ncol(z)),
                             z, nlevels=30, ...) {
  ma <- max(abs(z))
  lvls <- seq(-ma, ma, length.out = nlevels)
  cols <- colorRampPalette(c("blue","white","red")) (nlevels - 1)
  filled.contour(x, y, z, plot.axes={axis(1); axis(2)},
                 col=cols, levels=lvls, ...)
}

使用filled.contour 并再次使用水下火山dd：

myFilled.contour(z=d)

使用您的数据：

myFilled.contour(x,y,z)

警告：图例包括等高线图中未使用的级别。

希望对你有帮助，

亚历克斯