R语言terra包，focal焦点量
# Thu Dec 24 11:35:25 2020 -
# 字符编码：UTF-8
# R 版本：R x64 4.0.3 for window 10
# cgh163email@163.com
# 个人笔记不负责任
# —— 拎了个梨🍐
.rs.restartR()
require(terra)
# focal焦点量
rm(list=ls());gc()
# 使用权重矩阵（可能与函数结合）计算焦点细胞附近的焦点（“移动窗口”）值。
r <- rast(ncols=10, nrows=10, ext(0, 10, 0, 10))
values(r) <- 1:ncell(r)
f <- focal(r, w=3, fun=function(x, ...)quantile(x, .25, ...), na.rm=TRUE)
f <- focal(r, w=3, fun="mean")
plot(r)
plot(f)
# 以下两个语句是等效的
a <- focal(r, w=matrix(1/9, nc=3, nr=3))
b <- focal(r, w=3, fun=mean, na.rm=FALSE)
# but this is different
d <- focal(r, w=3, fun=mean, na.rm=TRUE)
a
b
d
