Improve threshold determination: (1) Choose automatic thresholds by fitting a
mixture of univariate Gaussians to the recovered amplitudes (currently uses
matlab's function, which chooses bandwidth automatically); (2) allow multiple
thresholds, so that similar waveforms with different amplitudes are handled
correctly - iterating on this will require the operator-slitting method mentioned
below,