[R-sig-Geo] Spacetime variogram - non square cells

[Edzer Pebesma]

Mark, what does

diff(dat at sp@grid at cellsize)

give, as output?

On 06/06/2014 06:56 PM, Mark Payne wrote:
> 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]]
> 
> 
> 
> _______________________________________________
> 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

-------------- next part --------------
A non-text attachment was scrubbed...
Name: signature.asc
Type: application/pgp-signature
Size: 490 bytes
Desc: OpenPGP digital signature
URL: <https://stat.ethz.ch/pipermail/r-sig-geo/attachments/20140607/eb472105/attachment.bin>