是否有相当于 GRASS GIS 的 r.neighbors 范围的 R？

Question

使用这些参数的QGIS GRASS r.neighbors函数中有几千个DEM地理iff出错

GRASS r.neighbors邻域运算=范围，邻域大小=3

出现以下错误：-

2021-04-23T15:51:51     WARNING    Duplicate parameter coordinates registered for alg v.net.visibility

QGIS GRASS r.neighbors发布

我想在R.中复制这个操作

看起来只有内置的focal((函数是mean((，var((而不是range((？与GRASS不同，输出光栅的确切范围与输入光栅范围不匹配？

在focal((函数或一般的r中，是否存在与GRASS r.neighbors邻域运算=范围，邻域大小=3等价的矩阵表达式？

R聚焦函数参考

Answer 1

我假设；范围"；你指的是焦点最小值和最大值之间的差异？你可以用(至少(两种方法计算

library(terra)
f <- system.file("ex/elev.tif", package="terra")
r <- rast(f)

xmn <- focal(r, 3, min)
xmx <- focal(r, 3, max)
rng <- xmx - xmn

或者像这个

x <- focal(r, 3, function(i) diff(range(i)))

在存在NA的情况下，结果略有不同。您可以使用na.rm=TRUE。