rasterVis (version 0.40)

xyplot-methods: xyplot for Raster objects

Description

Scatter plots of space-time Raster objects for directions defined by xyLayer

Usage

"xyplot"(x, data=NULL, dirXY=y, stat='mean', xlab='Time', ylab='', digits=0, par.settings=rasterTheme(),...)

Arguments

x

A RasterStackBrick object whose z slot is not NULL.

data

Not used.

dirXY

A direction as a function of the coordinates (see xyLayer).

stat

a function to be applied to summarize the values by zone. See zonal for details.

xlab, ylab

Labels of the axis.

par.settings

Customization of lattice. See xyplot for details.

digits

An integer, number of digits for zonal.

...

Additional arguments for the xyplot function.

See Also

zonal

Examples

## Not run: 
# ##Solar irradiation data from CMSAF http://dx.doi.org/10.5676/EUM_SAF_CM/RAD_MVIRI/V001
# old <- setwd(tempdir())
# download.file('https://raw.github.com/oscarperpinan/spacetime-vis/master/data/SISmm2008_CMSAF.zip',
#    'SISmm2008_CMSAF.zip', method='wget')
# unzip('SISmm2008_CMSAF.zip')
# 
# listFich <- dir(pattern='\.nc')
# stackSIS <- stack(listFich)
# stackSIS <- stackSIS*24 ##from irradiance (W/m2) to irradiation Wh/m2
# 
# idx <- seq(as.Date('2008-01-15'), as.Date('2008-12-15'), 'month')
# 
# SISmm <- setZ(stackSIS, idx)
# names(SISmm) <- month.abb
# 
# xyplot(SISmm)
# ## Formula interface 
# xyplot(Jan ~ Jul, data=SISmm)
# ## Different scatterplots for each latitude zone
# xyplot(Jan ~ Dec|cut(y, 4), data=SISmm)
# ## More sophisticated bands can be defined using the dirXY argument 
# xyplot(Jan ~ Dec|cut(dirXY, 4), dirXY=x^2 + y^2, data=SISmm)
# ## End(Not run)

## Not run: 
# 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)
# 
# spSST <- SpatialPointsDataFrame(sst.ll, sst.dat)
# gridded(spSST) <- TRUE
# proj4string(spSST) = "+proj=longlat +datum=WGS84"
# SST <- brick(spSST)
# 
# idx <- seq(as.Date('1970-01-01'), as.Date('2003-03-01'), by='month')
# idx <- as.yearmon(idx)
# SST <- setZ(SST, idx)
# names(SST) <- as.character(idx)
# 
# xyplot(SST)
# ## End(Not run)