Illustrations

Maps, colours, and areas of the MEASO regions

All figures on this page are made in R from the published files and the package data sets; expand “Code” to see how.

Colours

The colours of the MEASO Summary for Policymakers (Figure 5); see Definition.

Code 
colour_table <- function(x) {
cell <- function(name, fill) {
  if (is.na(fill)) return(sprintf("<td>%s</td>", name))
  rgb <- col2rgb(fill)
  fg <- if (sum(rgb * c(0.299, 0.587, 0.114)) > 140) "black" else "white"
  sprintf('<td style="background:%s;color:%s">%s<br><code style="background:none;color:%s">%s</code></td>',
          substr(fill, 1, 7), fg, name, fg, substr(fill, 1, 7))
}
zones <- c("Antarctic", "Subantarctic", "Northern", NA)
rows <- vapply(unique(x$sector), function(s) {
  d <- x[x$sector == s, ]
  cells <- vapply(zones, function(z) {
    i <- which(if (is.na(z)) is.na(d$zone) else d$zone %in% z)
    cell(d$name[i], d$fill[i])
  }, "")
  sprintf("<tr><th>%s</th>%s</tr>", s, paste(cells, collapse = ""))
}, "")
cat('<table class="table" style="text-align:center">',
    "<tr><th></th><th>Antarctic</th><th>Subantarctic</th><th>Northern</th><th>(outside)</th></tr>",
    rows, "</table>", sep = "\n")
}
colour_table(x)
Table 1: Fill colours by sector and zone (version 06). T regions have no colour.
Antarctic Subantarctic Northern (outside)
Atlantic AOA
#F9BDBC 		AOS
#F27071 		AON
#E11F26 		AOT
CentralIndian CIA
#F9D7C2 		CIS
#F4B187 		CIN
#EC8444 		CIT
EastIndian EIA
#D1D3EA 		EIS
#B0B4D9 		EIN
#878DC5 		EIT
WestPacific WPA
#E0B8D6 		WPS
#C579B2 		WPN
#843A75 		WPT
EastPacific EPA
#9FC5C7 		EPS
#609FA0 		EPN
#4B7E7F 		EPT

Version 05 colours

The colours of version 05 (measo_names in the R package), as released. Sectors that alternate around the continent share a palette. The first of the two West Pacific Antarctic (WPA) rows in measo_names has a different colour (an error); this table shows the second.

Code 
colour_table(measo_names[!duplicated(measo_names$name, fromLast = TRUE), ])
Table 2: Fill colours by sector and zone (version 05).
Antarctic Subantarctic Northern (outside)
EastPacific EPA
#FFFFFF 		EPS
#5F9EA0 		EPN
#4B7D7E 		EPT
Atlantic AOA
#BCECFE 		AOS
#1094AF 		AON
#054e70 		AOT
CentralIndian CIA
#FFFFFF 		CIS
#5F9EA0 		CIN
#4B7D7E 		CIT
EastIndian EIA
#EAFAFF 		EIS
#00AFD5 		EIN
#016074 		EIT
WestPacific WPA
#BCECFE 		WPS
#1094AF 		WPN
#054e70 		WPT

Polar map

Code 
cst <- measo_regions06_coastline
par(mar = rep(0, 4))
plot(st_union(st_geometry(x)), col = "grey80", border = NA)  ## continent shows through
plot(st_geometry(cst), col = cst$fill, border = "grey60", lwd = 0.3, add = TRUE)
pos <- st_coordinates(suppressWarnings(st_point_on_surface(st_geometry(x))))
text(pos, x$name, cex = 0.7)
Figure 1: MEASO regions in the polar equal area projection, with the continent cut out.

Longitude/latitude map

Code 
ll <- measo_regions06_ll_coastline
par(mar = c(2.5, 2.5, 0.5, 0.5))
plot(st_geometry(ll), col = ifelse(is.na(ll$fill), "grey97", ll$fill),
     border = "grey50", axes = TRUE, xlim = c(-180, 180), ylim = c(-85, -30))
## label positions from the polar centroids
cds <- st_coordinates(st_transform(st_centroid(st_geometry(measo_regions06)), 4326))
text(cds, measo_regions06$name, cex = 0.55)
Figure 2: MEASO regions in longitude/latitude. The West Pacific regions are in two parts, either side of the anti-meridian.

Sectors and zones

Code 
op <- par(mfrow = c(1, 2), mar = rep(0.2, 4))
## each sector's hue (its Northern zone colour)
sec_col <- with(x[x$zone %in% "Northern", ], setNames(fill, sector))
plot(st_geometry(x), col = sec_col[x$sector], border = "white")
legend("bottomleft", legend = names(sec_col), fill = sec_col, bty = "n", cex = 0.7)
zone_col <- c(Antarctic = "#BEE7FA", Subantarctic = "#0E94AF", Northern = "#0A4F70")
plot(st_geometry(x), col = ifelse(is.na(x$zone), "grey95", zone_col[x$zone]), border = "white")
legend("bottomleft", legend = c(names(zone_col), "outside MEASO"), fill = c(zone_col, "grey95"),
       bty = "n", cex = 0.7)
par(op)
Figure 3: Sectors (left, in their hues from Figure 5) and zones (right, in the zone colours of Figure 7) of the MEASO Summary for Policymakers.

Areas

Areas are calculated in the polar equal area projection, for the full regions and with the continent cut out.

Code 
tab <- data.frame(name = measo_regions06$name, label = measo_regions06$label,
                  full = round(as.numeric(st_area(measo_regions06)) / 1e12, 2),
                  without_continent = round(as.numeric(st_area(measo_regions06_coastline)) / 1e12, 2))
knitr::kable(tab, row.names = FALSE)
Table 3: Region areas (million km²).
name label full without_continent
AOA Antarctic Atlantic 12.11 9.89
AOS Subantarctic Atlantic 4.49 4.49
AON Northern Atlantic 4.63 4.63
AOT Temperate Atlantic 11.74 11.74
CIA Antarctic Central Indian 8.02 3.76
CIS Subantarctic Central Indian 8.12 8.12
CIN Northern Central Indian 3.90 3.90
CIT Temperate Central Indian 9.94 9.94
EIA Antarctic East Indian 4.05 1.69
EIS Subantarctic East Indian 3.70 3.70
EIN Northern East Indian 3.62 3.62
EIT Temperate East Indian 7.69 7.69
WPA Antarctic West Pacific 5.63 4.91
WPS Subantarctic West Pacific 2.84 2.84
WPN Northern West Pacific 7.28 7.28
WPT Temperate West Pacific 7.50 7.50
EPA Antarctic East Pacific 3.56 2.19
EPS Subantarctic East Pacific 3.28 3.28
EPN Northern East Pacific 6.26 6.26
EPT Temperate East Pacific 8.60 8.60

Plots on a map

The regions make a natural layout for small plots of data summarised by region. The map-patches vignette in the R package shows one way to do this with ggplot2 and patchwork.