findCinZ: find contour for a given quantile value, within an envelope...

Description Usage Arguments Details Value Examples

Description

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

Usage

1
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

Details

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

Value

a list with components

area

area of SpatialPolygons corresponding to contour

contourSp

SpatialPolygons corresponding to contour

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
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)

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