Sophisticated Study of Time, Frequency and Statistical Analysis for Gradient-Switching-Induced Potentials during MRI
Abstract
:1. Introduction
2. State of the Art of Stationaity Test
2.1. The Kwiatkowski–Phillips–Schmidt–Shin Stationarity Test
2.2. Stationarity Test with a Time-Frequency Approach
2.2.1. The Time-Frequency Approach
2.2.2. Surrogates
2.2.3. Distances
2.2.4. Stationarity Test
2.2.5. Index of Non-Stationarity
3. Signal Acquisition and Treatment
3.1. Recording of Induced Potentials
3.2. Pre-Treatment
3.2.1. Normalization
3.2.2. Puff Extraction
3.2.3. Time and Frequency Analysis of Induced Potentials
- -
- RMS values of the global signal and RMS values of the different bursts.
- -
- Estimation of the average curve of the chirps, and calculation of its RMS value.
- -
- Measurement of the similarity between the puffs by calculating the mean square error between each puff and the mean curve according to the following equation:
- -
- The calculation of the power spectral density (PSD) is performed by the Welche–WOSA method, and estimation of the characteristic parameters, the average frequency, the maximum amplitude frequency and the standard deviation of the overall spectrum and on the set of puffs are obtained by segmentation.
3.3. Stationarity Study
3.3.1. Kpss Method
- (a)
- the KPSS of the six 5 s recordings (FSE/axial/coronal/sagittal-CINE/axial/coronal/sagittal).
- (b)
- the KPSS of the 480 extracted puffs and evaluation of the variabilities by estimating the mean values and standard deviation of the obtained series of values. The values were also grouped and graphed to highlight the degree of stationarity or non-stationarity of the different studied segments of the induced potentials.
3.3.2. Surrogate-Based Method
- (1)
- Time-frequency representation: The choice was made for the multi-window Wigner–Ville spectrogram, having successive short-term windows of a Hermite function base. This allows the possibility to adapt the window sizes of our recordings to the MRI sequences.
- (2)
- Surrogate generation: A set of surrogates each having the same power spectral density as the original signal was created. This was achieved by keeping the Fourier transform modulus unchanged but replacing its phase with another randomly taken on .
- (3)
- The stationarity test is based on the distances between the local and global spectra. The distance calculation was carried out by combining the Kullback–Leiber divergence (KL) and log spectral deviation (LSD) methods.
- Number of substitutes: 5000;
- Number of windows: 5;
- Window size range: [0.03:0.04:0.005:0.07:0.075] adjusted for slice orientation.
4. Results and Discussion
4.1. Puff Analysis
4.2. Global and Local Power Spectral Density
4.3. Stationarity
5. Conclusions
Author Contributions
Funding
Institutional Review Board Statement
Informed Consent Statement
Data Availability Statement
Conflicts of Interest
References
- Cosoli, G.; Scalise, L.; De Leo, A.; Russo, P.; Tricarico, G.; Tomasini, E.P.; Cerri, G. Development of a Novel Medical Device for Mucositis and Peri-Implantitis Treatment. Bioengineering 2020, 7, 87. [Google Scholar] [CrossRef] [PubMed]
- Hauser, P.V.; Chang, H.M.; Nishikawa, M.; Kimura, H.; Yanagawa, N.; Hamon, M. Bioprinting scaffolds for vascular tissues and tissue vascularization. Bioengineering 2021, 8, 178. [Google Scholar] [CrossRef] [PubMed]
- Mustahsan, V.M.; Anugu, A.; Komatsu, D.E.; Kao, I.; Pentyala, S. Biocompatible Customized 3D Bone Scaffolds Treated with CRFP, an Osteogenic Peptide. Bioengineering 2021, 8, 199. [Google Scholar] [CrossRef]
- Felblinger, J.; Slotboom, J.; Kreis, R.; Jung, B.; Boesch, C. Restoration of electrophysiological signals distorted by inductive effects of magnetic field gradients during MR sequences. Magn. Reson. Med. Off. J. Int. Soc. Magn. Reson. Med. 1999, 41, 715–721. [Google Scholar] [CrossRef]
- Bresch, E.; Nielsen, J.; Nayak, K.; Narayanan, S. Synchronized and noise-robust audio recordings during realtime magnetic resonance imaging scans. J. Acoust. Soc. Am. 2006, 120, 1791–1794. [Google Scholar] [CrossRef]
- Price, D.L.; De Wilde, J.P.; Papadaki, A.M.; Curran, J.S.; Kitney, R.I. Investigation of acoustic noise on 15 MRI scanners from 0.2 T to 3 T. J. Magn. Reson. Imaging Off. J. Int. Soc. Magn. Reson. Med. 2001, 13, 288–293. [Google Scholar] [CrossRef]
- Casale, R.; De Angelis, R.; Coquelet, N.; Mokhtari, A.; Bali, M.A. The Impact of Edema on MRI Radiomics for the Prediction of Lung Metastasis in Soft Tissue Sarcoma. Diagnostics 2023, 13, 3134. [Google Scholar] [CrossRef]
- Feng, M.; Xu, J. Detection of ASD Children through Deep-Learning Application of fMRI. Children 2023, 10, 1654. [Google Scholar] [CrossRef]
- Von Smekal, A.; Seelos, K.; Küper, C.; Reiser, M. Patient monitoring and safety during MRI examinations. Eur. Radiol. 1995, 5, 302–305. [Google Scholar] [CrossRef]
- Shellock, F.G. Patient monitoring in the MRI environment. Magn. Reson. Proced. Health Eff. Saf. 2001, 217–241. Available online: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC3387746/ (accessed on 3 September 2023).
- Ramírez, W.A.; Gizzi, A.; Sack, K.L.; Filippi, S.; Guccione, J.M.; Hurtado, D.E. On the role of ionic modeling on the signature of cardiac arrhythmias for healthy and diseased hearts. Mathematics 2020, 8, 2242. [Google Scholar] [CrossRef]
- Liu, J.Z.; Zhang, L.; Yao, B.; Yue, G.H. Accessory hardware for neuromuscular measurements during functional MRI experiments. Magn. Reson. Mater. Physics Biol. Med. 2001, 13, 164–171. [Google Scholar] [CrossRef] [PubMed]
- van Duinen, H.; Zijdewind, I.; Hoogduin, H.; Maurits, N. Surface EMG measurements during fMRI at 3T: Accurate EMG recordings after artifact correction. NeuroImage 2005, 27, 240–246. [Google Scholar] [CrossRef] [PubMed]
- Ganesh, G.; Franklin, D.W.; Gassert, R.; Imamizu, H.; Kawato, M. Accurate real-time feedback of surface EMG during fMRI. J. Neurophysiol. 2007, 97, 912–920. [Google Scholar] [CrossRef] [PubMed]
- Van der Meer, J.; Tijssen, M.; Bour, L.; Van Rootselaar, A.; Nederveen, A. Robust EMG–fMRI artifact reduction for motion (FARM). Clin. Neurophysiol. 2010, 121, 766–776. [Google Scholar] [CrossRef] [PubMed]
- Lemieux, L.; Salek-Haddadi, A.; Hoffmann, A.; Gotman, J.; Fish, D.R. EEG-correlated functional MRI: Recent methodologic progress and current issues. Epilepsia 2002, 43, 64–68. [Google Scholar] [CrossRef]
- Allen, P.J.; Josephs, O.; Turner, R. A method for removing imaging artifact from continuous EEG recorded during functional MRI. Neuroimage 2000, 12, 230–239. [Google Scholar] [CrossRef]
- Abächerli, R.; Pasquier, C.; Odille, F.; Kraemer, M.; Schmid, J.J.; Felblinger, J. Suppression of MR gradient artefacts on electrophysiological signals based on an adaptive real-time filter with LMS coefficient updates. Magn. Reson. Mater. Physics. Biol. Med. 2005, 18, 41–50. [Google Scholar] [CrossRef]
- Lou, Q.; Wan, X.; Jia, B.; Song, D.; Qiu, L.; Yin, S. Application Study of Empirical Wavelet Transform in Time–Frequency Analysis of Electromagnetic Radiation Induced by Rock Fracture. Minerals 2022, 12, 1307. [Google Scholar] [CrossRef]
- Xi, Q.; Sahakian, A.V.; Ng, J.; Swiryn, S. Stationarity of surface ECG atrial fibrillatory wave characteristics in the time and frequency domains in clinically stable patients. In Proceedings of the Computers in Cardiology, Thessaloniki, Greece, 21–24 September 2003; pp. 133–136. [Google Scholar]
- Lenka, B. Time-frequency analysis of non-stationary electrocardiogram signals using Hilbert-Huang Transform. In Proceedings of the 2015 International Conference on Communications and Signal Processing (ICCSP), Melmaruvathur, India, 2–4 April 2015; pp. 1156–1159. [Google Scholar]
- Nazmi, N.; Abdul Rahman, M.A.; Yamamoto, S.i.; Ahmad, S.A.; Malarvili, M.; Mazlan, S.A.; Zamzuri, H. Assessment on stationarity of EMG signals with different windows size during isotonic contractions. Appl. Sci. 2017, 7, 1050. [Google Scholar] [CrossRef]
- Kwiatkowski, D.; Phillips, P.C.; Schmidt, P.; Shin, Y. Testing the null hypothesis of stationarity against the alternative of a unit root: How sure are we that economic time series have a unit root? J. Econom. 1992, 54, 159–178. [Google Scholar] [CrossRef]
- Xiao, J.; Borgnat, P.; Flandrin, P. Testing stationarity with time-frequency surrogates. In Proceedings of the 2007 15th European Signal Processing Conference, Poznan, Poland, 3–7 September 2007; pp. 2020–2024. [Google Scholar]
- Xiao, J.; Borgnat, P.; Flandrin, P.; Richard, C. Testing stationarity with surrogates-a one-class SVM approach. In Proceedings of the 2007 IEEE/SP 14th Workshop on Statistical Signal Processing, Madison, WI, USA, 26–29 August 2007; pp. 720–724. [Google Scholar]
- Flandrin, P. Time-Frequency/Time-Scale Analysis; Academic Press: Cambridge, MA, USA, 1998. [Google Scholar]
- Bayram, M.; Baraniuk, R.G. Multiple window time-varying spectrum estimation. Nonlinear Nonstationary Signal Process. 2000, 49, 292–316. [Google Scholar]
- Theiler, J.; Eubank, S.; Longtin, A.; Galdrikian, B.; Farmer, J.D. Testing for nonlinearity in time series: The method of surrogate data. Phys. Nonlinear Phenom. 1992, 58, 77–94. [Google Scholar] [CrossRef]
- Keylock, C. Constrained surrogate time series with preservation of the mean and variance structure. Phys. Rev. 2006, 73, 036707. [Google Scholar] [CrossRef] [PubMed]
- Fokapu, O.; El-Tatar, A. An Experimental Setup to Characterize MR Switched Gradient-Induced Potentials. IEEE Trans. Biomed. Circuits Syst. 2012, 7, 355–362. [Google Scholar] [CrossRef]
FSE | CINE | ||||||
---|---|---|---|---|---|---|---|
Axial | Coronal | Sagital | Axial | Coronal | Sagital | ||
RMS | Global | 1.9022 | 1.2566 | 1.7004 | 1.4715 | 2.2606 | 2.1468 |
[min–max] | [0.8770–1.0017] | [1.0947–1.3463] | [1.2018–1.3360] | [1.3902–1.4610] | [2.1046–2.1855] | [2.0844–2.1394] | |
Mean–stdev | 0.9393–0.0279 | 1.2684–0.0644 | 1.2770–0.0293 | 1.4154–0.0152 | 2.1463–0.0181 | 2.1028–0.0125 | |
MSE | [min–max] | [0.0001–0.0133] | [0.0004–0.0140] | [0.0003–0.0109] | [0.0003–0.0151] | [0.0003–0.0096] | [0.0004–0.0099] |
Mean–stdev | 0.0029–0.0036 | 0.0054–0.0027 | 0.0036–0.0031 | 0.0053–0.0037 | 0.0043–0.0025 | 0.0043–0.0027 |
FSE | CINE | |||||
---|---|---|---|---|---|---|
Frequency | Coronal | Axial | Sagital | Coronal | Axial | Sagital |
[246.163–254.390] | [87.66–97.451] | [121.364–161.663] | [217.309–231.005] | [254.220–295.720] | [230.406–237.216] | |
[280.681–287.959] | [9.765–9.765] | [9.548–10.184] | [234.375–234.375] | [234.375–234.375] | [234.375–234.375] | |
[0.0013–0.0013] | [0.004–0.005] | [0.0033–0.0049] | [0.0043–0.0051] | [0.0013–0.0015] | [0.0079–0.0084] |
FSE | CINE | ||||||
---|---|---|---|---|---|---|---|
Stationarity Test | Coronal | Axial | Sagital | Coronal | Axial | Sagital | |
KPSS test | Statistical value | 0.0678 | 0.0921 | 0.0403 | 0.0098 | 0.0019 | 0.0020 |
Surrogates | Theta | 0.0089 | 0.0048 | 0.0033 | 0.0072 | 0.0032 | 0.0055 |
Threshold | 0.2698 | 0.0994 | 0.4134 | 0.0311 | 0.0573 | 0.0320 | |
INS | 0.0776 | 0.0941 | 0.1208 | 0.0779 | 0.0377 | 0.0652 | |
INS threshold | 1.3457 | 1.3314 | 1.3341 | 1.6109 | 1.5754 | 1.5632 |
KPSS Test | Surrogate Test | |||||
---|---|---|---|---|---|---|
FSE | Statistical Value | Theta | Threshold | INS | INS Threshold | |
Coronal | Mean–stdev | 0.1280–0.1472 | 0.0051–0.0009 | 0.0058–0.005 | 1.2607–0.0984 | 1.5920–0.0413 |
Axial | Mean–stdev | 0.4376–0.1838 | 0.0041–0.0004 | 0.0067–0.0011 | 1.3963–0.4480 | 0.8594–0.3205 |
Sagittal | Mean–stdev | 0.3142–0.1475 | 0.0117–0.0011 | 0.0067–0.0011 | 1.9565–0.5376 | 1.2336–0.2032 |
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content. |
© 2023 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).
Share and Cite
Bouzrara, K.; Fokapu, O.; Fakhfakh, A.; Derbel, F. Sophisticated Study of Time, Frequency and Statistical Analysis for Gradient-Switching-Induced Potentials during MRI. Bioengineering 2023, 10, 1282. https://doi.org/10.3390/bioengineering10111282
Bouzrara K, Fokapu O, Fakhfakh A, Derbel F. Sophisticated Study of Time, Frequency and Statistical Analysis for Gradient-Switching-Induced Potentials during MRI. Bioengineering. 2023; 10(11):1282. https://doi.org/10.3390/bioengineering10111282
Chicago/Turabian StyleBouzrara, Karim, Odette Fokapu, Ahmed Fakhfakh, and Faouzi Derbel. 2023. "Sophisticated Study of Time, Frequency and Statistical Analysis for Gradient-Switching-Induced Potentials during MRI" Bioengineering 10, no. 11: 1282. https://doi.org/10.3390/bioengineering10111282
APA StyleBouzrara, K., Fokapu, O., Fakhfakh, A., & Derbel, F. (2023). Sophisticated Study of Time, Frequency and Statistical Analysis for Gradient-Switching-Induced Potentials during MRI. Bioengineering, 10(11), 1282. https://doi.org/10.3390/bioengineering10111282