Description Usage Arguments Details Value Author(s) See Also Examples
Plots several layer of information, overlaying different kind of information. This funtion make it easy to plot several shapefiles/data pair information over one single plot.
1 |
... |
a list of lists, where each element of the list describes a layer of
information. The list must contain a |
This function can also be embeded into filled.contour
function.
Example two produces the following graphic.
Function does not return any value.
A.M. Sajo-Castelli
vetools
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 | library(maptools)
library(vetools)
# Example 1 ####
ZU <- get.shape.state("ZU")
border <- list(FUN = plot, ZU, asp = 1, lwd = 2,
border = "blue", col = NA, add = TRUE)
r <- get.shape.range(ZU)
x <- seq(r[1], r[2], length.out = nrow(volcano))
y <- seq(r[3], r[4], length.out = ncol(volcano))
image(x, y, volcano, col = heat.colors(100),
axes = FALSE, xlab = NA, ylab = NA, asp = 1)
plotLayers(border)
plotArrow(ZU, cex = 0.666, offset.arrow = c(0.1, 0))
title(main = "Raster image combined with plotLayers")
# Example 2 ####
long=c(-72.5, -71.5, -71.0); lat=c(9.5, 8.75, 10.5);
mm = 1.5 * c(2.5, 3.8, 4.2)
data <- list(FUN = points, x = long, y = lat, pch = 21,
bg = rgb(0, 1, 0, 0.666), col = "blue",
cex = mm)
filled.contour(x, y, volcano, xlab = "Longitude",
ylab = "Latitude", asp = 1,
color.palette = heat.colors,
plot.axis = { plotLayers(border, data) },
main = "plotLayers & filled.contour example")
# Example 3 ####
pts <- cbind(r[1] + 2 * runif(10), r[3] + 3 * runif(10))
sts <- runif(10)
stations <- list(FUN = plot, x = pts[, 1], y = pts[, 2],
asp = 1, pch = 21, col = rgb(sts, 0, 0),
bg = 'white', cex = 2, lwd = 2,
xlim = r[1:2], ylim = r[3:4], axes = FALSE,
xlab = NA, ylab = NA)
labs <- list(FUN=text, x=pts[,1], y=pts[,2], labels=1:10,
cex=0.7)
type = 1 + round(2 * sts)
LABELS = c('optimal', 'normal', 'critical')
status <- list(FUN = text, x = pts[, 1], y = pts[, 2],
labels = LABELS[type], cex = 0.7,
pos = 4, col = rgb(sts, 0, 0))
arrow <- list(FUN = plotArrow, shape = ZU, cex = 0.7)
plotLayers(stations, border, labs, status, arrow)
title(main = "plotLayers example", sub = "Zulia state")
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.