filled.contour -描述z值级别的变化

Question

我必须生成在z值中具有某种对称性的filled.contour图，如下面的代码所示

x <- 1:5
y <- 1:5
z <- matrix(outer(x,y,"+"),nrow=5)
filled.contour(x,y,z)
filled.contour(x,y,z,color.palette=rainbow)
z2 <- z
z2[5,5] <- Inf
filled.contour(x,y,z2,col=rainbow(200),nlevels=200)

我希望能够描绘z值的级别，例如某些级别或增量(2,4,6,8,10)增量，以便获得以下图

​

​

我想保留col=rainbow(200)的级别，以便绘制平滑的颜色和颜色过渡，但我想描绘某些值/增量。

这个是可能的吗？我将感谢一个基本的R解决方案，但如果这是不可能的，任何解决方案将是好的(ggplot，网格)

Answer 1

您可以使用countourLines()添加行，使用text()添加文本：

txtlab <- c("2", "4", "6", "8", "Inf")
txtpos <- c(1.25, 2, 3, 4, 4.75)
ctlns <- contourLines(x, y, z2, levels=c(3, 5, 7, 9))
filled.contour(x,y,z2,col=rainbow(200),nlevels=200, 
     plot.axes={axis(1); axis(2); text(txtpos, txtpos, txtlab, cex=1.5);
     sapply(1:4, function(x) lines(ctlns[[x]][[2]], ctlns[[x]][[3]], lwd=2))
})

​

​

Answer 2

这里的每个整数都是黑色的(使用is.integer中的一个例子)，它可以很容易地修改为只返回偶数的black：

getcolor <- function(z, n) {
  is.wholenumber <-
    function(x, tol = .Machine$double.eps^0.5)  abs(x - round(x)) < tol
  pp <- pretty(range(z, finite=TRUE), n)
  ifelse(is.wholenumber(pp), 'black', rainbow(n))
}

filled.contour(x,y,z2,col=getcolor(z2, 200),nlevels=200)

​

​

Answer 3

我最近做了一些类似的事情，我使用了这种类型的代码：

x <- 1:5
y <- 1:5
z <- matrix(outer(x,y,"+"),nrow=5)
z2 <- z
z2[5,5] <- Inf
filled.contour(x,y,z2,col=rainbow(200),nlevels=200
plot.axis = {axis(1);axis(2);contour(x,y,z2,add=TRUE)}) #that's when you add the contours

Result (click on link to see)