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

In geozoning: Zoning Methods for Spatial Data

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

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

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.