geozoning (version 1.0.0)

getPoly: getPoly

Description

getPoly

Usage

getPoly(Z, iZ, k)

Arguments

Z

zoning geometry (list of SpatialPolygons)

iZ

current zone index

k

polygon number within current zone

Value

a polygon (object of class Polygon)

Details

get the kth polygon of the current zone in zoning Z

Examples

# NOT RUN {
ZK=initialZoning(qProb=c(0.4,0.2,0.7),mapTest)
Z=ZK$resZ$zonePolygone
P1=getPoly(Z,5,1)
P2=getPoly(Z,5,2) # second polygon is a hole
plot(P1@coords,type="l")
lines(P2@coords,type="l",col="blue")
# }