blind separation of fetal ecg from single mixture using ...changec/publications/fetalslide.pdf ·...
TRANSCRIPT
![Page 1: Blind Separation of Fetal ECG from Single Mixture using ...changec/publications/FetalSlide.pdf · One Approach 1. Find the occurrences of maternal heartbeat by identifying the peak](https://reader036.vdocument.in/reader036/viewer/2022081402/5f0ee1777e708231d44164a4/html5/thumbnails/1.jpg)
Blind Separation of Fetal ECG from Single Mixture using
SVD and ICA
Ping GaoDepartment of Computational Science
National University of Singapore
Ee-Chien Chang School of Computing
National University of Singapore
Lonce Wyse Institute for Infocomm Research
![Page 2: Blind Separation of Fetal ECG from Single Mixture using ...changec/publications/FetalSlide.pdf · One Approach 1. Find the occurrences of maternal heartbeat by identifying the peak](https://reader036.vdocument.in/reader036/viewer/2022081402/5f0ee1777e708231d44164a4/html5/thumbnails/2.jpg)
Mixed ECG’s for a period of roughly 4.6 seconds
time
![Page 3: Blind Separation of Fetal ECG from Single Mixture using ...changec/publications/FetalSlide.pdf · One Approach 1. Find the occurrences of maternal heartbeat by identifying the peak](https://reader036.vdocument.in/reader036/viewer/2022081402/5f0ee1777e708231d44164a4/html5/thumbnails/3.jpg)
Maternal heartbeat (R)
![Page 4: Blind Separation of Fetal ECG from Single Mixture using ...changec/publications/FetalSlide.pdf · One Approach 1. Find the occurrences of maternal heartbeat by identifying the peak](https://reader036.vdocument.in/reader036/viewer/2022081402/5f0ee1777e708231d44164a4/html5/thumbnails/4.jpg)
Fetal heartbeat
![Page 5: Blind Separation of Fetal ECG from Single Mixture using ...changec/publications/FetalSlide.pdf · One Approach 1. Find the occurrences of maternal heartbeat by identifying the peak](https://reader036.vdocument.in/reader036/viewer/2022081402/5f0ee1777e708231d44164a4/html5/thumbnails/5.jpg)
![Page 6: Blind Separation of Fetal ECG from Single Mixture using ...changec/publications/FetalSlide.pdf · One Approach 1. Find the occurrences of maternal heartbeat by identifying the peak](https://reader036.vdocument.in/reader036/viewer/2022081402/5f0ee1777e708231d44164a4/html5/thumbnails/6.jpg)
![Page 7: Blind Separation of Fetal ECG from Single Mixture using ...changec/publications/FetalSlide.pdf · One Approach 1. Find the occurrences of maternal heartbeat by identifying the peak](https://reader036.vdocument.in/reader036/viewer/2022081402/5f0ee1777e708231d44164a4/html5/thumbnails/7.jpg)
108.raw
Separation
Goal: separation from a single mixture
fetal ECG complex
maternal ECG complex
mixture
![Page 8: Blind Separation of Fetal ECG from Single Mixture using ...changec/publications/FetalSlide.pdf · One Approach 1. Find the occurrences of maternal heartbeat by identifying the peak](https://reader036.vdocument.in/reader036/viewer/2022081402/5f0ee1777e708231d44164a4/html5/thumbnails/8.jpg)
One Approach
1. Find the occurrences of maternal heartbeat by identifying the peak
2. Find the maternal ECG complex by “averaging”3. Subtract the maternal ECG complex from the
mixture.4. Repeat the above for fetal ECG.
Disadvantage: require significant fine-tuning in step1 and 3.
![Page 9: Blind Separation of Fetal ECG from Single Mixture using ...changec/publications/FetalSlide.pdf · One Approach 1. Find the occurrences of maternal heartbeat by identifying the peak](https://reader036.vdocument.in/reader036/viewer/2022081402/5f0ee1777e708231d44164a4/html5/thumbnails/9.jpg)
Main Idea
• Identify the heart beat in the spectrogram.• Observation: with the right window size used inthe spectrogram, a ECG complex in the spectrogramcan be viewed as a separable function.
f(s,t) = f(s) f(t)
![Page 10: Blind Separation of Fetal ECG from Single Mixture using ...changec/publications/FetalSlide.pdf · One Approach 1. Find the occurrences of maternal heartbeat by identifying the peak](https://reader036.vdocument.in/reader036/viewer/2022081402/5f0ee1777e708231d44164a4/html5/thumbnails/10.jpg)
Main Idea
• Identify the heart beat in the spectrogram.• Observation: with the right window size used inthe spectrogram, a ECG complex in the spectrogramcan be viewed as a separable function.
f(s,t) = f(s) f(t)
![Page 11: Blind Separation of Fetal ECG from Single Mixture using ...changec/publications/FetalSlide.pdf · One Approach 1. Find the occurrences of maternal heartbeat by identifying the peak](https://reader036.vdocument.in/reader036/viewer/2022081402/5f0ee1777e708231d44164a4/html5/thumbnails/11.jpg)
Main Idea
• Identify the heart beat in the spectrogram.• Observation: with the right window size used inthe spectrogram, a ECG complex in the spectrogramcan be viewed as a separable function.
f(s,t) = f(s) f(t) f(s)
f(t)
![Page 12: Blind Separation of Fetal ECG from Single Mixture using ...changec/publications/FetalSlide.pdf · One Approach 1. Find the occurrences of maternal heartbeat by identifying the peak](https://reader036.vdocument.in/reader036/viewer/2022081402/5f0ee1777e708231d44164a4/html5/thumbnails/12.jpg)
Main Idea
• Identify the heart beat in the spectrogram.• Observation: with the right window size used inthe spectrogram, a ECG complex in the spectrogramcan be viewed as a separable function.
f(s,t) = f(s) f(t)
![Page 13: Blind Separation of Fetal ECG from Single Mixture using ...changec/publications/FetalSlide.pdf · One Approach 1. Find the occurrences of maternal heartbeat by identifying the peak](https://reader036.vdocument.in/reader036/viewer/2022081402/5f0ee1777e708231d44164a4/html5/thumbnails/13.jpg)
Main Idea
• Identify the heart beat in the spectrogram.• Observation: with the right window size used inthe spectrogram, a ECG complex in the spectrogramcan be viewed as a separable function.
f(s,t) = f(s) f(t)
![Page 14: Blind Separation of Fetal ECG from Single Mixture using ...changec/publications/FetalSlide.pdf · One Approach 1. Find the occurrences of maternal heartbeat by identifying the peak](https://reader036.vdocument.in/reader036/viewer/2022081402/5f0ee1777e708231d44164a4/html5/thumbnails/14.jpg)
S = Um Vm + Uf Vf
= +
* *
we want to find the “best” Um, Vm , Uf and Vf s.t.
![Page 15: Blind Separation of Fetal ECG from Single Mixture using ...changec/publications/FetalSlide.pdf · One Approach 1. Find the occurrences of maternal heartbeat by identifying the peak](https://reader036.vdocument.in/reader036/viewer/2022081402/5f0ee1777e708231d44164a4/html5/thumbnails/15.jpg)
Suppose we want to find the solutions that minimized the noise in the sense that || N ||2 is minimized,
S = Um Vm + Uf Vf + N
then, we can employ the well-known SVD. However, experimental studies show that it gives unsatisfactory results.
We borrow idea of ICA (Independent Component Analysis).In the proposed algorithm, we attempt to find the solutions that are “statistically independent”and non-Gaussian like.
![Page 16: Blind Separation of Fetal ECG from Single Mixture using ...changec/publications/FetalSlide.pdf · One Approach 1. Find the occurrences of maternal heartbeat by identifying the peak](https://reader036.vdocument.in/reader036/viewer/2022081402/5f0ee1777e708231d44164a4/html5/thumbnails/16.jpg)
Experiments
292.raw
Maternal heartbeat trend Um using ica+svd
![Page 17: Blind Separation of Fetal ECG from Single Mixture using ...changec/publications/FetalSlide.pdf · One Approach 1. Find the occurrences of maternal heartbeat by identifying the peak](https://reader036.vdocument.in/reader036/viewer/2022081402/5f0ee1777e708231d44164a4/html5/thumbnails/17.jpg)
Experiments
![Page 18: Blind Separation of Fetal ECG from Single Mixture using ...changec/publications/FetalSlide.pdf · One Approach 1. Find the occurrences of maternal heartbeat by identifying the peak](https://reader036.vdocument.in/reader036/viewer/2022081402/5f0ee1777e708231d44164a4/html5/thumbnails/18.jpg)
Experiments
108.raw
Fetal heartbeat trend Uv using SVD
![Page 19: Blind Separation of Fetal ECG from Single Mixture using ...changec/publications/FetalSlide.pdf · One Approach 1. Find the occurrences of maternal heartbeat by identifying the peak](https://reader036.vdocument.in/reader036/viewer/2022081402/5f0ee1777e708231d44164a4/html5/thumbnails/19.jpg)
Experiments
Fetal heartbeat trend Uv using svd+ica
![Page 20: Blind Separation of Fetal ECG from Single Mixture using ...changec/publications/FetalSlide.pdf · One Approach 1. Find the occurrences of maternal heartbeat by identifying the peak](https://reader036.vdocument.in/reader036/viewer/2022081402/5f0ee1777e708231d44164a4/html5/thumbnails/20.jpg)
Independent Components Analysis
x1 (t) = a11 s1 (t) + a12 s2(t)x2 (t) = a21 s1 (t) + a22 s2(t)
We want to find the s1 and s2 which hasmaximum independence and minimum nongaussianity.
∑ s1(t) s2 (t) = ∑ s1(t) ∑ s2(t)Negentropy
unknowns > equations, under constrained.
In our experiment, we use FastICA
![Page 21: Blind Separation of Fetal ECG from Single Mixture using ...changec/publications/FetalSlide.pdf · One Approach 1. Find the occurrences of maternal heartbeat by identifying the peak](https://reader036.vdocument.in/reader036/viewer/2022081402/5f0ee1777e708231d44164a4/html5/thumbnails/21.jpg)
S = Um Vm + Uf Vf
= +
* *
we want to find the “independent” Um, Vm , Uf and Vf s.t.
![Page 22: Blind Separation of Fetal ECG from Single Mixture using ...changec/publications/FetalSlide.pdf · One Approach 1. Find the occurrences of maternal heartbeat by identifying the peak](https://reader036.vdocument.in/reader036/viewer/2022081402/5f0ee1777e708231d44164a4/html5/thumbnails/22.jpg)
Proposed method1. Compute Spectrogram S2. Perform SVD on S, S=U T Vt
3. Apply ICA on the k most significant spectral components, i.e. on V1, V2, V3, ..., Vk.
4. Update the U using the mixture obtained in step 3.
5. Apply ICA on the k most significant spectral of the updated U.
6. Choose two ``best’’ components in U and call them um and uf.
S = U T
Vt
![Page 23: Blind Separation of Fetal ECG from Single Mixture using ...changec/publications/FetalSlide.pdf · One Approach 1. Find the occurrences of maternal heartbeat by identifying the peak](https://reader036.vdocument.in/reader036/viewer/2022081402/5f0ee1777e708231d44164a4/html5/thumbnails/23.jpg)
Proposed method1. Compute Spectrogram S2. Perform SVD on S, S=U T Vt
3. Apply ICA on the k most significant spectral components, i.e. on V1, V2, V3, ..., Vk.
4. Update the U using the mixture obtained in step 3.
5. Apply ICA on the k most significant spectral of the updated U.
6. Choose two ``best’’ components in U and call them um and uf.
S = U T
VtICA
![Page 24: Blind Separation of Fetal ECG from Single Mixture using ...changec/publications/FetalSlide.pdf · One Approach 1. Find the occurrences of maternal heartbeat by identifying the peak](https://reader036.vdocument.in/reader036/viewer/2022081402/5f0ee1777e708231d44164a4/html5/thumbnails/24.jpg)
Proposed method1. Compute Spectrogram S2. Perform SVD on S, S=U T Vt
3. Apply ICA on the k most significant spectral components, i.e. on V1, V2, V3, ..., Vk.
4. Update the U using the mixture obtained in step 3.
5. Apply ICA on the k most significant spectral of the updated U.
6. Choose two ``best’’ components in U and call them um and uf.
S = U T
Vt~~ ~
![Page 25: Blind Separation of Fetal ECG from Single Mixture using ...changec/publications/FetalSlide.pdf · One Approach 1. Find the occurrences of maternal heartbeat by identifying the peak](https://reader036.vdocument.in/reader036/viewer/2022081402/5f0ee1777e708231d44164a4/html5/thumbnails/25.jpg)
Proposed method1. Compute Spectrogram S2. Perform SVD on S, S=U T Vt
3. Apply ICA on the k most significant spectral components, i.e. on V1, V2, V3, ..., Vk.
4. Update the U using the mixture obtained in step 3.
5. Apply ICA on the k most significant spectral of the updated U.
6. Choose two ``best’’ components in U and call them um and uf.
S = U T
Vt~~ ~
ICA
![Page 26: Blind Separation of Fetal ECG from Single Mixture using ...changec/publications/FetalSlide.pdf · One Approach 1. Find the occurrences of maternal heartbeat by identifying the peak](https://reader036.vdocument.in/reader036/viewer/2022081402/5f0ee1777e708231d44164a4/html5/thumbnails/26.jpg)
Components obtained after SVD
Components obtained after SVD + ICA
![Page 27: Blind Separation of Fetal ECG from Single Mixture using ...changec/publications/FetalSlide.pdf · One Approach 1. Find the occurrences of maternal heartbeat by identifying the peak](https://reader036.vdocument.in/reader036/viewer/2022081402/5f0ee1777e708231d44164a4/html5/thumbnails/27.jpg)
Remarks & Future works
1. The proposed method does not find the intendeduf and um as stated in the formulation. It is an approximation.
2. Extend the method to other applications.
3. More analysis.
![Page 28: Blind Separation of Fetal ECG from Single Mixture using ...changec/publications/FetalSlide.pdf · One Approach 1. Find the occurrences of maternal heartbeat by identifying the peak](https://reader036.vdocument.in/reader036/viewer/2022081402/5f0ee1777e708231d44164a4/html5/thumbnails/28.jpg)
Acknowledgement
We would like to thank Prof Ho Ting Fei, Department of Physiology, National University of Singapore, for providing the data-sets.