Skip to contents

Photon irradiance (i.e. quantum irradiance) for one or more wavebands of a light source spectrum.

Usage

q_irrad(
  spct,
  w.band,
  quantity,
  time.unit,
  scale.factor,
  wb.trim,
  use.cached.mult,
  use.hinges,
  allow.scaled,
  ...
)

# S3 method for default
q_irrad(
  spct,
  w.band,
  quantity,
  time.unit,
  scale.factor,
  wb.trim,
  use.cached.mult,
  use.hinges,
  allow.scaled,
  ...
)

# S3 method for source_spct
q_irrad(
  spct,
  w.band = NULL,
  quantity = "total",
  time.unit = NULL,
  scale.factor = 1,
  wb.trim = getOption("photobiology.waveband.trim", default = TRUE),
  use.cached.mult = getOption("photobiology.use.cached.mult", default = FALSE),
  use.hinges = NULL,
  allow.scaled = !quantity %in% c("average", "mean", "total"),
  naming = "default",
  return.tb = FALSE,
  ...
)

# S3 method for source_mspct
q_irrad(
  spct,
  w.band = NULL,
  quantity = "total",
  time.unit = NULL,
  scale.factor = 1,
  wb.trim = getOption("photobiology.waveband.trim", default = TRUE),
  use.cached.mult = getOption("photobiology.use.cached.mult", default = FALSE),
  use.hinges = NULL,
  allow.scaled = !quantity %in% c("average", "mean", "total"),
  naming = "default",
  ...,
  attr2tb = NULL,
  idx = "spct.idx",
  .parallel = FALSE,
  .paropts = NULL
)

Arguments

spct

an R object.

w.band

a list of waveband objects or a waveband object.

quantity

character string One of "total", "average" or "mean", "contribution", "contribution.pc", "relative" or "relative.pc".

time.unit

character or lubridate::duration object.

scale.factor

numeric vector of length 1, or length equal to that of w.band. Numeric multiplier applied to returned values.

wb.trim

logical if TRUE wavebands crossing spectral data boundaries are trimmed, if FALSE, they are discarded.

use.cached.mult

logical indicating whether multiplier values should be cached between calls.

use.hinges

logical Flag indicating whether to insert "hinges" into the spectral data before integration so as to reduce interpolation errors at the boundaries of the wavebands.

allow.scaled

logical indicating whether scaled or normalized spectra as argument to spct are flagged as an error.

...

other arguments (possibly ignored).

naming

character one of "long", "default", "short" or "none". Used to select the type of names to assign to returned value.

return.tb

logical Flag forcing a tibble to be always returned, even for a single spectrum as argumnet to spct. The default is FALSE for backwards compatibility.

attr2tb

character vector, see add_attr2tb for the syntax for attr2tb passed as is to formal parameter col.names.

idx

character Name of the column with the names of the members of the collection of spectra.

.parallel

if TRUE, apply function in parallel, using parallel backend provided by foreach

.paropts

a list of additional options passed into the foreach function when parallel computation is enabled. This is important if (for example) your code relies on external data or packages: use the .export and .packages arguments to supply them so that all cluster nodes have the correct environment set up for computing.

Value

A named numeric vector in the case of a _spct object containing a single spectrum and return.tb = FALSE. The vector has one member one value for each waveband passed to parameter

w.band. In all other cases a tibble, containing one column for each waveband object, an index column with the names of the spectra, and optionally additional columns with metadata values retrieved from the attributes of the member spectra.

By default values are only integrated, but depending on the argument passed to parameter quantity they can be re-expressed as relative fractions or percentages. In the case of vector output, names attribute is set to the name of the corresponding waveband unless a named list is supplied in which case the names of the list members are used. The time.unit attribute is copied from the spectrum object to the output. Units are as follows: If time.unit is second, [W m-2 nm-1] -> [mol s-1 m-2] If time.unit is day, [J d-1 m-2 nm-1] -> [mol d-1 m-2]

Methods (by class)

  • q_irrad(default): Default for generic function

  • q_irrad(source_spct): Calculates photon irradiance from a source_spct object.

  • q_irrad(source_mspct): Calculates photon (quantum) irradiance from a source_mspct object.

Note

The last two parameters control speed optimizations. The defaults should be suitable in most cases. If you will use repeatedly the same SWFs on many spectra measured at exactly the same wavelengths you may obtain some speed up by setting use.cached.mult=TRUE. However, be aware that you are responsible for ensuring that the wavelengths are the same in each call, as the only test done is for the length of the w.length vector.

See also

Other irradiance functions: e_fluence(), e_irrad(), fluence(), irrad(), q_fluence()

Examples

q_irrad(sun.spct, waveband(c(400,700)))
#> Q_range.400.700 
#>    0.0008941483 
#> attr(,"time.unit")
#> [1] "second"
#> attr(,"radiation.unit")
#> [1] "total photon irradiance"
q_irrad(sun.spct, split_bands(c(400,700), length.out = 3))
#>        Q_wb1        Q_wb2        Q_wb3 
#> 0.0002633524 0.0003141513 0.0003166446 
#> attr(,"time.unit")
#> [1] "second"
#> attr(,"radiation.unit")
#> [1] "total photon irradiance"
q_irrad(sun.spct, split_bands(c(400,700), length.out = 3), quantity = "total")
#>        Q_wb1        Q_wb2        Q_wb3 
#> 0.0002633524 0.0003141513 0.0003166446 
#> attr(,"time.unit")
#> [1] "second"
#> attr(,"radiation.unit")
#> [1] "total photon irradiance"
q_irrad(sun.spct, split_bands(c(400,700), length.out = 3), quantity = "average")
#>    Q(wl)_wb1    Q(wl)_wb2    Q(wl)_wb3 
#> 2.633524e-06 3.141513e-06 3.166446e-06 
#> attr(,"time.unit")
#> [1] "second"
#> attr(,"radiation.unit")
#> [1] "average photon irradiance"
q_irrad(sun.spct, split_bands(c(400,700), length.out = 3), quantity = "relative")
#> Q/Qsum_wb1 Q/Qsum_wb2 Q/Qsum_wb3 
#>  0.2945288  0.3513414  0.3541299 
#> attr(,"time.unit")
#> [1] "second"
#> attr(,"radiation.unit")
#> [1] "relative photon irradiance"
q_irrad(sun.spct, split_bands(c(400,700), length.out = 3), quantity = "relative.pc")
#> Q/Qsum_wb1 Q/Qsum_wb2 Q/Qsum_wb3 
#>   29.45288   35.13414   35.41299 
#> attr(,"time.unit")
#> [1] "second"
#> attr(,"radiation.unit")
#> [1] "relative.pc photon irradiance"
q_irrad(sun.spct, split_bands(c(400,700), length.out = 3), quantity = "contribution")
#> Q/Qtot_wb1 Q/Qtot_wb2 Q/Qtot_wb3 
#>  0.2097834  0.2502492  0.2522353 
#> attr(,"time.unit")
#> [1] "second"
#> attr(,"radiation.unit")
#> [1] "contribution photon irradiance"
q_irrad(sun.spct, split_bands(c(400,700), length.out = 3), quantity = "contribution.pc")
#> Q/Qtot_wb1 Q/Qtot_wb2 Q/Qtot_wb3 
#>   20.97834   25.02492   25.22353 
#> attr(,"time.unit")
#> [1] "second"
#> attr(,"radiation.unit")
#> [1] "contribution.pc photon irradiance"