[SciPy-User] Fitting Gaussian in spectra

Joe Philip Ninan indiajoe@gmail....
Fri Sep 28 13:45:21 CDT 2012


Hi,
I have a spectra with multiple gaussian emission lines over a noisy
continuum.
My primary objective is to find areas under all the gaussian peaks.
For that, the following is the algorithm i have in mind.
1) fit the continuum and subtract it.
2) find the peaks
3) do least square fit of gaussian at the peaks to find the area under each
gaussian peaks.
I am basically stuck at the first step itself. Simple 2nd or 3rd order
polynomial fit is not working because the contribution from peaks are
significant. Any tool exist to fit continuum ignoring the peaks?
For finding peaks, i tried find_peaks_cwt in signal module of scipy. But it
seems to be quite sensitive of the width of peak and was picking up
non-existing peaks also.
The wavelet used was default mexican hat. Is there any better wavelet i
should try?

Or is there any other module in python/scipy which i should give a try?
Thanking you.
-cheers
joe
-- 
/---------------------------------------------------------------
"GNU/Linux: because a PC is a terrible thing to waste" -  GNU Generation
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mail.scipy.org/pipermail/scipy-user/attachments/20120929/dc8784c7/attachment-0001.html 


More information about the SciPy-User mailing list