gpt4 book ai didi

R:从光栅对象计算基台、范围和 block 金

转载 作者:行者123 更新时间:2023-12-04 11:25:23 25 4
gpt4 key购买 nike

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

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

我想从拟合到这些数据层的半变异函数参数中获取基石、金块和范围,以创建与此类似的图:enter image description here

原始数据层可用 here作为多频段 tiff。这是来自 this 的图论文进一步说明了这个概念。

enter image description here

最佳答案

使用 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)来验证情节。

关于R:从光栅对象计算基台、范围和 block 金,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/36033864/

25 4 0
Copyright 2021 - 2024 cfsdn All Rights Reserved 蜀ICP备2022000587号
广告合作:1813099741@qq.com 6ren.com