Transition_Zone_Near_Boundary: Transition_Zone_Near_Boundary

In geozoning: Zoning Methods for Spatial Data

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.

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.

Value

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

Examples

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

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