inst/doc/proxy.R

In CopernicusMarine: Search Download and Handle Data from Copernicus Marine Service Information

## ----include = FALSE----------------------------------------------------------
knitr::opts_chunk$set(
  collapse = TRUE,
  comment = "#>"
)
eval_chunks <-
  CopernicusMarine:::has_blosc &&
  curl::has_internet() &&
  CopernicusMarine::cms_get_password() != "" &&
  sf::st_drivers("raster", "^HDF5$")$vsi

## ----create_proxy, eval=eval_chunks-------------------------------------------
library(CopernicusMarine)
library(stars, quietly = TRUE)

my_proxy_gc <- cms_zarr_proxy(
  product       = "GLOBAL_ANALYSISFORECAST_PHY_001_024",
  layer         = "cmems_mod_glo_phy-thetao_anfc_0.083deg_P1D-m",
  asset         = "geoChunked")

my_proxy_tc <- cms_zarr_proxy(
  product       = "GLOBAL_ANALYSISFORECAST_PHY_001_024",
  layer         = "cmems_mod_glo_phy-thetao_anfc_0.083deg_P1D-m",
  asset         = "timeChunked")

print(my_proxy_tc)

## ----slice_time, eval=eval_chunks---------------------------------------------
time_slice <- my_proxy_gc[,2000, 1000, 48, 1:400]
show(time_slice)

## ----plot_time_slice, eval=eval_chunks----------------------------------------
time_slice <- st_as_stars(time_slice)
plot(st_get_dimension_values(time_slice, "time"), time_slice$thetao,
     xlab = "date", ylab = "temperature", type = "l")

## ----slice_area, eval=eval_chunks---------------------------------------------
geo_slice <- my_proxy_tc[,2000:2500, 1500:1750, 50, 500]
plot(geo_slice, col = hcl.colors(10))