R: クリギングふたたび

ひさしぶりにRでクリギングをしてみる。かなりわすれている。library(geoR)のdata(ca20)をサンプルにしてみる。

library(geoR)
data(ca20)

png("Rplot001.png", width = 900, height = 900,
    pointsize = 14)
plot(ca20)
dev.off()

出力

バリオグラム。こういうやり方でよかったのか。あとで教科書を確認しよう。

v <- variog(ca20)
vf <- variofit(v, ini.cov.pars = c(100, 100),
               cov.model = "gaussian")
summary(vf)

png("Rplot002.png", width = 500, height = 500,
    pointsize = 14)
plot(v)
lines(vf)
dev.off()

出力

$pmethod
[1] "WLS (weighted least squares)"

$cov.model
[1] "gaussian"

$spatial.component
  sigmasq       phi 
 90.90173 309.20958 

$spatial.component.extra
kappa 
  0.5 

$nugget.component
   tausq 
60.39267 

$fix.nugget
[1] FALSE

$fix.kappa
[1] TRUE

$practicalRange
[1] 535.1856

$sum.of.squares
   value 
244168.1 

$estimated.pars
    tausq   sigmasq       phi 
 60.39267  90.90173 309.20958 

$weights
[1] "npairs"

$call
variofit(vario = v, ini.cov.pars = c(100, 100), cov.model = "gaussian")

attr(,"class")
[1] "summary.variomodel"

で、クリギング。

loc <- expand.grid(seq(4900, 6000, 25),
                   seq(4800, 5800, 25))
k <- krige.conv(ca20, locations = loc,
                krige = krige.control(cov.pars = vf$cov.pars))

png("Rplot003.png", width = 900, height = 900,
    pointsize = 14)
image(k)
dev.off()

出力

タグ：R