data-raw/add_landseamask.R

library(tidyverse)
landmask <- readr::read_csv("data-raw/landmask_global_0p2deg.csv",
                col_names = c("longitude", "latitude", "land"),
                skip = 1) %>%
  dplyr::mutate(land = factor(land))

landmask <- landmask %>%
  arrange(latitude) %>%
  mutate(latitude = format(latitude, digits = 1, trim = TRUE)) %>%
  pivot_wider(names_from = "latitude", values_from = land) %>%
  arrange(longitude) %>%
  mutate(longitude = format(longitude, digits = 1, trim = TRUE)) %>%
  column_to_rownames("longitude") %>%
  as.matrix()

# ggplot(landmask, aes(x = longitude, y = latitude, color = land)) +
#   geom_point()

usethis::use_data(landmask, overwrite = TRUE)
geanders/stormwindmodel documentation built on Sept. 27, 2022, 6:47 a.m.