inst/doc/geospatial_data.R

In ColOpenData: Download Colombian Demographic, Climate and Geospatial Data

## ----include = FALSE----------------------------------------------------------
knitr::opts_chunk$set(
  collapse = TRUE,
  comment = "#>",
  eval = identical(tolower(Sys.getenv("NOT_CRAN")), "true")
)

## ----available datasets, echo = FALSE-----------------------------------------
#  code <- c(
#    "DPTO", "MPIO",
#    "MPIOCL", "MZN",
#    "SECR", "SECU",
#    "SETR", "SETU",
#    "ZU"
#  )
#  level <- c(
#    "Department", "Municipality",
#    "Municipality including Class", "Block", "Rural Sector", "Urban Sector",
#    "Rural Section", "Urban Section", "Urban Zone"
#  )
#  dictionary_key <- c(
#    "DANE_MGN_2018_DPTO", "DANE_MGN_2018_MPIO",
#    "DANE_MGN_2018_MPIOCL", "DANE_MGN_2018_MZN",
#    "DANE_MGN_2018_SECR", "DANE_MGN_2018_SECU",
#    "DANE_MGN_2018_SETR", "DANE_MGN_2018_SETU",
#    "DANE_MGN_2018_ZU"
#  )
#  
#  mgncnpv <- data.frame(
#    Code = code, Level = level, Name = dictionary_key,
#    stringsAsFactors = FALSE
#  )
#  knitr::kable(mgncnpv)

## ----library imports, results = "hide", warning = FALSE, message = FALSE------
#  library(ColOpenData)
#  library(dplyr)
#  library(sf)
#  library(ggplot2)
#  library(leaflet)

## ----download data------------------------------------------------------------
#  dpto <- download_geospatial(
#    spatial_level = "dpto",
#    simplified = TRUE,
#    include_geom = TRUE,
#    include_cnpv = TRUE
#  )
#  
#  head(dpto)

## ----dictionary for urban sections--------------------------------------------
#  dict <- geospatial_dictionary(spatial_level = "dpto", language = "EN")
#  
#  head(dict)

## -----------------------------------------------------------------------------
#  internet_cov <- dpto %>% mutate(internet = round(viv_internet / viviendas, 2))

## ----ggplot2------------------------------------------------------------------
#  ggplot(data = internet_cov) +
#    geom_sf(mapping = aes(fill = internet), color = NA) +
#    theme_minimal() +
#    theme(
#      plot.background = element_rect(fill = "white", colour = "white"),
#      panel.background = element_rect(fill = "white", colour = "white"),
#      panel.grid = element_blank(),
#      axis.text = element_blank(),
#      axis.ticks = element_blank()
#    ) +
#    scale_fill_gradient("Percentage", low = "#10bed2", high = "#deff00") +
#    ggtitle(
#      label = "Internet coverage",
#      subtitle = "Colombia"
#    )

## ----define color palette-----------------------------------------------------
#  colfunc <- colorRampPalette(c("#10bed2", "#deff00"))
#  pal <- colorNumeric(
#    palette = colfunc(100),
#    domain = internet_cov[["internet"]]
#  )

## ----leaflet------------------------------------------------------------------
#  leaflet(internet_cov) %>%
#    addProviderTiles(providers$CartoDB.Positron) %>%
#    addPolygons(
#      stroke = TRUE,
#      weight = 0,
#      color = NA,
#      fillColor = ~ pal(internet_cov[["internet"]]),
#      fillOpacity = 1,
#      popup = paste0(internet_cov[["internet"]])
#    ) %>%
#    addLegend(
#      position = "bottomright",
#      pal = pal,
#      values = ~ internet_cov[["internet"]],
#      opacity = 1,
#      title = "Internet Coverage"
#    )