rasterblaster: Plots a kriging result in color!

Description Usage Arguments Note Author(s) Examples

Description

Pass in a kriging result, and watch the world explode into technicolor, just like for Dorothy.

Usage

1
2
3
rasterblaster(krigePlot, type = "predict", z = NULL, ramp = "rainbow",
  invert = T, alpha = 0.8, dark = 0.85, contour = T, xlab = "X coord",
  ylab = "Y coord", nlevels = 10, ...)

Arguments

krigePlot

The kriging object to use, returned from krige.conv

type

The type of plot to produce. Specifying type="predict" will produce a map of kriging predictions, specifying type="se" will produce a map of kriging standard errors, and specifying type="lower95" or type="upper95" will produce a map of lower or upper 95 percent CIs of kriging predictions. Defaults to "predict".

z

Alternately, a vector of quantities can be used, with the same length as the prediction locations. Defaults to NULL.

ramp

The color ramp to use for plotting. Currently, allowed ramps are "grey", "red", "rainbow", "stoplight", and "heat". Defaults to "rainbow".

invert

Whether to invert the color ramp. Try it both ways, see if it looks better.

alpha

Color transparency to use for plotting. This ranges from 0 (completely transparent) to 1 (true color). Defaults to 0.8.

dark

Color saturation to use for plotting. This ranges from 0 (black) to 1 (true color). Defaults to 0.85.

contour

Whether to add contour lines. Defaults to TRUE.

nlevels

The number of contour levels to use. Defaults to 10.

...

Additional plotting parameters

Note

Can also produce a plot from an object returned from pois.krige.bayes or binom.krige.bayes, but adding contour lines is awfully slow.

Author(s)

Matt Tyers

Examples

1
2
3
4
5
6
7
8
9
data(krigePlot,locations)

rasterblaster(krigePlot)
rasterblaster(krigePlot,invert=T)
rasterblaster(krigePlot,alpha=.4)
rasterblaster(krigePlot,alpha=.4,dark=.7)
rasterblaster(krigePlot,ramp="grey")

rasterblaster(krigePlot,type="se",ramp="red")

mbtyers/rasterblaster documentation built on May 22, 2019, 12:58 p.m.