R:从栅格对象计算窗台,范围和块

Question

我需要从栅格图层计算窗台,范围和块金.我已经探索了gstat,usdm包,其中可以创建变异函数但是我找不到给定栅格图层估计这些参数的函数.在大多数函数中,这些参数必须被定义,例如.krigging.

我有不同高度的栅格数据层,看起来类似于

我想从半变异函数的参数到这些数据层得到窗台,金块和范围,以创建类似于此的图:

原始数据层在此处可用作多频带tiff.下面是从一个数字此纸张这进一步说明了这一概念.

Answer 1

这只是一个猜测。这就是我估计半方差的方法

其中n是均值小于总均值的层数。m是所有层的总平均值。r是低于总平均值的每一层的平均值。

s <- stack("old_gap_.tif")
m <- cellStats(mean(s), stat="mean", na.rm=T) # 0.5620522
r <- m[m < 0.5620522]
sem <- 1/53 * (0.5620522 - r)^2
plot(sem, r)