data-raw/rawdata.R

library(rgdal)
library(raster)
library(dplyr)

file <- "data-raw/marmap_coord_-180;-90;180;90_res_10.csv"
nam <- c("lon", "lat", "z")
agg <- 20

d.bath <- read.csv(file) %>%  tbl_df %>% setNames(nam)
r <- raster(extent(-180, 180, -90, 90), res=1/6)
projection(r) <- "+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs"
r <- setValues(r, d.bath$z)

bathymetry <- aggregate(r, agg) %>% rasterToPoints %>% tbl_df %>% setNames(nam)
borders <- ggplot2::map_data("world")
names(borders)[1] <- "lon"
devtools::use_data(bathymetry, overwrite=TRUE)
devtools::use_data(borders, overwrite=TRUE)
leonawicz/mapmate documentation built on May 21, 2019, 5:09 a.m.