Bad map – r.iresmi.net

Author

Michaël

Published

2023-11-04

Modified

2023-11-19

A photo of The Alte Nahebrücke, a medieval stone arch bridge in Bad Kreuznach, in western Germany, dating from around 1300 — Bad Kreuznach, Germany – CC-BY-NC by Billy Wilson

Day 4 of 30DayMapChallenge: « Bad map » (previously).

A map of all populated places containing the word “Bad” from Geonames.

library(tidyverse)
library(leaflet)
library(httr)
library(fs)
library(glue)

gn_file <- "~/data/geonames/allCountries.zip"

if (!file_exists(gn_file)) {
  GET("http://download.geonames.org/export/dump/allCountries.zip",
      write_disk(gn_file))
}

gn <- read_delim(
  gn_file, 
  delim = "\t", 
  col_names = c("geonameid",     
                "name",          
                "asciiname",     
                "alternatenames",
                "latitude",      
                "longitude",     
                "feature_class", 
                "feature_code",  
                "country_code",  
                "cc2",           
                "admin1_code",   
                "admin2_code",   
                "admin3_code",   
                "admin4_code",   
                "population",    
                "elevation",     
                "dem",           
                "timezone",      
                "modification_date"))

bad <- gn |> 
  filter(str_detect(feature_code, "^PPL"),
         str_detect(asciiname, "\\b[Bb]ad\\b"))

bad |> 
  leaflet() |> 
  addCircleMarkers(popup = ~ glue("{name} ({asciiname})<br />
                                  <b>{country_code}</b><br />
                                  {population} inhab.")) |> 
  addTiles(attribution = "Geonames and © OpenStreetMap, ODbL")

Figure 1: Bad map

So we have 471 “Bad” populated places…

---
title: "Bad map"
description: "Day 4 of 30DayMapChallenge"
author: "Michaël"
date: "2023-11-04T12:00:00"
date-modified: last-modified
categories:
  - R
  - 30DayMapChallenge
  - spatial
  - datavisualization
draft: false
freeze: true
editor_options:
  chunk_output_type: console
chunk_output_type: console
execute:
  eval: true
---

[![Bad Kreuznach, Germany -- CC-BY-NC by Billy Wilson](images/49007861517_19dcc9082c_c.jpg){fig-alt="A photo of The Alte Nahebrücke, a medieval stone arch bridge in Bad Kreuznach, in western Germany, dating from around 1300"}](https://flic.kr/p/2hEE6PB)

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

A map of all populated places containing the word "Bad" from [Geonames](https://www.geonames.org/).

```{r}
library(tidyverse)
library(leaflet)
library(httr)
library(fs)
library(glue)

gn_file <- "~/data/geonames/allCountries.zip"

if (!file_exists(gn_file)) {
  GET("http://download.geonames.org/export/dump/allCountries.zip",
      write_disk(gn_file))
}

gn <- read_delim(
  gn_file, 
  delim = "\t", 
  col_names = c("geonameid",     
                "name",          
                "asciiname",     
                "alternatenames",
                "latitude",      
                "longitude",     
                "feature_class", 
                "feature_code",  
                "country_code",  
                "cc2",           
                "admin1_code",   
                "admin2_code",   
                "admin3_code",   
                "admin4_code",   
                "population",    
                "elevation",     
                "dem",           
                "timezone",      
                "modification_date"))
```

```{r}
#| label: fig-bad-map
#| fig-cap: Bad map
#| fig-alt: A map of populated places containing the word Bad
bad <- gn |> 
  filter(str_detect(feature_code, "^PPL"),
         str_detect(asciiname, "\\b[Bb]ad\\b"))

bad |> 
  leaflet() |> 
  addCircleMarkers(popup = ~ glue("{name} ({asciiname})<br />
                                  <b>{country_code}</b><br />
                                  {population} inhab.")) |> 
  addTiles(attribution = "Geonames and © OpenStreetMap, ODbL")
```

So we have `r nrow(bad)` "Bad" populated places...