R/burnFromAge.R

In tati-micheletti/usefulFuns: Useful functions for my modules and packages

burnFromAge <- function(ageMap,
                        oldBurnTime,
                        newBurnName,
                        oldBurnName){
  
  newBurn <- oldBurn <- ageMap
  newAge <- oldAge <- getValues(ageMap)
  
  newAge[newAge > oldBurnTime] <- -1
  newAge[!is.na(newAge) & newAge >= 0] <- 1
  newBurn <- setValues(newBurn, newAge)
  newBurn[newBurn == -1] <- 0
  names(newBurn) <- newBurnName
  
  oldAge[oldAge <= oldBurnTime | oldAge > (20 + oldBurnTime)] <- -1 # Educated guess that we can only id fires for 60 years. Then it is not old fire anymore
  oldAge[!is.na(oldAge) & oldAge > 0] <- 1
  oldBurn <- setValues(oldBurn, oldAge)
  oldBurn[oldBurn == -1] <- 0
  names(oldBurn) <- oldBurnName
  
  return(list(newBurn, oldBurn))
}