geozoning (version 1.0.0)

findCinZ: find contour for a given quantile value, within an envelope and englobing current zone

Description

find contour for a given quantile value, within an envelope and englobing current zone

Usage

findCinZ(iC, Z, K, map, vRef, envel)

Arguments

iC

zone number

Z

zoning geometry (list of SpatialPolygons)

K

zoning object (such as returned by calNei function)

map

object returned by genMap function

vRef

quantile value

envel

SpatialPolygons within which the contour must be contained

Value

a list with components

area

area of SpatialPolygons corresponding to contour

contourSp

SpatialPolygons corresponding to contour

Details

withing a zoning, find contour for a given vRef quantile value, contour contains current zone and is included in envel (spatial Polygon)

Examples

# NOT RUN {
data(mapTest)
qProb=c(0.3,0.5)
criti = correctionTree(qProb,mapTest)
best = criti$zk[[2]][[1]]
Z=best$zonePolygone
plotZ(Z)
iC=4
envel=calFrame(iC,Z,best$zoneNModif)
sp::plot(envel,col="blue",add=TRUE)
vRef=quantile(mapTest$krigGrid,0.6)
resp=geozoning:::findCinZ(iC,Z,best,mapTest,vRef,envel)
sp::plot(resp$contourSp,col="red",add=TRUE)
# }