Язык R, проблемы с SpatialPixelsDataFrame
Следующие два сценария сгенерируют объект "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(термины (формула), as(data, "data.frame"), na.action = na.fail): недопустимый тип (закрытие) для переменной 'v ar'
Я буду очень благодарен, если кто-то проявит любезность и скажет мне, как это исправить, предоставив мне хитрость, чтобы обойти проблему памяти / возможностей в первом случае (предпочтительно), или исправив ошибку, сгенерированную во втором случае.
ЗАРАНЕЕ БОЛЬШОЕ СПАСИБО!!!
Перель