interCB: interCB

In geozoning: Zoning Methods for Spatial Data

Description

interCB

Usage

interCB(co, step, bd = list(x = c(0, 0, 1, 1, 0), y = c(0, 1, 1, 0, 0)),
  envel, disp = 0)

Arguments

co	contour line
step	map grid resolution
bd	map boundary
envel	envelope
disp	info level (0-no info, 1- add lines to plot)

Details

generates SpatialPolygons object corresponding to intersection of contour with boundary, must be within SpatialPolygons given in envel argument

Value

a SpatialPolygons

Examples

data(mapTest)
pG=polyToSp2(sp::Polygon(mapTest$boundary)) #SpatialPolygons corresponding to map boundary
cL=contourAuto(list(),mapTest$step,mapTest$xsize,mapTest$ysize,
   mapTest$krigGrid,c(5,7),mapTest$boundary)
ps = geozoning:::interCB(cL[[8]],mapTest$step,mapTest$boundary,pG)#envelope is the whole map
sp::plot(pG)
sp::plot(ps,col="red",add=TRUE)

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