0

我对 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
4

1 回答 1

0

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

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

于 2014-06-24T12:41:44.733 回答