[R-sig-Geo] Running spacetime eof on large raster stack

[Shiva Khanal]

Hi,

I am trying to run the spatial empirical orthogonal function on a raster stack using eof function in spacetime 1.2-1 package.

My raster stack is

> stk
class       : RasterStack
dimensions  : 2423, 2470, 5984810, 690  (nrow, ncol, ncell, nlayers)
resolution  : 499.9339, 499.9718  (x, y)
extent      : -85163.44, 1149673, 5656775, 6868207  (xmin, xmax, ymin, ymax)
coord. ref. : NA
names       : X2001.01.01, X2001.01.09, X2001.01.17, X2001.01.25, X2001.02.02, X2001.02.10, X2001.02.18, X2001.02.26, X2001.03.06, X2001.03.14, X2001.03.22, X2001.03.30, X2001.04.07, X2001.04.15, X2001.04.23, ...
min values  :           0,           0,           0,           0,           0,           0,           0,           0,           0,           0,           0,           0,           0,           0,           0, ...
max values  :         100,         100,         100,         100,         100,         100,         100,         100,         100,         100,         100,         100,         100,         100,         100, ...

>object.size(stk)
8699416 bytes

The size of raster stack is ~8 mb but when I try to convert to STFDF it become several GBs until I receive error
about memory.

My system is
sysname        release        version        machine
"Windows"     ">= 8 x64"    "build 9200"      "x86-64"

>memory.limit()
[1] 32659

My R version is
"R version 3.5.0 (2018-04-23)"

I would appreciate any advise on possible solutions.

Thank you.

Shiva


	[[alternative HTML version deleted]]