[R-sig-Geo] question about wireframe

[Edzer Pebesma]

On 02/21/2013 06:03 PM, Hodgess, Erin wrote:
> Hello!
>
> I have a krigeST object:
>
>> str(kr1)
> Formal class 'STFDF' [package "spacetime"] with 4 slots
>    ..@ data   :'data.frame':     8250 obs. of  1 variable:
>    .. ..$ var1.pred: num [1:8250] 95 95.1 95.3 95.4 95.6 ...
>    ..@ sp     :Formal class 'SpatialPoints' [package "sp"] with 3 slots
>    .. .. ..@ coords     : num [1:825, 1:2] -108 -108 -107 -106 -106 ...
>    .. .. .. ..- attr(*, "dimnames")=List of 2
>    .. .. .. .. ..$ : NULL
>    .. .. .. .. ..$ : chr [1:2] "lon" "lat"
>    .. .. ..@ bbox       : num [1:2, 1:2] -108 24 -92 36
>    .. .. .. ..- attr(*, "dimnames")=List of 2
>    .. .. .. .. ..$ : chr [1:2] "lon" "lat"
>    .. .. .. .. ..$ : chr [1:2] "min" "max"
>    .. .. ..@ proj4string:Formal class 'CRS' [package "sp"] with 1 slots
>    .. .. .. .. ..@ projargs: chr " +proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0"
>    ..@ time   :An ‘xts’ object on 2001-01-01/2010-01-01 containing:
>    Data: int [1:10, 1] 1 2 3 4 5 6 7 8 9 10
>   - attr(*, "dimnames")=List of 2
>    ..$ : NULL
>    ..$ : chr "timeIndex"
>    Indexed by objects of class: [Date] TZ: UTC
>    xts Attributes:
>   NULL
>    ..@ endTime: POSIXct[1:10], format: "2002-01-01" "2003-01-01" ...
>
>
> I would like to draw a wireframe, but I keep getting the following:
>
>> str(pred.grd)
> Formal class 'STF' [package "spacetime"] with 3 slots
>    ..@ sp     :Formal class 'SpatialPoints' [package "sp"] with 3 slots
>    .. .. ..@ coords     : num [1:825, 1:2] -108 -108 -107 -106 -106 ...
>    .. .. .. ..- attr(*, "dimnames")=List of 2
>    .. .. .. .. ..$ : NULL
>    .. .. .. .. ..$ : chr [1:2] "lon" "lat"
>    .. .. ..@ bbox       : num [1:2, 1:2] -108 24 -92 36
>    .. .. .. ..- attr(*, "dimnames")=List of 2
>    .. .. .. .. ..$ : chr [1:2] "lon" "lat"
>    .. .. .. .. ..$ : chr [1:2] "min" "max"
>    .. .. ..@ proj4string:Formal class 'CRS' [package "sp"] with 1 slots
>    .. .. .. .. ..@ projargs: chr " +proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0"
>    ..@ time   :An ‘xts’ object on 2001-01-01/2010-01-01 containing:
>    Data: int [1:10, 1] 1 2 3 4 5 6 7 8 9 10
>   - attr(*, "dimnames")=List of 2
>    ..$ : NULL
>    ..$ : chr "timeIndex"
>    Indexed by objects of class: [Date] TZ: UTC
>    xts Attributes:
>   NULL
>    ..@ endTime: POSIXct[1:10], format: "2002-01-01" "2003-01-01" ...
>> wireframe(pred.grd at coords[,1],pred.grd at coords[,2],pred.grd at time,data=kr1$var1.pred)
> Error in wireframe(pred.grd at coords[, 1], pred.grd at coords[, 2], pred.grd at time,  :
>    no slot of name "coords" for this object of class "STF"
>>

the error message is pretty clear; showClass("STF") would have told you 
where to look, as does your output above: assuming kr1 and pred.grd 
refer to the same geometry, you should have used pred.grd at sp@coords and 
index(pred.grd at time).

Simpler would be to use as(kr1, "data.frame") and specify your wireframe 
from there, using the formula var1.pred~lon+lat .

I don't understand what you expect lattice::wireframe to do when you 
pass it four vectors.

>
>
> Thanks,
> Erin
> PS this is with the latest version of gstat and spacetime (2/19)
>
>
> 	[[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
Weseler Straße 253, 48151 Münster, Germany. Phone: +49 251
8333081, Fax: +49 251 8339763  http://ifgi.uni-muenster.de
http://www.52north.org/geostatistics      e.pebesma at wwu.de