将曲线添加到栅格地图？

Question

是否可以将等值线或类似“filled.contour”的内容添加到栅格地图中？

E.g。到Bioclim变量？ （RasterLayer）

# Download Bioclim data 
library(dismo)
tmin_06<-getData("worldclim", var="tmin", res=0.5, lon=10, lat=70)
# plot tmin 
plot(tmin$tmin1_06)

And then add e.g. filled.contour {graphics}

我似乎无法从RasterLayer中提取正确的值并得到错误按摩： “不能强迫'S4'类型为'double'类型的矢量”

或者您会用什么来说明地图中的变化方向，如温度或高程地图中的等距离？

谢谢！

丹尼尔

Answer 1

以下是两种方法：

# example data
f <- system.file("external/test.grd", package="raster")
r <- raster(f)

plot(r)
contour(r, add=TRUE)

# or
filledContour(r)