260 likes | 400 Views
The Best Method of Noise Filtering. Yuri Kalambet, Sergey Maltsev, Ampersand Ltd., Moscow, Russia; Yuri Kozmin , Shemyakin Institute of Bioorganic Chemistry, Moscow, Russia kalambet@ampersand.ru. History: Adaptive peak approximation. Rough slope width estimate.
E N D
The Best Method of Noise Filtering Yuri Kalambet, Sergey Maltsev, Ampersand Ltd., Moscow, Russia; Yuri Kozmin, ShemyakinInstitute of Bioorganic Chemistry, Moscow, Russia kalambet@ampersand.ru
Rough slope width estimate • Evaluate baseline using default gap (minimum peak width Integration parameter) • Evaluate peak height using default gap • Count all points from peak apex to slope end with height bigger than half-height of the peak. Count obtained is an estimate of the slope width.
Properties of adaptive peak approximation • Good noise suppression at each slope • Minimal peak shape disturbances • All peak parameters are resistant to oversampling • Baseline approximation may be poor – either noisy (small gap) or disturbed (large gap). • No approximation outside of peaks • Does not improve formal signal/noise ratio • Baseline position is one of the most important sources of error
Confidence interval estimate where n - number of data points used for polynomial approximation (gap of the filter); p - power of the polynomial; X - matrix of x power values on independent axis (time); Y - vector of detector response values; - Student’s coefficient for confidence probability (1-δ) and m degrees of freedom x* - position at which smoothed (approximated) value is estimated.
G 1 x * G confidence interval 2 Approximation using confidence intervals
Algorithm of simple Confidence filter approximation • Evaluate points and confidence intervals for new (shifted) window
Algorithm of simple Confidence filter approximation • Evaluate points and confidence intervals for new (shifted) window • Compare new confidence interval with that for previously evaluated point. If the new one is smaller than previous, replace approximated point and its confidence interval.
Algorithm of simple Confidence filter approximation • Evaluate points and confidence intervals for new (shifted) window • Compare new confidence interval with that for previously evaluated point. If the new one is smaller than previous, replace approximated point and its confidence interval. • Computational complexity of Confidence filter is comparable to that of simple convolution, (e.g. Savitzky-Golay) and linearly depends on the product gap∙ (degree of the polynomial).
Bonus #1: Correct handling of baseline steps and array boundaries dotted – raw data; thick line – Confidence Filter; thin line – Savitzky-Golay filter
Confidence filter algorithm improvement: Adaptive gap of the polynomial • Repeat confidential filter algorithm for approximations with different windows (gaps) • Computational complexity: degree∙gap∙(gap-1)/2 • Logarithmic step: next gap is k times smaller, than previous, e.g. gap2 = gap1/k, k>1; Computational complexity: degree∙gap∙k/(k-1)
Confidence interval estimate where n - number of data points used for polynomial approximation (gap of the filter); p - power of the polynomial; X - matrix of x power values on independent axis (time); Y - vector of detector response values; - Student’s coefficient for confidence probability (1-δ) and m degrees of freedom x* - position at which smoothed (approximated) value is estimated.
Confidence interval profiles for differentslits (degree = 3)
σ evaluation problems: • Small gaps: accidental perfect fit • Large gaps: treating small peaks as a noise due to large number of degrees of freedom • Is pump pulsation a noise or a signal? • Small gaps: confidence interval depends on confidence level σ evaluation solutions: • Evaluate in advance using the whole data array • Use the estimate for evaluation of confidence intervals
Automatic selection of degree and gap of approximating polynomial
Conclusions: • Confidence filter introduces a measure of approximation quality • Confidence filter helps to select the best set of functions that approximate the data set • Confidence filter is metrologically the best noise filtering method and can be used in the fight with legal metrology Patent pending