是否有相當於 GRASS GIS 的 r.neighbors 范圍的 R ？

Question

使用這些參數在 QGIS GRASS r.neighbors function 中出現了數千個 DEM geotiffs 錯誤

GRASS r.neighbors 鄰域操作=范圍，鄰域大小=3

出現以下錯誤：-

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

QGIS GRASS r.neighbors 問題

我想在 R 中復制這個操作。

看起來唯一內置的焦點（）函數是mean（），var（）而不是range（）？ 並且 output 柵格的確切范圍與 GRASS 不同的輸入柵格范圍不匹配？

在焦點（） function 或 ZE1E1D3D40573127E9EE0480CAF1283D 中是否存在相當於 GRASS r.neighbors 鄰域操作 = 范圍，鄰域大小 = 3 的矩陣表達式？

R 焦距 function 參考

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 。