geozoning (version 1.0.0)

Extreme_Zone: Extreme_Zone

Description

Extreme_Zone

Usage

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

Value

approximated values of the values in zone (numZ).

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

Examples

# NOT RUN {
 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
# }