[R-sig-Geo] Spacetime variogram - non square cells
Edzer Pebesma
edzer.pebesma at uni-muenster.de
Sat Jun 7 19:25:37 CEST 2014
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>
More information about the R-sig-Geo
mailing list