[R-sig-Geo] spplot for raster not drawn

[Els Ducheyne]

dear Oscar
thanks, very helpful

Best

On 13 Sep 2011, at 13:01, Oscar Perpiñan Lamigueiro wrote:

> Hi,
> 
> Another approach is to use the layout argument of levelplot. 
> 
> As an example, I use the data from Cressie and Wikle, following the
> superb Edzer Pebesma's notes in demo/CressieWikle.R in the 
> spacetime package.
> 
> Best, 
> 
> Oscar.
> 
> library(raster)
> library(rasterVis)
> library(zoo)
> 
> url <- "ftp://ftp.wiley.com/public/sci_tech_med/spatio_temporal_data/"
> 
> sst.dat = read.table(paste(url, "SST011970_032003.dat", sep=''), header
> 	= FALSE) 
> sst.ll = read.table(paste(url, "SSTlonlat.dat", sep=''),
> 	header = FALSE)
> 
> sp <- SpatialPointsDataFrame(sst.ll, sst.dat)
> gridded(sp) <- TRUE
> proj4string(sp) = "+proj=longlat +datum=WGS84"
> s <- brick(sp)
> 
> idx <- seq(as.Date('1970-01-01'), as.Date('2003-03-01'), by='month')
> idx <- as.yearmon(idx)
> s <- setZ(s, idx)
> layerNames(s) <- as.character(idx)
> 
> ##One year per page with layout=c(3,4)
> trellis.device('pdf', file='cressie_level.pdf')
> levelplot(s, layout=c(3, 4))
> dev.off()
> 
> ## layout inside histogram works correctly 
> ## only with rasterVis 0.10-5 (at R-Forge) 
> trellis.device('pdf', file='cressie_hist.pdf')
> histogram(s, layout=c(3, 4))  
> dev.off()
> 
> 
> 
> -------------------
> Oscar Perpiñán Lamigueiro
> Dpto. Ingeniería Eléctrica
> EUITI-UPM
> 
> http://procomun.wordpress.com
> 
> -----------------------------------------
>> Use print on your levelplot command. This is in the R FAQ (7.22).
>> 
>> Paul
>> 
>> On 09/13/2011 09:10 AM, Els Ducheyne wrote:
>>> dear list
>>> I want to display a time series of data using either spplot or
>>> levelplot to verify consistency in the raw images prior to
>>> processing (visual quality check) to do this I am use spplot (or
>>> levelplot). This works up to 16 images, however if I want to see
>>> more I have to split the inputfiles and put them in  a loop. When
>>> the loop is executed, I can see the quartz window but the content is
>>> not shown. Does anyone have an idea how to solve this ? See below
>>> for information on OS and versions + code
>>> 
>>> Thanks for help
>>> Els
>>> 
>>> 
>>> require(raster) #raster package = 1.9-5
>>> require(rasterVis) #rasterVis = 0.10-4
>>> setwd(outputdir)
>>> #get all the inputdata
>>> par(ask=T)
>>> inputlist <- list.files(pattern="LSTDay")
>>> datadisplay <- seq(1,length(inputlist),16)
>>> datadisplay <- c(datadisplay,length(inputlist))
>>> n <- length(datadisplay)-1
>>> for (i in 1:n){
>>> 	input <- inputlist[datadisplay[i]:datadisplay[i+1]-1]
>>> 	images <- stack(inputlist)
>>> 	NAvalue(images)<-0
>>> 	#spplot(images,col.regions=terrain.colors(80)) # quartz is
>>> shown but remains empty
>>> levelplot(images,col.regions=terrain.colors(80)) # quartz is shown
>>> but remains empty hist(images,xlim=c(220*50,340*50),n=20)
>>> #histograms are displayed in quartz }
>>> 
>>> 
>>> $platform
>>> [1] "i386-apple-darwin9.8.0"
>>> 
>>> $arch
>>> [1] "i386"
>>> 
>>> $os
>>> [1] "darwin9.8.0"
>>> 
>>> $system
>>> [1] "i386, darwin9.8.0"
>>> 
>>> $status
>>> [1] ""
>>> 
>>> $major
>>> [1] "2"
>>> 
>>> $minor
>>> [1] "13.1"
>>> 
>>> $year
>>> [1] "2011"
>>> 
>>> $month
>>> [1] "07"
>>> 
>>> $day
>>> [1] "08"
>>> 
>>> $`svn rev`
>>> [1] "56322"
>>> 
>>> $language
>>> [1] "R"
>>> 
>>> $version.string
>>> [1] "R version 2.13.1 (2011-07-08)"
>>> _______________________________________________
>>> R-sig-Geo mailing list
>>> R-sig-Geo at r-project.org
>>> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>> 
>> 
>> -- 
>> Paul Hiemstra, Ph.D.
>> Global Climate Division
>> Royal Netherlands Meteorological Institute (KNMI)
>> Wilhelminalaan 10 | 3732 GK | De Bilt | Kamer B 3.39
>> P.O. Box 201 | 3730 AE | De Bilt
>> tel: +31 30 2206 494
>> 
>> http://intamap.geo.uu.nl/~paul
>> http://nl.linkedin.com/pub/paul-hiemstra/20/30b/770