KrigeLegend {StatDA} | R Documentation

## 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.

### Author(s)

P.Filzmoser <P.Filzmoser@tuwien.ac.at>

### References

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

### 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)

[Package StatDA version 1.0 Index]