我已经从包含R:
中景观高度值的栅格创建了以下过渡对象wd <- "C:/Users/LG/Dropbox/Random Walk"
setwd(wd)
library(gdistance)
library(raster)
r <- raster("altitude.tif")
altDiff <- function(x){x[2] - x[1]}
hd <- transition(r, altDiff, 8, symm=FALSE)
slope <- geoCorrection(hd)
正如您所看到的,我用来创建所述目标的函数是两个单元格之间的高度差异,然后通过“地理定位”功能将其校正为单元格之间的斜率。如何才能访问我选择的两个单元格之间的斜率?为了完成我的脚本,我需要获得我选择的单元格和与其相邻的所有八个单元格之间的斜率值。所以基本上我想知道如何访问我选择的两个单元格的转换值? 如果有人能帮助我,我会非常感激。 提前谢谢。
答案 0 :(得分:1)
我建议你阅读gdistance包的插图,这完全解释了这一点。
要获得单元格1和单元格2之间的斜率,您可以执行以下操作:
slope[1,2]
查看栅格包中的adjacent
函数,以获取8个相邻单元格的完整列表。