対象地域の DEM によって提供される陰影起伏の平均気温をどのように求めることができるlevelplot
でしょうか? raster
以下の例を考えてみましょう (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 ではなく温度です。
これについて何か考えはありますか?