我正在尝试使用Raster包中的hillShade函数,并获得“找不到函数”覆盖“”错误。
我的程序读入shapefile。 将其转换为栅格。 使用terrain函数来获取斜率和方面。 然后当它应该创建hillShade时会发生以下错误: “hillShade中的错误(fieldSlope,fieldAspect,angle = 45,direction = 180 +: 找不到功能“叠加”“
如果有什么关于hillShade函数和stackoverflow上的覆盖我错过了它。所以请随意指出我正确的方向。
该计划如下:
library(rgdal)
library(raster)
setwd("C:/Users/Etching Memories/Documents/Devenshire/Shape/2015/Harvest")
field<-readOGR(".","DEVENSHIRE_Bouge_Big_Harvesting_20151010")
ras<-raster(ncols=40,nrows=40,xmn=min(field$coords.x1),xmx=max(field$coords.x1),ymn= min(field$coords.x2),ymx=max(field$coords.x2))
rasfield<-rasterize(field,ras,field="Elevation",fun=mean)
fieldSlope<-terrain(rasfield,opt="slope")
fieldAspect<-terrain(rasfield,opt="aspect")
hillShade(fieldSlope,fieldAspect,angle=45,direction=180+45)
压缩的形状文件位于: https://www.dropbox.com/s/06pfbusfrxun1fr/bouge%20file.zip?dl=0
非常感谢任何帮助。感谢
答案 0 :(得分:0)
正如Jbaums所说,重新安装光栅包是可行的,并且是正确的答案。