Hi,

I have a space-time array derived from (patchy) satellite data that I am
trying to do space-time kriging on. However, I'm not having much success
with trying to calculate the variogram. In particular, I get the following
error:

>vg <- variogramST(logchl~1,location=dat)

|
|   0%
Error: gstat: value not allowed for: cannot deal with non-square cells
>

As far as I can see, my cells are "square" - I have a regular lat-lon grid,
and regularly spaced time intervals. Or have I mis-understood the meaning
of this error message?

A summary of my STFDF is included below, and package versions as well.

Best wishes,

Mark

> summary(dat)
Object of class STFDF
 with Dimensions (s, t, attr): (11520, 301, 1)
[[Spatial:]]
Object of class SpatialPixels
Coordinates:
         min       max
s1 -4.979167 -1.020833
s2 43.020835 47.979168
Is projected: FALSE
proj4string : [+proj=longlat +datum=WGS84]
Number of points: 11520
Grid attributes:
   cellcentre.offset   cellsize cells.dim
s1         -4.979167 0.04166666        96
s2         43.020835 0.04166666       120
[[Temporal:]]
     Index              timeIndex
 Min.   :2002-05-28   Min.   :  1
 1st Qu.:2002-08-11   1st Qu.: 76
 Median :2002-10-25   Median :151
 Mean   :2002-10-25   Mean   :151
 3rd Qu.:2003-01-08   3rd Qu.:226
 Max.   :2003-03-24   Max.   :301
[[Data attributes:]]
     logchl
 Min.   :-6.8
 1st Qu.:-1.4
 Median :-0.9
 Mean   :-0.8
 3rd Qu.:-0.3
 Max.   : 4.0
 NA's   :2461741
> packageVersion("gstat")
[1] ‘1.0.19’
> packageVersion("spacetime")
[1] ‘1.1.0’
>

	[[alternative HTML version deleted]]

