以下两个脚本将生成一个" 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)网格时,除了我的计算资源的能力之外,第二个生成了一个我无法用于进一步分析的对象。例如,当我将它用于krige目的时
x <- krige(V3~var, points, elev.grid)
我得到以下内容:
model.frame.default中的错误(条款(公式),如(数据,&#34; data.frame&#34;),na.action = na.fail):无效的类型(闭包)变量&#39; v ar&#39;
如果有人能够告诉我如何修复它,是否在第一种情况下(最好是)为我提供绕过内存/能力问题的技巧,或者修复第二种情况产生的错误,我将非常感激。
感谢很多!!!“
perep