data-raw/prepare_world_map_data.R

In sidora-tools/sidora.cli: Basic CLI for PANDORA (MPI-SHH DAG)

library(magrittr)

#### input info ####

# dimensions of the ascii plot
plot_width <- 55
plot_height <- 12

#### prepare data ####

# world map landmass as a raster
world_map_raster <- raster::rasterize(
  rnaturalearthdata::countries110, 
  raster::raster(xmn = -180, xmx = 180, ymn = -90, ymx = 90, nrow = plot_height, ncol = plot_width)
)

# extract worldmass positions
world_map_matrix <- t(raster::as.matrix(world_map_raster))
landmass_cells <- raster::xyFromCell(
  world_map_raster, 
  cell = which(!is.na(world_map_matrix))
) %>% 
  tibble::as_tibble()

# remapping of point positions to the ascii plot grid
x_breaks <- tibble::tibble(
  x_plot_pos = 1:plot_width
) %>%
  dplyr::mutate(
    x_coord_start = -180 + 360/plot_width * (x_plot_pos - 1),
    x_coord_stop = -180 + 360/plot_width * x_plot_pos
  )

y_breaks <- tibble::tibble(
  y_plot_pos = 1:plot_height
) %>%
  dplyr::mutate(
    y_coord_start = -90 + 180/plot_height * (y_plot_pos - 1),
    y_coord_stop = -90 + 180/plot_height * y_plot_pos
  )

landmass_cells_plot <- landmass_cells %>%
  dplyr::mutate(
    x = as.integer(cut(x, breaks = c(x_breaks$x_coord_start[1], x_breaks$x_coord_stop), labels = 1:plot_width)),
    y = as.integer(cut(y, breaks =  c(y_breaks$y_coord_start[1], y_breaks$y_coord_stop), labels = 1:plot_height)),
    land = TRUE
  )

# compiling plot grid
plot_grid <- expand.grid(
  x = 1:plot_width,
  y = 1:plot_height
) %>% 
  dplyr::left_join(
    landmass_cells_plot, by = c("x", "y")
  ) %>%
  dplyr::mutate(
    land =  tidyr::replace_na(land, FALSE)
  )

# store data
world_map_data <- list(
  plot_width = plot_width, 
  plot_height = plot_height, 
  x_breaks = x_breaks, 
  y_breaks = y_breaks, 
  plot_grid = plot_grid
)
usethis::use_data(world_map_data, overwrite = T)