以下两个脚本将生成一个“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”的类型无效(闭包)
如果有人能告诉我如何修复它,我将非常感激,无论是在第一种情况下(最好)为我提供绕过内存/功能问题的技巧,还是修复第二种情况产生的错误。
非常感谢!
佩雷普