geoR (version 1.3-2)

image.kriging: Image or Perspective Plot with Kriging Results

Description

Plots image or perspective plots with results of the kriging calculations.

Usage

## S3 method for class 'kriging':
image(x, locations, borders, values = obj$predict,
              coords.data, xlim, ylim, x.leg, y.leg, cex.leg = 0.75,
              vertical = FALSE, \dots)
persp.kriging(x, locations, borders, values = obj$predict, ...)

Arguments

x

an object of the class kriging, typically with the output of the functions krige.conv or ksline.

locations

an $n \times 2$ matrix with the coordinates of the prediction locations, which should define a regular grid in order to be plotted by image or persp

borders

an $n \times 2$ matrix with the coordinates defining the borders of a region inside the grid defined by locations. Elements in the argument values are assigned to locations internal to the borders and NA

values

a vector with values to be plotted. Defaults to obj$predict.

coords.data

optional. If an $n \times 2$ matrix with the data coordinates is provided, points indicating the data locations are included in the plot.

xlim, ylim

x.leg, y.leg

limits for the legend in the horizontal and vertical directions.

cex.leg

controls the size of the numerical values in the legend.

vertical

logical, defining whether the legend is drawn in horizontal (default) or vertical position.

...

further arguments to be passed to the functions image or persp.

Value

• An image or perspective plot is produced o the current graphics device. No values are returned.

References

Further information about geoR can be found at: http://www.maths.lancs.ac.uk/~ribeiro/geoR.

See Also

krige.conv and ksline for kriging calculations. Documentation for image and persp contain basic information on the plotting functions.

Examples

data(s100) 
loci <- expand.grid(seq(0,1,l=31), seq(0,1,l=31))
kc <- krige.conv(s100, loc=loci,
                 krige=krige.control(cov.pars=c(1, .25)))

op <- par(no.readonly = TRUE)
par(mfrow=c(1,2), mar=c(3,3,0,0), mgp=c(1.5, .8,0))
image(kc, loc=loci)
image(kc, loc=loci, val=kc$krige.var)

# two different ways to add legends 
image(kc, loc=loci, ylim=c(-0.2, 1), x.leg=c(0,1), y.leg=c(-0.2, -0.1))
image(kc, loc=loci, val=kc$krige.var, ylim=c(-0.2, 1))
legend.krige(y.leg=c(-0.2,-0.1), x.leg=c(0,1), val=kc$krige.var)

par(op)