R/outerProbability.R

In move: Visualizing and Analyzing Animal Track Data

# if (!isGeneric('outerProbability')){
setGeneric("outerProbability", function(raster, border = 0.1, ...) {standardGeneric("outerProbability")})
# }

setMethod(f = "outerProbability", 
          signature = c(raster = "RasterLayer"), 
          definition = function(raster, border, ...) {
            rowRange <- ceiling(nrow(raster) * border)
            colRange <- ceiling(ncol(raster) * border)
            innerProbability <- sum(getValuesBlock(x = raster, row = rowRange, nrows = nrow(raster) - 
              rowRange, col = colRange, ncols = ncol(raster) - colRange))
            outerProbability <- cellStats(raster, stat = sum) - innerProbability
            return(outerProbability/cellStats(raster, stat = sum))
          })
setMethod(f = "outerProbability", 
          signature = c(raster = "DBBMMStack"), 
          definition = function(raster, border, ...) {
            lapply(lapply(split(raster), as, "RasterLayer"), outerProbability, border=border)
          })