R语言terra包做栅格算术运算报错 raster包正常如何解决

• terra包的global()函数返回值是data.frame类型，行名对应栅格图层名称，第一列为统计计算得到的数值，并非可直接参与栅格算术运算的数值标量。原代码直接将该data.frame代入四则运算时，对象类型不匹配，触发了重复行名的报错。
• raster包的minValue()、maxValue()函数直接返回数值型向量，因此相同逻辑的运算可以正常运行。

只需要提取global()返回结果中的数值部分再参与运算即可，推荐两种常用写法，同时避免了原代码重复三次调用global()计算统计值的性能浪费：

基础修正写法（逻辑和原代码完全一致） #

library(terra)

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

# 提前提取全局最值为数值对象
r_min <- global(r, fun = "min", na.rm=TRUE)[[1]]
r_max <- global(r, fun = "max", na.rm=TRUE)[[1]]

r_n <- (r - r_min) / (r_max - r_min)

更高效的内置函数写法 #

terra提供了栅格标准化的内置函数，大栅格运算速度更快，多波段场景下也支持自动逐波段归一化：

# 用scale函数直接做线性缩放
r_n_scale <- scale(r, center = r_min, scale = r_max - r_min)

运算完成后可查看最值确认归一化结果范围在0-1之间：

global(r_n, c("min", "max"), na.rm=TRUE)

内容的提问来源于stack exchange，提问作者UseR10085