barak.spec.combine¶

barak.spec.combine(spectra, cliphi=None, cliplo=None, verbose=False)[source]¶

Combine spectra pixel by pixel, weighting by the inverse variance of each pixel. Clip high sigma values by sigma times clip values Returns the combined spectrum.

If the wavelength scales of the input spectra differ, combine() will rebin the spectra to a common linear (not log-linear) wavelength scale, with pixel width equal to the largest pixel width in the input spectra. If this is not what you want, rebin the spectra by hand with rebin() before using combine().