Я пытаюсь использовать функцию autoKrige() в пакете automap для простого применения универсального кригинга. У меня есть неравномерная сетка измерений, и я хочу интерполировать между ними в точном пространственном масштабе. Пример кода:
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 точно такими же, как в исходных данных, но я все еще получаю эту ошибку. Я новичок в геостатистических методах, так что извините, если я делаю основную ошибку. Может ли кто-нибудь посоветовать, где я ошибаюсь? Спасибо.
pred <-cbind(pred[,1], pred[,2])
для меня можно опустить,autoKrige
работает и без него. - person Paul Hiemstra   schedule 08.09.2015