Понимание вариограммы Гстата - PullRequest
Новая
       71

Понимание вариограммы Гстата

0 голосов
/ 03 октября 2019

Я новичок в использовании gstat и вариограмм и надеюсь, что кто-то может помочь мне объяснить некоторые детали, которых нет в документации.

Я использовал данные по умолчанию и примерс одной из страниц справки для размещения сферической вариограммы https://gisgeography.com/semi-variogram-nugget-range-sill/ и многих других веб-сайтов со следующими определениями:

  • SILL: значение, при котором модель сначала выравнивается.
  • RANGE: Расстояние, на котором модель сначала выравнивается.
  • NUGGET: Значение, при котором полувариограмма (почти) пересекает значение y.

Что имеет смысл. В данных примера встроенная вариограмма дает следующие данные. Выровняйте подоконник, кажется, ближе к 0,45. Кто-нибудь может объяснить эти числа? 

library(gstat)
library(sp)
library(lattice)

data(meuse)
coordinates(meuse) = ~x+y
v = variogram(log(zinc)~x+y, meuse)
v.fit = fit.variogram(v, vgm(1, "Sph", 700, 1))
v.fit
#>   model      psill    range
#> 1   Nug 0.08234213    0.000
#> 2   Sph 0.38866509 1098.571

testplot<-  plot(v, v.fit,
          panel = function(x, ...) {
            panel = vgm.panel.xyplot(x, ...);
            panel.abline(h = v.fit$psill[1:2], v = v.fit$range[2])})
print(testplot)

Создано в 2019-10-03 пакетом Представить (v0.3.0)

1 Ответ

0 голосов
/ 03 октября 2019

Хорошо, я только что понял, что два значения psill складываются в верхнее значение. Это может быть очевидным для других, но я отвечу здесь вместо удаления, если другие люди запутаются, как я, в будущем. 

    plot(thisvgm, fitvgm2, main = paste0(Reduce(paste, deparse(FORMULAS[[MOD]] )), " range = 15"), ylim = c(0,1500),
         panel = function(x, ...) {
           panel = vgm.panel.xyplot(x, ...);
           panel.abline(h = c(fitvgm2$psill[1],sum(fitvgm2$psill)), v = fitvgm2$range[2])})
...