如何使 ggplot2 等高线图类似于 lattice:filled.contour()？

Question

我一直在学习 ggplot2，并希望将它用于我所有的 R 绘图。但是，我还没有找到一种方法来制作看起来类似于传统等高线图的等高线图，就像可以使用 lattice:filled.contour() 获得的一样。例如：

#define data
x<-seq(1,11,1)
y<-seq(1,11,1)
xyz.func<-function(x,y) {-10.4+6.53*x+6.53*y-0.167*x^2-0.167*y^2+0.0500*x*y}

#contour plot using lattice graphics and R Color Brewer
library(lattice) #for filled.contour()
library(RColorBrewer) #for brewer.pal()
z.lattice<-outer(x,y,xyz.func)
filled.contour(x,y,z.lattice,nlevels=6,col=brewer.pal(6,"YlOrRd"))

这给了我一个很好的等高线图。

现在，让我们在 ggplot2 中尝试同样的事情。根据我读过的所有内容（特别是Drawing labels on flat section of contour lines in ggplot2），我能想到的最好的是：

#contour plot using ggplot2
library(ggplot2)
library(reshape2) #for melt()
z.molten<-melt(z.lattice) 
names(z.molten) <- c("x", "y", "z")
v<-ggplot(z.molten, aes(x,y,z=z))+
    geom_tile(aes(fill=z))+
    stat_contour(bins=6,aes(x,y,z=z), color="black", size=0.6)+
    scale_fill_gradientn(colours=brewer.pal(6,"YlOrRd"))
v

这个图和fill.contour()有同样的基本思想，但是彩色的瓦片不太符合轮廓。

我也没有成功更改图块的大小。

关于如何使 ggplot2 的输出更接近 fill.contour() 的输出有什么建议吗？

Answer 1

您的问题的本质似乎是如何在 ggplot 中生成具有 离散 填充轮廓的等高线图，而不是使用传统的 geom_tile(...) 方法获得的连续等高线。这是一种方法。

x<-seq(1,11,.03)                    # note finer grid
y<-seq(1,11,.03)
xyz.func<-function(x,y) {-10.4+6.53*x+6.53*y-0.167*x^2-0.167*y^2+0.0500*x*y}
gg <- expand.grid(x=x,y=y)
gg$z <- with(gg,xyz.func(x,y))      # need long format for ggplot
library(ggplot2)
library(RColorBrewer)               #for brewer.pal()
brks <- cut(gg$z,breaks=seq(0,100,len=6))
brks <- gsub(","," - ",brks,fixed=TRUE)
gg$brks <- gsub("\\(|\\]","",brks)  # reformat guide labels
ggplot(gg,aes(x,y)) + 
  geom_tile(aes(fill=brks))+
  scale_fill_manual("Z",values=brewer.pal(6,"YlOrRd"))+
  scale_x_continuous(expand=c(0,0))+
  scale_y_continuous(expand=c(0,0))+
  coord_fixed()

使用，例如，scale_x_continuos(...) 只是为了摆脱多余的空间 ggplot 围绕轴限制；大多数事情都很好，但在等高线图中会分散注意力。 coord_fixed(...) 的使用只是将纵横比设置为 1:1。这些是可选的。