R通用克里金與autoKrige（）

我想在自動映射包中使用autoKrige（）函數進行通用克里金簡單應用。我有一個不規則間距的測量網格，我想在精細的空間尺度上插入它們之間。示例代碼：R通用克里金與autoKrige（）

library('automap') 

    # create an irregularly spaced grid 
    y <-x <-c(-5,-4,-2,-1,-0.5,0,0.5,1,2,4,5) 
    grid <-expand.grid(x,y) 
    names(grid) <-c('x', 'y') 

    # create some measurements, greatest in the centre, with some noise 
    vals <-apply(grid,1, function(x) {12/(0.1+sqrt(x[1]^2 + x[2]^2))+rnorm(1,2,1.5)}) 

    # get data into sp format 
    s <-SpatialPointsDataFrame(grid, data.frame(vals)) 

    # make some prediction locations and get them into sp format 
    pred <-expand.grid(seq(-5,5,by=0.5), seq(-5,5,by=0.5)) 
    pred <-cbind(pred[,1], pred[,2]) # this seems to be needed, not sure why 
    pred <-SpatialPoints(pred) 

    # try universal kriging 
    surf <-autoKrige(vals~x+y, s, new_data=pred)

這導致錯誤：

Error in gstat.formula.predict(d$formula, newdata, na.action = na.action, : 
     NROW(locs) != NROW(X): this should not occur

我試圖使NEW_DATA具有相同的行數與原始數據，並且甚至試圖使在NEW_DATA的座標與原始數據完全相同，但仍然出現此錯誤。如果我犯了一個基本錯誤，我對地質統計學技術很陌生，所以很抱歉。任何人都可以告訴我哪裏出錯了？謝謝。

來源

2015-09-07 Jonathan Denniss

的'預解碼<-cbind（PRED [，1]，預解碼[，2]）'一步可以省略對我來說，'autoKrige'工作，即使沒有它。 –

問題是您的autoKrige函數的語法錯誤。公式輸入autoKrige指定要使用線性模式，例如：

log(zinc) ~ dist

從meuse數據集。在這種情況下，您使用線性模型對log(zinc)與dist進行建模，並使用變差函數對該模型的殘差進行插值。基本上，通用克里金是具有空間相關殘差的線性迴歸。

在你的情況，你指定：

所以autoKrige（gstat實際上）將首先嚐試模式vals與x和y（多元迴歸）的線性模型，並利用變差函數插值殘差模型。但是，x和y變量不存在於SpatialPointsDataFrame中。

我認爲你想要做的只是插值空間使用變差函數模型。在這種情況下，線性模型是很簡單的，實際上只是裝配一個平均值：

vals ~ 1

其中vals平均值被確定，並且殘差使用變差函數模型內插。這實際上被稱爲普通克里格。您對autoKrige調用將是這樣的：

surf <-autoKrige(vals ~ 1, s, new_data=pred)

來源

2015-09-08 13:57:34

R通用克里金與autoKrige（）

回答

相關問題