Creates a grid of colored or grayscale rectangles with colors
corresponding to the values in z
. This can be used to display
threedimensional or spatial data aka images.
This is a generic function.
The functions heat.colors
, terrain.colors
and topo.colors
create heatspectrum (red to white) and
topographical color schemes suitable for displaying ordered data, with
n
giving the number of colors desired.
1 2 3 4 5 6 
x, y 
locations of grid lines at which the values in 
z 
a numeric or logical matrix containing the values to be plotted
( 
zlim 
the minimum and maximum 
xlim, ylim 
ranges for the plotted 
col 
a list of colors such as that generated by

add 
logical; if 
xaxs, yaxs 
style of x and y axis. The default 
xlab, ylab 
each a character string giving the labels for the x and
y axis. Default to the ‘call names’ of 
breaks 
a set of finite numeric breakpoints for the colours: must have one more breakpoint than colour and be in increasing order. Unsorted vectors will be sorted, with a warning. 
oldstyle 
logical. If true the midpoints of the colour intervals
are equally spaced, and 
useRaster 
logical; if 
... 
graphical parameters for 
The length of x
should be equal to the nrow(z)+1
or
nrow(z)
. In the first case x
specifies the boundaries
between the cells: in the second case x
specifies the midpoints
of the cells. Similar reasoning applies to y
. It probably
only makes sense to specify the midpoints of an equallyspaced
grid. If you specify just one row or column and a lengthone x
or y
, the whole user area in the corresponding direction is
filled. For logarithmic x
or y
axes the boundaries between
cells must be specified.
Rectangles corresponding to missing values are not plotted (and so are
transparent and (unless add = TRUE
) the default background
painted in par("bg")
will show through and if that is
transparent, the canvas colour will be seen).
If breaks
is specified then zlim
is unused and the
algorithm used follows cut
, so intervals are closed on
the right and open on the left except for the lowest interval which is
closed at both ends.
The axes (where plotted) make use of the classes of xlim
and
ylim
(and hence by default the classes of x
and
y
): this will mean that for example dates are labelled as
such.
Notice that image
interprets the z
matrix as a table of
f(x[i], y[j])
values, so that the x axis corresponds to row
number and the y axis to column number, with column 1 at the bottom,
i.e. a 90 degree counterclockwise rotation of the conventional
printed layout of a matrix.
Images for large z
on a regular grid are rendered more
efficiently with useRaster = TRUE
and can prevent rare
antialiasing artifacts, but may not be supported by all graphics
devices. Some devices (such as postscript
and X11(type =
"Xlib")
) which do not support semitransparent colours may emit
missing values as white rather than transparent, and there may be
limitations on the size of a raster image. (Problems with the
rendering of raster images have been reported by users of
windows()
devices under Remote Desktop, at least under its
default settings.)
The graphics files in PDF and PostScript can be much smaller under this option.
If useRaster
is not specified, raster images are used when the
getOption("preferRaster")
is true, the grid is regular
and either dev.capabilities("rasterImage")$rasterImage
is "yes"
or it is "nonmissing"
and there are no missing
values.
Originally based on a function by Thomas Lumley.
filled.contour
or heatmap
which can
look nicer (but are less modular),
contour
;
The lattice equivalent of image
is
levelplot
.
heat.colors
, topo.colors
,
terrain.colors
, rainbow
,
hsv
, par
.
dev.capabilities
to see if useRaster = TRUE
is
supported on the current device.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22  require(grDevices) # for colours
x < y < seq(4*pi, 4*pi, len = 27)
r < sqrt(outer(x^2, y^2, "+"))
image(z = z < cos(r^2)*exp(r/6), col = gray((0:32)/32))
image(z, axes = FALSE, main = "Math can be beautiful ...",
xlab = expression(cos(r^2) * e^{r/6}))
contour(z, add = TRUE, drawlabels = FALSE)
# Volcano data visualized as matrix. Need to transpose and flip
# matrix horizontally.
image(t(volcano)[ncol(volcano):1,])
# A prettier display of the volcano
x < 10*(1:nrow(volcano))
y < 10*(1:ncol(volcano))
image(x, y, volcano, col = terrain.colors(100), axes = FALSE)
contour(x, y, volcano, levels = seq(90, 200, by = 5),
add = TRUE, col = "peru")
axis(1, at = seq(100, 800, by = 100))
axis(2, at = seq(100, 600, by = 100))
box()
title(main = "Maunga Whau Volcano", font.main = 4)

Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.