R：在轮廓图上添加一条线

Question

我想在R中使用filled.contour（）绘制函数f（x，y）并添加一条线/曲线以标识该函数的值为0的点。为了使示例可重写，我的函数的值是火山数据集中的值，而不是寻找f（x，y）= 0，我们想添加一条线/曲线来标识火山的值500。如何执行此操作？ 下面的代码正确地在X = 500和Y = 500处添加了一个点，但是如何添加一条线以使只有volcano = 500的点通过一条线连接？ 我只想使用基本图形。

x <- 10*1:nrow(volcano)
y <- 10*1:ncol(volcano)
filled.contour(x, y, volcano, color = terrain.colors, plot.axes = { axis(1); axis(2); points(500, 500)})

Answer 1

一个不完美的解决方案：

#create a new plot
contour(x,y,volcano)
#use .filled.colour instead of filled.contour
.filled.contour(x, y, volcano,levels=seq(90,200,1),col=terrain.colors(109))
# add a contour plot with specific levels on top of the filled contour
contour(x,y,volcano, level=130,add=T)

您可以通过将vectoor传递到参数level来指定几条轮廓线：level = c（130,150）将绘制高度为130和150的轮廓

问题在于，fill.contour广告的图例栏使简单地调用fill.contour然后轮廓使两张图片不匹配（请参阅help（filled.contour）中的注释）。

但是，使用此解决方案，您可以松开图例栏。