ModelEsoSkyCalc

class popurri.telluricutils.ModelEsoSkyCalc(filin, dirout='./', diroutplot='./', tag='', wk='w', fk='ftransmission')

Bases: object

warray

Wavelength [A].

ftransmission, femissionarray

Transmission and emission flux.

extra: Dictionary with keys ‘trans_o3’, ‘trans_ma’,’trans_rs’, ‘trans_ms’

Methods Summary

cut_w(wmin, wmax)

Cut model to a specific wavelength range.

fig_spec([figsize, lisfk, filout, sh, sv, svext])

fig_spec_lines([figsize, filout, sh, sv, ...])

kwargs : dict

fig_spec_transmission([figsize, wmin, wmax, ...])

Plot telluric model transmission spectrum, total and different components.

find_lines()

Find minima in the telluric spectrum.

make_mask([fcut, dw, sv, filout, verbose, svobj])

Make telluric mask.

plot_lines([ax, ylim, fimin, fimaxl, fimaxr])

Plot identified lines.

plot_spec([ax, xunit, color, xlabel, ...])

plot_spec_transmission([ax, xunit, color, ...])

rmv_lines_depth([depthmin])

Find and remove lines with depth less than depthmin (i.e. shallow lines).

rmv_lines_npix([npixmin])

Find and remove lines with less than npixmin pixels.

Methods Documentation

cut_w(wmin, wmax)

Cut model to a specific wavelength range.

fig_spec(figsize=(16, 4), lisfk=None, filout=None, sh=False, sv=True, svext=['pdf'], **kwargs)
fig_spec_lines(figsize=(16, 4), filout=None, sh=False, sv=True, svext=['pdf'], fimin={'color': '0.5', 'zorder': -10}, fimaxl={'alpha': 0.5, 'color': 'b', 'linestyle': 'dashed', 'zorder': -10}, fimaxr={'alpha': 0.5, 'color': 'r', 'linestyle': 'dotted', 'zorder': -10}, **kwargs)
kwargsdict

Only passed to plot_spec().

fig_spec_transmission(figsize=(16, 10), wmin=6270, wmax=6300, color='k', cmap='viridis', alpha=0.8, lwz=3, xunit='A', filout=None, sh=False, sv=True, svext=['pdf'], **kwargs)

Plot telluric model transmission spectrum, total and different components. Left panels are full wavelength range, right panels are a zoom on a smaller range. Top panels show the spectrum, middle panels the difference between the total and each components, and the bottom panels the ratio between the total and each component.

find_lines()

Find minima in the telluric spectrum.

make_mask(fcut=0.999, dw=0.001, sv=True, filout=None, verbose=True, svobj=False)

Make telluric mask.

svojbbool

Save mask in object dictionary (as a Mask object), with fcut as key. Useful if making several masks from the same model and want to compare them.

plot_lines(ax=None, ylim=None, fimin={'color': '0.5', 'zorder': -10}, fimaxl={'alpha': 0.5, 'color': 'b', 'linestyle': 'dashed', 'zorder': -10}, fimaxr={'alpha': 0.5, 'color': 'r', 'linestyle': 'dotted', 'zorder': -10})

Plot identified lines.

It is recommended to do plot_spec() first, get ylimits of tge plots and pass them here with ylim.

plot_spec(ax=None, xunit='A', color='k', xlabel=None, ylabel='Flux', leglabel='', title='', zorder=100, **kwargs)
plot_spec_transmission(ax=None, xunit='A', color='k', alpha=0.8, cmap='viridis', xlabel=None, ylabel='Flux', title='', zorder=100, **kwargs)
rmv_lines_depth(depthmin=0.99995)

Find and remove lines with depth less than depthmin (i.e. shallow lines).

rmv_lines_npix(npixmin=3)

Find and remove lines with less than npixmin pixels.