KR101380964B1 - Apparatus of brain-computer interface, and classification method thereof - Google Patents

Apparatus of brain-computer interface, and classification method thereof Download PDF

Info

Publication number
KR101380964B1
KR101380964B1 KR1020120039497A KR20120039497A KR101380964B1 KR 101380964 B1 KR101380964 B1 KR 101380964B1 KR 1020120039497 A KR1020120039497 A KR 1020120039497A KR 20120039497 A KR20120039497 A KR 20120039497A KR 101380964 B1 KR101380964 B1 KR 101380964B1
Authority
KR
South Korea
Prior art keywords
class
signal
dictionary
solution
test signal
Prior art date
Application number
KR1020120039497A
Other languages
Korean (ko)
Other versions
KR20120125948A (en
Inventor
이흥노
신영학
이승찬
Original Assignee
광주과학기술원
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by 광주과학기술원 filed Critical 광주과학기술원
Priority to KR1020120039497A priority Critical patent/KR101380964B1/en
Priority to PCT/KR2012/003572 priority patent/WO2012153965A2/en
Publication of KR20120125948A publication Critical patent/KR20120125948A/en
Application granted granted Critical
Publication of KR101380964B1 publication Critical patent/KR101380964B1/en

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235Details of waveform analysis
    • A61B5/7264Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/16Devices for psychotechnics; Testing reaction times ; Devices for evaluating the psychological state
    • A61B5/18Devices for psychotechnics; Testing reaction times ; Devices for evaluating the psychological state for vehicle drivers or machine operators
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/24Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/24Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/316Modalities, i.e. specific diagnostic methods
    • A61B5/369Electroencephalography [EEG]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F3/00Input arrangements for transferring data to be processed into a form capable of being handled by the computer; Output arrangements for transferring data from processing unit to output unit, e.g. interface arrangements
    • G06F3/01Input arrangements or combined input and output arrangements for interaction between user and computer
    • G06F3/011Arrangements for interaction with the human body, e.g. for user immersion in virtual reality
    • G06F3/015Input arrangements based on nervous system activity detection, e.g. brain waves [EEG] detection, electromyograms [EMG] detection, electrodermal response detection
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2218/00Aspects of pattern recognition specially adapted for signal processing
    • G06F2218/12Classification; Matching
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/20ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for computer-aided diagnosis, e.g. based on medical expert systems

Landscapes

  • Health & Medical Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Biomedical Technology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Molecular Biology (AREA)
  • Pathology (AREA)
  • Veterinary Medicine (AREA)
  • Animal Behavior & Ethology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Surgery (AREA)
  • Biophysics (AREA)
  • Medical Informatics (AREA)
  • Psychiatry (AREA)
  • General Engineering & Computer Science (AREA)
  • Psychology (AREA)
  • Theoretical Computer Science (AREA)
  • Artificial Intelligence (AREA)
  • Hospice & Palliative Care (AREA)
  • Human Computer Interaction (AREA)
  • Social Psychology (AREA)
  • Educational Technology (AREA)
  • Developmental Disabilities (AREA)
  • Child & Adolescent Psychology (AREA)
  • General Physics & Mathematics (AREA)
  • Dermatology (AREA)
  • Neurosurgery (AREA)
  • Neurology (AREA)
  • Evolutionary Computation (AREA)
  • Fuzzy Systems (AREA)
  • Mathematical Physics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Physiology (AREA)
  • Signal Processing (AREA)
  • Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
  • Measuring And Recording Apparatus For Diagnosis (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)

Abstract

뇌파 신호인 테스트 신호를 복수의 클래스 중 어느 하나로 분류하는 뇌-컴퓨터 접속 장치로서, 테스트 신호를 분류하기 위한 사전을 각 클래스의 훈련 신호를 기초로 설계하는 사전부, 사전을 기초로 테스트 신호의 해를 구하는 풀이부, 그리고 해의 각 클래스에 나타난 0이 아닌 계수의 분포를 기초로 테스트 신호를 분류하는 분류부를 포함한다.A brain-computer connection device for classifying test signals, which are brain wave signals, into any one of a plurality of classes, comprising: a dictionary unit for designing a dictionary for classifying test signals based on training signals of each class, and a solution of test signals based on a dictionary And a classifier for classifying the test signal based on the distribution of nonzero coefficients shown in each class of the solution.

Description

뇌-컴퓨터 접속 장치, 그리고 그의 분류 방법{APPARATUS OF BRAIN-COMPUTER INTERFACE, AND CLASSIFICATION METHOD THEREOF}Brain-computer connection device and its classification method {APPARATUS OF BRAIN-COMPUTER INTERFACE, AND CLASSIFICATION METHOD THEREOF}

본 발명은 뇌-컴퓨터 접속 장치, 그리고 그의 분류 방법에 관한 것이다.The present invention relates to a brain-computer connection device and a classification method thereof.

지금까지의 연구들에서 두피로부터 측정되는 뇌파(electroencephalogram, EEG)는 자발적인 움직임이 거의 없더라도 의사전달이 가능하다는 것을 보여주었다. 이것을 일반적으로 뇌파기반 뇌-컴퓨터 접속(brain-computer interface, BCI) 시스템이라고 부른다. 이 BCI시스템은 뇌파활동의 특정한 특징 신호를 추출하고 이 신호를 이용하여 컨트롤 신호를 만들어낸다. 움직임 상상 기반의 BCI(motor imagery-based BCI) 시스템은 두피의 감각운동 피질 부분의 Mu(8-14Hz)와 Beta(15-30Hz) 리듬을 포함하는 감각운동 리듬(sensorimotor rhythms, SMRs)을 사용한다. 피실험자가 왼쪽 또는 오른쪽 손 움직임을 상상했을 때, Mu 리듬 주파수 대역의 감쇄와 ERD(event related desynchronisation)과 같은 뇌파신호의 활동이 상상하는 손의 반대쪽 감각운동피질 부분에서 나타난다.To date, studies have shown that electroencephalograms (EEGs) measured from the scalp can communicate even with little spontaneous movement. This is commonly called an EEG-based brain-computer interface (BCI) system. The BCI system extracts specific characteristic signals of EEG activity and uses them to generate control signals. Motion imagery-based BCI (BCI) systems use sensorimotor rhythms (SMRs) that include Mu (8-14 Hz) and Beta (15-30 Hz) rhythms in the sensorimotor cortex of the scalp. . When the subject imagines left or right hand movements, the activity of EEG signals such as attenuation of the Mu rhythm frequency band and event related desynchronisation (ERD) occurs in the part of the sensory cortex opposite the imaginary hand.

BCI 시스템은 피실험자로부터 감지된 신호특징을 명령신호로 변환(translation)하는데, 이 과정은 분류 알고리즘을 통하여 수행된다. 그러나, BCI 시스템이 이용하는 움직임 상상 신호는 잡음이 매우 심하고 피실험자마다 다르기 때문에 이 신호를 분류하는 것은 쉽지 않다.The BCI system translates the signal feature detected from the test subject into a command signal, which is performed through a classification algorithm. However, it is not easy to classify the motion imaginary signal used by the BCI system because it is very noisy and varies from subject to subject.

본 발명이 해결하고자 하는 과제는 희소 표현 기반으로 움직임 상상 뇌파 신호를 분류하는 뇌-컴퓨터 접속 장치, 그리고 그의 분류 방법을 제공하는 것이다.The problem to be solved by the present invention is to provide a brain-computer connection device for classifying a motion imaginary brain wave signal based on the sparse expression, and a classification method thereof.

본 발명의 한 실시예에 따른 뇌파 신호인 테스트 신호를 복수의 클래스 중 어느 하나로 분류하는 뇌-컴퓨터 접속 장치로서, 상기 테스트 신호를 분류하기 위한 사전(dictionary)을 각 클래스의 훈련 신호를 기초로 설계하는 사전부, 상기 사전을 기초로 상기 테스트 신호의 해를 구하는 풀이부, 그리고 상기 해의 각 클래스에 나타난 0이 아닌 계수의 분포를 기초로 상기 테스트 신호를 분류하는 분류부를 포함한다.A brain-computer connection device for classifying a test signal, which is an EEG signal, according to an embodiment of the present invention into any one of a plurality of classes, wherein a dictionary for classifying the test signal is designed based on a training signal of each class And a classification unit for classifying the test signal based on a distribution of non-zero coefficients shown in each class of the solution.

상기 뇌-컴퓨터 접속 장치는 상기 훈련 신호를 하나의 클래스에 대해서는 최대화하고, 다른 클래스에 대해서는 최소화하는 CSP 필터링을 하여 CSP 필터링된 신호를 출력하는 CSP 필터링부를 더 포함하고, 상기 사전부는 상기 CSP 필터링된 신호를 기초로 상기 사전을 설계할 수 있다.The brain-computer connection device may further include a CSP filtering unit configured to output a CSP filtered signal by performing CSP filtering to maximize the training signal for one class and minimize the other class, and the dictionary unit is configured to filter the CSP. The dictionary can be designed based on the signal.

상기 사전부는 상기 복수의 클래스 각각에 대응하는 복수의 행렬로 구성되고, 행의 수는 CSP 필터의 수이며, 각 클래스의 열의 수는 각 클래스의 훈련 신호의 수인 상기 사전을 설계할 수 있다.The dictionary unit may be configured of a plurality of matrices corresponding to each of the plurality of classes, the number of rows is the number of CSP filters, and the number of columns of each class is the number of training signals of each class.

상기 뇌-컴퓨터 접속 장치는 상기 CSP 필터링된 신호 중에서 특정 주파수 대역의 신호를 특징 신호로 추출하는 특징신호추출부를 더 포함하고, 상기 사전부는 상기 특징 신호를 기초로 상기 사전을 설계할 수 있다.The brain-computer connection apparatus may further include a feature signal extracting unit extracting a signal of a specific frequency band from the CSP filtered signal as a feature signal, and the dictionary unit may design the dictionary based on the feature signal.

상기 특정 주파수 대역은 두피의 감각운동 피질 부분의 감각운동 리듬에 해당하는 주파수 대역일 수 있다.The specific frequency band may be a frequency band corresponding to the sensory rhythm of the sensorimotor cortex portion of the scalp.

상기 풀이부는 L1 최소화 방법으로 상기 해를 구할 수 있다.The solver can obtain the solution by L1 minimization.

상기 분류부는 상기 해의 각 클래스 중에서 0이 아닌 계수가 가장 많이 나타난 클래스로 상기 테스트 신호를 분류할 수 있다.The classifying unit may classify the test signal into a class in which a non-zero coefficient is highest among each class of the solution.

상기 분류부는 상기 해를 기초로 클래스별 나머지를 계산하여 나머지가 최소인 클래스로 상기 테스트 신호를 분류하고, 상기 클래스별 나머지는 상기 해 중에서 자신 이외의 클래스의 원소를 0으로 채워 생성한 벡터를 상기 사전에 투영한 제1값과 상기 테스트 신호에 해당하는 제2값의 차이로 계산되는 값일 수 있다.The classification unit calculates a residual for each class based on the solution to classify the test signal into a class having a minimum residual, and the residual for each class includes a vector generated by filling an element of a class other than itself with zero in the solution. It may be a value calculated as a difference between a first value projected in advance and a second value corresponding to the test signal.

상기 훈련 신호와 상기 테스트 신호는 움직임 상상 기반 뇌파 신호일 수 있다.The training signal and the test signal may be a motion imagination based EEG signal.

상기 복수의 클래스는 왼쪽 움직임 상상에 해당하는 왼쪽 클래스와 오른쪽 움직임 상상에 해당하는 오른쪽 클래스를 포함할 수 있다.The plurality of classes may include a left class corresponding to a left movement imagination and a right class corresponding to a right movement imagination.

본 발명의 다른 실시예에 따른 뇌-컴퓨터 접속 장치가 뇌파 신호인 테스트 신호를 복수의 클래스 중 어느 하나로 분류하는 방법으로서, 각 클래스에 해당하는 훈련 신호를 획득하는 단계, 상기 훈련 신호를 하나의 클래스에 대해서는 최대화하고, 다른 클래스에 대해서는 최소화하는 CSP 필터링하여 CSP 필터링한 신호를 생성하는 단계, 상기 CSP 필터링한 신호를 기초로 상기 테스트 신호의 해를 풀기 위한 사전(dictionary)을 설계하는 단계, 상기 사전을 기초로 상기 테스트 신호의 해를 구하는 단계, 그리고 상기 해의 각 클래스에 나타난 0이 아닌 계수의 분포를 기초로 상기 테스트 신호를 분류하는 단계를 포함한다.A method of classifying a test signal as an EEG signal into any one of a plurality of classes by a brain-computer connection device according to another embodiment of the present invention, the method comprising: acquiring a training signal corresponding to each class, and classifying the training signal into one class Generating a CSP filtered signal by maximizing the CSP and minimizing the other classes by using the CSP filtering signal; designing a dictionary for solving the test signal based on the CSP filtered signal; Obtaining a solution of the test signal on the basis of and classifying the test signal based on a distribution of nonzero coefficients represented in each class of the solution.

상기 훈련 신호와 상기 테스트 신호는 움직임 상상 기반 뇌파 신호일 수 있다.The training signal and the test signal may be a motion imagination based EEG signal.

상기 복수의 클래스는 왼쪽 움직임 상상에 해당하는 왼쪽 클래스와 오른쪽 움직임 상상에 해당하는 오른쪽 클래스를 포함할 수 있다.The plurality of classes may include a left class corresponding to a left movement imagination and a right class corresponding to a right movement imagination.

상기 사전을 설계하는 단계는 상기 복수의 클래스 각각에 대응하는 복수의 행렬로 구성되고, 행의 수는 CSP 필터의 수이며, 각 클래스의 열의 수는 각 클래스의 훈련 신호의 수인 상기 사전을 설계할 수 있다.The designing the dictionary may include a plurality of matrices corresponding to each of the plurality of classes, the number of rows is the number of CSP filters, and the number of columns of each class is the number of training signals of each class. Can be.

상기 분류 방법은 상기 CSP 필터링된 신호 중에서 특정 주파수 대역의 신호를 특징 신호로 추출하는 단계를 더 포함하고, 상기 사전을 설계하는 단계는 상기 특징 신호를 기초로 상기 사전을 설계할 수 있다.The classification method may further include extracting a signal of a specific frequency band from the CSP filtered signal as a feature signal, and designing the dictionary may design the dictionary based on the feature signal.

상기 특정 주파수 대역은 두피의 감각운동 피질 부분의 감각운동 리듬에 해당하는 주파수 대역일 수 있다.The specific frequency band may be a frequency band corresponding to the sensory rhythm of the sensorimotor cortex portion of the scalp.

상기 해를 구하는 단계는 L1 최소화 방법으로 상기 해를 구할 수 있다.Obtaining the solution may be obtained by the L1 minimization method.

상기 해를 구하는 단계는 0이 아닌 계수의 개수가 각 클래스의 훈련 신호의 개수보다 희소하다고 판단될 수 있는 정도로 적은 상기 해를 구할 수 있다.The obtaining of the solution may obtain the solution so that the number of nonzero coefficients may be determined to be less than the number of training signals of each class.

상기 테스트 신호를 분류하는 단계는 상기 해의 각 클래스 중에서 0이 아닌 계수가 가장 많이 나타난 클래스로 상기 테스트 신호를 분류할 수 있다.In the classifying of the test signal, the test signal may be classified into a class in which a non-zero coefficient is highest among each class of the solution.

상기 테스트 신호를 분류하는 단계는 상기 해를 기초로 클래스별 나머지를 계산하여 나머지가 최소인 클래스로 상기 테스트 신호를 분류하고, 상기 클래스별 나머지는 상기 해 중에서 자신 이외의 클래스의 원소를 0으로 채워 생성한 벡터를 상기 사전에 투영한 제1값과 상기 테스트 신호에 해당하는 제2값의 차이로 계산되는 값일 수 있다.The classifying the test signal may include classifying the test signal into a class having a minimum residual, by calculating a residual for each class based on the solution, and filling the element of a class other than itself among the solutions with zeros. The generated vector may be a value calculated by a difference between a first value projected beforehand and a second value corresponding to the test signal.

본 발명의 실시예에 따르면 잡음이 심하고 시간에 따라 변화하는 움직임 상상 뇌파 신호를 기존의 분류 방법보다 정확하게 분류할 수 있다. 특히, 판별이 어려운 테스트 신호를 분류할 때 기존의 선형 판별법보다 정확히 분류할 수 있다. According to an exemplary embodiment of the present invention, a motion imaginary EEG signal which is noisy and changes over time can be classified more accurately than a conventional classification method. In particular, when classifying a test signal that is difficult to discriminate, it may be classified more accurately than a conventional linear discrimination method.

도 1은 본 발명의 한 실시예에 따른 뇌-컴퓨터 접속 장치의 개략적인 블록도이다.
도 2는 본 발명의 다른 실시예에 따른 뇌-컴퓨터 접속 장치의 개략적인 블록도이다.
도 3은 본 발명의 한 실시예에 따른 설계된 사전을 나타내는 도면이다.
도 4는 본 발명의 한 실시예에 따른 선형방정식을 나타내는 도면이다.
도 5는 본 발명의 한 실시예에 따른 CSP 필터링 전의 훈련 신호 관계를 나타내는 도면이다.
도 6은 본 발명의 한 실시예에 따른 CSP 필터링 후의 훈련 신호 관계를 나타내는 도면이다.
도 7은 본 발명의 한 실시예에 따른 분류 방법의 흐름도이다.
도 8은 선형 판별법에 의한 쉬운 분류 결과를 나타내는 도면이다.
도 9는 본 발명의 한 실시예에 따른 분류 결과를 나타내는 도면이다.
도 10은 선형 판별법에 의한 어려운 분류 결과를 나타내는 도면이다.
도 11은 본 발명의 다른 실시예에 따른 분류 결과를 나타내는 도면이다.
1 is a schematic block diagram of a brain-computer connection device according to an embodiment of the present invention.
2 is a schematic block diagram of a brain-computer connection device according to another embodiment of the present invention.
3 is a diagram illustrating a dictionary designed according to an embodiment of the present invention.
4 is a diagram illustrating a linear equation according to an embodiment of the present invention.
5 is a diagram illustrating a training signal relationship before CSP filtering according to an embodiment of the present invention.
6 is a diagram illustrating a training signal relationship after CSP filtering according to an embodiment of the present invention.
7 is a flowchart of a classification method according to an embodiment of the present invention.
8 is a diagram showing an easy classification result by the linear discriminant method.
9 is a diagram illustrating a classification result according to an embodiment of the present invention.
10 is a diagram showing a difficult classification result by the linear discriminant method.
11 is a view showing a classification result according to another embodiment of the present invention.

아래에서는 첨부한 도면을 참고로 하여 본 발명의 실시예에 대하여 본 발명이 속하는 기술 분야에서 통상의 지식을 가진 자가 용이하게 실시할 수 있도록 상세히 설명한다. 그러나 본 발명은 여러 가지 상이한 형태로 구현될 수 있으며 여기에서 설명하는 실시예에 한정되지 않는다. 그리고 도면에서 본 발명을 명확하게 설명하기 위해서 설명과 관계없는 부분은 생략하였으며, 명세서 전체를 통하여 유사한 부분에 대해서는 유사한 도면 부호를 붙였다.Hereinafter, embodiments of the present invention will be described in detail with reference to the accompanying drawings so that those skilled in the art can easily carry out the present invention. The present invention may, however, be embodied in many different forms and should not be construed as limited to the embodiments set forth herein. In order to clearly illustrate the present invention, parts not related to the description are omitted, and similar parts are denoted by like reference characters throughout the specification.

명세서 전체에서, 어떤 부분이 어떤 구성요소를 "포함"한다고 할 때, 이는 특별히 반대되는 기재가 없는 한 다른 구성요소를 제외하는 것이 아니라 다른 구성요소를 더 포함할 수 있는 것을 의미한다.Throughout the specification, when an element is referred to as "comprising ", it means that it can include other elements as well, without excluding other elements unless specifically stated otherwise.

이제 도면을 참고하여 본 발명의 실시예에 따른 뇌-컴퓨터 접속 장치, 그리고 그의 분류 방법에 대해 설명한다.A brain-computer connection device and a classification method thereof according to an embodiment of the present invention will now be described with reference to the drawings.

도 1은 본 발명의 한 실시예에 따른 뇌-컴퓨터 접속 장치의 개략적인 블록도이다. 1 is a schematic block diagram of a brain-computer connection device according to an embodiment of the present invention.

도 1을 참고하면, 뇌-컴퓨터 접속 장치(100)는 희소 표현 기반 분류(Sparse Representation-based Classification, SRC)를 통해 테스트 신호(

Figure 112012030230487-pat00001
)를 어느 하나의 클래스(class)로 분류한다. 여기서, 클래스는 테스트 신호(
Figure 112012030230487-pat00002
)를 분류하기 위해 나누어진 단위로서, 복수 개일 수 있다. 간단히, 클래스는 왼쪽 클래스와 오른쪽 클래스로 가정하고 설명한다. 여기서 희소 표현(sparse representation)이란 다양한 자연신호가 특정 기저에서는 희소(sparse)하기 때문에, 압축센싱(compressive sensing, CS)을 통해 자연신호들을 희소하게 표현(Sparse Representation)할 수 있다는 이론이다. 그리고, 이렇게 희소하게 표현 가능한 신호의 각 샘플은 모든 신호의 전체적인 관점을 담고 있기 때문에 간단한 선형적 투영 연산을 통하여 홀리스틱(holistic) 샘플로 압축될 수 있다. 이런 홀리스틱 샘플들은 양호한 신호 복구에 필요한 압축된 샘플들의 개수가 섀넌-나이퀴스트(Shannon-Nyquist) 샘플 개수보다 상당히 적다.Referring to FIG. 1, the brain-computer connection apparatus 100 may test a test signal through a sparse representation-based classification (SRC).
Figure 112012030230487-pat00001
) Into one class. Where class represents the test signal (
Figure 112012030230487-pat00002
) Are divided into units, which may be plural. For simplicity, we assume that classes are left and right classes. Here, the sparse representation is a theory that sparse representation of natural signals can be sparse through compression sensing (CS) because various natural signals are sparse at a specific basis. Each sample of such a sparse signal can be compressed into a holistic sample through a simple linear projection operation because it contains an overall view of all signals. These holistic samples have a significantly smaller number of compressed samples for good signal recovery than Shannon-Nyquist samples.

뇌-컴퓨터 접속 장치(100)는 분류하고자 하는 신호, 즉 테스트 신호(

Figure 112012030230487-pat00003
)를 사전(dictionary,
Figure 112012030230487-pat00004
)의 희소 표현을 기초로 해(
Figure 112012030230487-pat00005
)를 구한다. 그리고 뇌-컴퓨터 접속 장치(100)는 해(
Figure 112012030230487-pat00006
)를 기초로 테스트 신호(
Figure 112012030230487-pat00007
)를 분류한다. 이를 위해, 뇌-컴퓨터 접속 장치(100)는 사전부(dictionary unit)(110), 풀이부(solving unit)(130), 그리고 분류부(classification unit)(150)를 포함한다. The brain-computer connection device 100 is a signal to be classified, that is, a test signal (
Figure 112012030230487-pat00003
) Into a dictionary,
Figure 112012030230487-pat00004
Based on the rare expression of)
Figure 112012030230487-pat00005
). And brain-computer connection device 100
Figure 112012030230487-pat00006
Based on the test signal (
Figure 112012030230487-pat00007
Classify To this end, the brain-computer connection device 100 includes a dictionary unit 110, a solving unit 130, and a classification unit 150.

사전부(110)는 훈련 신호를 기초로 사전(

Figure 112012030230487-pat00008
)을 설계한다. 사전(
Figure 112012030230487-pat00009
)은 훈련 신호(training signal)를 기초로 희소하게 표현된 행렬(
Figure 112012030230487-pat00010
)로서, 테스트 신호(
Figure 112012030230487-pat00011
)의 해를 찾을 수 있는 기능을 한다. 사전(
Figure 112012030230487-pat00012
)은 수학식 1과 같이 나타낼 수 있다. 그리고, 수학식 1의 왼쪽 클래스에 대한 행렬(
Figure 112012030230487-pat00013
)과 오른쪽 클래스에 대한 행렬(
Figure 112012030230487-pat00014
)은 수학식 2와 같이 나타낼 수 있다. 여기서, 사전(
Figure 112012030230487-pat00015
)은
Figure 112012030230487-pat00016
차원이고,
Figure 112012030230487-pat00017
는 훈련 신호의 개수이다.The dictionary unit 110 is a dictionary (based on the training signal)
Figure 112012030230487-pat00008
). dictionary(
Figure 112012030230487-pat00009
) Is a sparse representation of a matrix based on a training signal (
Figure 112012030230487-pat00010
, The test signal (
Figure 112012030230487-pat00011
Function to find the solution. dictionary(
Figure 112012030230487-pat00012
) Can be expressed as in Equation 1. Then, the matrix for the left class of Equation 1
Figure 112012030230487-pat00013
) And the matrix for the right class (
Figure 112012030230487-pat00014
) May be expressed as in Equation 2. Where dictionary (
Figure 112012030230487-pat00015
)silver
Figure 112012030230487-pat00016
Dimension,
Figure 112012030230487-pat00017
Is the number of training signals.

Figure 112012030230487-pat00018
Figure 112012030230487-pat00018

Figure 112012030230487-pat00019
Figure 112012030230487-pat00019

사전부(110)는 훈련 신호를 열에 축적하여 사전(

Figure 112012030230487-pat00020
)을 만들 수 있다. 또는 뇌파 신호의 분류 성능을 높이기 위해, 사전부(110)는 뇌파 신호를 다양한 필터링 방법으로 필터링하여 사전(
Figure 112012030230487-pat00021
)을 만들 수 있다. 예를 들면, 사전부(110)는 CSP(common spatial pattern) 필터링 방법, PCA(Principal component analysis) 필터링 방법, 또는 ICA(Independent component analysis) 필터링 방법 등으로 훈련 신호를 필터링할 수 있다. 여기서는 CSP 필터링 방법 위주로 설명한다. The dictionary unit 110 accumulates a training signal in a column and advances the dictionary (
Figure 112012030230487-pat00020
). Alternatively, in order to increase the classification performance of the EEG signals, the dictionary unit 110 filters the EEG signals by various filtering methods.
Figure 112012030230487-pat00021
). For example, the dictionary unit 110 may filter the training signal using a common spatial pattern (CSP) filtering method, a principal component analysis (PCA) filtering method, or an independent component analysis (ICA) filtering method. This section focuses on the CSP filtering method.

여기서 훈련 신호는 움직임 상상 실험을 통해 획득한 뇌파(electroencephalogram, EGG) 신호로서, 분류하고자 하는 복수의 클래스, 예를 들면 왼쪽 클래스와 오른쪽 클래스 각각에 해당하는 신호들이다. 훈련 신호를 획득하는 방법은 다양하다. 예를 들면, 피실험자들은 왼손과 오른손의 움직임 상상(motor imagery)에 대해 BCI 실험을 수행한다. 피실험자들은 두피에 부착된 금도금 전극을 부착하고, 정해진 뇌파 채널 수에 대해 측정한다. 실험이 시작되면 ‘왼손’또는 ‘오른손’의 문자지시가 무작위로 보이면, 피실험자는 지시에 맞게 왼쪽이나 오른쪽 손 움직임, 예를 들면 주먹을 꽉 쥐는 행동을 상상한다. 이렇게 수집된 신호를 명령(Cue)이 나타난 후 특정 시점부터 특정 구간, 예를 들면 256샘플을 추출한다. 그리고, 피실험자 별로 최적의 시작 시점이 다르기 때문에, 각 피실험자에 대해 분류 정확도가 최대가 되는 시작시점을 선택할 수 있다.Here, the training signal is an electroencephalogram (EGG) signal obtained through a motion imagination experiment, and is a signal corresponding to a plurality of classes to be classified, for example, left and right classes. There are various ways to obtain training signals. For example, subjects perform BCI experiments on motor imagery of left and right hands. Subjects attached gold-plated electrodes attached to the scalp and measured for a defined number of EEG channels. When the experiment begins, if the letter indications of the left hand or the right hand appear random, the subject imagines the left or right hand movements, such as clenching fists, as directed. The collected signal is extracted from a specific time point, for example, 256 samples, after a command appears. Since the optimum start time differs for each test subject, a start time for maximizing classification accuracy can be selected for each test subject.

풀이부(130)는 사전(

Figure 112012030230487-pat00022
)을 통해서 테스트 신호(
Figure 112012030230487-pat00023
)의 해를 푼다. 즉, 풀이부(130)는 수학식 3의 선형방정식을 풀어 해(
Figure 112012030230487-pat00024
)를 구한다.The pool part 130 is a dictionary (
Figure 112012030230487-pat00022
Through the test signal (
Figure 112012030230487-pat00023
Solve the year. That is, the solver 130 solves the linear equation of equation (3) (
Figure 112012030230487-pat00024
).

Figure 112012030230487-pat00025
Figure 112012030230487-pat00025

해(

Figure 112012030230487-pat00026
)는 테스트 신호(
Figure 112012030230487-pat00027
)와 관계된 몇 개의 원소만 0이 아니고, 나머지 원소는 0인 희소 벡터(sparse vector)이다. 이때, 테스트 신호(
Figure 112012030230487-pat00028
)의 희소 표현은 해(
Figure 112012030230487-pat00029
)의 0이 아닌 원소의 개수가
Figure 112012030230487-pat00030
보다 훨씬 적을 때 만들어진다.year(
Figure 112012030230487-pat00026
) Is the test signal (
Figure 112012030230487-pat00027
Are not sparse, only a few elements are associated with the zero, and the rest are zero sparse vectors. At this time, the test signal (
Figure 112012030230487-pat00028
The rare expression of)
Figure 112012030230487-pat00029
), The number of nonzero elements
Figure 112012030230487-pat00030
Created when much less than.

풀이부(130)가 수학식 3을 풀 때, 사전(

Figure 112012030230487-pat00031
)의 열의 크기(
Figure 112012030230487-pat00032
)가 행(
Figure 112012030230487-pat00033
)보다 큰 언더디터민드(underdetermined) 상태이다. 따라서, 풀이부(130)는 L1 최소화 방법(L1 norm minimization)으로 해(
Figure 112012030230487-pat00034
)를 구하며, 수학식 4와 같이 표현된다.When the solver 130 solves Equation 3, the dictionary (
Figure 112012030230487-pat00031
Column's size (
Figure 112012030230487-pat00032
) Rows (
Figure 112012030230487-pat00033
Is underdetermined greater than). Therefore, the solver 130 is a solution L1 minimization (L1 norm minimization)
Figure 112012030230487-pat00034
) And is expressed as Equation 4.

Figure 112012030230487-pat00035
Figure 112012030230487-pat00035

예를 들어, 오른쪽 움직임 상상 신호에 해당하는 테스트 신호(

Figure 112012030230487-pat00036
)가 입력되는 경우, 해(
Figure 112012030230487-pat00037
)는 수학식 5와 같이 구해진다. 이때, 테스트 신호(
Figure 112012030230487-pat00038
)는 오른쪽 움직임 상상 신호이므로, 수학식 6과 같이 0이 아닌 계수(coefficient)가 오른쪽 클래스에 해당하는 해(
Figure 112012030230487-pat00039
)에 더 많이 나타난다. For example, the test signal corresponding to the right motion imaginary signal (
Figure 112012030230487-pat00036
) Is entered, the solution (
Figure 112012030230487-pat00037
) Is obtained as in Equation 5. At this time, the test signal (
Figure 112012030230487-pat00038
) Is a right motion imaginary signal, so a solution whose nonzero coefficient corresponds to the right class
Figure 112012030230487-pat00039
Appears more).

Figure 112012030230487-pat00040
Figure 112012030230487-pat00040

Figure 112012030230487-pat00041
Figure 112012030230487-pat00041

분류부(150)는 해(

Figure 112012030230487-pat00042
)를 기초로 테스트 신호(
Figure 112012030230487-pat00043
)를 분류한다. 이때, 최소화 문제를 풀어 획득한 해(
Figure 112012030230487-pat00044
)는 0이 아닌 원소들을 가지는데, 이 0이 아닌 원소들은 어느 하나의 클래스에 해당해야 한다. 하지만, 뇌파 신호는 잡음이 매우 심하기 때문에, 수학식 6과 같이 0이 아닌 원소들이 다른 클래스에도 나타날 수 있다. 따라서, 분류부(150)는 특성함수(
Figure 112012030230487-pat00045
)를 이용하여, 해(
Figure 112012030230487-pat00046
)가 어느 클래스에 해당하는지 분류한다. 여기서, 특성함수(
Figure 112012030230487-pat00047
)는 해(
Figure 112012030230487-pat00048
)에서, 다른 클래스와 관련된 모든 원소들을 0으로 채운(nullify) 벡터이다. 예를 들어, 특성함수(
Figure 112012030230487-pat00049
)는 해(
Figure 112012030230487-pat00050
)에서,
Figure 112012030230487-pat00051
의 원소들을 0으로 채운 벡터(
Figure 112012030230487-pat00052
)이고, 특성함수(
Figure 112012030230487-pat00053
)는 해(
Figure 112012030230487-pat00054
)에서,
Figure 112012030230487-pat00055
의 원소들을 0으로 채운 벡터(
Figure 112012030230487-pat00056
)이다. The classification unit 150 is the solution (
Figure 112012030230487-pat00042
Based on the test signal (
Figure 112012030230487-pat00043
Classify In this case, the solution obtained by solving the minimization problem (
Figure 112012030230487-pat00044
) Has nonzero elements, and these nonzero elements must belong to either class. However, since the EEG signal is very noisy, non-zero elements may appear in other classes as shown in Equation (6). Therefore, the classification unit 150 is a characteristic function (
Figure 112012030230487-pat00045
), The solution (
Figure 112012030230487-pat00046
) Classify which class it belongs to. Where the characteristic function (
Figure 112012030230487-pat00047
) Is the year (
Figure 112012030230487-pat00048
), All the elements associated with the other class It is a null filled vector. For example, the characteristic function (
Figure 112012030230487-pat00049
) Is the year (
Figure 112012030230487-pat00050
)in,
Figure 112012030230487-pat00051
A vector of zero elements
Figure 112012030230487-pat00052
) And the characteristic function (
Figure 112012030230487-pat00053
) Is the year (
Figure 112012030230487-pat00054
)in,
Figure 112012030230487-pat00055
A vector of zero elements
Figure 112012030230487-pat00056
)to be.

분류부(150)는 수학식 7과 같이 특성함수(

Figure 112012030230487-pat00057
)를 사전(
Figure 112012030230487-pat00058
)에 투영한 값[
Figure 112012030230487-pat00059
]과 테스트 신호(
Figure 112012030230487-pat00060
)의 차이를 기초로 클래스별 나머지[
Figure 112012030230487-pat00061
]를 계산한다. 그리고, 분류부(150)는 수학식 8과 같이 나머지[
Figure 112012030230487-pat00062
]가 최소가 되는 클래스로 테스트 신호(
Figure 112012030230487-pat00063
)를 분류한다.The classification unit 150 is a characteristic function (Equation 7)
Figure 112012030230487-pat00057
)
Figure 112012030230487-pat00058
Projected to) [
Figure 112012030230487-pat00059
] And test signal (
Figure 112012030230487-pat00060
Rest of each class based on the difference between
Figure 112012030230487-pat00061
]. And, the classification unit 150 is the rest [
Figure 112012030230487-pat00062
] Is the class with the minimum test signal (
Figure 112012030230487-pat00063
Classify).

Figure 112012030230487-pat00064
Figure 112012030230487-pat00064

Figure 112012030230487-pat00065
Figure 112012030230487-pat00065

도 2는 본 발명의 다른 실시예에 따른 뇌-컴퓨터 접속 장치의 개략적인 블록도이고, 도 3은 본 발명의 한 실시예에 따른 설계된 사전을 나타내는 도면이며, 도 4는 본 발명의 한 실시예에 따른 선형방정식을 나타내는 도면이다.2 is a schematic block diagram of a brain-computer connection device according to another embodiment of the present invention, FIG. 3 is a diagram showing a dictionary designed according to an embodiment of the present invention, and FIG. 4 is an embodiment of the present invention. It is a figure which shows the linear equation by.

먼저, 도 2를 참고하면, 뇌-컴퓨터 접속 장치(100)는 사전부(110), 풀이부(130), 분류부(150), CSP(common spatial pattern) 필터링부(170) 그리고 특징 신호 추출부(feature extraction unit)(190)를 포함한다. 사전부(110)는 CSP 필터링부(170)와 특징 신호 추출부(190)를 거친 훈련 신호를 이용하여 사전(

Figure 112012030230487-pat00066
)을 설계할 수 있다. First, referring to FIG. 2, the brain-computer connection device 100 includes a dictionary unit 110, a solver 130, a classifier 150, a CSP filtering unit 170, and feature signal extraction. A feature extraction unit 190. The dictionary unit 110 uses a training signal that has passed through the CSP filtering unit 170 and the feature signal extractor 190.
Figure 112012030230487-pat00066
) Can be designed.

CSP 필터링부(170)는 입력 받은 훈련 신호를 공간적으로 필터링하여 하나의 클래스에 대해서는 최대화하고, 다른 클래스에 대해서는 최소화한다. 즉, CSP 필터링부(170)는 훈련 신호를 두 클래스, 즉 왼쪽 클래스와 오른쪽 클래스 사이의 연관(correlation)이 적은 신호로 만든다. The CSP filtering unit 170 spatially filters the received training signal to maximize one class and minimize the other class. That is, the CSP filtering unit 170 makes the training signal a signal having a low correlation between two classes, that is, a left class and a right class.

CSP 필터링부(170)가 수행하는 CSP 필터링은 다음과 같다. 여기서,

Figure 112012030230487-pat00067
를 뇌파 신호의 세그먼트로 정의한다. 이때
Figure 112012030230487-pat00068
는 뇌파 채널 수이고, T는 한 명의 피실험자에 대한 모든 시행의 샘플 시점(time point)이다. 뇌파 훈련 신호는 왼손 움직임 상상과 오른손 움직임 상상에 각각 해당하는 두 클래스(
Figure 112012030230487-pat00069
,
Figure 112012030230487-pat00070
)를 포함한다.CSP filtering performed by the CSP filtering unit 170 is as follows. here,
Figure 112012030230487-pat00067
Is defined as the segment of the EEG signal. At this time
Figure 112012030230487-pat00068
Is the number of EEG channels and T is the sample time point of all trials for one subject. EEG training signals are divided into two classes, one for imagining the left hand movement and one for imagining the right hand movement.
Figure 112012030230487-pat00069
,
Figure 112012030230487-pat00070
) .

CSP 필터링부(170)는 CSP 알고리즘을 기초로 CSP 필터들의 셋트(set)인 필터 행렬(

Figure 112012030230487-pat00071
)을 찾는다. 이때, 필터 행렬(
Figure 112012030230487-pat00072
)의 각 열벡터(
Figure 112012030230487-pat00073
)는
Figure 112012030230487-pat00074
이다. CSP 필터링부(170)는 필터 행렬(
Figure 112012030230487-pat00075
)의 앞쪽 n개의 CSP 필터와 뒤쪽 n개의 CSP 필터를 사용하여 2n CPS 필터(
Figure 112012030230487-pat00076
)를 만든다. CSP 필터 행렬(
Figure 112012030230487-pat00077
)은 수학식 9와 같이 나타내진다.The CSP filtering unit 170 may filter the filter matrix, which is a set of CSP filters based on the CSP algorithm.
Figure 112012030230487-pat00071
). In this case, the filter matrix (
Figure 112012030230487-pat00072
Each column vector ()
Figure 112012030230487-pat00073
)
Figure 112012030230487-pat00074
to be. The CSP filtering unit 170 is a filter matrix (
Figure 112012030230487-pat00075
2n CPS filter (using front n CSP filters and rear n CSP filters)
Figure 112012030230487-pat00076
). CSP filter matrix (
Figure 112012030230487-pat00077
) Is represented by Equation (9).

Figure 112012030230487-pat00078
Figure 112012030230487-pat00078

왼손과 오른손 움직임 상상에 해당하는 훈련 신호(

Figure 112012030230487-pat00079
,
Figure 112012030230487-pat00080
)는 수학식 10과 같이 CSP 필터링된 신호로 정의될 수 있다.Training signals corresponding to imagining left and right hand movements (
Figure 112012030230487-pat00079
,
Figure 112012030230487-pat00080
) May be defined as a CSP filtered signal as shown in Equation 10.

Figure 112012030230487-pat00081
Figure 112012030230487-pat00081

특징신호추출부(190)는 CSP 필터링한 신호 중에서 특정 주파수 대역의 신호를 특징 신호로 추출한다. 여기서 특징 신호는 두피의 감각운동 피질 부분의 감각운동 리듬(sensorimotor rhythms, SMRs)에 해당하는 신호일 수 있다. 특징신호추출부(190)는 다양한 방법으로 특징 신호를 추출할 수 있다. 특징신호추출부(190)는 특정 주파수 대역의 밴드파워(band power)를 계산하여 특징 신호를 추출한다. 또는 특징신호추출부(190)는 각각의 CSP 채널에 대해 FFT(fast Fourier Transform)로 획득한 특정 주파수 대역의 계수를 특징 신호로 추출할 수 있다. 특정 주파수 대역은 감각운동 리듬(SMRs)의 Mu(8-14Hz) 또는 Beta(15-30Hz) 대역의 일부 대역일 수 있다.The feature signal extractor 190 extracts a signal of a specific frequency band from the CSP filtered signal as a feature signal. The feature signal may be a signal corresponding to sensorimotor rhythms (SMRs) of the sensorimotor cortex of the scalp. The feature signal extractor 190 may extract the feature signal in various ways. The feature signal extractor 190 extracts a feature signal by calculating a band power of a specific frequency band. Alternatively, the feature signal extractor 190 may extract a coefficient of a specific frequency band obtained by FFT (fast Fourier transform) for each CSP channel as a feature signal. The specific frequency band may be part of the Mu (8-14 Hz) or Beta (15-30 Hz) band of sensorimotor rhythms (SMRs).

사전부(110)는 CSP 필터링부(170)와 특징 신호 추출부(190)를 거친 훈련 신호를 이용하여 사전(

Figure 112012030230487-pat00082
)을 설계한다. 도 3을 참고하면, 사전(
Figure 112012030230487-pat00083
)은 왼쪽 클래스에 대한 행렬(
Figure 112012030230487-pat00084
)과 오른쪽 클래스에 대한 행렬(
Figure 112012030230487-pat00085
)로 구성된다. 그리고 사전(
Figure 112012030230487-pat00086
)의 행의 수는 CSP 필터의 수이고, 각 클래스의 열의 수는 각 클래스의 훈련 신호의 수이다.The dictionary unit 110 uses a training signal that has passed through the CSP filtering unit 170 and the feature signal extractor 190.
Figure 112012030230487-pat00082
). Referring to Figure 3, the dictionary (
Figure 112012030230487-pat00083
) Is the matrix (for the left class
Figure 112012030230487-pat00084
) And the matrix for the right class (
Figure 112012030230487-pat00085
). And the dictionary (
Figure 112012030230487-pat00086
The number of rows) is the number of CSP filters, and the number of columns in each class is the number of training signals in each class.

다시 도 2를 참고하면, 풀이부(130)는 사전(

Figure 112012030230487-pat00087
)을 통해서 테스트 신호(
Figure 112012030230487-pat00088
)의 해를 푼다. 이때, 풀이부(130)는 L1 최소화 방법(L1 norm minimization)으로 해(
Figure 112012030230487-pat00089
)를 구한다. 이때, 테스트 신호(
Figure 112012030230487-pat00090
)는 사전(
Figure 112012030230487-pat00091
)의 열을 만드는 방식과 같이 CSP 필터링부(170)와 특징 신호 추출부(190)를 거친 신호이며, 벡터화를 통해 m 차원의 벡터(
Figure 112012030230487-pat00092
)로 변환된다. 도 4를 참고하면, 풀이부(130)는 수학식 4의 선형방정식(
Figure 112012030230487-pat00093
)을 푼다. 만약, 오른쪽 클래스에 해당하는 테스트 신호(
Figure 112012030230487-pat00094
)가 입력되는 경우, 해(
Figure 112012030230487-pat00095
) 중에서 오른쪽 클래스(
Figure 112012030230487-pat00096
)에서 0이 아닌 계수(scalar coefficient)들이 다른 클래스보다 더 나타난다.Referring back to FIG. 2, the pooling unit 130 is a dictionary (
Figure 112012030230487-pat00087
Through the test signal (
Figure 112012030230487-pat00088
Solve the year. At this time, the pool 130 is a solution L1 minimization (L1 norm minimization)
Figure 112012030230487-pat00089
). At this time, the test signal (
Figure 112012030230487-pat00090
) Is a dictionary (
Figure 112012030230487-pat00091
Like the way to create a column of) is a signal passed through the CSP filtering unit 170 and the feature signal extraction unit 190, through the vectorization of the m-dimensional (
Figure 112012030230487-pat00092
Is converted to). 4, the solver 130 is a linear equation (4)
Figure 112012030230487-pat00093
) If the test signal for the right class (
Figure 112012030230487-pat00094
) Is entered, the solution (
Figure 112012030230487-pat00095
Of the right class (
Figure 112012030230487-pat00096
), Non-scalar coefficients appear more than other classes.

다시 도 2를 참고하면, 분류부(150)는 특성함수(

Figure 112012030230487-pat00097
)를 기초로 클래스별 나머지(
Figure 112012030230487-pat00098
)를 계산하고, 나머지가 최소가 되는 클래스로 테스트 신호(
Figure 112012030230487-pat00099
)를 분류한다.Referring back to FIG. 2, the classifying unit 150 has a characteristic function (
Figure 112012030230487-pat00097
Based on the rest of the class (
Figure 112012030230487-pat00098
), And the test signal (
Figure 112012030230487-pat00099
Classify

다음에서 특징신호추출부(190)가 특징 신호를 추출하는 다양한 방법에 대해 설명한다.Next, various methods for extracting the feature signal by the feature signal extractor 190 will be described.

특징신호추출부(190)는 밴드파워 계산(band power computation)을 통해 특징 신호를 추출할 수 있다. 이 경우, 특징신호추출부(190)는 CSP 필터링한 신호의 세컨드 모멘트(second moment)를 계산하여 특정 주파수 대역의 밴드파워를 계산한다. 사전(

Figure 112012030230487-pat00100
)은 밴드파워 계산으로 추출한 각 성분(element)을 이용하여 설계된다. 즉, 사전(
Figure 112012030230487-pat00101
)의 열은 2n개의 밴드파워로 구성될 수 있다.The feature signal extractor 190 may extract the feature signal through band power computation. In this case, the feature signal extractor 190 calculates a band moment of a specific frequency band by calculating a second moment of the CSP filtered signal. dictionary(
Figure 112012030230487-pat00100
) Is designed using each element extracted by band power calculation. That is, the dictionary (
Figure 112012030230487-pat00101
) May be composed of 2n band powers.

또는, 특징신호추출부(190)는 FFT를 통해 특징 신호를 추출할 수 있다. 이 경우, 특징신호추출부(190)는 특정 주파수 대역에 해당하는 계수만을 남긴다. 특징신호추출부(190)에서 출력된 신호는 수학식 11과 같이 나타낼 수 있다.Alternatively, the feature signal extractor 190 may extract the feature signal through the FFT. In this case, the feature signal extractor 190 leaves only coefficients corresponding to a specific frequency band. The signal output from the feature signal extractor 190 may be represented by Equation (11).

Figure 112012030230487-pat00102
Figure 112012030230487-pat00102

여기서,

Figure 112012030230487-pat00103
는 특징 신호로 추출된 열의 개수이다. 예를 들어, 샘플링 주파수가 256 samples/sec이고, 주파수 해상도가 1Hz이며, Mu(8~14Hz)와 Beta(15~30Hz)에 해당하는 밴드파워를 추출하는 경우,
Figure 112012030230487-pat00104
는 최대 24이다. here,
Figure 112012030230487-pat00103
Is the number of columns extracted as the feature signal. For example, if the sampling frequency is 256 samples / sec, the frequency resolution is 1 Hz, and the band power corresponding to Mu (8-14 Hz) and Beta (15-30 Hz) is extracted,
Figure 112012030230487-pat00104
Is up to 24.

이 경우, 사전부(110)는

Figure 112012030230487-pat00105
의 2n개의 행들을 연결한 뒤 트랜스포즈(transpose)하여 벡터화하여 사전(
Figure 112012030230487-pat00106
)의 왼쪽 행렬(
Figure 112012030230487-pat00107
)을 설계한다. 마찬가지로, 사전부(110)는
Figure 112012030230487-pat00108
의 2n개의 행들을 연결한 뒤 트랜스포즈하여 벡터화하여 사전(
Figure 112012030230487-pat00109
)의 오른쪽 행렬(
Figure 112012030230487-pat00110
)을 설계한다. In this case, the dictionary unit 110
Figure 112012030230487-pat00105
Concatenate 2n rows of and transpose and vectorize
Figure 112012030230487-pat00106
Left matrix ()
Figure 112012030230487-pat00107
). Similarly, the dictionary unit 110
Figure 112012030230487-pat00108
Concatenate 2n rows of and transpose and vectorize
Figure 112012030230487-pat00109
Matrix on right
Figure 112012030230487-pat00110
).

도 5는 본 발명의 한 실시예에 따른 CSP 필터링 전의 훈련 신호 관계를 나타내는 도면이고, 도 6은 본 발명의 한 실시예에 따른 CSP 필터링 후의 훈련 신호 관계를 나타내는 도면이다.5 is a diagram illustrating a training signal relationship before CSP filtering according to an embodiment of the present invention, and FIG. 6 is a diagram showing a training signal relationship after CSP filtering according to an embodiment of the present invention.

먼저 도 5를 참고하면, CSP 필터링 전의 훈련 신호는 왼쪽과 오른쪽 샘플이 연관(correlation)되어 있다.First, referring to FIG. 5, left and right samples are correlated with a training signal before CSP filtering.

다음 도 6을 참고하면, CSP 필터링 후의 훈련 신호는 CSP 필터링 전보다 연관이 적다. 즉, CSP 필터링이란 입력받은 훈련 신호를 공간적으로 필터링하여 하나의 클래스에 대해서는 최대화하고, 다른 클래스에 대해서는 최소화하기 때문에, CSP 필터링 후의 왼쪽과 오른쪽 샘플들은 CSP 필터링 전의 왼쪽과 오른쪽 샘플들에 비해 연관이 적어 잘 구분된다Referring to FIG. 6, training signals after CSP filtering are less relevant than before CSP filtering. In other words, since CSP filtering spatially filters the input training signal and maximizes one class and minimizes the other class, the left and right samples after CSP filtering are more related than the left and right samples after CSP filtering. Write well

따라서, 뇌-컴퓨터 접속 장치(100)는 훈련 신호를 CSP 필터링하여 사전(

Figure 112012030230487-pat00111
)을 만들면, 불규칙(incoherent)한 사전(
Figure 112012030230487-pat00112
)을 설계할 수 있다. 그리고, 압축센싱 시, 사전(
Figure 112012030230487-pat00113
)이 불규칙할수록 희소 표현 기반 분류의 정확도를 높일 수 있다. 결국, 뇌-컴퓨터 접속 장치(100)는 잡음이 매우 심하고 시간에 따라 변화하여 분류가 쉽지 않은 움직임 상상 뇌파 신호를 더 정확히 분류할 수 있다.Therefore, the brain-computer connection device 100 performs CSP filtering on the training signal.
Figure 112012030230487-pat00111
), An incoherent dictionary (
Figure 112012030230487-pat00112
) Can be designed. And during compression sensing, the dictionary (
Figure 112012030230487-pat00113
), The more irregular the higher the accuracy of sparse representation-based classification. As a result, the brain-computer connection device 100 may classify the motion imaginary EEG signal which is very noisy and changes over time, which is not easy to classify.

도 7은 본 발명의 한 실시예에 따른 분류 방법의 흐름도이다.7 is a flowchart of a classification method according to an embodiment of the present invention.

도 7을 참고하면, 뇌-컴퓨터 접속 장치(100)는 각 클래스에 해당하는 훈련 신호를 획득한다(S710). 훈련 신호는 움직임 상상 실험을 통해 획득한 뇌파 신호이다.Referring to FIG. 7, the brain-computer connection device 100 obtains a training signal corresponding to each class (S710). The training signal is an EEG signal obtained through a motion imagination experiment.

뇌-컴퓨터 접속 장치(100)는 훈련 신호를 하나의 클래스에 대해서는 최대화하고, 다른 클래스에 대해서는 최소화하는 CSP 필터링을 한다(S720). 뇌-컴퓨터 접속 장치(100)는 CSP 필터링을 통해 훈련 신호에서 왼쪽 클래스와 오른쪽 클래스의 연관을 줄여 구분이 잘 되도록 처리한다.The brain-computer connection device 100 performs CSP filtering to maximize the training signal for one class and minimize the other class (S720). The brain-computer connection device 100 processes the CSP filtering so that the left class and the right class are reduced in the training signal so as to be distinguished.

뇌-컴퓨터 접속 장치(100)는 CSP 필터링한 신호를 기초로 테스트 신호(

Figure 112012030230487-pat00114
)의 해(
Figure 112012030230487-pat00115
)를 풀기 위한 사전(
Figure 112012030230487-pat00116
)을 설계한다(S730). 뇌-컴퓨터 접속 장치(100)는 훈련 신호의 정보가 희소 표현된 사전(
Figure 112012030230487-pat00117
)을 만든다. 특히, 뇌-컴퓨터 접속 장치(100)는 CSP 필터링한 신호를 기초로 불규칙(incoherent)한 사전(
Figure 112012030230487-pat00118
)을 설계할 수 있다. 이때, 뇌-컴퓨터 접속 장치(100)는 CSP 필터링한 신호 중에서 특정 주파수 대역의 신호를 특징 신호로 추출하고, 이 특징 신호를 기초로 사전(
Figure 112012030230487-pat00119
)을 설계할 수 있다. 즉, 뇌-컴퓨터 접속 장치(100)는 감각운동 리듬에 해당하는 신호를 특징 신호로 추출할 수 있다.Brain-computer connection device 100 is based on the CSP filtered signal test signal (
Figure 112012030230487-pat00114
)due to(
Figure 112012030230487-pat00115
Dictionary () to solve
Figure 112012030230487-pat00116
) Is designed (S730). The brain-computer connection device 100 is a dictionary (in which the information of the training signal is scarcely expressed)
Figure 112012030230487-pat00117
). In particular, the brain-computer connection device 100 is an incoherent dictionary (based on the CSP filtered signal).
Figure 112012030230487-pat00118
) Can be designed. At this time, the brain-computer connection device 100 extracts a signal of a specific frequency band from the CSP filtered signal as a feature signal, and based on the feature signal,
Figure 112012030230487-pat00119
) Can be designed. That is, the brain-computer connection device 100 may extract a signal corresponding to the sensory motor rhythm as a feature signal.

뇌-컴퓨터 접속 장치(100)는 테스트 신호(

Figure 112012030230487-pat00120
)를 입력받아 선형방정식(
Figure 112012030230487-pat00121
)의 해(
Figure 112012030230487-pat00122
)를 구한다(S740). 이때, 선형방정식은 언더디터민드(underdetermined) 상태이므로, 뇌-컴퓨터 접속 장치(100)는 L1 최소화 방법으로 해(
Figure 112012030230487-pat00123
)를 구한다.Brain-computer connection device 100 is a test signal (
Figure 112012030230487-pat00120
), The linear equation (
Figure 112012030230487-pat00121
)due to(
Figure 112012030230487-pat00122
To obtain (S740). At this time, since the linear equation is underdetermined, the brain-computer connection device 100 uses the L1 minimization method (
Figure 112012030230487-pat00123
).

뇌-컴퓨터 접속 장치(100)는 클래스가 구분된 해(

Figure 112012030230487-pat00124
) 중에서 0이 아닌 계수의 분포를 기초로 테스트 신호(
Figure 112012030230487-pat00125
)를 분류한다(S750). 해(
Figure 112012030230487-pat00126
)는 테스트 신호(
Figure 112012030230487-pat00127
)와 관계된 몇 개의 원소만 0이 아니고, 나머지 원소는 0인 희소 벡터이다. 따라서, 뇌-컴퓨터 접속 장치(100)는 왼쪽 클래스에 해당하는 해(
Figure 112012030230487-pat00128
)의 계수들이 다른 클래스보다 더 나타나면 테스트 신호(
Figure 112012030230487-pat00129
)를 왼쪽으로 분류하고, 오른쪽 클래스에 해당하는 해(
Figure 112012030230487-pat00130
)의 계수들이 다른 클래스보다 더 나타나면 테스트 신호(
Figure 112012030230487-pat00131
)를 오른쪽으로 분류한다. 이를 위해, 뇌-컴퓨터 접속 장치(100)는 해(
Figure 112012030230487-pat00132
)에서, 다른 클래스와 관련된 모든 원소들을 0으로 채운 벡터인 특성함수(
Figure 112012030230487-pat00133
)를 이용하여 수학식 7과 같이 클래스별 나머지를 계산한다. 그리고, 뇌-컴퓨터 접속 장치(100)는 수학식 8과 같이 나머지가 최소가 되는 클래스로 테스트 신호(
Figure 112012030230487-pat00134
)를 분류할 수 있다.The brain-computer connection device 100 is a classed solution (
Figure 112012030230487-pat00124
Of test signals based on the distribution of nonzero coefficients
Figure 112012030230487-pat00125
) Are classified (S750). year(
Figure 112012030230487-pat00126
) Is the test signal (
Figure 112012030230487-pat00127
Only a few elements related to) are nonzero, and the rest are zero sparse vectors. Therefore, the brain-computer connection device 100 has a solution corresponding to the left class (
Figure 112012030230487-pat00128
), If more coefficients than other classes, the test signal (
Figure 112012030230487-pat00129
) To the left, and the solution for the right class (
Figure 112012030230487-pat00130
), If more coefficients than other classes, the test signal (
Figure 112012030230487-pat00131
) To the right. To this end, the brain-computer connection device 100 can solve (
Figure 112012030230487-pat00132
), All the elements associated with the other class A characteristic function that is a vector filled with zeros (
Figure 112012030230487-pat00133
) To calculate the remainder for each class as shown in equation (7). In addition, the brain-computer connection device 100 has a test signal in a class in which the rest is minimum, as shown in Equation (8).
Figure 112012030230487-pat00134
) Can be classified.

도 8은 선형 판별법에 의한 쉬운 분류 결과를 나타내는 도면이고, 도 9는 본 발명의 한 실시예에 따른 분류 결과를 나타내는 도면이다.FIG. 8 is a diagram illustrating an easy classification result by a linear discrimination method, and FIG. 9 is a diagram illustrating a classification result according to an embodiment of the present invention.

도 8을 참고하면, 선형 판별법(Linear Discriminant Analysis, LDA)은 오른쪽 훈련 신호(o 표시)와 왼쪽 훈련 신호(x 표시)를 각각 매핑한다. 그리고, 오른쪽 훈련 신호와 왼쪽 훈련 신호를 가르는 판별선(decision line)을 기초로 테스트 신호를 분류한다. 즉, 선형 판별법은 테스트 신호가 판별선 아래로 매핑되면 오른쪽 상상 신호로 분류하고, 판별선 위에 매핑되면 왼쪽 상상 신호로 분류한다. 이때, 판별이 쉬운 왼쪽 테스트 신호가 입력되면, 왼쪽 테스트 신호는 판별선에서 떨어진 위치에 매핑되어, 선형 판별법은 쉽게 테스트 신호를 왼쪽 상상 신호로 분별한다.Referring to FIG. 8, a linear discriminant analysis (LDA) maps a right training signal (o mark) and a left training signal (x mark), respectively. The test signal is classified based on a decision line separating the right training signal and the left training signal. That is, in the linear discrimination method, when the test signal is mapped below the discrimination line, it is classified as the right imaginary signal, and when it is mapped over the discrimination line, it is classified as the left imaginary signal. At this time, if the left test signal that is easy to discriminate is input, the left test signal is mapped to a position away from the discrimination line, and the linear discrimination method easily distinguishes the test signal into the left imaginary signal.

도 9를 참고하면, 뇌-컴퓨터 접속 장치(100) 역시 판별이 쉬운 테스트 신호가 입력되면, 0이 아닌 계수들이 거의 왼쪽 해에 나타나므로, 테스트 신호를 왼쪽 상상 신호로 쉽게 분별한다.Referring to FIG. 9, when the brain-computer connection device 100 also inputs a test signal that is easy to discriminate, non-zero coefficients appear in a nearly left solution, so that the test signal is easily classified into a left imaginary signal.

도 10은 선형 판별법에 의한 어려운 분류 결과를 나타내는 도면이고, 도 11은 본 발명의 다른 실시예에 따른 분류 결과를 나타내는 도면이다.10 is a diagram illustrating a difficult classification result by a linear discrimination method, and FIG. 11 is a diagram showing a classification result according to another embodiment of the present invention.

도 10을 참고하면, 선형 판별법은 판별이 어려운 테스트 신호가 입력되면, 판별선 근처에 매핑된다. 따라서, 왼쪽 테스트 신호는 왼쪽 신호임에도 불구하고, 판별선 아래에 매핑될 수 있다. 이 경우, 선형 판별법은 왼쪽 테스트 신호를 오른쪽 상상 신호로 잘못 분류할 수 있다. 이렇게 선형 판별법은 판별이 어려운 지역(grey region)이 존재하기 때문에, 테스트 신호를 정확히 분류하기에는 한계가 있다.Referring to FIG. 10, when a test signal that is difficult to discriminate is input, the linear discrimination method is mapped near the discrimination line. Thus, although the left test signal is a left signal, it may be mapped below the discrimination line. In this case, the linear discrimination method may misclassify the left test signal as the right imaginary signal. Since there is a gray region in which the linear discrimination method is difficult to discriminate, there is a limit in classifying the test signal accurately.

도 11을 참고하면, 뇌-컴퓨터 접속 장치(100)는 판별이 어려운 테스트 신호가 입력되면, 0이 아닌 계수들이 왼쪽 해뿐만 아니라 오른쪽 해에도 일부 나타난다. 그러나, 뇌-컴퓨터 접속 장치(100)는 테스트 신호의 희소 표현을 기초로 희소 벡터를 계산하므로, 왼쪽에 나타난 계수들이 오른쪽에 나타난 계수들보다 더 표현된다. 따라서, 뇌-컴퓨터 접속 장치(100)는 판별이 어려운 테스트 신호가 입력되더라도 정확히 분류할 가능성이 높다.Referring to FIG. 11, when a test signal that is difficult to discriminate is input to the brain-computer connection device 100, coefficients other than zero appear in the right solution as well as the left solution. However, since the brain-computer connection device 100 calculates the sparse vector based on the sparse representation of the test signal, the coefficients shown on the left are more represented than the coefficients shown on the right. Therefore, the brain-computer connection device 100 is more likely to classify correctly even if a difficult test signal is input.

이상에서 설명한 본 발명의 실시예는 장치 및 방법을 통해서만 구현이 되는 것은 아니며, 본 발명의 실시예의 구성에 대응하는 기능을 실현하는 프로그램 또는 그 프로그램이 기록된 기록 매체를 통해 구현될 수도 있다.The embodiments of the present invention described above are not implemented only by the apparatus and method, but may be implemented through a program for realizing the function corresponding to the configuration of the embodiment of the present invention or a recording medium on which the program is recorded.

이상에서 본 발명의 실시예에 대하여 상세하게 설명하였지만 본 발명의 권리범위는 이에 한정되는 것은 아니고 다음의 청구범위에서 정의하고 있는 본 발명의 기본 개념을 이용한 당업자의 여러 변형 및 개량 형태 또한 본 발명의 권리범위에 속하는 것이다.While the present invention has been particularly shown and described with reference to exemplary embodiments thereof, it is to be understood that the invention is not limited to the disclosed exemplary embodiments, It belongs to the scope of right.

Claims (20)

뇌파 신호인 테스트 신호를 왼쪽 클래스와 오른쪽 클래스 중 어느 하나로 분류하는 뇌-컴퓨터 접속 장치로서,
왼쪽 클래스의 훈련 신호와 오른쪽 클래스의 훈련 신호를 입력받고, 각 클래스의 훈련 신호를 필터링하여, 상기 왼쪽 클래스에 해당하는 제1행렬과 상기 오른쪽 클래스에 해당하는 제2행렬을 추출하고, 상기 제1행렬과 상기 제2행렬을 축적하여 사전을 설계하는 사전부,
상기 사전을 기초로 테스트 신호의 해를 구하는 풀이부, 그리고
상기 해의 각 클래스에 나타난 0이 아닌 계수의 분포를 기초로 상기 테스트 신호를 상기 왼쪽 클래스와 상기 오른쪽 클래스 중 어느 하나로 분류하는 분류부
를 포함하고,
상기 제1행렬과 상기 제2행렬 각각은 상기 필터링에 의해 희소하게 표현된 행렬인, 뇌-컴퓨터 접속 장치.
A brain-computer connection device that classifies a test signal, an EEG signal, into one of a left class and a right class,
The training signal of the left class and the training signal of the right class are received, and the training signal of each class is filtered to extract a first matrix corresponding to the left class and a second matrix corresponding to the right class. A dictionary unit for designing a dictionary by accumulating a matrix and the second matrix,
A solution unit for solving a test signal based on the dictionary, and
A classification unit classifying the test signal into one of the left class and the right class based on a distribution of nonzero coefficients represented in each class of the solution.
Lt; / RTI >
Wherein each of the first matrix and the second matrix is a matrix sparsely represented by the filtering.
제1항에서,
상기 훈련 신호를 하나의 클래스에 대해서는 최대화하고, 다른 클래스에 대해서는 최소화하는 CSP(common spatial pattern) 필터링을 하여 CSP 필터링된 신호를 출력하는 CSP 필터링부를 더 포함하고,
상기 사전부는 상기 CSP 필터링된 신호를 기초로 상기 사전을 설계하는 뇌-컴퓨터 접속 장치.
In claim 1,
The apparatus further includes a CSP filtering unit configured to output a CSP filtered signal by performing CSP (common spatial pattern) filtering to maximize the training signal for one class and minimize the other class.
And the dictionary unit designs the dictionary based on the CSP filtered signal.
제2항에서,
상기 사전부는
행의 수는 CSP 필터의 수이며, 각 클래스의 열의 수는 각 클래스의 훈련 신호의 수인 상기 사전을 설계하는 뇌-컴퓨터 접속 장치.
3. The method of claim 2,
The dictionary part
The number of rows is the number of CSP filters, and the number of columns of each class is the number of training signals of each class.
제2항에서,
상기 CSP 필터링된 신호 중에서 특정 주파수 대역의 신호를 특징 신호로 추출하는 특징신호추출부를 더 포함하고,
상기 사전부는 상기 특징 신호를 기초로 상기 사전을 설계하는 뇌-컴퓨터 접속 장치.
3. The method of claim 2,
Further comprising a feature signal extraction unit for extracting a signal of a specific frequency band from the CSP filtered signal as a feature signal,
And the dictionary unit designs the dictionary based on the feature signal.
제4항에서,
상기 특정 주파수 대역은 두피의 감각운동 피질 부분의 감각운동 리듬에 해당하는 주파수 대역인 뇌-컴퓨터 접속 장치.
5. The method of claim 4,
Wherein said specific frequency band is a frequency band corresponding to the sensory rhythm of the sensorimotor cortex portion of the scalp.
제1항에서,
상기 풀이부는
L1 최소화 방법으로 상기 해를 구하는 뇌-컴퓨터 접속 장치.
In claim 1,
The glue part
A brain-computer connection device that solves the solution by L1 minimization.
제1항에서,
상기 분류부는
상기 해의 각 클래스 중에서 0이 아닌 계수가 가장 많이 나타난 클래스로 상기 테스트 신호를 분류하는 뇌-컴퓨터 접속 장치.
In claim 1,
The classifier
A brain-computer connection device for classifying the test signal into a class in which the most non-zero coefficients appear among each class of the solution.
제1항에서,
상기 분류부는
상기 해를 기초로 클래스별 나머지를 계산하여 나머지가 최소인 클래스로 상기 테스트 신호를 분류하고,
상기 클래스별 나머지는 상기 해 중에서 자신 이외의 클래스의 원소를 0으로 채워 생성한 벡터를 상기 사전에 투영한 제1값과 상기 테스트 신호에 해당하는 제2값의 차이로 계산되는 값인 뇌-컴퓨터 접속 장치.
In claim 1,
The classifier
Calculate the remainder for each class based on the solution and classify the test signal into a class having the least remainder,
The remainder for each class is a value calculated by the difference between a first value projecting a vector generated by filling an element of a class other than itself with zero in the solution and a second value corresponding to the test signal. Device.
제1항에서,
상기 훈련 신호와 상기 테스트 신호는 움직임 상상 기반 뇌파 신호인 뇌-컴퓨터 접속 장치.
In claim 1,
And the training signal and the test signal are motion imaginary based brain wave signals.
제1항에서,
상기 왼쪽 클래스는 왼쪽 움직임 상상에 해당하는 클래스이고, 상기 오른쪽 클래스는 오른쪽 움직임 상상에 해당하는 클래스인 뇌-컴퓨터 접속 장치.
In claim 1,
The left class is a class corresponding to the left movement imagination, and the right class is a class corresponding to the right movement imagination.
뇌-컴퓨터 접속 장치가 뇌파 신호인 테스트 신호를 왼쪽 클래스와 오른쪽 클래스 중 어느 하나로 분류하는 방법으로서,
왼쪽 클래스에 해당하는 훈련 신호와 오른쪽 클래스에 해당하는 훈련 신호를 획득하는 단계,
각 클래스의 훈련 신호를 해당 클래스에 대해서는 최대화하고 다른 클래스에 대해서는 최소화하는 필터링하여, 상기 왼쪽 클래스에 해당하는 제1행렬과 상기 오른쪽 클래스에 해당하는 제2행렬을 추출하는 단계,
상기 제1행렬과 상기 제2행렬을 축적하여 사전(dictionary)을 설계하는 단계,
상기 사전을 기초로 상기 테스트 신호의 해를 구하는 단계, 그리고
상기 해의 각 클래스에 나타난 0이 아닌 계수의 분포를 기초로 상기 테스트 신호를 상기 왼쪽 클래스와 상기 오른쪽 클래스 중 어느 하나로 분류하는 단계
를 포함하고,
상기 제1행렬과 상기 제2행렬은 상기 필터링에 의해 희소하게 표현된 행렬인, 분류 방법.
The brain-computer connection device classifies a test signal, which is an EEG signal, into one of a left class and a right class,
Acquiring a training signal corresponding to the left class and a training signal corresponding to the right class;
Filtering the training signals of each class to maximize the corresponding class and minimize the other class, extracting a first matrix corresponding to the left class and a second matrix corresponding to the right class;
Designing a dictionary by accumulating the first matrix and the second matrix;
Solving the test signal based on the dictionary, and
Classifying the test signal into one of the left class and the right class based on a distribution of nonzero coefficients represented in each class of the solution.
Lt; / RTI >
And the first matrix and the second matrix are matrices sparsely represented by the filtering.
제11항에서,
상기 훈련 신호와 상기 테스트 신호는 움직임 상상 기반 뇌파 신호인 분류 방법.
12. The method of claim 11,
And the training signal and the test signal are motion imaginary based brain wave signals.
제11항에서,
상기 왼쪽 클래스는 왼쪽 움직임 상상에 해당하는 클래스이고, 상기 오른쪽 클래스는 오른쪽 움직임 상상에 해당하는 클래스인 분류 방법.
12. The method of claim 11,
And the left class is a class corresponding to a left movement imagination, and the right class is a class corresponding to a right movement imagination.
제11항에서,
상기 필터링은 CSP(Common Spatial Pattern) 필터링이고,
상기 사전을 설계하는 단계는
행의 수는 CSP 필터의 수이며, 각 클래스의 열의 수는 각 클래스의 훈련 신호의 수인 상기 사전을 설계하는 분류 방법.
12. The method of claim 11,
The filtering is Common Spatial Pattern (CSP) filtering,
The step of designing the dictionary
And the number of rows is the number of CSP filters and the number of columns of each class is the number of training signals of each class.
제14항에서,
상기 CSP 필터링된 신호 중에서 특정 주파수 대역의 신호를 특징 신호로 추출하는 단계를 더 포함하고,
상기 사전을 설계하는 단계는 상기 특징 신호를 기초로 상기 사전을 설계하는 분류 방법.
The method of claim 14,
Extracting a signal of a specific frequency band from the CSP filtered signal as a feature signal;
Designing the dictionary comprises designing the dictionary based on the feature signal.
제15항에서,
상기 특정 주파수 대역은 두피의 감각운동 피질 부분의 감각운동 리듬에 해당하는 주파수 대역인 분류 방법.
16. The method of claim 15,
The specific frequency band is a frequency band corresponding to the sensory motor rhythm of the sensorimotor cortex portion of the scalp.
제11항에서,
상기 해를 구하는 단계는
L1 최소화 방법으로 상기 해를 구하는 분류 방법.
12. The method of claim 11,
The solution step is
A classification method for solving the solution by L1 minimization.
제11항에서,
상기 해를 구하는 단계는
0이 아닌 계수의 개수가 각 클래스의 훈련 신호의 개수보다 희소하다고 판단될 수 있는 정도로 적은 상기 해를 구하는 분류 방법.
12. The method of claim 11,
The solution step is
A classification method for obtaining the solution that is so small that the number of nonzero coefficients can be determined to be less than the number of training signals of each class.
제11항에서,
상기 분류하는 단계는
상기 해의 각 클래스 중에서 0이 아닌 계수가 가장 많이 나타난 클래스로 상기 테스트 신호를 분류하는 분류 방법.
12. The method of claim 11,
The classifying step
A classification method for classifying the test signal into classes in which the most non-zero coefficients appear among the classes of the solution.
제11항에서,
상기 분류하는 단계는
상기 해를 기초로 클래스별 나머지를 계산하여 나머지가 최소인 클래스로 상기 테스트 신호를 분류하고, 상기 클래스별 나머지는 상기 해 중에서 자신 이외의 클래스의 원소를 0으로 채워 생성한 벡터를 상기 사전에 투영한 제1값과 상기 테스트 신호에 해당하는 제2값의 차이로 계산되는 값인 분류 방법.
12. The method of claim 11,
The classifying step
The test signal is classified into a class having a least residual by calculating a residual for each class based on the solution, and the residual for each class projects a vector generated by filling an element of a class other than itself with zero in the solution. And a value calculated by a difference between one first value and a second value corresponding to the test signal.
KR1020120039497A 2011-05-09 2012-04-17 Apparatus of brain-computer interface, and classification method thereof KR101380964B1 (en)

Priority Applications (2)

Application Number Priority Date Filing Date Title
KR1020120039497A KR101380964B1 (en) 2011-05-09 2012-04-17 Apparatus of brain-computer interface, and classification method thereof
PCT/KR2012/003572 WO2012153965A2 (en) 2011-05-09 2012-05-07 Brain-computer interface device and classification method therefor

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
KR1020110043567 2011-05-09
KR20110043567 2011-05-09
KR1020120039497A KR101380964B1 (en) 2011-05-09 2012-04-17 Apparatus of brain-computer interface, and classification method thereof

Publications (2)

Publication Number Publication Date
KR20120125948A KR20120125948A (en) 2012-11-19
KR101380964B1 true KR101380964B1 (en) 2014-04-04

Family

ID=47139796

Family Applications (1)

Application Number Title Priority Date Filing Date
KR1020120039497A KR101380964B1 (en) 2011-05-09 2012-04-17 Apparatus of brain-computer interface, and classification method thereof

Country Status (2)

Country Link
KR (1) KR101380964B1 (en)
WO (1) WO2012153965A2 (en)

Families Citing this family (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101524585B1 (en) * 2013-08-08 2015-06-04 중앙대학교 산학협력단 Apparatus and method for classifying brain wave
CN104473635B (en) * 2014-12-16 2016-07-06 重庆邮电大学 Right-hand man's Mental imagery EEG feature extraction method of hybrid wavelet and common space pattern
CN104771163B (en) * 2015-01-30 2018-01-02 杭州电子科技大学 EEG feature extraction method based on CSP and R CSP algorithms
CN105962932B (en) * 2016-04-20 2018-09-11 西安电子科技大学 Spike potential detection method based on Orthogonal Subspaces vector
CN107423668B (en) * 2017-04-14 2022-09-27 山东建筑大学 Electroencephalogram signal classification system and method based on wavelet transformation and sparse expression
CN107292329B (en) * 2017-04-19 2020-10-16 西北工业大学 Event idea classification method based on CI-CSP algorithm
US10582316B2 (en) 2017-11-30 2020-03-03 Starkey Laboratories, Inc. Ear-worn electronic device incorporating motor brain-computer interface
KR102248796B1 (en) * 2018-10-31 2021-05-07 울산과학기술원 Method and device of predicting user's decision
KR102175997B1 (en) * 2018-12-13 2020-11-06 성균관대학교산학협력단 Methods and apparatuses for recognizing user intention based on brain wave
US20200229730A1 (en) * 2019-01-17 2020-07-23 Mitsubishi Electric Research Laboratories, Inc. Methods and Systems for Source Connectivity Estimation With Hierarchical Structural Vector Autoregressive Models
CN110321856B (en) * 2019-07-08 2023-01-10 山东大学 Time-frequency multi-scale divergence CSP brain-computer interface method and device
CN110516711B (en) * 2019-07-25 2022-10-21 安徽大学 Training set quality evaluation method of MI-BCI system and optimization method of single training sample
KR102340544B1 (en) 2019-10-30 2021-12-20 고려대학교 산학협력단 Device and method for motor imagery classification using time domain parameter and correlation coefficients
KR102347688B1 (en) 2019-11-01 2022-01-06 고려대학교 산학협력단 Device method for correlation coefficient based cognitive task classification using demixing eeg channel pair
CN112971813A (en) * 2021-01-30 2021-06-18 上海麓联智能科技有限公司 Peak potential classification system and method for neural signals
CN113780392B (en) * 2021-08-31 2024-03-19 杭州电子科技大学 Channel selection method for motor imagery electroencephalogram identification

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20110064302A1 (en) * 2008-01-31 2011-03-17 Yi Ma Recognition via high-dimensional data classification
KR20110037726A (en) * 2009-10-07 2011-04-13 포항공과대학교 산학협력단 Method of analysing composite common spatial pattern for brain computer interface and method of analysing electroencephalogram using the same

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20050085744A1 (en) * 2003-10-20 2005-04-21 Stmicroelectronics S.R.I. Man-machine interfaces system and method, for instance applications in the area of rehabilitation
US7120486B2 (en) * 2003-12-12 2006-10-10 Washington University Brain computer interface

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20110064302A1 (en) * 2008-01-31 2011-03-17 Yi Ma Recognition via high-dimensional data classification
KR20110037726A (en) * 2009-10-07 2011-04-13 포항공과대학교 산학협력단 Method of analysing composite common spatial pattern for brain computer interface and method of analysing electroencephalogram using the same

Also Published As

Publication number Publication date
WO2012153965A3 (en) 2013-02-21
WO2012153965A2 (en) 2012-11-15
KR20120125948A (en) 2012-11-19

Similar Documents

Publication Publication Date Title
KR101380964B1 (en) Apparatus of brain-computer interface, and classification method thereof
Cui et al. Automatic sleep stage classification based on convolutional neural network and fine-grained segments
García-Salinas et al. Transfer learning in imagined speech EEG-based BCIs
Shin et al. Sparse representation-based classification scheme for motor imagery-based brain–computer interface systems
Zhang et al. Bayesian learning for spatial filtering in an EEG-based brain–computer interface
Ameri et al. Projective dictionary pair learning for EEG signal classification in brain computer interface applications
KR101068017B1 (en) Method of analysing composite common spatial pattern for brain computer interface and method of analysing electroencephalogram using the same
Hettiarachchi et al. Motor imagery data classification for BCI application using wavelet packet feature extraction
Uehara et al. Robust averaging of covariances for EEG recordings classification in motor imagery brain-computer interfaces
Lotte et al. EEG feature extraction
Siuly et al. Classification of EEG signals using sampling techniques and least square support vector machines
Taherisadr et al. EEG signal dimensionality reduction and classification using tensor decomposition and deep convolutional neural networks
Ghaheri et al. Extracting common spatial patterns from EEG time segments for classifying motor imagery classes in a Brain Computer Interface (BCI)
Abbas et al. FBCSP-based multi-class motor imagery classification using BP and TDP features
Nakra et al. Feature Extraction and Dimensionality Reduction Techniques with Their Advantages and Disadvantages for EEG-Based BCI System: A Review.
Upadhyay et al. Feature extraction and classification of imagined motor movement electroencephalogram signals
KR101345216B1 (en) Mehtod and apparatus for analyzing multi-subject electroencephalograph
KR102340544B1 (en) Device and method for motor imagery classification using time domain parameter and correlation coefficients
Ghaheri et al. Temporal windowing in CSP method for multi-class Motor Imagery Classification
Puri et al. Wavelet packet sub-band based classification of alcoholic and controlled state EEG signals
Gupta et al. Improved classification of motor imagery datasets for BCI by using approximate entropy and WOSF features
Gu et al. Heterogeneous classifier ensembles for EEG-based motor imaginary detection
Higashi et al. Regularization using geometric information between sensors capturing features from brain signals
Mo et al. EEG classification based on sparse representation
Thang et al. Investigation of regularization theory for four-class classification in brain-computer interface

Legal Events

Date Code Title Description
A201 Request for examination
E902 Notification of reason for refusal
E701 Decision to grant or registration of patent right
GRNT Written decision to grant
FPAY Annual fee payment

Payment date: 20161219

Year of fee payment: 4

FPAY Annual fee payment

Payment date: 20180226

Year of fee payment: 5

LAPS Lapse due to unpaid annual fee