planck: Calculate the spectral radiance of the Absolutely Black Body

In omega1x/spectrotest: Manipulate with DRIFT-spectroscopy data

Description

Calculate the spectral radiance (power) emitted per unit projected area of the Absolutely Black Body (ABB) as a function of the given wavenumbers according to the Planck's law.

Usage

planck(
  x,
  temperature = y0^(1/3) * 560.323665873058,
  x0 = temperature * 1.9610115857552,
  y0 = (x0 * 0.001)^3 * 0.75377608961114
)

Arguments

x	numeric vector of wavenumbers (1/cm) which the spectral radiance should be calculated at.
temperature	ABB temperature in kelvins given as a numeric vector of length 1, or alternatively.
x0	abscissa of the radiance peak (in 1/cm) given as a numeric vector of length 1, or alternatively.
y0	ordinate of the radiance peak (SI units of power: Wm^(-2)sr^(-1)(1/cm)^(-1)) given as a numeric vector of length 1.

Details

Only one of the next ABB radiance parameters should be provided:

• ABB temperature (temperature) or

• abscissa of the radiance peak (x0) or

• ordinate of the radiance peak (y0).

Value

Calculated spectral radiance (SI units of power: Wm^(-2)sr^(-1)(1/cm)^(-1)) as a numeric vector.

References

Calculating Blackbody Radiance / www.spectralcalc.com - High-resolution spectral modelling: GATS, Inc.

See Also

Other math: pinv()

Examples

 # Typical middle infrared range exposed by FTIR-spectrometers:
 wave_numbers <- seq(from = 349.115696, by = 1.928816, length.out = 3709)

 # radiance parameterized with absolute temperature in kelvins:
 planck(wave_numbers, temperature = 944.2387)

 # radiance parameterized with peak position in 1/cm:
 planck(wave_numbers, x0 = wave_numbers[780])

 # radiance parameterized with peak value in SI-units of radiance:
 planck(wave_numbers, y0 = 4.785513)

omega1x/spectrotest documentation built on Oct. 1, 2020, 4 p.m.