[R-sig-Geo] Spatio-temporal geostatistics with Raster stack

[Edzer Pebesma]

library(spacetime)
EVI.st = as(EVI, "STFDF")
# This requires that the layers in your EVI have time stamps.

library(gstat)
vignette("st")

Having said this, you may well run into problems with memory or
computation time, or both -- I haven't run tests with data this large,
and would be interested to hear from you also if it doesn't work.

On 11/26/2013 01:39 PM, julio cesar oliveira wrote:
> Hi,
> 
> I have a raster stack with 250 layers, and i would like to Fitt a
> spatio-temporal variogram model.
> 
> How to convert the data for analysis?
> 
> 
>> EVI
> class       : RasterBrick
> dimensions  : 476, 336, 159936, 250  (nrow, ncol, ncell, nlayers)
> resolution  : 250, 250  (x, y)
> extent      : 607868.3, 691868.3, 8421531, 8540531  (xmin, xmax, ymin, ymax)
> coord. ref. : +proj=utm +zone=21 +south +datum=WGS84 +units=m +no_defs
> +ellps=WGS84 +towgs84=0,0,0
> data source : in memory
> names       : MOD13Q1.A//6_days_EVI, MOD13Q1.A//6_days_EVI,
> MOD13Q1.A//6_days_EVI, MOD13Q1.A//6_days_EVI, MOD13Q1.A//6_days_EVI,
> MOD13Q1.A//6_days_EVI, MOD13Q1.A//6_days_EVI, MOD13Q1.A//6_days_EVI,
> MOD13Q1.A//6_days_EVI, MOD13Q1.A//6_days_EVI, MOD13Q1.A//6_days_EVI,
> MOD13Q1.A//6_days_EVI, MOD13Q1.A//6_days_EVI, MOD13Q1.A//6_days_EVI,
> MOD13Q1.A//6_days_EVI, ...
> 
> 
> 
> Regards,
> 
> Julio
> 
> 	[[alternative HTML version deleted]]
> 
> _______________________________________________
> R-sig-Geo mailing list
> R-sig-Geo at r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
> 

-- 
Edzer Pebesma
Institute for Geoinformatics (ifgi), University of Münster
Heisenbergstraße 2, 48149 Münster, Germany. Phone: +49 251
83 33081 http://ifgi.uni-muenster.de GPG key ID 0xAC227795