make_data.Rmd

---
title: "Making polygons"
csl: the-american-naturalist.csl
output:
  html_document:
    theme: cerulean
    toc: yes
  pdf_document:
    toc: yes
<!-- bibliography: references.bib -->
editor_options: 
  chunk_output_type: console
---

<!--
IMAGES:
Insert them with: ![alt text](image.png)
You can also resize them if needed: convert image.png -resize 50% image.png
If you want to center the image, go through HTML code:
<div style="text-align:center"><img src ="image.png"/></div>

REFERENCES:
For references: Put all the bibTeX references in the file "references.bib"
in the current folder and cite the references as @key or [@key] in the text.
Uncomment the bibliography field in the above header and put a "References"
title wherever you want to display the reference list.
-->

<style type="text/css">
.main-container {
  max-width: 1370px;
  margin-left: auto;
  margin-right: auto;
}
</style>

```{r general_options, include = FALSE}
knitr::knit_hooks$set(
  margin = function(before, options, envir) {
    if (before) par(mgp = c(1.5, .5, 0), bty = "n", plt = c(.105, .97, .13, .97))
    else NULL
  },
  prompt = function(before, options, envir) {
    options(prompt = if (options$engine %in% c("sh", "bash")) "$ " else "> ")
  })

knitr::opts_chunk$set(margin = TRUE, prompt = TRUE, comment = "",
                      collapse = TRUE, cache = FALSE, autodep = TRUE,
                      dev.args = list(pointsize = 11), fig.height = 3.5,
                      fig.width = 4.24725, fig.retina = 2, fig.align = "center")

options(width = 137)
```

```{r}
library(dplyr)
library(sf)
```

```{r}
decide <- "../../cleaned_data/GIS/Layers/Lao Decide/data-827086/data.shp" %>%
  st_read() %>% 
  select(uscne, uucne) %>% 
  transmute(district = sub(" *District *$", "", uscne),
            village = sub("^ *B\\. *", "", uucne))
```

This is the one sent by Olivier through Slack on June 21st:

```{r}
olivier <- "../../cleaned_data/GIS/Layers/Modified Villages/VIllage1806.shp" %>% 
  st_read(crs = 4326) %>% 
  transmute(village = Name)
```

```{r}
pacs <- readr::read_csv("../census/data/census.csv",
                        col_types = paste(c("ddccciiiiilllc", rep("d", 246)), collapse = ""))
```

```{r}
sort(setdiff(olivier$village, pacs$Village_Name))
```

```{r}
sort(setdiff(pacs$Village_Name, olivier$village))
```

```{r}
raw_qgis <- st_read("../../raw_data/Lao Statistics Bureau/QGIS/Modified Village/Vientiane.shp")
```