raster (version 1.3-4)

plot: Plot a Raster* object

Description

Plot (that is, make a map of) the values of a RasterLayer or RasterStack, or make a scatterplot of their values. Points, lines, and polygons can be drawn on top of a map using plot(..., add=TRUE), or functions like points, lines, polygons

Methods

x=RasterLayer plot(x, col=rev(terrain.colors(25)), maxpixels=100000, axes = TRUE, xlab="", ylab="", ...) rll{ x a RasterLayer object col A color palette, i.e. a vector of n contiguous colors such as rainbow, heat.colors, and topo.colors, on one or your own making maxpixels Maximum number of pixels used for the map ... Any argument that can be passed to image.plot and to plot } If x has a 'Color table', the default all to plot is: plot(x, maxpixels=500000, extent=NULL, interpolate=FALSE, axes=FALSE, xlab="", ylab="", ...) rll{ extent An extent object to zoom in a region interpolate Logical. Should the image be interpolated? axes Logical. Should axes be drawn? ... Graphical parameters } x=RasterStack or RasterBrick, y=numeric or missing plot(x, col=rev(terrain.colors(25)), subsample=TRUE, maxpixels=100000, axes=TRUE, xlab="", ylab="", ...) rll{ x a RasterStack or RasterBrick object y The index of a RasterLayer in a RasterStack to be plotted. Either a single number between 1 and nlayers(rstack), a vector of numbers, or missing. If y is missing, all RasterLayers in the RasterStack will be plotted. col A color palette, i.e. a vector of n contiguous colors such as rainbow, heat.colors, and topo.colors, on one or your own making maxpixels Maximum number of pixels used for the map axes see ?image xlab see ?image ylab see ?image ... Any argument that can be passed to image.plot and to plot } x=RasterLayer, y=RasterLayer Produces a scatter plot of the values of x against those of y. plot(x, y, maxpixels=100000, cex=0.1, ...) rll{ x a RasterLayer object y a RasterLayer object maxpixels Maximum number of pixels in the map display cex point size ... Any argument that can be passed to plot } x=Extent, y="ANY" Draws the (rectangular) bounding box of a raster object. plot(ext) rll{ ext a Extent object }

Details

Most of the code for this function was taken from image.plot (fields package)

See Also

plotRGB, plot3D, persp, contour, pairs, hist

Examples

# RasterLayer
r <- raster(nrows=36, ncols=18)
r <- setValues(r, runif(ncell(r)))
plot(r)
plot(extent(r))


# 2 RasterLayers
r2 <- sqrt(r)
plot(r, r2)

# RasterStack
s <- stack(r, r2, r/r)
plot(s, 2)
plot(s)

# plotting of a file with a color table:
data(drg)
e <- c(-220185.3, -219058.2, -30911.26, -29966.9)
plot(drg, extent=e)
#plot(drg, maxpixels=ncell(drg), interpolate=TRUE)