KrigeLegend              package:StatDA              R Documentation

_K_r_i_g_e

_D_e_s_c_r_i_p_t_i_o_n:

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

_U_s_a_g_e:

     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)

_A_r_g_u_m_e_n_t_s:

       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 

_D_e_t_a_i_l_s:

     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.

_A_u_t_h_o_r(_s):

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

_R_e_f_e_r_e_n_c_e_s:

     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.

_E_x_a_m_p_l_e_s:

     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)

