R/resolucao.R

In thaisdoria/habitaR: Tools for mapping the species’ area of habitat (AOH) from geospatial data

resolucao <- function(x, y, type = 0, lbl = 'blank'){
  if(class(y) %in% c('RasterBrick', 'RasterStack')){
    y = y[[1]]
  }
  #RasterStack or #RasterBrick
  if (class(x) %in% c('RasterBrick', 'RasterStack')) {
    new.x <- list()
    for (i in 1:nlayers(x)) {
      new.x[[i]] <- x[[i]]
      if(lbl %in% 'climSuit' | type == 1){
        new.x[[i]] <- resample(new.x[[i]], y)
      }
      if(lbl %in% c('lc', 'alt')| type == 2){
        factor <- trunc(xres(y) / xres(new.x[[i]]))
        if(factor > 1){
          new.x[[i]] <- aggregate(new.x[[i]], fact = factor, fun = sum)
          new.x[[i]] <- new.x[[i]] / (factor^2)
        }
        new.x[[i]] <- resample(new.x[[i]], y, method = 'ngb')
      }
    }
    new.x <- stack(new.x)
    return(new.x)
  } else {
    #RasterLayer
    if(lbl %in% 'climSuit'| type == 1){
      x <- resample(x, y)
    }
    if(lbl %in% c('lc', 'alt')| type == 2){
      factor <- trunc(xres(y) / xres(x))
      if(factor > 1){
        x <- aggregate(x, fact = factor, fun = sum)
        x <- x / (factor^2)
      }
      x <- resample(x, y, method = 'ngb')
    }
    return(x)
  }
}