bc_spline function

Fit and subtract a smoothing spline to each spectrum in a dataset.