r.iresmi.net - Scales

Author

Michaël

Published

2022-11-12

Modified

2024-04-21

Photo of sale rulers — Scale rulers – © by Frank

Day 12 of 30DayMapChallenge: « scales » (previously).

French « départements » at the same size but tidily ordered by area, using data from Adminexpress made in Simplifying polygons layers.

library(tidyverse)
library(sf)
library(units)
library(tmap)

"~/data/adminexpress/adminexpress_cog_simpl_000_2022.gpkg" |> 
  read_sf(layer = "departement") |> 
  filter(!is.na(nom)) |> 
  mutate(surf_km2 = set_units(st_area(geom), km^2), 
         rank = dense_rank(surf_km2), 
         nom_dep = fct_rev(fct_reorder(as_factor(nom), rank))) |> 
  tm_shape() +
  tm_polygons(fill = "lightgrey") +
  tm_facets(by = "nom_dep", free.scales = TRUE)

A small multiple map of french départements ordered by area — Figure 1: Départements

---
title: "Scales"
description: "Day 12 of 30DayMapChallenge"
author: "Michaël"
date: "2022-11-12"
date-modified: last-modified
categories:
  - R
  - 30DayMapChallenge
  - datavisualization
  - spatial 
draft: false
freeze: true
editor_options: 
  chunk_output_type: console
---

[![Scale rulers -- © by Frank](images/38274549456_f1d36db345_w.jpg){fig-alt="Photo of sale rulers"}](https://flic.kr/p/21jc567)

Day 12 of [30DayMapChallenge](https://30daymapchallenge.com/): « scales » ([previously](/#category=30DayMapChallenge)).

French « départements » at the same size but tidily ordered by area, using data from [Adminexpress](https://geoservices.ign.fr/adminexpress) made in [Simplifying polygons layers](/posts/simplifying_polygon_layers/).

```{r}
#| label: fig-scale
#| fig-cap: Départements
#| fig-alt: A small multiple map of french départements ordered by area
library(tidyverse)
library(sf)
library(units)
library(tmap)

"~/data/adminexpress/adminexpress_cog_simpl_000_2022.gpkg" |> 
  read_sf(layer = "departement") |> 
  filter(!is.na(nom)) |> 
  mutate(surf_km2 = set_units(st_area(geom), km^2), 
         rank = dense_rank(surf_km2), 
         nom_dep = fct_rev(fct_reorder(as_factor(nom), rank))) |> 
  tm_shape() +
  tm_polygons(fill = "lightgrey") +
  tm_facets(by = "nom_dep", free.scales = TRUE)
```