R/get_fires.R

Defines functions get_fires

Documented in get_fires

get_fires <- function(level = NULL, panel = FALSE){
  # The code below help to avoid the note "no visible binding for global variable [variable name]"
  CD_GEOCMU <-  cod.state <- sigla.state <- year <- country <- fires <- v1 <- NULL


  temp <- tempfile()
  if(is.null(level) == TRUE){
    message('Pleace, wait for the data to download!')
    dt <- openxlsx::read.xlsx('https://trello-attachments.s3.amazonaws.com/5ea83f462064047eed09f846/5ea83fc7d5707865983b4f14/16378a8eeeb8d1651493fe4e9606cf51/INCENDIOS_2001_2019.xlsx')
  }else{
    if(level == 'municipality'){
      message('Pleace, wait for the data to download!')
      dt <- openxlsx::read.xlsx('https://trello-attachments.s3.amazonaws.com/5ea83f462064047eed09f846/5ea83fc7d5707865983b4f14/16378a8eeeb8d1651493fe4e9606cf51/INCENDIOS_2001_2019.xlsx')
    }
    if(level == 'state'){
      message('Pleace, wait for the data to download!')
      dt <- openxlsx::read.xlsx('https://trello-attachments.s3.amazonaws.com/5ea83f462064047eed09f846/5ea83fc7d5707865983b4f14/16378a8eeeb8d1651493fe4e9606cf51/INCENDIOS_2001_2019.xlsx')
      dt <- dplyr::mutate(dt, cod.state = substr(CD_GEOCMU, 1, 2),
                          sigla.state = ifelse(cod.state == '11', 'RO' , 'DF'),
                          sigla.state = ifelse(cod.state == '12', 'AC' , sigla.state), sigla.state = ifelse(cod.state == '13', 'AM' , sigla.state),
                          sigla.state = ifelse(cod.state == '14', 'RR' , sigla.state), sigla.state = ifelse(cod.state == '15', 'PA' , sigla.state),
                          sigla.state = ifelse(cod.state == '16', 'AP' , sigla.state), sigla.state = ifelse(cod.state == '17', 'TO' , sigla.state),
                          sigla.state = ifelse(cod.state == '21', 'MA' , sigla.state), sigla.state = ifelse(cod.state == '22', 'PI' , sigla.state),
                          sigla.state = ifelse(cod.state == '23', 'CE' , sigla.state), sigla.state = ifelse(cod.state == '24', 'RN' , sigla.state),
                          sigla.state = ifelse(cod.state == '25', 'PB' , sigla.state), sigla.state = ifelse(cod.state == '26', 'PE' , sigla.state),
                          sigla.state = ifelse(cod.state == '27', 'AL' , sigla.state), sigla.state = ifelse(cod.state == '28', 'SE' , sigla.state),
                          sigla.state = ifelse(cod.state == '29', 'BA' , sigla.state), sigla.state = ifelse(cod.state == '31', 'MG' , sigla.state),
                          sigla.state = ifelse(cod.state == '32', 'ES' , sigla.state), sigla.state = ifelse(cod.state == '33', 'RJ' , sigla.state),
                          sigla.state = ifelse(cod.state == '35', 'SP' , sigla.state), sigla.state = ifelse(cod.state == '41', 'PR' , sigla.state),
                          sigla.state = ifelse(cod.state == '42', 'SC' , sigla.state), sigla.state = ifelse(cod.state == '43', 'RS' , sigla.state),
                          sigla.state = ifelse(cod.state == '50', 'MS' , sigla.state), sigla.state = ifelse(cod.state == '51', 'MT' , sigla.state),
                          sigla.state = ifelse(cod.state == '52', 'GO' , sigla.state))
      dt <- dplyr::group_by(dt, sigla.state)
      dt<- dplyr::summarise_at(dt, 2:20, sum)
    }
    if(level == 'country'){
      message('Pleace, wait for the data to download!')
      dt <- openxlsx::read.xlsx('https://trello-attachments.s3.amazonaws.com/5ea83f462064047eed09f846/5ea83fc7d5707865983b4f14/16378a8eeeb8d1651493fe4e9606cf51/INCENDIOS_2001_2019.xlsx')
      dt<- dplyr::mutate(dt, country = 'Brazil')
      dt<-dplyr::group_by(dt, country)
      dt <- dplyr::summarise_at(dt, 2:20, sum)
    }
  }
  if(panel == TRUE){
    dt<-tidyr::gather(dt, v1, fires, 'FIRE_2001':'FIRE_2019')
    dt <- dplyr::mutate(dt, year = substr(v1, 6, 9))
    dt <- dplyr::select(dt, -v1)
  }
  warning('\n
  -------------------------------------------------------------------------------------- \n
  Please, cite https://earthdata.nasa.gov/earth-observation-data/near-real-time/firms \n
  -------------------------------------------------------------------------------------- ')
  return(dt)
}
Helson-Gomes/BReco documentation built on May 17, 2020, 4:37 p.m.