radiance2Kelvin
index
../src3/radiance2Kelvin.py

radiance2Kelvin
convert radiance (intensity) to equivalent blackbody temperatures according to Planck
 
usage:
  radiance2Kelvin [options] file[s]
 
options:
  -c      char     comment character(s) used in files (default '#', several characters allowed)
  -C      int      number of y column to be processed (default: all y columns)
  -h               help
  -o      file     output (print) file (default: standard out)
  -P               plot:  quicklook using matplotlib
  -x      string   unit of x-axis (wavenumber, frequency or wavelength): 'cm-1' (default), 'm-1', 'Hz', 'MHz', 'GHz', 'mue'
  -y      string   unit of y-axis (radiance): 'erg/s' for erg/s/(cm^2 sr cm^-1) (default) or 'W' for W(cm^2 sr cm^-1)
 -xFormat string   format to be used for x values (default:'% 12.6f')
 -yFormat string   format to be used for y values (default:'% 10.4f')

 
Modules
       
numpy

 
Functions
       
ergs_to_Kelvin(x, y)
Convert radiance in erg/s/(cm^2 sr cm^-1) to BlackBody temperature via inverse Planck.
ergs_to_Kelvin_frequency(x, y, xUnit='Hz')
Convert radiance in erg/s/(cm^2 sr Hz) to BlackBody temperature via inverse Planck.
ergs_to_Kelvin_wavelength(x, y, xUnit)
Convert radiance in erg/s/(cm^2 sr micrometer) to BlackBody temperature via inverse Planck.
parse_radiance_file_header(comments, commentChar='#')
Read the header of the radiance file and try to find x and y units.
radiance2Kelvin(xGrid, yValues, xUnit='cm-1', yUnit='erg/s/(cm2 sr cm-1)', flux=False)
Convert radiance spectra to equivalent Black Body Temperature by inversion of Planck's Law.
 
ARGUMENTS:
----------
xGrid     wavenumber|frequency|wavelength
yValues   radiance value(s)
xUnit     physical unit of the xGrid data  [default cm-1]
yUnit     physical unit of the radiances   [default erg/s/(cm2 sr cm-1)]
flux      flag, assume yValues is flux, so divide by pi
 
RETURNS:
--------
temperature(s)  [Kelvin]

 
Data
        C2 = 1.4387751601679206
c = 29979245800.0
h = 6.62606896e-27
hcc2 = 1.1910427584934559e-05
k = 1.3806504e-16
pi = 3.141592653589793
xUnitsValid = ['cm-1', 'Hz', 'kHz', 'THz', 'GHz', 'MHz', 'mue', 'nm']