geozoning (version 1.0.0)

genMap: Wrapper for randKmap, generate 2D map

Description

Wrapper for randKmap, generate 2D map

Usage

genMap(DataObj = NULL, seed = 80, krig = 2, Vpsill = 5, Vrange = 0.2,
  Vnugget = 0.2, Vmean = 8, typeMod = "Exp", nPointsK = 1000,
  boundary = list(x = c(0, 0, 1, 1, 0), y = c(0, 1, 1, 0, 0)), disp = 0,
  FULL = FALSE)

Arguments

DataObj

=NULL: simulated data with seed or a data frame with real data

seed

numeric, seed used to randomly generate data points

krig

numeric, 1: kriging with vgm model, 2: inverse distance kriging

Vpsill

numeric parameter of the variogram model,

Vrange

numeric parameter of the variogram model,

Vnugget

numeric parameter of the variogram model,

Vmean

numeric parameter of the variogram model,

typeMod

type of variogram model (see vgm) "Gau", "Sph", "Exp"

nPointsK

number of generated points after kriging

boundary

list, contains x and y coordinates of map boundaries

disp

numeric,

FULL

logical, if TRUE the returned list is complete

Value

a map object as a list with components

rawData

simulated or real raw data within the boundary

step

grid step

krigData

kriged data as a SpatialPointsDataFrame

krigGrid

kriged data in form of a grid-useful for image plots.

krigN

list of neighbours of each kriged data point

krigSurfVoronoi

list of areas of Voronoi polygons in the tesselation of kriged data

modelGen

random fields model

VGMmodel

vgm model

boundary

(x,y) list of boundary points

ratio

ratio used to normalize x data

Details

Wrapper for randKmap, generates a 2D map with a Gaussian field, either by simulating data or by reading data in a data frame. Kriged data are normalized so that x-coordinates are between 0 and 1. y-coordinates are normalized with the same ratio used for x-coordinates. Kriging is either done with inverse distance interpolation, or with a variogram model. It creates a structure that contains the data and parameters necessary to perform a zoning. The structure is identical wether the data are simulated or not.

Examples

# NOT RUN {
m=genMap(seed=1,krig=2,disp=1,nPointsK = 200) #generates a map and plots data
mean(m$krigGrid) # mean of generated kriged data
# }