代码之家 › 专栏 › 技术社区 › Rose

是否有GRASS GIS的R.neighbors范围的R等效值?

grass qgis raster r

Rose · 技术社区 · 3 年前

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

GRASS r.neighbors邻域操作=范围,邻域大小=3

出现以下错误:-

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

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

看起来唯一内置的focal()函数是mean()、var(),而不是range()?与GRASS不同,输出光栅的确切范围与输入光栅范围不匹配?

在focal()函数或一般的r中,是否有与GRASS r.neighbors邻域操作=范围,邻域大小=3等效的矩阵表达式?

0 回复 | 直到 3 年前

Robert Hijmans 3 年前

我假设“范围”是指焦点最小值和最大值之间的差异?你可以用(至少)两种方式计算

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 .

推荐文章

Marc B. · 使用ggplot2创建条形图时“缺少值”

1 年前

deschen · tidyverse与外部向量发生突变,该外部向量的元素是数据帧中的列值

1 年前

Laura · 在Shiny中使用可排序的包拖放名称,这些名称将成为图表

1 年前

Mallikarjun M · 如何使用随机森林进行时间序列预测?

1 年前

ly li · 模型摘要:当表格形状改变时,拟合优度消失

1 年前

C.Robin · 将marginaffects::predictions()的结果连接回main df?

1 年前

monotonic · 如何将格式为“col1+col3+col4”的数据帧的行名转换为一列数字向量“c(1,3,4)”?

2 年前

Shawn Hemelstrand · 为什么我的自定义errorbar函数不能在R中工作?

2 年前

RoyBatty · 统计每个字符在整个数据集中出现的次数

2 年前

stats_noob · R: 记录某个“行为”发生的循环的索引?

2 年前