library(raster)
library(terra)
f <- system.file("external/test.grd", package="raster")
r <- raster(f)
s <- stack(r, r-500, r+500)
bwplot(s)
bwplot(s, violin = FALSE, box.ratio = 1)
## Modify colours
myTheme <- bwTheme(
box.rectangle = list(col = 'green', fill = 'lightgreen'),
plot.polygon = list(col = 'blue'),
plot.symbol = list(col = 'gray', cex = 0.8, alpha = 0.1)
)
## Display raw points
bwplot(s, draw.points = TRUE,
par.settings = myTheme)
if (FALSE) {
dataURL <- "https://raw.github.com/oscarperpinan/bookvis/master/data/"
##Solar irradiation data from CMSAF http://dx.doi.org/10.5676/EUM_SAF_CM/RAD_MVIRI/V001
old <- setwd(tempdir())
download.file(paste0(dataURL, "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
setwd(old)
idx <- seq(as.Date('2008-01-15'), as.Date('2008-12-15'), 'month')
SISmm <- setZ(stackSIS, idx)
names(SISmm) <- month.abb
bwplot(SISmm)
##FUN applies to z if not NULL
library(zoo)
bwplot(SISmm, FUN=as.yearqtr)
}
if (FALSE) {
##http://neo.sci.gsfc.nasa.gov/Search.html?group=64
old <- setwd(tempdir())
download.file(paste0(dataURL, "875430rgb-167772161.0.FLOAT.TIFF"),
"875430rgb-167772161.0.FLOAT.TIFF", method='wget')
pop <- raster("875430rgb-167772161.0.FLOAT.TIFF")
pop[pop==99999] <- NA
levelplot(pop, zscaleLog=10, par.settings=BTCTheme,
panel=panel.levelplot.raster, interpolate=TRUE)
##http://neo.sci.gsfc.nasa.gov/Search.html?group=20
download.file(paste0(dataURL, "241243rgb-167772161.0.TIFF"),
"241243rgb-167772161.0.TIFF", method='wget')
landClass <- raster("241243rgb-167772161.0.TIFF")
landClass[landClass==254] <- NA
s <- stack(pop, landClass)
names(s) <- c('pop', 'landClass')
bwplot(asinh(pop) ~ landClass|cut(y, 3), data = s,
layout = c(3, 1))
bwplot(asinh(pop) ~ cut(y, 5)|landClass, data = s,
scales = list(x=list(rot=45)), layout = c(4, 5),
strip = strip.custom(strip.levels = TRUE))
## Modify colours
bwplot(asinh(pop) ~ cut(y, 5)|landClass, data = s,
scales = list(x=list(rot=45)), layout = c(4, 5),
strip = strip.custom(strip.levels = TRUE),
par.settings = bwTheme(plot.polygon = list(col = 'lightgray'),
box.rectangle = list(fill = 'lightgreen')))
}
Run the code above in your browser using DataLab