interCB: interCB

Description Usage Arguments Details Value Examples

Description

interCB

Usage

1
2
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

1
2
3
4
5
6
7
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)

hazaeljones/geozoning documentation built on May 30, 2019, 3:06 p.m.