R - ggplot - stat_contour无法生成轮廓线

Question

我试图通过stat_contour()我的ggplot/ ggplot2-plot 添加轮廓线.不幸的是,我无法提供应该评估点值的真实数据.但是,另一个易于复制的示例表现相同:

testPts <- data.frame(x=rep(seq(7.08, 7.14, by=0.005), 200))
testPts$y <- runif(length(testPts$x), 50.93, 50.96)
testPts$z <- sin(testPts$y * 500)

ggplot(data=testPts, aes(x=x, y=y, z=z)) + geom_point(aes(colour=z))
       + stat_contour()

这会导致以下错误消息:

if(nrow(layer_data)== 0)return()中的错误:参数长度为零另外:警告消息:无法生成轮廓数据

该示例与stackoverflow或官方手册/教程中发布的其他示例没有什么不同,如果我提供更多规范,seeminlgy无关紧要stat_contour.看来,该函数不会像错误消息所指示的那样传递数据(-layer).

感谢您的想法和建议!

Answer 1

使用stat_density2d而不是stat_contour使用不规则间隔的数据.

library(ggplot2)

testPts <- data.frame(x=rep(seq(7.08, 7.14, by=0.005), 200))
testPts$y <- runif(length(testPts$x), 50.93, 50.96)
testPts$z <- sin(testPts$y * 500)

(ggplot(data=testPts, aes(x=x, y=y, z=z))
+ geom_point(aes(colour=z))
+ stat_density2d()
)

Answer 2

该问题的一个解决方案是生成规则网格以及关于该网格的点值的插值.以下是我为多个数据字段中的一个执行此操作的方法:

pts.grid <- interp(as.data.frame(pts)$coords.x1, as.data.frame(pts)$coords.x2, as.data.frame(pts)$GWLEVEL_TI)
pts.grid2 <- expand.grid(x=pts.grid$x, y=pts.grid$y)
pts.grid2$z <- as.vector(pts.grid$z)

这导致数据帧可以在ggplot中使用,stat_contour()当在该函数的data-parameter中定义时:

(ggplot(as.data.frame(pts), aes(x=coords.x1, y=coords.x2, z=GWLEVEL_TI))
#+ geom_tile(data=na.omit(pts.grid2), aes(x=x, y=y, z=z, fill=z))
+ stat_contour(data=na.omit(pts.grid2), binwidth=2, colour="red", aes(x=x, y=y, z=z))
+ geom_point()
)

这个解决方案很可能包括不必要的转换,因为我还不知道更好.此外,我必须为每个数据字段单独生成相同的网格生成,然后再将它们组合在一个数据框中 - 这对于更大的数据集而言并不像我希望的那样高效.