0

以下两个脚本将生成一个“SpatialPixelDataFrame”对象:

# FIRST
library(rgdal)
elev.grid <- readGDAL("whatever.asc")
elev.grid <- as(elev.grid, "SpatialPixelsDataFrame")

# SECOND
library(raster)
library(SDMTools)
library(adehabitat)
elev.grid <- raster("whatever.asc")
elev.grid.asc <- asc.from.raster(elev.grid) 
elev.grid.SPDF <- asc2spixdf(elev.grid.asc)

但是,第一个在将其应用于大 (15000 x 16000) 网格时超出了我的计算资源的能力,而第二个生成了一个我不能用于我的一些进一步分析的对象。例如,当我将它用于克里格目的时

x <- krige(V3~var, points, elev.grid) 

我得到以下信息:

model.frame.default(terms(formula), as(data, "data.frame"), na.action = na.fail) 中的错误:变量“v ar”的类型无效(闭包)

如果有人能告诉我如何修复它,我将非常感激,无论是在第一种情况下(最好)为我提供绕过内存/功能问题的技巧,还是修复第二种情况产生的错误。

非常感谢!

佩雷普

4

0 回答 0