Automatic background subtraction for multiple spectra.

Hi guys and girls.

First time poster on here, and a total novice at Matlab in general. I have better part of 1300 spectra that I need to do background subtraction on. Now, doing that by hand is going to suck royally, so I wanted to make a Matlab script that would enable me to take all the text files for the spectra and dump them in a folder and have Matlab grab them one by one, do a wavelet denoising to smooth things out nicely, and then subtract the background and save the resulting "tidy" spectrum under the same filename. If someone could help me work out how to do this for a single spectrum, I'm pretty certain that I can smash together the code to just make it into a loop that burns through the entire folder.

I've had a crack at it myself already, but I just don't know enough about these things to do a decent job of it (I just can't get it to be a beautiful curve beneath the data. It starts picking up bits of peaks, etc). If you'd like some of my actual data to see, I've uploaded that, too.

Thanks in advance to any and all kind folks who take a crack at proposing a solution!

Answer 1

[Image Analyst]

See the envelope() function, if you have R2015b or later, sgolayfilt() otherwise.