data-raw/04_ubs_pontos.R

In bruno-pinheiro/asha: Analisys of Spatial Health Acessibility

library(dplyr)
library(sf)
library(ggmap)

# IMPORTAR OS DADOS -----------------------------------------------------------
file <- system.file("extdata", "ubs_pontos-raw.rda", package = "asha")
load(file)
data("ubs_malha")

# LIMPAR OS DADOS -------------------------------------------------------------
ubs_pontos <- ubs_pontos %>%
  dplyr::select(cnes = eq_cnes) %>%
  sf::st_set_crs(31983)

## Limpar UBS fechadas da base de pontos
diferencas <- list(
  novas = dplyr::filter(ubs_malha, !(cnes %in% ubs_pontos$cnes)) %>%
    dplyr::select(cnes, endereco),
  fechadas = dplyr::filter(ubs_pontos, !(cnes %in% ubs_malha$cnes))
)

## Extrair as que não foram fechadas
diferencas$antigas <- filter(ubs_pontos, !(cnes %in% diferencas$fechadas$cnes))

## GEOCOLOCALIZAR NOVAS UBS ---------------------------------------------------
diferencas$novas_geocode <- geocode(diferencas$novas$endereco)

diferencas$novas <- diferencas$novas_geocode %>%
  st_as_sf(coords = c("lon", "lat")) %>%
  st_set_crs(4326) %>%
  st_transform(31983) %>%
  bind_cols(diferencas$novas %>% as_tibble()) %>%
  mutate(cnes = as.character(cnes)) %>%
  select(cnes)

## UNIR PONTOS DE UBS NOVAS E ANTIGAS -----------------------------------------
ubs_pontos <- rbind(diferencas$antigas, diferencas$novas)

## CORRIGIR UBS COM ERRO DE LOCALIZAÇÃO ---------------------------------------
erros <- list(
  cnes = as.character(
    c(2787954, 2788306, 2788608, 2815273, 3219771, 3561186, 3708411,
      3762831, 4050223, 5731143, 6018912, 6194990, 6329527, 6332420, 6332471)
  ),
  endereco =
    c(" UBS Jardim Vitoria R. Lajedo, 123 - Cidade Tiradentes, São Paulo - SP",
      "Av. Edu Chaves, 1197 - Parque Edu Chaves, São Paulo - SP",
      "R. Frei Fidelis Mota, 1001 – Parque Cruzeiro do Sul, São Paulo - SP",
      "R. Humberto de Almeida, 279, Chacara Santana, São Paulo - SP",
      "Tv. Somos Todos Iguais, 915 - Jardim da Conquista, São Paulo - SP",
      "UBS Vera Poty, São Paulo",
      "R. Rolando Curti, 701 - Vila Clara, São Paulo - SP, 04413-070",
      "UBS Vila Calú - Estrada do M'Boi Mirim - Vila Calu, São Paulo - SP",
      "UBS Parque do Largo - Estrada da Baronesa. 1000, Parque do Lago, São Paulo - SP",
      "Rua Catarina Gabrielli, 236 - Americanópolis, São Paulo - SP",
      "R. Benedito Schunck, 8 - Emburá, São Paulo - SP, 04893-050",
      "AMA CHÁCARA CRUZEIRO DO SUL - Rua Mercedes Lopes - Vila Santana, São Paulo - SP",
      "R. Francisco de Sáles, 10 - Jardim das Palmas, São Paulo - SP, 05749-280",
      "UBS Barragem - Estrada Evangelista de Souza - Barragem, São Paulo - SP",
      "UBS MARSILAC - Estrada Engenheiro Marsilac - Emburá, São Paulo - SP")
)

erros$geocode <- ggmap::geocode(erros$endereco)
erros$corrigidos <- erros$geocode %>%
  sf::st_as_sf(coords = c("lon", "lat")) %>%
  sf::st_set_crs(4326) %>%
  sf::st_transform(31983) %>%
  dplyr::bind_cols(cnes = erros$cnes) %>%
  dplyr::mutate(cnes = as.character(cnes)) %>%
  dplyr::select(cnes)
erros$check <- dplyr::filter(ubs_pontos, cnes %in% erros$cnes)

ubs_pontos <- ubs_pontos %>%
  dplyr::filter(!(cnes %in% erros$corrigidos$cnes)) %>%
  rbind(erros$corrigidos)

# EXPORTAR OS DADOS -----------------------------------------------------------
usethis::use_data(ubs_pontos, overwrite = TRUE)

rm(list = ls())