Fast wavelet estimation of weak bio-signals using novel algorithms for generating multiple additional data frames
View Patent ↗A method and apparatus for de-noising weak bio-signals having a relatively low signal to noise ratio utilizes an iterative process of wavelet de-noising a data set comprised of a new set of frames of wavelet coefficients partially generated through a cyclic shift algorithm. The method preferably operates on a data set having 2 N frames, and the iteration is performed N−1 times. The resultant wavelet coefficients are then linearly averaged and an inverse discrete wavelet transform is performed to arrive at the de-noised original signal. The method is preferably carried out in a digital processor.
1. A signal processor for de-noising an input signal, said signal processor including a memory for storing said input signal as an array A1 of a plurality of component frames, and a processor configured to:
(a) perform a wavelet transform on each of said component frames to thereby create a corresponding array of wavelet coefficients;
(b) iteratively for at least one repetition (i) combine multiple numbers of said wavelet coefficients into a plurality of resultant wavelet coefficients forming an array A2 of a like number 2 N of wavelet coefficients, each of said resultant wavelet coefficients being formed from a different combination of wavelet coefficients determined from the original set of wavelet coefficients or the previous iteration, as appropriate; and (ii) wavelet de-noising each resultant wavelet coefficient of array A2; and
(c) combine the resultant wavelet coefficients resulting from the last iteration to form a combined resultant wavelet coefficient.
2. The signal processor of claim 1 wherein said input signal is characterized by a relatively low signal to noise ratio, and wherein said processor is configured to combine different wavelet coefficients to create each resultant wavelet coefficient.
3. The signal processor of claim 2 wherein said processor is configured to use a different threshold for discrete wavelet de-noising in each iteration.
4. The signal processor of claim 2 wherein said processor is configured to combine adjacent pairs of wavelet coefficients in a cyclic manner to create at least some of the resultant wavelet coefficients.
5. The signal processor of claim 2 wherein the processor is configured to combine the final array of iterated wavelet coefficients by averaging them.
6. The signal processor of claim 2 wherein the processor is configured to combine said wavelet coefficients by Euler-Fermat reordering.
7. The signal processor of claim 2 wherein the processor is configured to combine said wavelet coefficients by creating all mathematical permutations thereof to produce N! frames of wavelet coefficients.
8. A method for de-noising an input signal, said method utilizing a signal processor including a memory for storing said input signal as an array A1 of 2 N component frames, said method comprising the steps of:
(a) performing a discrete wavelet transform on each of said component frames to thereby determine a corresponding array of wavelet coefficients;
(b) iteratively for N−1 repetitions (i) combining multiple numbers of the wavelet coefficients determined from the previous step or iteration, as appropriate, into a plurality of resultant wavelet coefficients forming an array A2 of a like number 2 N of wavelet coefficients, each of said resultant wavelet coefficients being formed from a different combination of wavelet coefficients; and (ii) wavelet de-noising each wavelet coefficient of array A2; and
(c) combining the resultant wavelet coefficients resulting from the iteration to form a combined resultant wavelet coefficient; and
(d) performing an inverse discrete wavelet transform on the combined resultant wavelet coefficient to thereby produce the de-noised input signal.
9. The method of claim 8 wherein said input signal is characterized by a relatively low signal to noise ratio, and further comprising the step of combining different wavelet coefficients to create each resultant wavelet coefficient.
10. The method of claim 9 further comprising the step of using a different threshold for discrete wavelet de-noising in each iteration.
11. The method of claim 9 further comprising the step of combining adjacent pairs of wavelet coefficients in a cyclic manner to create at least some of the resultant wavelet coefficients.
12. The method of claim 9 further comprising the step of combining the final array of iterated wavelet coefficients by averaging them.
13. The method of claim 9 further comprising the step of combining said wavelet coefficients by Euler-Fermat reordering.
14. The method of claim 9 further comprising the step of combining said wavelet coefficients by creating all mathematical permutations thereof to produce N! wavelet coefficients.