# What is a country ?

This vignette shows how rnaturalearth allows mapping countries using different definitions of what a country is. What a country is can be more complicated than you might expect.

For example, from my own parochial perspective, it allows mapping the UK as a whole or separating out England, Scotland, Wales and Northern Ireland. It also allows you to exclude far away places like the Falkland Islands, or not. Mapping France it allows the inclusion or exclusion of French Guiana and islands in the South Pacific.

rnaturalearth is an R package to hold and facilitate interaction with natural earth vector map data.

Natural Earth is a public domain map dataset including vector country boundaries.

This vignette uses `plot` as a simple, quick way to plot the data obtained using `rnaturalearth`. `rnaturalearth` data can also be used to make more elaborate maps with `ggplot2`, `tmap` and other options.

## Load required packages

    library(rnaturalearth)

## Country types : countries, map_units and sovereignty

Natural Earth data are classified by `countries`, `map_units` and `sovereignty`. Below you will see that specifying `united kingdom` for

- `countries` gives the UK undivided
- `map_units` gives England, Scotland, Wales and Northern Ireland
- `sovereignty` includes the Falkland Islands

Filtering by `geounit` can give finer control, e.g. to plot Scotland alone, or France without French Guiana.

    # countries, UK undivided
    plot(ne_countries(country = "united kingdom", type = "countries")["geometry"])

    # map_units, UK divided into England, Scotland, Wales and Northern Ireland
    plot(ne_countries(country = "united kingdom", type = "map_units")["geometry"])

    # map_units, select by geounit to plot Scotland alone
    plot(ne_countries(geounit = "scotland", type = "map_units")["geometry"])

    # sovereignty, Falkland Islands included in UK
    plot(
      ne_countries(country = "united kingdom", type = "sovereignty")["geometry"],
      col = "red"
    )

    plot(ne_coastline(scale = 110L)["geometry"], lty = 3L, add = TRUE)

    # France, country includes French Guiana
    plot(ne_countries(country = "france")["geometry"])

    # France map_units includes French Guiana too
    plot(ne_countries(country = "france", type = "map_units")["geometry"])

    # France filter map_units by geounit to exclude French Guiana
    plot(ne_countries(geounit = "france", type = "map_units")["geometry"])

    # France sovereignty includes South Pacicic islands
    plot(
      ne_countries(country = "france", type = "sovereignty")["geometry"],
      col = "red"
    )

    plot(ne_coastline(scale = 110L)["geometry"], lty = 3L, add = TRUE)

## Country scales : small, medium and large

The different definitions of a country outlined above are available at different scales.

    plot(
      ne_countries(country = "united kingdom", scale = "large")["geometry"],
      main = "large scale"
    )

    # countries, medium scale
    plot(
      ne_countries(country = "united kingdom", scale = "medium")["geometry"],
      main = "medium scale"
    )

    # countries, small scale
    plot(
      ne_countries(country = "united kingdom", scale = "small")["geometry"],
      main = "small scale"
    )

## States, admin level 1, select by country or geounit

    # states country='united kingdom'
    plot(ne_states(country = "united kingdom")["geometry"])

    # states geounit='england'
    plot(ne_states(geounit = "england")["geometry"])

    # states country='france'
    plot(ne_states(country = "france")["geometry"])

    # states geounit='france'
    plot(ne_states(geounit = "france")["geometry"])