【R】afrilearndata

1. はじめに

afrilearndataは、アフリカのデータをまとめたパッケージで、以下の6つのデータから構成されています。

1. africountries polygons, country boundaries
2. africontinent polygons, continent outline including madagascar
3. afrihighway lines, trans african highway network
4. africapitals points, capital cities
5. afripop2020 raster grid, population density 2020 from WorldPop aggregated to 20km squares
6. afripop2000 raster grid, population density 2000 from WorldPop aggregated to 20km squares

これらのデータを可視化してみたいと思います。

2. インストール

インストールは、Githubで。

remotes::install_github("afrimapr/afrilearndata")

3. つかってみる

まずは、人工毎のコロプレス図を描きます。

library(afrilearndata)
library(ggplot2)
library(ggrepel)
library(sf)
library(tidyverse)
library(kableExtra)

dat <- africountries %>%
  mutate( 
    centroid = st_centroid(geometry),
    x = st_coordinates(centroid)[, 1],
    y = st_coordinates(centroid)[, 2]
  ) %>%
  arrange(desc(pop_est))

dat %>%
  ggplot() +
  geom_sf(aes(fill = pop_est)) +  
  coord_sf(datum = NA) + 
  scale_fill_viridis_c(alpha = 0.6) +
  theme_void()+ 
  geom_text_repel(aes(x = x, y = y, label = name), col="black", size = 3)


table_df<-data.frame(Country=dat$name, Population=dat$pop_est)

library(clipr)
table_df %>% 
  kable(align = "c", row.names=FALSE) %>%
  kable_styling(full_width = F) %>%
  column_spec(1, bold = T) %>%
  collapse_rows(columns = 1, valign = "middle") %>% 
  write_clip

例にあるように、すべてのデータを使ってtmapで表現。

library(tmap)

tm_shape(afripop2020) +
  tm_raster(palette = rev(viridisLite::magma(5)), breaks=c(0,2,20,200,2000,25000)) +
  tm_shape(africountries) +
  tm_borders("red", lwd = .7) +
  tm_shape(afrihighway) +
  tm_lines(col = "green", alpha=0.6) + 
  tm_shape(africapitals) +
  tm_symbols(col = "blue", alpha=0.6, scale = .4 )+
  tm_legend(show = FALSE)

4. さいごに

こうやって、データをまとめてくれると非常に使いやすいです（笑）。恥ずかしながら不勉強でアフリカのことをよく知らないのですが、これを機に勉強してみたいと思います。