CEDA ocean colour placeholder

[mdsumner]

@raymondben

these are merged across sensor, available patchily for daily, 5day, 8day, monthly, clim and for L3bin and mapped, and there is

• Kd
• chlor_a
• iop
• RRS

EDIT: release 6.0 goes up to end of 2022 but only for 5d, monthly, annual - it's a bit inconsistent wha'ts available in the different releases - , just I used 5.0 for my examples originally

https://dap.ceda.ac.uk/neodc/esacci/ocean_colour/data/v6.0-release/geographic/netcdf/

I download because gdal is slow via url netcdf/vsicurl etc prefix needed to avoid HDF5 driver and opendap ....

library(whatarelief)
u <- "https://dap.ceda.ac.uk/neodc/esacci/ocean_colour/data/v5.0-release/geographic/netcdf/chlor_a/5day/v5.0/2020/ESACCI-OC-L3S-CHLOR_A-MERGED-5D_DAILY_4km_GEO_PML_OCx-20201231-fv5.0.nc?download=1"

curl::curl_download(u, tf <- tempfile(fileext = basename(u)))
vsicurl <- function(x) sprintf("/vsicurl/%s", x)
netcdf <- function(x) sprintf("NETCDF:%s", x)
dsn <- vapour::vapour_vrt(vapour::vapour_sds_names(netcdf(tf))[6], projection = "OGC:CRS84")


ex <- c(-1, 1, -1, 1) * 1e7
prj <- "+proj=laea +lon_0=147 +lat_0=-55"
prj <- "OGC:CRS84"
ex <- c(-180, 180, -90, 90)
library(whatarelief)
s <- elevation(source = dsn, extent = ex, dimension = c(512, 512), projection = prj, resample = "near")
pal <- palr::chl_pal(palette = TRUE)

ximage::ximage(log(s), extent = ex, col = pal$cols)
maps::map(add = TRUE)

thanks to Philippe Massicotte and this gist

https://gist.github.com/mdsumner/dbade9bb2a52bc53ce5413bf20c274e2