使用levelplot R在DEM上显示栅格

时间:2016-09-08 17:51:48

标签: r raster terrain levelplot

一个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。

对此有何想法?

non-SOAP HTTP request library for node

0 个答案:

没有答案