Make a Red-Green-Blue plot based on three layers (in a RasterBrick or RasterStack). Three layers (sometimes referred to as "bands" because they may represent different bandwidths in the electromagnetic spectrum) are combined such that they represent the red, green and blue channel. This function can be used to make 'true (or false) color images' from Landsat and other multi-band satellite images.
# S4 method for RasterStackBrick
plotRGB(x, r=1, g=2, b=3, scale, maxpixels=500000, stretch=NULL,
ext=NULL, interpolate=FALSE, colNA='white', alpha, bgalpha, addfun=NULL, zlim=NULL,
zlimcol=NULL, axes=FALSE, xlab='', ylab='', asp=NULL, add=FALSE, ...)RasterBrick or RasterStack
integer. Index of the Red channel, between 1 and nlayers(x)
integer. Index of the Green channel, between 1 and nlayers(x)
integer. Index of the Blue channel, between 1 and nlayers(x)
integer. Maximum (possible) value in the three channels. Defaults to 255 or to the maximum value of x if that is known and larger than 255
integer > 0. Maximum number of pixels to use
character. Option to stretch the values to increase the contrast of the image: "lin" or "hist"
An Extent object to zoom in to a region of interest (see drawExtent)
logical. If TRUE, interpolate the image when drawing
color for the background (NA values)
transparency. Integer between 0 (transparent) and 255 (opaque)
Background transparency. Integer between 0 (transparent) and 255 (opaque)
Function to add additional items such as points or polygons to the plot (map). See plot
numeric vector of length 2. Range of values to plot (optional)
If NULL the values outside the range of zlim get the color of the extremes of the range. If zlimcol has any other value, the values outside the zlim range get the color of NA values (see colNA)
logical. If TRUE axes are drawn (and arguments such as main="title" will be honored)
character. Label of x-axis
character. Label of y-axis
numeric. Aspect (ratio of x and y. If NULL, and appropriate value is computed to match data for the longitude/latitude coordinate reference system, and 1 for planar coordinate reference systems
logical. If TRUE add values to current plot
graphical parameters as in plot or rasterImage
# NOT RUN {
b <- brick(system.file("external/rlogo.grd", package="raster"))
plotRGB(b)
plotRGB(b, 3, 2, 1)
plotRGB(b, 3, 2, 1, stretch='hist')
# }
Run the code above in your browser using DataLab