data-raw/DATASET.R

In nick-gauthier/anthromes-12k: Tools for Analyzing Long-Term Population and Land Use

## code to prepare `DATASET` dataset goes here

library(anthromes)
library(stars)
library(tidyverse)

bbox <- st_bbox(c(xmin = 29, xmax = 39, ymin = 30, ymax = 40), crs = 4326)

hyde_med <- get_hyde(dir = 'data-raw/') %>%
  st_crop(bbox) %>%
  # current bug in stars loses offset info when filtered
  #filter(time %in% time_steps_millennia) %>%
  .[,,,c(8,9,10,11,21,58)] %>%  # so subset instead
  st_as_stars() %>%
  # change the names to something more readable
  setNames(c('crops', 'grazing', 'rice', 'pop', 'irrigation', 'urban'))

inputs_med <- get_hyde('inputs', dir = 'data-raw/') %>%
  st_crop(bbox)

dgg_med <- read_sf('dgg_land.gpkg')[st_as_sfc(bbox),, op = st_covered_by]

anthromes_med <- get_anthromes(hyde_med, inputs_med)

usethis::use_data(hyde_med,
                  inputs_med,
                  dgg_med,
                  anthromes_med,
                  overwrite = TRUE)