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）中的注釋）。

但是，使用此解決方案，您可以松開圖例欄。