【发布时间】:2017-03-18 12:05:25
【问题描述】:
我正在尝试使用 Raster 包中的 hillShade 函数,并得到“找不到函数“覆盖””错误。
我的程序读入一个 shapefile。 将其转换为栅格。 使用地形函数获取坡度和坡向。 然后,当它应该创建 hillShade 时,会发生此错误: “hillShade 中的错误(fieldSlope,fieldAspect,角度 = 45,方向 = 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
非常感谢任何帮助。谢谢
【问题讨论】:
-
我无法重现这个 - 对我来说很好。您是否尝试过重新安装
raster? -
我在另一台计算机上尝试并按照您的建议更新了软件包,并且对我来说也很好。抱歉,我一定没有尝试在我的电脑上更新软件包。感谢您的帮助,它为我节省了数小时的额外痛苦。