R/plot.bsims_landscape.R

In psolymos/bSims: Agent-based Bird Point Count Simulator

plot.bsims_landscape <-
function(x, col_H, col_E, col_R, xlim=NULL, ylim=NULL, ...) {

  .bsims_theme <- getOption("bsims_theme")
  if (missing(col_H))
    col_H <- .bsims_theme$col_H
  if (missing(col_E))
    col_E <- .bsims_theme$col_E
  if (missing(col_R))
    col_R <- .bsims_theme$col_R

  A <- diff(x$strata) * diff(range(x$strata))
  A <- c(h=A[1]+A[5], e=A[2]+A[4], r=A[3])
  names(A) <- c("H", "E", "R")
  plot(0, type="n",
    xlim=if (is.null(xlim)) range(x$box[,"x"]) else xlim,
    ylim=if (is.null(ylim)) range(x$box[,"y"]) else ylim,
    xlab="", ylab="", axes=FALSE, asp=1, ...)
  if (A[1] > 0)
    polygon(x$box, col=col_H, border=NA)
  if (A[2] > 0)
    polygon(x$strata[c("he", "he", "eh", "eh")],
      x$strata[c("+h", "h+", "h+", "+h")],
      col=col_E, border=NA)
  if (A[3] > 0)
    polygon(x$strata[c("er", "er", "re", "re")],
      x$strata[c("+h", "h+", "h+", "+h")],
      col=col_R, border=NA)
  invisible(x)
}