Язык 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'

Я буду очень благодарен, если кто-то проявит любезность и скажет мне, как это исправить, предоставив мне хитрость, чтобы обойти проблему памяти / возможностей в первом случае (предпочтительно), или исправив ошибку, сгенерированную во втором случае.

ЗАРАНЕЕ БОЛЬШОЕ СПАСИБО!!!

Перель

0 ответов

Другие вопросы по тегам