maschette / somap Goto Github PK

View Code? Open in Web Editor NEW

8.0 8.0 6.0 49.02 MB

Southern Ocean round maps

R 100.00%

somap's People

Contributors

Stargazers

Watchers

Forkers

australianantarcticdivision

somap's Issues

Why not use coord_proj("+proj=laea")?

A rhetorical question, but should be interesting to explore different approaches.

How would you make this map?

The chlorophyll-a data:

u <- "https://oceandata.sci.gsfc.nasa.gov/cgi/getfile/A20022132017243.L3m_MC_CHL_chl_ocx_9km.nc"
download.file(u, basename(u), mode = "wb")

The chlorophyll-a palette (breaks and colours):

pal <- palr::chlPal(palette = TRUE)

The projection and extents:

outer boundary is -20 (20S)
central longitude is 180, central latitude -90
inner graticule is xx <- c(0, 90, 180, 270, 360); yy <- c(-60, -40)

"+proj=laea +lat_0=-90 +lon_0=180 +datum=WGS84 +ellps=WGS84 +no_defs +towgs84=0, 0, 0"

The fronts are from orsifronts::orsifronts:

subset(orsifronts::orsifronts, front %in% c("pf", "saf"))

to dos

extend to coastline if given, control what the extents is defined by
increase used of sf
ggplot converter for output object
allow custom prj input
input types (e.g. spatial points)
palettes

SOleg subticks

Insert options for subticks with SOleg.

Legends for default plot.

Currently no option for legends in default_somap

croptograt

Try

build graticule of tiles, or just one big grat
union them all together (rgeos::gUnionCascaded, or sf::st_union)
mask the raster with the grat (or find the cell numbers that overlap, and set the rest to NA)
intersect the coastline with the grat
add the border

SOcool

Prototype:

library(sospatial)  ## devtools::install_github("AustralianAntarcticDivision/sospatial")
library(SOmap)
SOcool <- function(xs, ys, icedate = "2017-09-01") {
  ice <-  tibble::tibble(lon = rep_len(seq(-180, 179), length.out = length(packed_lats)), 
                lat = packed_lats/10,
                day = as.Date(rep(ice_dates, each = 360)))
  icedate <- as.Date(icedate)
  f <- "%Y-%m-%d"
  if (icedate < min(ice$day) || icedate > max(ice$day)) {
    stop(sprintf("sorry, available range of icedate is %s : %s", format(min(icedate), f), format(max(icedate), f)))
  }
                                                                
                                                                
  ind <- which(icedate == ice$day)
  x <- default_somap(xs, ys)
  ice <- dplyr::slice(ice, ind)
  ice[c("x", "y")] <- rgdal::project(as.matrix(ice[c("lon", "lat")]), projection(x$bathy))
  ix <- which(ice$x >= xmin(x$bathy) & ice$x <= xmax(x$bathy) &
                ice$y >= ymin(x$bathy) & ice$y <= ymax(x$bathy))
  ice <- dplyr::slice(ice, ix)
  lines(ice[c("x", "y")], lty = 2)
  x$ice <- ice
  invisible(x)
}

SOcool(xs, ys)

I expect sospatial will get incorporated into this package, though first we need a way to keep the ice data up to date - it's about a year out atm.

Why not just coord_proj("+proj=laea")?

A rhetorical question, but might be interesting to explore. How would you make this map?

The chlorophyll-a data:

u <- "https://oceandata.sci.gsfc.nasa.gov/cgi/getfile/A20022132017243.L3m_MC_CHL_chl_ocx_9km.nc"
download.file(u, basename(u), mode = "wb")

The chlorophyll-a palette (breaks and colours):

pal <- palr::chlPal(palette = TRUE)

The projection and extents:

outer boundary is -20 (20S)
central longitude is 180, central latitude -90
inner graticule is xx <- c(0, 90, 180, 270, 360); yy <- c(-60, -40)

"+proj=laea +lat_0=-90 +lon_0=180 +datum=WGS84 +ellps=WGS84 +no_defs +towgs84=0, 0, 0"

The fronts are from orsifronts::orsifronts:

subset(orsifronts::orsifronts, front %in% c("pf", "saf"))

waht

zzz

cast problem with albatross data

Why?

data("albatross", package = "adehabitatLT")
track <- rgdal::project(as.matrix(purrr::map_df(albatross, ~rbind(.x[c("x", "y")], NA))), "+proj=utm +zone=42 +south +datum=WGS84", inv = TRUE)
SOauto_map(track[,1], track[,2])
 Error in st_cast_sfc_default(x) : list item(s) not of class sfg

enable pkgdown, travis, appveyor, covr

why not just

Why not just coord_proj("+proj=laea")?

A rhetorical question, but should be interesting to explore different approaches.

How would you make this map?

The chlorophyll-a data:

u <- "https://oceandata.sci.gsfc.nasa.gov/cgi/getfile/A20022132017243.L3m_MC_CHL_chl_ocx_9km.nc"
download.file(u, basename(u), mode = "wb")

The chlorophyll-a palette (breaks and colours):

pal <- palr::chlPal(palette = TRUE)

The projection and extents:

outer boundary is -20 (20S)
central longitude is 180, central latitude -90
inner graticule is xx <- c(0, 90, 180, 270, 360); yy <- c(-60, -40)

"+proj=laea +lat_0=-90 +lon_0=180 +datum=WGS84 +ellps=WGS84 +no_defs +towgs84=0, 0, 0"

The fronts are from orsifronts::orsifronts:

subset(orsifronts::orsifronts, front %in% c("pf", "saf"))

what

aaa

maschette / somap Goto Github PK

somap's People

Contributors

Stargazers

Watchers

Forkers

somap's Issues

Recommend Projects

Recommend Topics

Recommend Org