autoKrige - eval（expr，envir，enclos）中的错误：找不到对象'....'

Question

我对 Automap 包有疑问。

我已经尝试了普通的克里金法，我的数据没有问题并且工作正常。但是当我尝试使用 Universal Kriging 时，收到一条我无法理解的错误消息。

我认为可能很简单，但我不知道为什么该函数没有找到“silt_clay”。

谢谢。

在我的代码上方：

library(automap)
data_interpol<-read.table("merge_guaiba_ceco2.txt", header=T)
head(data_interpol)
   long     lat depth silt_clay    avgran pcsand    om
1 471699 6680184   1.9    51.626 0.6683944 48.313 11.50
2 473777 6679430   1.8    57.810 0.6669483 42.186  8.80
3 471931 6679087   2.4    54.625 0.6718675 45.370  9.30
4 471337 6678456   1.9     2.887 0.4128039 97.113  0.00
5 472936 6678421   2.5    49.136 0.6344626 50.864  9.90
6 473945 6678482   1.9    99.614 0.9984355  0.386  9.45

coordinates(data_interpol)<- ~long+lat
kriging_om = autoKrige(om~ 1, data_interpol, grid4)# works fine
kriging_om_Universal= autoKrige(om~ silt_clay + avgran + pcsand, data_interpol, grid4)
Error in eval(expr, envir, enclos) : object 'silt_clay' not found # my error

####code to obtain grid4:
    grid <- raster( )
    grid <- raster(ncol=1544, nrow=1056, xmn=468318.3, xmx=499998.3, ymn=6634918, ymx=6681238) projection (grid) <- '+proj=utm +zone=22 +south +ellps=aust_SA +units=m +no_defs ‘
    res(grid) <- 30
    grid4 <- as(grid, ‘SpatialPoints’)

编辑讨论错误。更改网格后错误仍然存​​在：

d_margin_xyz 是我的新网格。

str(d_margin_xyz)
Formal class 'SpatialPointsDataFrame' [package "sp"] with 5 slots
  ..@ data       :'data.frame': 1628352 obs. of  1 variable:
  .. ..$ DIST_MARGIN: num [1:1628352] 0 0 0 0 0 0 0 0 0 0 ...
  ..@ coords.nrs : int [1:2] 1 2
  ..@ coords     : num [1:1628352, 1:2] 468333 468363 468393 468423 468453 ...
  .. ..- attr(*, "dimnames")=List of 2
  .. .. ..$ : NULL
  .. .. ..$ : chr [1:2] "x" "y"
  ..@ bbox       : num [1:2, 1:2] 468333 6634963 499983 6681193
  .. ..- attr(*, "dimnames")=List of 2
  .. .. ..$ : chr [1:2] "x" "y"
  .. .. ..$ : chr [1:2] "min" "max"
  ..@ proj4string:Formal class 'CRS' [package "sp"] with 1 slots
  .. .. ..@ projargs: chr NA

使用包含通用克里金法预测器的新网格再次尝试，错误再次出现。有什么线索吗？

kriging_om_Universal= autoKrige(om~ DIST_MARGIN, data_interpol, d_margin_xyz)
Error in eval(expr, envir, enclos) : object 'DIST_MARGIN' not found

Answer 1

问题在于您的预测对象grid4 只是一个没有属性的SpatialGrid 对象。为了能够使用通用克里金法进行预测，您需要预测位置上所有预测变量（即silt_clay、avgran 等）的值，即作为grid4 的属性。您的错误是由silt_clay 不是属性引起的。

要解决此问题，您需要将属性添加到grid4。为此，您需要有某种区域覆盖silt_clay 和其他预测变量的地图。您可以将此地图转换为具有属性的SpatialGrid 对象，或者从该地图中提取grid4 中位置的值。