chla-is-back.R

## https://github.com/AustralianAntarcticDivision/raadtools/commit/f26cafbc10efce72f7cd24abbe70f5acd65ba270
library(raadtools)

## Loading required package: raster

## Loading required package: sp

dates <- seq(as.Date("2016-12-15"), length = 30, by = "1 day")
print(range(dates))

## [1] "2016-12-15" "2017-01-13"

xylim <- raster::extent(100, 160, -70, -30)
system.time({
johnson <- readchla(dates, product = "MODISA", xylim = xylim)
nasa    <- readchla(dates, product = "MODISA", xylim = xylim, 
                    algorithm = "nasa")
})

##    user  system elapsed 
##  47.988   2.952  54.910

pal <- palr::chlPal(palette = TRUE)
par(mar = rep(0, 4))
image(johnson, col = pal$cols[-1], breaks = pal$breaks, asp = "", useRaster = TRUE)

image(nasa, col = pal$cols[-1], breaks = pal$breaks, asp = "", useRaster = TRUE)

print(johnson)

## class       : RasterLayer 
## dimensions  : 944, 1440, 1359360  (nrow, ncol, ncell)
## resolution  : 0.04166667, 0.04166667  (x, y)
## extent      : 100, 160, -69.33333, -30  (xmin, xmax, ymin, ymax)
## coord. ref. : +init=epsg:4326 +proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0 
## data source : in memory
## names       : layer 
## values      : 0.0250624, 16.9833  (min, max)