geozoning (version 1.0.0)

Transition_Zone_Near_Boundary: Transition_Zone_Near_Boundary

Description

Transition_Zone_Near_Boundary

Usage

Transition_Zone_Near_Boundary(map, Z, numZ)

Arguments

map

object returned by function genMap

Z

list of zones.

numZ

number of the zone whose values will be approximated.

Value

approximated values of the zone (numZ) given as parameter.

Details

funtion that approximates the value in a transition zone (which has commun boundary with the map) by the solution of the Laplace's equation. The numerical resolution of the Laplace's equation will be based on the discretisation of the data on the grid (map$krigGrid). The domaine of study is a transition zone which have a commun border with the map.

Examples

# NOT RUN {
seed=21
map=genMap(DataObj=NULL,seed=seed,krig=2,typeMod="Exp",nPointsK=800)
ZK=initialZoning(qProb=c(0.55,0.85),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 3 is a transition zone that has commun boundary with the map
numZ = 3
Estimation = Transition_Zone_Near_Boundary(map = map, Z = Z, numZ = numZ)
# 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 3 is a zone with gradient
# }