Add tests

Author: danielvartan

.gitattributes

@@ -2,6 +2,6 @@
 
 *.qmd linguist-language=R
 
-images/* linguist-documentation
-docs/* linguist-documentation
-ttf/* linguist-documentation
+images/** linguist-documentation
+docs/** linguist-documentation
+ttf/** linguist-documentation

R/loop-geobr.R

@@ -43,37 +43,36 @@ year <- 2024
 #     sort()
 # )
 
-code <- c(
-  double_quote("all"),
-  c(
-    2507507, # João Pessoa (PB)
-    3106200, # Belo Horizonte (MG)
-    3550308, # São Paulo (SP)
-    4314902, # Porto Alegre (RS)
-    1501402 # Belém (PA)
-  )
-)
-
-# code <-
-#   read_state(
-#     "all",
-#     year = 2020
-#   ) |>
-#   arrange(1) |>
-#   pull(1) %>%
-#   c(double_quote("all"), .)
-
-# code <-
-#   here("pmtiles", "geobr", "read_municipality") |>
-#   dir_ls(type = "file") |>
-#   basename() |>
-#   str_extract("\\d+") |>
-#   as.integer() |>
-#   unique() |>
-#   sort() %>%
-#   setdiff(code, .)
+# code <- c(
+#   double_quote("all"),
+#   c(
+#     2507507, # João Pessoa (PB)
+#     3106200, # Belo Horizonte (MG)
+#     3550308, # São Paulo (SP)
+#     4314902, # Porto Alegre (RS)
+#     1501402 # Belém (PA)
+#   )
+# )
 
-# length(code)
+code <-
+  read_municipality(year = year) |>
+  arrange(1) |>
+  pull(1) # %>%
+# c(double_quote("all"), .)
+
+code <-
+  here("pmtiles", "geobr", "read_municipality") |>
+  dir_ls(type = "file") |>
+  basename() |>
+  str_subset(paste0("year-", year)) |>
+  str_extract("\\d+") |>
+  as.integer() |>
+  unique() |>
+  sort() %>%
+  setdiff(code, .) |>
+  sort()
+
+length(code)
 
 # Replace Fixed Values -----
 

qmd/geobr.qmd

@@ -80,7 +80,7 @@ for (i in c(data_dir)) {
 }
 ```
 
-### Set Initial Variables
+### Set Initial Parameters
 
 ::: {.callout-note}
 The `code` variable specifies the IBGE code of the Brazilian territorial unit to retrieve. Use `"all"` to retrieve all units. For the full list of IBGE codes and their corresponding territorial units, see <https://sidra.ibge.gov.br/territorio>.
@@ -92,7 +92,7 @@ The `code` variable specifies the IBGE code of the Brazilian territorial unit to
 # brazil_municipality_code("São Paulo")
 # "all" to select all units
 
-code <- 1501402
+code <- 5106000
 ```
 
 ::: {.callout-note}
@@ -128,7 +128,7 @@ The `h3jsr` variable specifies whether to process the spatial data with the [`h3
 :::
 
 ```{r}
-h3jsr <- TRUE
+h3jsr <- FALSE
 ```
 
 ::: {.callout-note}

qmd/geodata.qmd

@@ -80,7 +80,7 @@ for (i in c(data_dir, raw_data_dir)) {
 }
 ```
 
-### Set Initial Variables
+### Set Initial Parameters
 
 ::: {.callout-note}
 The `country` variable specifies which country to retrieve. Use a three-letter country code that follows the [ISO 3166-1 alpha-3](https://en.wikipedia.org/wiki/ISO_3166-1_alpha-3) standard (e.g., `BRA`, `ITA`, or `USA`).

qmd/tests.qmd

@@ -0,0 +1,110 @@
+```{r}
+#| label: setup
+#| include: false
+
+library(here)
+
+here("R", ".setup.R") |> source()
+```
+
+## Overview
+
+This document describes the testing process used after the creation of the [PMTiles](https://docs.protomaps.com/pmtiles/) files.
+
+## Methods
+
+### Code Style
+
+The Tidyverse [tidy tools manifesto](https://tidyverse.tidyverse.org/articles/manifesto.html) [@wickham2023c], [code style guide](https://style.tidyverse.org/) [@wickhama] and [design principles](https://design.tidyverse.org/) [@wickhamc] were followed to ensure consistency and enhance readability.
+
+### Reproducibility
+
+The pipeline is fully reproducible and can be run again at any time. To ensure consistent results, the [`renv`](https://rstudio.github.io/renv/) package [@ushey2025] is used to manage and restore the R environment. See the [README](https://github.com/cem-usp/pmtilesbr/blob/main/README.md) file in the code repository to learn how to run it.
+
+## Set the Environment
+
+### Load Packages
+
+```{r}
+#| label: Set the Environment
+#| output: false
+
+library(cli)
+library(fs)
+library(here)
+library(pmtiles) # github.com/walkerke/pmtiles
+library(stringr)
+```
+
+### Set Data Directories
+
+```{r}
+raw_data_dir <- here("data-raw")
+data_dir <- here("pmtiles")
+```
+
+```{r}
+for (i in c(data_dir, raw_data_dir)) {
+  if (!dir_exists(i)) {
+    dir_create(i, recurse = TRUE)
+  }
+}
+```
+
+### Set Initial Parameters
+
+::: {.callout-note}
+The `year` variable specifies the year of the data.
+:::
+
+```{r}
+year <- 2024
+```
+
+```{r}
+pmtiles_dir <- data_dir
+```
+
+## Test and Delete Corrupted PMTiles
+
+```{r}
+#| label: Test and Delete Corrupted PMTiles
+
+files <-
+  pmtiles_dir |>
+  dir_ls(
+    type = "file",
+    recurse = TRUE,
+    glob = "*.pmtiles"
+  ) |>
+  str_subset(
+    paste0("year-", year)
+  )
+```
+
+```{r}
+cli_progress_bar(
+  name = "Testing PMTiles",
+  total = length(files)
+)
+
+for (i in files) {
+  test <-
+    i |>
+    pm_show(tilejson = TRUE) |>
+    try(silent = TRUE)
+
+  if (inherits(test, "try-error")) {
+    file_delete(i)
+  }
+
+  cli_progress_update()
+}
+
+cli_progress_done()
+```
+
+## References {.unnumbered}
+
+::: {#refs}
+:::