我需要从栅格图层计算窗台,范围和块金。我已经探索了gstat,usdm包,其中可以创建变异函数但是我找不到给定栅格图层估计这些参数的函数。在大多数函数中,这些参数必须被定义,例如。 krigging。
我想从适合这些数据层的半变异函数参数得到窗台,金块和范围,以创建类似于此的图:
答案 0 :(得分:3)
使用gstat,这是一个例子:
library(raster)
library(gstat)
demo(meuse, ask = FALSE, echo = FALSE)
set.seed(131) # make random numbers reproducible
# add some noise with .1 variance
meuse.grid$dist = meuse.grid$dist + rnorm(nrow(meuse.grid), sd=sqrt(.1))
r = raster(meuse.grid["dist"])
v = variogram(dist~1, as(r, "SpatialPixelsDataFrame"))
(f = fit.variogram(v, vgm("Sph")))
# model psill range
# 1 Nug 0.09035948 0.000
# 2 Sph 0.06709838 1216.737
f$psill[2] # sill
# [1] 0.06709838
f$range[2] # range
# [1] 1216.737
f$psill[1] # nugget
# [1] 0.09035948
为r
插入自己的栅格,它应该可以使用。更改Sph
以适应另一个变异函数模型,尝试plot(v,f)
验证图。
答案 1 :(得分:1)