KrigeLegend: Krige

In StatDA: Statistical Analysis for Environmental Data

Krige

Description

Plots Krige maps and Legend based on continuous or percentile scale.

Usage

KrigeLegend(X, Y, z, resol = 100, vario, type = "percentile",
whichcol = "gray", qutiles = c(0, 0.05, 0.25, 0.5, 0.75, 0.9, 0.95, 1),borders=NULL,
leg.xpos.min = 780000, leg.xpos.max = 8e+05, leg.ypos.min = 7760000,
leg.ypos.max = 7870000, leg.title = "mg/kg", leg.title.cex = 0.7,
leg.numb.cex = 0.7, leg.round = 2, leg.numb.xshift = 70000, leg.perc.xshift = 40000,
leg.perc.yshift = 20000, tit.xshift = 35000)

Arguments

X	X-coordinates
Y	Y-coordinates
z	values on the coordinates
resol	resolution of blocks for Kriging
vario	variogram model
type	"percentile" for percentile legend, "contin" for continous grey-scale or colour map
whichcol	type of colour scheme to use: "gray", "rainbow", "rainbow.trunc", "rainbow.inv", "terrain", "topo"
qutiles	considered quantiles if type="percentile" is used
borders	either NULL or character string with the name of the list with list elements x and y for x- and y-coordinates of map borders
leg.xpos.min	minimum value of x-position of the legend
leg.xpos.max	maximum value of x-position of the legend
leg.ypos.min	minimum value of y-position of the legend
leg.ypos.max	maximum value of y-position of the legend
leg.title	title for legend
leg.title.cex	cex for legend title
leg.numb.cex	cex for legend number
leg.round	round legend to specified digits "pretty"
leg.numb.xshift	x-shift of numbers in legend relative to leg.xpos.max
leg.perc.xshift	x-shift of "Percentile" in legend relative to leg.xpos.min
leg.perc.yshift	y-shift of numbers in legend relative to leg.ypos.max
tit.xshift	x-shift of title in legend relative to leg.xpos.max

Details

Based on a variogram model a interpolation of the spatial data is computed. The variogram has to be provided by the user and based on this model the spatial prediction is made. To distinguish between different values every predicted value is plotted in his own scale of the choosen colour.

Value

No return value, creates a plot.

Author(s)

Peter Filzmoser <[e-mail]> http://cstat.tuwien.ac.at/filz/

References

C. Reimann, P. Filzmoser, R.G. Garrett, and R. Dutter: Statistical Data Analysis Explained. Applied Environmental Statistics with R. John Wiley and Sons, Chichester, 2008.

Examples

data(chorizon)
data(kola.background)
X=chorizon[,"XCOO"]
Y=chorizon[,"YCOO"]
#el=chorizon[,"As"]
#vario.b <- variog(coords=cbind(X,Y), data=el, lambda=0, max.dist=300000)
#data(res.eyefit.As_C_m) #need the data 
#v5=variofit(vario.b,res.eyefit.As_C_m,cov.model="spherical",max.dist=300000)

plot(X,Y,frame.plot=FALSE,xaxt="n",yaxt="n",xlab="",ylab="",type="n")

# to inclrease the resolution, set e.g. resol=100
#data(bordersKola) # x and y coordinates of project boundary
#KrigeLegend(X,Y,el,resol=25,vario=v5,type="percentile",whichcol="gray",
#    qutiles=c(0,0.05,0.25,0.50,0.75,0.90,0.95,1),borders="bordersKola",
#    leg.xpos.min=7.8e5,leg.xpos.max=8.0e5,leg.ypos.min=77.6e5,leg.ypos.max=78.7e5,
#    leg.title="mg/kg", leg.title.cex=0.7, leg.numb.cex=0.7, leg.round=2,
#    leg.numb.xshift=0.7e5,leg.perc.xshift=0.4e5,leg.perc.yshift=0.2e5,tit.xshift=0.35e5)
#
#plotbg(map.col=c("gray","gray","gray","gray"),map.lwd=c(1,1,1,1),add.plot=TRUE)

StatDA documentation built on June 7, 2023, 6:26 p.m.