脳波信号解析に注目したノイズ除去、特徴抽出、実験観測応用を最適化する数理基盤に関する研究

Electroencephalography (EEG) data inevitably contains a large amount of noise particularly from ocular potentials in tasks with eye-movements and eye-blink, known as electrooculography (EOG) artifact, which has been a crucial issue in the braincomputer- interface (BCI) study. The eye-movements and eye-blinks have different time-frequency properties mixing together in EEGs of interest. This time-frequency characteristic has been substantially dealt with past proposed denoising algorithms relying on the consistent assumption based on the single noise component model. However, the traditional model is not simply applicable for biomedical signals consist of multiple signal components, such as weak EEG signals easily recognized as a noise because of the signal amplitude with respect to the EOG signal. In consideration of the realistic signal contamination, we newly designed the EEG-EOG signal contamination model for quantitative validations of the artifact removal from EEGs, and then proposed the two-stage wavelet shrinkage method with the undecimated wavelet decomposition (UDWT), which is suitable for the signal structure. The features of EEG-EOG signal has been extracted with existing decomposition methods known as Principal Component Analysis (PCA), Independent Component Analysis (ICA) based on a consistent assumption of the orthogonality of signal vectors or statistical independence of signal components. In the viewpoint of the signal morphology such as spiking, waves and signal pattern transitions, A systematic decomposition method is proposed to identify the type of signal components or morphology on the basis of sparsity in time-frequency domain. Morphological Component Analysis (MCA) is extended the traditional concept of signal decomposition including Fourier and wavelet transforms and provided a way of reconstruction that guarantees accuracy in reconstruction by using multiple bases being independent of each other and uniqueness representation, called the concept of “dictionary”. MCA is applied to decompose the real EEG signal and clarified the best combination of dictionaries for the purpose. In this proposed semi-realistic biological signal analysis, target EEG data was prepared as mixture signals of artificial eye movements and blinks and iEEG recorded from electrodes embedded into the brain intracranially and then those signals were successfully decomposed into original types by a linear expansion of waveforms such as redundant transforms: UDWT, DCT,LDCT, DST and DIRAC. The result demonstrated that the most suitable combination for EEG data analysis was UDWT, DST and DIRAC to represent the baseline envelop, multi frequency wave forms and spiking activities individually as representative types of EEG morphologies. MCA proposed method is used in negative-going Bereitschaftspotential (BP). It is associated with the preparation and execution of voluntary movement. Thus far, the BP for simple movements involving either the upper or lower body segment has been studied. However, the BP has not yet been recorded during sit-to-stand movements, which use the upper and lower body segments. Electroencephalograms were recorded during movement. To detect the movement of the upper body segment, a gyro sensor was placed on the back, and to detect the movement of the lower body segment, an electromyogram (EMG) electrode was placed on the surface of the hamstrings and quadriceps. Our study revealed that a negative-going BP was evoked around -3 to -2 seconds before the onset of the upper body movement in the sit-to-stand movement in response to the start cue. The BP had a negative peak before the onset of the movement. The potential was followed by premotor positivity, a motor-related potential, and a reafferent potential. The BP for the sit-to-stand movement had a steeper negative slope (-0.8 to -0.001 seconds) just before the onset of the upper body movement. The slope correlated with the gyro peak and the max amplitude of hamstrings EMG. A BP negative peak value was correlated with the max amplitude of the hamstring EMG. These results suggested that the observed BP is involved in the preparation/execution for a sit-to-stand movement using the upper and lower body. In summary, this thesis is help to pave the practical approach of real time analysis of desired EEG signal of interest toward the implementation of rehabilitation device which may be used for motor disabled people. We also pointed out the EEG-EOG contamination model that helps in removal of the artifacts and explicit dictionaries are representing the EEG morphologies.九州工業大学博士学位論文 学位記番号:生工博甲第290号 学位授与年月日:平成29年3月24日1 Introduction|2 Research Background and Preliminaries|3 Introduction of Morphological Component Analysis|4 Two-Stage Undecimated Wavelet Shrinkage Method|5 Morphologically Decomposition of EEG Signals|6 Bereitschaftspotential for Rise to Stand-Up Behavior九州工業大学平成28年