data-raw/cetesb.R

In curso-r/basesCursoR: Bases Usadas nos Cursos da Curso-R

## code to prepare `cetesb` dataset goes here

estacoes <- c(73, 65, 92, 95, 96, 120, 83, 85, 99, 72, 86, 64, 98, 63)

tab_estacoes <- Rpollution::cetesb_station_ids %>%
  dplyr::filter(id %in% estacoes)

base_CO <- purrr::map_dfr(
  na.omit(tab_estacoes$id[tab_estacoes$CO == "yes"]),
  Rpollution::scraper_cetesb,
  parameter = 16,
  start = "01/01/2017",
  end = "31/05/2020",
  login = "user",
  password = "senha",
  safe = TRUE
)

base_NO <- purrr::map_dfr(
  na.omit(tab_estacoes$id[tab_estacoes$NO == "yes"]),
  Rpollution::scraper_cetesb,
  parameter = 17,
  start = "01/01/2017",
  end = "31/05/2020",
  login = "user",
  password = "senha",
  safe = TRUE
)

base_NO2 <- purrr::map_dfr(
  na.omit(tab_estacoes$id[tab_estacoes$NO2 == "yes"]),
  Rpollution::scraper_cetesb,
  parameter = 15,
  start = "01/01/2017",
  end = "31/05/2020",
  login = "user",
  password = "senha",
  safe = TRUE
)

base_MP10 <- purrr::map_dfr(
  na.omit(tab_estacoes$id[tab_estacoes$MP10 == "yes"]),
  Rpollution::scraper_cetesb,
  parameter = 12,
  start = "01/01/2017",
  end = "31/05/2020",
  login = "user",
  password = "senha",
  safe = TRUE
)

base_o3 <- purrr::map_dfr(
  na.omit(tab_estacoes$id[tab_estacoes$O3 == "yes"]),
  Rpollution::scraper_cetesb,
  parameter = 63,
  start = "01/01/2017",
  end = "31/05/2020",
  login = "user",
  password = "senha",
  safe = TRUE
)

cetesb <- dplyr::bind_rows(
  base_CO,
  base_NO,
  base_NO2,
  base_MP10,
  base_o3
) %>%
  dplyr::rename(
    estacao_cetesb = stationname,
    poluente = parameter,
    data = date,
    hora = hour,
    concentracao = conc
  )

cetesb <- cetesb %>%
  dplyr::mutate(estacao_cetesb = dplyr::case_when(
    estacao_cetesb == "Cid.Universitária-USP-Ipen" ~ "IPEN-USP",
    estacao_cetesb == "N.Senhora do Ó" ~ "Nossa Senhora do Ó",
    TRUE ~ estacao_cetesb
  )) %>% #dplyr::distinct(estacao_cetesb)
  dplyr::left_join(
    readxl::read_excel("data-raw/excel/cetesb_station_geoposition.xlsx"),
    by = c("estacao_cetesb" = "stationname")
  )

readr::write_rds(cetesb, "inst/extdata/cetesb.rds", compress = "xz")