map_crpd.R

crpd_colors <- c("Signed Convention" = "#FFCC67", "Signed Convention & Protocol" = "lightblue", 
                 "Ratified Convention" = "#346699", "Ratified Convention & Protocol" = "#FF4001")

denmark_crpd <- dat_all |> 
  filter(country == "Denmark") |> 
  pull(crpd)
dat_all <- dat_all |> 
  add_row(
    crpd = denmark_crpd, country = "Greenland"
  )

library(rnaturalearth)
world <- ne_countries(returnclass = "sf")

dat_all <- dat_all |> 
  mutate(
    crpd_label = case_match(
      crpd,
      "1" ~ "Signed Convention", 
      "2" ~ "Signed Convention & Protocol", 
      "3" ~ "Ratified Convention",
      "4" ~ "Ratified Convention & Protocol"
      ) |> 
      factor(levels = names(crpd_colors))
    )

# -------------------------------------------------------------------------

## Plotly is a headache..

# sfdat <- world |> 
#   left_join(
#     dat_all, by = join_by(iso_a3_eh == iso3c)
#   ) |> 
#   st_as_sf() |> 
#   filter(!is.na(geometry))

# plot_ly() |>
#   # add_sf(
#   #   data = sfdat, 
#   #   type = "scatter",
#   #   text = ~name,
#   #   ## Warning about markers...?
#   #   # mode = "markers",
#   #   split = ~name,
#   #   hoveron = "fills",
#   #   hoverinfo = "text",
#   #   color = ~crpd_label,
#   #   colors = crpd_colors,
#   #   stroke = NA,
#   #   alpha = 1,
#   #   legendgroup = ~name
#   #   # ,
#   #   # colorscale = crpd_colors
#   # ) |> 
#   add_sf(
#     data = sfdat,
#     hovertext = "none",
#     hoverinfo = "none",
#     color = ~crpd_label,
#     colors = crpd_colors,
#     stroke = NA,
#     alpha = 1
#     ) |>
#   # add_sf(
#   #   data = sfdat,
#   #   text = ~name,
#   #   split = ~name,
#   #   hovertext = "text",
#   #   hoverinfo = "text",
#   #   visible = FALSE
#   # ) |> 
#   layout(
#     title = "CRPD and Optional Protocol Signatures and Ratifications",
#     geo = list(
#       projection = list(type = "orthographic"),
#       showframe = TRUE,
#       showcoastlines = TRUE,
#       coastlinecolor = "RebeccaPurple",
#       lataxis = list(range = c(-59, 90))
#       ),
#     legend = list(
#       orientation = "h",
#       x = 0.5,
#       xanchor = "center",
#       y = 1.1,  # Positions the legend above the top of the map
#       yanchor = "top"
#     )
#   ) 
# # |>
#   add_trace(
#     text = ~name,  # Assuming 'name' is the column with country names
#     mode = 'text',
#     textposition = 'middle center',
#     textfont = list(color = 'black'),
#     lon = ~lon,
#     lat = ~lat
#   )
# # |> 
# #   colorbar(
# #     title = "CRPD Levels",
# #     orientation = 'h',  # Set the colorbar horizontally
# #     x = 0.5,  # Center the colorbar
# #     xanchor = 'center',  # Anchor the center at x
# #     y = 1.05,  # Position it above the top
# #     yanchor = "bottom"
# #   )

# GGplot ------------------------------------------------------------------

dat_all <- dat_all |> 
  mutate(
    crpd_label = case_match(
      crpd,
      "1" ~ "Signed Convention", 
      "2" ~ "Signed Convention & Protocol", 
      "3" ~ "Ratified Convention",
      "4" ~ "Ratified Convention & Protocol"
    ) |> 
      factor(levels = names(crpd_colors))
  )

# reproject from longlat to robinson
world <- st_transform(world, "+proj=robin")

world <- world |> 
  mutate(center = st_centroid(geometry),
         lon = st_coordinates(center)[,1],
         lat = st_coordinates(center)[,2])

world <- world |> 
  filter(!is.na(name) & !name %in% c("Antarctica", "Fr. S. Antarctic Lands"))

denmark_greenland <- world |> 
  filter(name %in% c("Greenland", "Denmark"))

# Combine their geometries
combined_geometry <- st_union(denmark_greenland$geometry)

world <- world |>
  mutate(geometry = if_else(name == "Denmark", combined_geometry, geometry)) |>
  filter(name != "Greenland") |>
  select(geometry, name, iso_a3_eh, lon, lat)

sfdat <- world |> 
  right_join(
    dat_all, by = join_by(iso_a3_eh == iso3c)
  ) |> 
  st_as_sf() |>
  filter(!is.na(geometry))

graticules <- sf::st_graticule(crs = "+proj=robin")

p <- ggplot(data = sfdat) +
  # geom_sf(data = graticules, fill = "grey90", color = "white", alpha = 0.5) +
  geom_sf(aes(fill = crpd_label),
          # text = ~name, 
          # color = I("black"),
          alpha = 1) +
  scale_fill_manual(values = crpd_colors) + 
  labs(title = "CRPD and Optional Protocol Signatures and Ratifications", fill = "") +
  theme(
    legend.position = "top", 
    # Increase title text size significantly (adjust the size value as needed)
    plot.title = element_text(hjust = 0.5, face = "bold", color = "#346699", size = 22),
    # Increase legend (category) text size
    legend.title = element_text(size = 16),
    legend.text = element_text(size = 14),
    panel.grid = element_line(color = "lightblue", linetype = 1),
    axis.title.x = element_blank(),
    axis.text.x = element_blank(),
    axis.title.y = element_blank(),
    axis.ticks.x = element_blank()
  )

p <- p + 
  geom_label_repel(
    size = 2, 
    aes(x = lon, y = lat, label = name),
    label.padding = 0.15,
    box.padding = 0.15,
    label.size = 0.05,
    max.overlaps = 15
  )

ggsave("www/map_crpd.png", p, width = 30, height = 15, dpi = 300)