Extreme_Zone: Extreme_Zone

Description Usage Arguments Details Value Examples

Description

Extreme_Zone

Usage

1
Extreme_Zone(map, Z, numZ, label.is.min = TRUE)

Arguments

map

object returned by function genMap

Z

list of zones.

numZ

number of the zone whose values will be approximated.

label.is.min

boolean value that is TRUE if the label of the zone is minimum and FALSE if the label is maximum

Details

funtion that approximates the value in a extreme zone (zone with label maximum or minimum, zones which have only one neighbour) by the solution of the Laplace's equation. The iso contours plotted on the approximate data will take the form of concentric circles as we supposed the extreme value of the zone is at the zone center (furthest point from the zone boundary.)

Value

approximated values of the values in zone (numZ).

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
 seed=35
 map=genMap(DataObj=NULL,seed=seed,krig=2,typeMod="Exp",nPointsK=500)
 ZK=initialZoning(qProb=c(0.8),map)
 Z=ZK$resZ$zonePolygone # list of zones
 lab = ZK$resZ$lab # label of zones
 plotM(map = map,Z = Z,lab = lab, byLab = FALSE)
 # zone 2 is a zone with maximum label
 numZ = 2
 Estimation = Extreme_Zone(map = map, Z = Z, numZ = numZ, label.is.min = FALSE)
 # compute the cost
 cL = Cost_By_Laplace(map = map, Z = Z, numZ = numZ, Estimation = Estimation)
 cM = Cost_By_Mean(map = map, Z = Z, numZ = numZ)
 print(cL$cost_Laplace)
 print(cM$cost_Mean)
 # zone 2 is homogeneous

geozoning documentation built on May 2, 2019, 9:43 a.m.