r - 使用 levelplot R 在 DEM 上显示栅格

Question

levelplot一个人怎么能raster说由感兴趣区域的 DEM 提供的山体阴影的平均温度？考虑下面的例子（来自 Robert H）：

library(raster)
library(rasterVis)
dem <- getData('alt', country='CHE')# get DEM for region of interest
slas <- slopeAspect(dem)#get slope and aspect
hill <- hillShade(slas[[1]], slas[[2]], 40, 270)# compute hillshade 

plot(hill, col=grey(0:100/100), legend=FALSE, main='Switzerland')

• 我喜欢山体阴影，因为它清楚地显示了山丘和山谷等。
• 现在给定瑞士的平均温度栅格，我怎样才能将它叠加在山上plot(dem, col=rainbow(25, alpha=0.35), add=TRUE)呢？我正在levelplot为我的映射使用函数。
• 温度数据集可以从获得climate <- getData('worldclim', var='bio', res=2.5)，然后plot(climate$bio1, main="Annual Mean Temperature")

所需的输出应如下图所示，但为温度而不是 DEM。

对此有什么想法吗？