我想对r中的光栅数据进行模糊求和,以形成累积电阻层进行研究。我找到了用矢量数据进行模糊求和的包和函数,想知道是否有人可以共享资源,专门将光栅层与模糊逻辑相结合。
谢谢
您可以将spatialEco::fuzzySum
用于矢量数据和光栅数据。
例如,对于三个terra光栅:rast1
、rast2
、rast3
,其工作方式如下:
rFuzzySum <- spatialEco::fuzzySum(c(rast1, rast2, rast3))
如果你把它写开,它会是:
rFuzzySum <- (1 - ( (1 - rast1) *
(1 - rast2) *
(1 - rast3) ) )
下面是如何使用MattKummu 的建议来实现这一点的示例
示例数据
library(terra)
x <- rast(system.file("ex/logo.tif", package="terra"))
x <- x / max(minmax(x))
的两种方法
a <- 1 - prod(1 - x)
b <- spatialEco::fuzzySum(x)