JP6859622B2 - EEG signal processing system, EEG signal processing method and program - Google Patents
EEG signal processing system, EEG signal processing method and program Download PDFInfo
- Publication number
- JP6859622B2 JP6859622B2 JP2016147559A JP2016147559A JP6859622B2 JP 6859622 B2 JP6859622 B2 JP 6859622B2 JP 2016147559 A JP2016147559 A JP 2016147559A JP 2016147559 A JP2016147559 A JP 2016147559A JP 6859622 B2 JP6859622 B2 JP 6859622B2
- Authority
- JP
- Japan
- Prior art keywords
- feature amount
- stimulus
- electroencephalogram
- features
- brain wave
- Prior art date
- Legal status (The legal status 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 status listed.)
- Active
Links
Images
Landscapes
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
- Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
- Measuring And Recording Apparatus For Diagnosis (AREA)
Description
本発明は、ユーザーの脳波を観測することによってユーザーの心理状態あるいは生理状態を知る際に用いて好適な脳波信号処理システム、脳波信号処理方法及びプログラムに関する。 The present invention relates to an electroencephalogram signal processing system, an electroencephalogram signal processing method and a program suitable for use in knowing the psychological state or physiological state of the user by observing the electroencephalogram of the user.
脳内活動の分布や変化を解析することによって、ユーザーの心理状態や生理状態を推定することができる。例えば、ブレイン・マシン・インタフェース(Brain−Machine Interface)は、ユーザーの脳内活動を解析することで当該ユーザーの意思を推定して装置の操作入力とするものである(非特許文献1)。脳内活動からユーザーの感情を推定しようとする試みも様々行われている(非特許文献2)。 By analyzing the distribution and changes in brain activity, the psychological and physiological states of the user can be estimated. For example, the Brain-Machine Interface estimates the user's intention by analyzing the user's brain activity and uses it as an operation input of the device (Non-Patent Document 1). Various attempts have been made to estimate the user's emotions from the activity in the brain (Non-Patent Document 2).
脳内活動の状況を推定するために広く用いられている計測手段の一つは脳波計側である。頭皮上に複数個の脳波電極を装着し、複数チャンネルの増幅が可能な脳波アンプを用いることで、脳内活動の分布を反映した脳波信号が得られる。フーリエ変換などの周波数解析を用いて求められる脳波のパワースペクトルの変化や、電極間のパワースペクトルの差異を、脳内信号源の挙動を示す指標として、ユーザーの状態推定のために用いる方法が広く用いられている。脳の異なる部位の活動に対応する異なる電極位置で観測される複数チャンネルの脳波信号の相互関係性から、ユーザーの状態を推定する方法(非特許文献3)もある。脳内信号源の相互関係性は脳波チャンネル間の相互相関係数などで数値化することができる。 One of the widely used measuring means for estimating the state of intracerebral activity is the electroencephalograph side. By mounting a plurality of electroencephalogram electrodes on the scalp and using an electroencephalogram amplifier capable of amplifying a plurality of channels, an electroencephalogram signal reflecting the distribution of intracerebral activity can be obtained. The method of using the change in the power spectrum of the brain wave obtained by using frequency analysis such as Fourier transform and the difference in the power spectrum between the electrodes as an index showing the behavior of the signal source in the brain for estimating the user's state is widespread. It is used. There is also a method of estimating the user's state from the interrelationship of the electroencephalogram signals of a plurality of channels observed at different electrode positions corresponding to the activities of different parts of the brain (Non-Patent Document 3). The interrelationship of signal sources in the brain can be quantified by the mutual correlation coefficient between EEG channels.
ただし、頭皮上に装着された電極を用いて計測される脳波信号(以下、観測信号と称する)は、脳内や脳外に散在する複数の信号源からの電気信号が混合されたものであり、脳内の特定部位の挙動を直接的に反映するものではない。複数の信号源が存在すると仮定される場合に、特定の事象に関係する信号源を分離して解析する方法として、信号源が発する原信号の統計的独立性を仮定した独立成分分析(特許文献1、非特許文献4)などのアルゴリズムを用いた信号源分離手法が知られている。 However, the electroencephalogram signal (hereinafter referred to as the observation signal) measured using the electrodes mounted on the scalp is a mixture of electrical signals from a plurality of signal sources scattered in and out of the brain. , It does not directly reflect the behavior of a specific part in the brain. Independent component analysis assuming statistical independence of the original signal emitted by the signal source as a method of separating and analyzing the signal sources related to a specific event when it is assumed that there are a plurality of signal sources (Patent Documents). 1. A signal source separation method using an algorithm such as Non-Patent Document 4) is known.
生理心理的影響が既知の外部刺激をユーザーに与え、この外部刺激を受容しているユーザの脳波を測定し、測定した脳波を解析して得られる特徴量(以下、脳波特徴量と記す)を求める。そして、ユーザに与えた外部刺激によってユーザに生理心理的状態の変化が生じるため、このユーザの生理心理的状態の変化に対応して連動して変化する脳波特徴量を探索することで、ユーザに与えた外部刺激における脳波特徴量と生理心理的状態の変化がとの相互関係性を求めることができる。
さらに、この相互関係性を既知の情報として用いることにより、未知の状態にあるユーザの脳波から求めた脳波特徴量に基づき、ユーザの生理心理的状態を推定することが可能となる。
An external stimulus with a known physiological and psychological effect is given to the user, the brain wave of the user receiving this external stimulus is measured, and the feature amount obtained by analyzing the measured brain wave (hereinafter referred to as the brain wave feature amount) is described. Ask. Then, since the external stimulus given to the user causes a change in the physiological and psychological state of the user, the user is asked to search for an electroencephalogram feature amount that changes in conjunction with the change in the physiological and psychological state of the user. It is possible to determine the interrelationship between the EEG features and changes in the physiological and psychological state of the given external stimulus.
Furthermore, by using this interrelationship as known information, it becomes possible to estimate the physiological and psychological state of the user based on the electroencephalogram feature amount obtained from the electroencephalogram of the user in an unknown state.
また、ユーザに対して感情や場面や運動などの状況を想起するように指示し、この状況の想起を内部刺激として用い、ユーザに生理心理的状態を誘導させ、この内部刺激を与えた際の脳波を測定する。そして、この測定結果からユーザの生理心理的状態の変化と連動して変化する脳波特徴量を求めることもできる。 In addition, when the user is instructed to recall situations such as emotions, scenes, and movements, and the recall of this situation is used as an internal stimulus, the user is induced to have a physiological and psychological state, and when this internal stimulus is given. Measure brain waves. Then, from this measurement result, it is possible to obtain an electroencephalogram feature amount that changes in conjunction with a change in the user's physiological and psychological state.
ところで、脳波は、ニューロンネットワークを流れたインパルスを外部から観測した電気信号である。このため、脳波は、主要な振動信号の周波数によって、0.4〜4Hzのデルタ波、4〜8Hzのシータ波、8〜14Hzのアルファ波および14〜26Hzのベータ波に分類することができる。 By the way, an electroencephalogram is an electrical signal obtained by observing an impulse flowing through a neuron network from the outside. Therefore, brain waves can be classified into 0.4 to 4 Hz delta waves, 4 to 8 Hz theta waves, 8 to 14 Hz alpha waves, and 14 to 26 Hz beta waves, depending on the frequency of the main vibration signal.
電気信号である脳波特徴量を、周波数帯域の各々における信号成分の振幅とした場合、取得する脳波特徴量の数は、以下のように計算される。すなわち、脳波を測定するためにユーザの頭部に取り付ける電極の数をN個とし、振幅を抽出するための周波数帯域の数をMとしたい場合、脳波特徴量の数はN×M個となる。
さらに、頭部に取り付けた電極間の相互関係性、例えば位相差または相関係数を脳波特徴量として加えると、電極の数をN個とし、周波数帯域の数をMとした場合、脳波特徴量の数はN×(N−1)/2×M個となる。
When the electroencephalogram feature amount which is an electric signal is the amplitude of the signal component in each frequency band, the number of electroencephalogram feature amount to be acquired is calculated as follows. That is, if the number of electrodes attached to the user's head for measuring brain waves is N and the number of frequency bands for extracting amplitude is M, the number of brain wave features is N × M. ..
Furthermore, when the interrelationship between the electrodes attached to the head, for example, the phase difference or the correlation coefficient, is added as the electroencephalogram feature amount, the number of electrodes is N, and the number of frequency bands is M, the electroencephalogram feature amount. The number of is N × (N-1) / 2 × M.
上述した位相差または相関係数だけでなく、位相差や時系列的な振幅値や波形の変化を脳波特徴量として用いることもできる。
さらに、脳波の各々の周波帯域の信号により求められる関数や、得られた脳波特徴量の四則演算を行なうことで二次的な脳波特徴量を導出することもできる。
この結果、脳波特徴量の数は膨大となり、得られた脳波特徴量の中から、ユーザの生理心理的状態の変化に対応した脳波特徴量、例えば「快適さ」や「楽しさ」の生理心理的状態の変化に関係する脳波特徴量を探索するためは、探索における解析に多大な計算負荷を要する。
In addition to the above-mentioned phase difference or correlation coefficient, the phase difference, time-series amplitude value, and waveform change can be used as the electroencephalogram feature amount.
Furthermore, secondary EEG features can be derived by performing the functions obtained from the signals of each frequency band of the EEG and the four arithmetic operations of the obtained EEG features.
As a result, the number of electroencephalogram features becomes enormous, and from the obtained electroencephalogram features, the electroencephalogram features corresponding to changes in the user's physiological and psychological state, such as "comfort" and "enjoyment", are physiological psychology. In order to search for EEG features related to changes in the target state, a large computational load is required for the analysis in the search.
本発明は、このような状況に鑑みてなされたものであり、多くの脳波特徴量の中から、ユーザーの生理心理的状態変化に関係している脳波特徴量を、実験者(分析者)に対して負荷を与えずに容易に探索することができる脳波信号処理システム、脳波信号処理方法及びプログラムを提供することを目的とする。 The present invention has been made in view of such a situation, and among many electroencephalogram features, the electroencephalogram features related to the change in the physiological and psychological state of the user are given to the experimenter (analyst). An object of the present invention is to provide an electroencephalogram signal processing system, an electroencephalogram signal processing method, and a program that can be easily searched without applying a load.
上述した課題を解決するために、本発明の脳波信号処理システムは、与えられた刺激に対応して発生する脳波信号の特徴量である脳波特徴量、前記刺激の特徴量である刺激特徴量の関連性を求める脳波信号処理システムであって、前記刺激に対応した前記脳波信号の複数の前記脳波特徴量を求める特徴量算出部と、前記刺激の1以上の前記刺激特徴量の各々と、複数の前記脳波特徴量それぞれとの組合わせ毎に、当該組合わせの関連の度合を示す相互関連度を算出する相互関連度算出部と、前記相互関連度が所定の閾値を超える、前記刺激特徴量と前記脳波特徴量との組合わせを抽出し、前記刺激特徴量に対して関連性の高い前記脳波特徴量を選択する特徴量選択部と、前記脳波信号を複数系列の脳波の特徴を有する成分信号に分離する信号分離部とを備え、前記特徴量算出部が、前記成分信号から前記脳波特徴量を算出することを特徴とする。
本発明の脳波信号処理システムは、与えられた刺激に対応して発生する脳波信号の特徴量である脳波特徴量、前記刺激の特徴量である刺激特徴量の関連性を求める脳波信号処理システムであって、前記刺激の各々に対応して、脳内信号源それぞれが発した前記脳波信号の周波数に基づく複数の前記脳波特徴量を求める特徴量算出部と、前記刺激の1以上の前記刺激特徴量の各々と、複数の前記脳波特徴量それぞれとの組合わせ毎に、当該組合わせの関連の度合を示す相互関連度を算出する相互関連度算出部と、前記相互関連度が所定の閾値を超える、前記刺激特徴量と前記脳波特徴量との組合わせを抽出し、前記刺激特徴量に対して関連性の高い前記脳波特徴量を選択する特徴量選択部とを備えることを特徴とする。
本発明の脳波信号処理システムは、与えられた刺激に対応して発生する脳波信号の特徴量である脳波特徴量のうち、前記刺激の特徴量である刺激特徴量の各々に対して関連性の高い前記脳波特徴量を選択する脳波信号処理システムであって、前記刺激に対応した前記脳波信号の複数の前記脳波特徴量を求める特徴量算出部と、前記刺激の1以上の前記刺激特徴量の各々と、複数の前記脳波特徴量それぞれとの組合わせ毎に、当該組合わせの関連の度合を示す相互関連度を算出する相互関連度算出部と、前記相互関連度が所定の閾値を超える、前記刺激特徴量と前記脳波特徴量との組合わせを抽出し、前記刺激特徴量に対して関連性の高い前記脳波特徴量を選択する特徴量選択部とを備え、前記複数の前記脳波特徴量のうち、少なくとも2以上の前記脳波特徴量は互いに異なる種類の単位であることを特徴とする。
In order to solve the above-mentioned problems, the electroencephalogram signal processing system of the present invention has an electroencephalogram feature amount which is a feature amount of an electroencephalogram signal generated in response to a given stimulus and a stimulus feature amount which is a feature amount of the stimulus. a brain wave signal processing system for determining the relevance, the feature amount calculating unit for obtaining a plurality of the EEG feature value of the brain wave signal corresponding to the stimulus, and each of the one or more of the stimulation characteristic quantity of the stimulation, a plurality For each combination with each of the electroencephalogram feature amounts of the above, an interrelationship degree calculation unit that calculates the degree of interrelationship indicating the degree of association of the combination, and the stimulus feature amount that the degree of interrelationship exceeds a predetermined threshold. A feature amount selection unit that extracts a combination of the electroencephalogram feature amount and the electroencephalogram feature amount and selects the electroencephalogram feature amount that is highly related to the stimulus feature amount, and a component having the electroencephalogram signal having a plurality of series of electroencephalogram features. A signal separation unit for separating into signals is provided, and the feature amount calculation unit calculates the electroencephalogram feature amount from the component signal.
The electroencephalogram signal processing system of the present invention is an electroencephalogram signal processing system that obtains the relationship between the electroencephalogram feature amount, which is the feature amount of the electroencephalogram signal generated in response to a given stimulus, and the stimulus feature amount, which is the feature amount of the stimulus. Therefore, a feature amount calculation unit for obtaining a plurality of the electroencephalogram feature amounts based on the frequency of the electroencephalogram signal emitted by each of the signal sources in the brain corresponding to each of the stimuli, and one or more of the stimulus features of the stimulus. For each combination of each amount and each of the plurality of electroencephalogram feature amounts, a degree of interrelationship calculation unit that calculates the degree of interrelationship indicating the degree of association of the combination, and the degree of interrelationship set a predetermined threshold value. It is characterized by including a feature amount selection unit that extracts a combination of the stimulus feature amount and the brain wave feature amount that exceeds the stimulus feature amount and selects the brain wave feature amount that is highly related to the stimulus feature amount.
The electroencephalogram signal processing system of the present invention is related to each of the stimulus feature quantities, which is the feature quantity of the stimulus, among the electroencephalogram feature quantities, which are the feature quantities of the brain wave signals generated in response to the given stimulus. An electroencephalogram signal processing system that selects a high electroencephalogram feature amount, which comprises a feature amount calculation unit for obtaining a plurality of the electroencephalogram feature amounts of the electroencephalogram signal corresponding to the stimulus, and one or more of the stimulus feature amounts of the stimulus. For each combination of each and the plurality of electroencephalogram feature amounts, a degree of interrelationship calculation unit that calculates the degree of relevance indicating the degree of relevance of the combination, and the degree of interrelationship exceeding a predetermined threshold. A combination of the stimulus feature amount and the brain wave feature amount is extracted, and the feature amount selection unit for selecting the brain wave feature amount having a high relevance to the stimulus feature amount is provided, and the plurality of the brain wave feature amounts are provided. of, characterized in that at least two or more of the EEG feature amount Ru different types of unit der each other.
本発明の脳波信号処理システムは、前記相互関連度として、相互情報量を用いることを特徴とする。 The electroencephalogram signal processing system of the present invention is characterized by using a mutual information amount as the degree of interrelationship.
本発明の脳波信号処理システムは、前記相互関連度として、相関係数を用いることを特徴とする。 The electroencephalogram signal processing system of the present invention is characterized in that a correlation coefficient is used as the degree of interrelationship.
本発明の脳波信号処理システムは、ユーザーに対して前記刺激を提示する刺激提示部をさらに有し、前記刺激提示部が、種類の異なる前記刺激特徴量を有する前記刺激をユーザーに対して提示することを特徴とする。 The electroencephalogram signal processing system of the present invention further includes a stimulus presenting unit that presents the stimulus to the user, and the stimulus presenting unit presents the stimulus having a different type of stimulus feature amount to the user. It is characterized by that.
本発明の脳波信号処理システムは、脳波推定データベースをさらに備え、前記特徴量選択部が、前記刺激特徴量に対して選択した、当該刺激特徴量に対して関連性の高い前記脳波特徴量を、前記刺激特徴量と組として、前記脳波推定データベースに書き込んで記憶させることを特徴とする。 The electroencephalogram signal processing system of the present invention further includes an electroencephalogram estimation database, and the electroencephalogram feature amount that is highly related to the stimulus feature amount selected by the feature amount selection unit for the stimulus feature amount. It is characterized in that it is written and stored in the electroencephalogram estimation database as a set with the stimulus feature amount.
本発明の脳波信号処理システムは、前記脳波信号の前記脳波特徴量が供給されると、当該脳波特徴量に対応する前記刺激特徴量を前記脳波推定データベースから抽出して、前記刺激の刺激特徴量を推定するユーザー状態推定部をさらに備えることを特徴とする。 In the electroencephalogram signal processing system of the present invention, when the electroencephalogram feature amount of the electroencephalogram signal is supplied, the stimulus feature amount corresponding to the electroencephalogram feature amount is extracted from the electroencephalogram estimation database, and the stimulus feature amount of the stimulus is extracted. It is characterized by further including a user state estimation unit for estimating.
本発明の脳波信号処理方法は、与えられた刺激に対応して発生する脳波信号の特徴量である脳波特徴量、前記刺激の特徴量である刺激特徴量の関係性を求める脳波信号処理方法であって、与えられた刺激に対応して発生する脳波信号の特徴量である脳波特徴量、前記刺激の特徴量である刺激特徴量の関連性を求める脳波信号処理方法であって、特徴量算出部が、前記刺激に対応した複数の前記脳波信号の前記脳波特徴量を求める特徴量算出過程と、相互関連度算出部が、前記刺激の1以上の前記刺激特徴量の各々と、複数の前記脳波特徴量それぞれとの組合わせ毎に、当該組合わせの関連の度合を示す相互関連度を算出する相互関連度算出過程と、特徴量選択部が、前記相互関連度が所定の閾値を超える、前記刺激特徴量と前記脳波特徴量との組合わせを抽出し、前記刺激特徴量に対して関連性の高い前記脳波特徴量を選択する特徴量選択過程と、信号分離部が、前記脳波信号を複数系列の脳波の特徴を有する成分信号に分離する信号分離過程とを含み、前記特徴量算出部が、前記成分信号から前記脳波特徴量を算出することを特徴とする。
本発明の脳波信号処理方法は、与えられた刺激に対応して発生する脳波信号の特徴量である脳波特徴量、前記刺激の特徴量である刺激特徴量の関連性を求める脳波信号処理方法であって、特徴量算出部が、前記刺激の各々に対応して、脳内信号源それぞれが発した前記脳波信号の周波数に基づく複数の前記脳波特徴量を求める特徴量算出過程と、相互関連度算出部が、前記刺激の1以上の前記刺激特徴量の各々と、複数の前記脳波特徴量それぞれとの組合わせ毎に、当該組合わせの関連の度合を示す相互関連度を算出する相互関連度算出過程と、特徴量選択部が、前記相互関連度が所定の閾値を超える、前記刺激特徴量と前記脳波特徴量との組合わせを抽出し、前記刺激特徴量に対して関連性の高い前記脳波特徴量を選択する特徴量選択過程とを含むことを特徴とする。
本発明の脳波信号処理方法は、与えられた刺激に対応して発生する脳波信号の特徴量である脳波特徴量のうち、前記刺激の特徴量である刺激特徴量の各々に対して関連性の高い前記脳波特徴量を選択する脳波信号処理方法であって、特徴量算出部が、前記刺激に対応した前記脳波信号の複数の前記脳波特徴量を求める特徴量算出過程と、相互関連度算出部が、前記刺激の1以上の前記刺激特徴量の各々と、複数の前記脳波特徴量それぞれとの組合わせ毎に、当該組合わせの関連の度合を示す相互関連度を算出する相互関連度算出過程と、特徴量選択部が、前記相互関連度が所定の閾値を超える、前記刺激特徴量と前記脳波特徴量との組合わせを抽出し、前記刺激特徴量に対して関連性の高い前記脳波特徴量を選択する特徴量選択過程とを含み、前記複数の前記脳波特徴量のうち、少なくとも2以上の前記脳波特徴量は互いに異なる種類の単位である。
The electroencephalogram signal processing method of the present invention is an electroencephalogram signal processing method for obtaining the relationship between the electroencephalogram feature amount, which is the feature amount of the electroencephalogram signal generated in response to a given stimulus, and the stimulus feature amount, which is the feature amount of the stimulus. Therefore, it is an electroencephalogram signal processing method for obtaining the relationship between the electroencephalogram feature amount, which is the feature amount of the electroencephalogram signal generated in response to a given stimulus, and the stimulus feature amount, which is the feature amount of the stimulus, and the feature amount is calculated. The feature amount calculation process for obtaining the brain wave feature amount of the plurality of electroencephalogram signals corresponding to the stimulus, and the interrelationship degree calculation unit are each of one or more of the stimulus feature amounts of the stimulus and a plurality of the above. For each combination with each electroencephalogram feature amount, the interrelationship degree calculation process for calculating the degree of interrelationship indicating the degree of association of the combination and the feature amount selection unit indicate that the degree of interrelationship exceeds a predetermined threshold. The feature amount selection process of extracting the combination of the stimulus feature amount and the brain wave feature amount and selecting the brain wave feature amount having a high relevance to the stimulus feature amount, and the signal separation unit generate the brain wave signal. It includes a signal separation process of separating into component signals having a plurality of series of electroencephalogram features, and the feature amount calculation unit calculates the brain wave feature amount from the component signals.
EEG signal processing method of the present invention, EEG feature quantity is a feature quantity of EEG signal generated in response to a given stimulus, the stimulus of the feature at which stimulation characteristic of relevance determined electroencephalogram signal processing method Therefore, the feature amount calculation unit obtains a plurality of the electroencephalogram feature amounts based on the frequencies of the electroencephalogram signals emitted by each of the signal sources in the brain in response to each of the stimuli, and the degree of interrelationship with the feature amount calculation process. The calculation unit calculates the degree of interrelationship indicating the degree of association of the combination for each combination of one or more of the stimulus feature amounts of the stimulus and each of the plurality of the electroencephalogram feature amounts. The calculation process and the feature amount selection unit extract a combination of the stimulus feature amount and the electroencephalogram feature amount in which the degree of interrelationship exceeds a predetermined threshold, and the feature amount is highly related to the stimulus feature amount. It is characterized by including a feature amount selection process for selecting an electroencephalogram feature amount.
The brain wave signal processing method of the present invention is related to each of the stimulus features, which are the stimulus features, among the brain wave features, which are the features of the brain wave signals generated in response to the given stimulus. A brain wave signal processing method for selecting a high brain wave feature amount, wherein the feature amount calculation unit obtains a plurality of the brain wave feature amounts of the brain wave signal corresponding to the stimulus, and the interrelationship degree calculation unit. However, for each combination of one or more of the stimulus features of the stimulus and each of the plurality of brain wave features, a process of calculating the degree of interrelationship indicating the degree of relevance of the combination is calculated. The feature amount selection unit extracts a combination of the stimulus feature amount and the brain wave feature amount whose degree of interrelationship exceeds a predetermined threshold, and the brain wave feature having a high relevance to the stimulus feature amount. The feature amount selection process of selecting the amount is included, and at least two or more of the brain wave feature amounts among the plurality of the brain wave feature amounts are units of different types from each other.
本発明のプログラムは、与えられた刺激に対応して発生する脳波信号の特徴量である脳波特徴量、前記刺激の特徴量である刺激特徴量の関連性を求めるためのプログラムであって、コンピュータを、前記刺激に対応した前記脳波信号の複数の前記脳波特徴量を求める特徴量算出手段、前記刺激の1以上の前記刺激特徴量の各々と、複数の前記脳波特徴量それぞれとの組合わせ毎に、当該組合わせの関連の度合を示す相互関連度を算出する相互関連度算出手段、前記相互関連度が所定の閾値を超える、前記刺激特徴量と前記脳波特徴量との組合わせを抽出し、前記刺激特徴量に対して関連性の高い前記脳波特徴量を選択する特徴量選択手段、前記脳波信号を複数系列の脳波の特徴を有する成分信号に分離する信号分離手段として動作させ、前記信号分離手段が、前記脳波信号を複数系列の脳波の特徴を有する成分信号に分離するプログラムである。
本発明のプログラムは、与えられた刺激に対応して発生する脳波信号の特徴量である脳波特徴量、前記刺激の特徴量である刺激特徴量の関連性を求めるためのプログラムであって、コンピュータを、前記刺激の各々に対応して、脳内信号源それぞれが発した前記脳波信号の周波数に基づく複数の前記脳波特徴量を求める特徴量算出手段、前記刺激の1以上の前記刺激特徴量の各々と、複数の前記脳波特徴量それぞれとの組合わせ毎に、当該組合わせの関連の度合を示す相互関連度を算出する相互関連度算出手段、前記相互関連度が所定の閾値を超える、前記刺激特徴量と前記脳波特徴量との組合わせを抽出し、前記刺激特徴量に対して関連性の高い前記脳波特徴量を選択する特徴量選択手段として動作させるためのプログラムである。
本発明のプログラムは、与えられた刺激に対応して発生する脳波信号の特徴量である脳波特徴量のうち、前記刺激の特徴量である刺激特徴量の各々に対して関連性の高い前記脳波特徴量を選択するためのプログラムであって、コンピュータを、前記刺激に対応した前記脳波信号の複数の前記脳波特徴量を求める特徴量算出手段、前記刺激の1以上の前記刺激特徴量の各々と、複数の前記脳波特徴量それぞれとの組合わせ毎に、当該組合わせの関連の度合を示す相互関連度を算出する相互関連度算出手段、前記相互関連度が所定の閾値を超える、前記刺激特徴量と前記脳波特徴量との組合わせを抽出し、前記刺激特徴量に対して関連性の高い前記脳波特徴量を選択する特徴量選択手段として動作させ、前記複数の前記脳波特徴量のうち、少なくとも2以上の前記脳波特徴量は互いに異なる種類の単位であるプログラムである。
Program of the present invention is a program for determining a characteristic quantity of EEG signal generated in response to a given stimulus EEG feature amount, the relevant stimulus feature quantity is a feature quantity of the stimulation, the computer A feature amount calculating means for obtaining a plurality of the brain wave feature amounts of the brain wave signal corresponding to the stimulus, each combination of one or more of the stimulus feature amounts of the stimulus and each of the plurality of the brain wave feature amounts. In addition, a means for calculating the degree of interrelationship, which indicates the degree of relevance of the combination, and a combination of the stimulus feature amount and the electroencephalogram feature amount in which the degree of interrelationship exceeds a predetermined threshold are extracted. , The feature amount selection means for selecting the electroencephalogram feature amount having a high relevance to the stimulus feature amount, and the signal separation means for separating the electroencephalogram signal into component signals having a plurality of series of electroencephalogram features. The separation means is a program that separates the electroencephalogram signal into component signals having the characteristics of a plurality of series of electroencephalograms.
Program of the present invention is a program for determining a characteristic quantity of EEG signal generated in response to a given stimulus EEG feature amount, the relevant stimulus feature quantity is a feature quantity of the stimulation, the computer A feature amount calculating means for obtaining a plurality of the electroencephalogram feature amounts based on the frequencies of the electroencephalogram signals emitted by each of the signal sources in the brain corresponding to each of the stimuli, and one or more of the stimulus feature amounts of the stimulus. Correlation degree calculation means for calculating the degree of interrelationship indicating the degree of association of the combination for each combination of each of the plurality of electroencephalogram feature amounts, and the degree of interrelationship exceeds a predetermined threshold. This is a program for extracting a combination of a stimulus feature amount and the electroencephalogram feature amount and operating it as a feature amount selection means for selecting the brain wave feature amount having a high relevance to the stimulus feature amount.
In the program of the present invention, among the brain wave features which are the features of the brain wave signal generated in response to a given stimulus, the brain waves which are highly related to each of the stimulus features which are the features of the stimulus. It is a program for selecting a feature amount, and a computer is used with a feature amount calculation means for obtaining a plurality of the brain wave feature amounts of the brain wave signal corresponding to the stimulus, and one or more of the stimulus feature amounts of the stimulus. , The interrelationship degree calculation means for calculating the degree of interrelationship indicating the degree of association of the combination for each combination with each of the plurality of the brain wave features, and the stimulus feature in which the degree of interrelationship exceeds a predetermined threshold. A combination of the amount and the brain wave feature amount is extracted and operated as a feature amount selection means for selecting the brain wave feature amount having a high relevance to the stimulus feature amount, and among the plurality of the brain wave feature amounts. At least two or more of the brain wave features are programs that are different types of units.
以上説明したように、本発明によれば、多くの脳波特徴量の中から、ユーザーの生理心理的状態変化に関係している脳波特徴量を、実験者(分析者)に対して負荷を与えずに容易に探索することができる脳波信号処理システム、脳波信号処理方法及びプログラムを提供することができる。 As described above, according to the present invention, the EEG feature amount related to the change in the physiological and psychological state of the user is given a load to the experimenter (analyst) from among many EEG feature amounts. It is possible to provide an electroencephalogram signal processing system, an electroencephalogram signal processing method, and a program that can be easily searched without.
以下、本発明の脳波信号処理システムの実施形態について図面を用いて説明する。なお、本実施形態は、一例として、独立成分分析による脳内信号源の分離を行って脳内信号源を推定するとともに、脳内信号源が発した特定周波数成分の信号の振幅と特定周波数成分の各々の位相差の変化を用いたユーザー状態推定を行う構成を備えている。
図1は、本発明の脳波信号処理システムの一実施形態であるユーザー状態推定装置の概略構成を示すブロック図である。ユーザー状態推定装置1は、コンピュータ10と、脳波計測部19と、表示装置18aと、スピーカ11aを備える。コンピュータ10は、例えば、パーソナルコンピュータ、スマートフォン、タブレット端末等の情報処理装置である。図1に示した例は、コンピュータ10の外部に表示装置18aとスピーカ11aを設けているが、コンピュータ10は表示装置18aおよびスピーカ11aを内蔵していてもよい。
また、図1に示した例において、刺激の提示としては、聴覚からの刺激として音を用いているが、静止画、動画、文章などの視覚からの刺激、あるいは触覚、嗅覚、味覚による刺激など音以外の刺激を用いても良い。
また、音、画像、臭いなどの外部刺激を用いるのではなく、ユーザーに所定の感情を想起したり、所定の場面を想像させるよう指示を与え、脳内における内部的な刺激状態に誘導してもよい。
図1に示した例は、脳波特徴量と刺激特徴量との相互関係度の指標として、相互情報量を用いているが、相関係数などの2変数間の関係性(関連性)の高さを示す他の指標を用いても良い。
Hereinafter, embodiments of the electroencephalogram signal processing system of the present invention will be described with reference to the drawings. In this embodiment, as an example, the signal source in the brain is estimated by separating the signal source in the brain by independent component analysis, and the amplitude and the specific frequency component of the signal of the specific frequency component emitted by the signal source in the brain are estimated. It has a configuration for estimating the user state using the change of each phase difference of.
FIG. 1 is a block diagram showing a schematic configuration of a user state estimation device according to an embodiment of the electroencephalogram signal processing system of the present invention. The user
Further, in the example shown in FIG. 1, sound is used as a stimulus from the auditory sense, but a visual stimulus such as a still image, a moving image, or a sentence, or a stimulus by the sense of touch, smell, or taste, etc. Stimuli other than sound may be used.
In addition, instead of using external stimuli such as sounds, images, and odors, the user is instructed to recall a predetermined emotion or imagine a predetermined scene, and induces an internal stimulus state in the brain. May be good.
In the example shown in FIG. 1, the mutual information amount is used as an index of the degree of interrelationship between the EEG feature amount and the stimulus feature amount, but the relationship (relationship) between two variables such as the correlation coefficient is high. Other indicators indicating the above may be used.
図1に示したコンピュータ10は、刺激提示部11と、信号分離部12と、特徴量算出部13と、刺激特徴量記憶部14と、相互情報量算出部15と、特徴量選択部16と、ユーザー状態推定部17と、表示部18とを備える。刺激提示部11、信号分離部12、特徴量算出部13、刺激特徴量記憶部14、相互情報量算出部15、特徴量選択部16およびユーザー状態推定部17は、コンピュータ10を構成するCPU(中央処理装置)、記憶装置等のハードウェアと、CPUが実行するソフトウエアとから構成されている。
The
図1に示した構成例において、刺激提示部11は、感情により変化する脳波を測定する際に、ユーザーの感情(喜怒哀楽、好き嫌いなど)を時系列的に変化させるため、時系列的に変化する刺激を、ユーザーに対して提示する。本実施形態においては、例えば、スピーカ11aを介してユーザーに一定時間間隔で、ユーザーの感情が大きく変化するように、曲調及び曲想の各々が異なる楽曲を切り替えて提示する。
In the configuration example shown in FIG. 1, the
次に、図2は、図1に示した脳波計測部19の概略構成を示す。脳波計測部19は接地電極21と計測チャンネルと同数以上の脳波電極22を含む。また、脳波計測部19は、増幅器24とAD(アナログ/デジタル)変換器25とを含む。接地電極21および脳波電極22は、ユーザーの頭部7にバンド22a、キャップ、ジェル剤等を用いて装着される。増幅器24は、微弱な電気信号である脳波を増幅し、AD変換器25の入力とする。AD変換器25は、増幅器24で増幅されたアナログ脳波信号をデジタル変換し、コンピュータ10の入力となるデジタル化された観測信号を生成する。ここで、観測信号は、複数チャンネルのデジタル脳波信号を含む。デジタル脳波信号は、脳波と雑音とを含むアナログ脳波信号をデジタル化した信号である。
Next, FIG. 2 shows a schematic configuration of the
AD変換器25は、各チャンネルのアナログ脳波信号を、例えば、4000回/秒のレートでサンプリングして、サンプリングした1サンプルを1チャンネル当たり16ビットのデジタルデータに変換する。この場合、2kHzまでの脳波信号が計測可能である。
The
図1に戻り、信号分離部12は、観測信号を、周波数分割や独立成分分析などにより複数系列の成分信号に分離する。すなわち、信号分離部12は、複数チャンネルの観測信号に対して独立成分分析を実行し、脳内信号源に対応すると推定される相互に独立した信号である独立成分に分離する。また、信号分離部12は、分離した独立成分を、内部に設けられた周波数フィルタを通してアルファ波成分とベータ波成分とに分離し成分信号を得る。この周波数フィルタは、分離された全ての周波数成分を含む信号から、解析に用いる周波数成分を抽出する。すなわち、周波数フィルタは、分離した各信号成分から解析に用いる目的信号として、周波数帯域のアルファ波の信号成分(以下、アルファ波成分とする)とベータ波の信号成分(以下、ベータ波成分とする)との各々を抽出して出力する。これにより、信号分離部12は、観測信号である脳波信号を複数系列の脳波の特徴を有する成分信号に分離する。
Returning to FIG. 1, the
特徴量算出部13は、信号分離部12によって分離された成分信号の脳波特徴量の時間区間毎の平均値を求める。本実施形態においては、脳波特徴量として、独立成分毎に抽出されたアルファ波成分とベータ波成分との各々の振幅、独立成分間のアルファ波成分とベータ波成分の各々の位相差を組み合わせて用いている。ここで、時間区間は、例えば、刺激提示部11によって提示される刺激の切り替わりを区切りとして用いる。また、この時間区間は、任意の時間間隔として設定してもよいし、時間区間が互いに重複していても良い。
The feature
刺激特徴量記憶部14には、刺激提示部11によって提示される各刺激に対応して、後述する快適度、高揚度、好感度や音量などの固有の刺激特徴量が書き込まれて保持されている。この刺激特徴量は、刺激提示部11によって提示される各刺激について、この各刺激を受けたユーザーが、各刺激に応じて感じた評価値として快適度や音量などの固有の特徴量である。
Unique stimulus features such as comfort, uplifting, likability, and volume, which will be described later, are written and held in the stimulus
図3は、刺激特徴量記憶部14に記憶させる刺激と刺激特徴量との対応関係を求めるための質問用紙の構成例を示す図である。例えば、この質問用紙を用いることにより、刺激提示部11が提示するそれぞれの刺激に対してユーザーが感じる印象の程度をユーザー自身に評価させる主観評価によって得られる心理量を、上記刺激特徴量として得る。評価結果は、例えば、「不快」の0点から「快適」の4点までの5段階評価、「退屈」の0点から「楽しい」の4点までの5段階評価、「嫌い」の0点から「好き」の4点までの5段階評価とし、それぞれの評価値を当該刺激の刺激特徴量としての「快適度」、「高揚度」、「好感度」とする。
FIG. 3 is a diagram showing a configuration example of a questionnaire for obtaining a correspondence relationship between a stimulus stored in the stimulus feature
刺激特徴量の主観評価は脳波計測と同時に行っても良いし、脳波計測と同時に行なわずに、刺激を提示してその回答(主観)を取得する処理を、脳波計測の事前または事後に行っても良い。
同一の刺激に対する複数の脳波計測の対象ユーザーの各々の評価値を平均し、この平均した平均評価値を刺激特徴量として用いてもよい。また、脳波計測の対象ユーザー以外の、他のユーザーによる刺激に対する評価値を刺激特徴量として用いてもよい。
The subjective evaluation of the stimulus feature amount may be performed at the same time as the EEG measurement, or the process of presenting the stimulus and obtaining the answer (subjective) is performed before or after the EEG measurement without performing the EEG measurement at the same time. Is also good.
The evaluation values of each of the target users of a plurality of electroencephalogram measurements for the same stimulus may be averaged, and the average evaluation value may be used as the stimulus feature amount. Further, the evaluation value for the stimulus by a user other than the target user of the electroencephalogram measurement may be used as the stimulus feature amount.
また、上述したような楽曲の刺激によりユーザーが感じた主観評価のみでなく、聴覚刺激の音量や視覚刺激の輝度などの物理量や、人が感じる音量感の評価指標であるラウドネスレベルなどの物理量の強度などから導出される心理量を刺激特徴量として用いても良い。ここで、刺激特徴量は、刺激の物理量や心理量と線形的な対応関係を有する必要はなく、刺激の物理量の強度に応じて心理量が変化すれば良い。
例えば、ユーザーに与えた刺激の各々に対して、それぞれ任意の番号を割り当てて、刺激を識別するようにしても良い。このとき、音楽のジャンル、画像の被写体の種類など、何らかの基準によって刺激をグループに分類し、分類したグループの各々に、それぞれのグループを識別するための番号を割り当て、同一グループの刺激の刺激特徴量は同一値としても良い。
Moreover, not only the subjective evaluation felt by the user due to the stimulus of the music as described above, but also the physical quantity such as the volume of the auditory stimulus and the brightness of the visual stimulus, and the physical quantity such as the loudness level which is an evaluation index of the volume feeling felt by the person. A psychological quantity derived from intensity or the like may be used as a stimulus feature quantity. Here, the stimulus feature quantity does not have to have a linear correspondence with the physical quantity and the psychological quantity of the stimulus, and the psychological quantity may change according to the intensity of the physical quantity of the stimulus.
For example, an arbitrary number may be assigned to each of the stimuli given to the user to identify the stimulus. At this time, the stimuli are classified into groups according to some criteria such as the genre of music and the type of the subject of the image, and a number for identifying each group is assigned to each of the classified groups, and the stimulus features of the stimuli of the same group are assigned. The amounts may be the same value.
相互情報量算出部15は、特徴量算出部13によって求められた、各成分信号の各時間区間における脳波特徴量と、刺激特徴量記憶部14に記憶された前記時間区間に対応する刺激特徴量、との組合わせの関連の度合を示す相互情報量を求める。相互情報量に代えて、変数間の関係度を示す他の指標、例えば、相関係数を用いることもできる。この相互情報量は2つの確率変数の相互依存性の強さを示す指標であり、2つの確率変数XとYとの相互情報量は、以下の(1)式により求められる。以下の(1)式において、p(x)とp(y)との各々は、それぞれX、Yそれぞれの周辺確率関数である。また、p(x,y)は、XとYとの各々の同時確率関数である。
The mutual information
特徴量選択部16は、相互情報量算出部15が求めた相互情報量に基づいて、特徴量算出部13が算出した脳波特徴量のなかから、刺激特徴量の各々に対して関係性の強い脳波特徴量それぞれを特定し選択する。
刺激特徴量と脳波特徴量との間に線形的な関係が得られる場合には、相互情報量に代えて、刺激特徴量と脳波特徴量との相関係数を用いて、この相関係数の大きさにより、刺激特徴量の各々に対して関係性の強い脳波特徴量それぞれを特定し選択する構成としてもよい。
The feature
When a linear relationship is obtained between the stimulus feature and the EEG feature, the correlation coefficient between the stimulus feature and the EEG feature is used instead of the mutual information, and this correlation coefficient is used. Depending on the size, each of the brain wave features having a strong relationship with each of the stimulus features may be specified and selected.
ユーザー状態推定部17は、特徴量選択部16において選択された脳波特徴量と刺激特徴量との関係性を求める。すなわち、ユーザー状態推定部17は、特徴量選択部16が刺激特徴量の各々に対応して、それぞれ選択した脳波特徴量の関係性を、刺激特徴量の示す感情を推定する際の脳波特徴量として用いる。
これにより、ユーザー状態推定部17は、刺激を受けた際のユーザーの脳波特徴量から、刺激よって生じるユーザーの心理状態を推定することができる。
The user
As a result, the user
表示部18は、ユーザー状態推定部17が推定したユーザーの生理心理状態に関する情報を、数値、グラフ、マップなど、ユーザーや他の実験者が認識可能な形式で表示する画像を生成し、液晶ディスプレイなどの表示装置18aに表示する。例えば、表示部18は、刺激特徴量における「快適度」、「高揚度」、「好感度」などを、表示装置18aに対してグラフ表示することにより、受けた刺激に対応した推定される心理的な状態を示す。
また、表示部18は、脳波計測部19が測定した脳波波形、信号分離結果、特徴量算出部13が算出した脳波特徴量、相互情報量算出部15が算出した相互情報量、特徴量選択部16が選択した脳波特徴量に関する情報などの、原データや中間情報を表示してもよい。
The
Further, the
次に、本実施形態における信号分離処理の流れを図4を用いて説明する。図4は、図1における各部の処理及び各部の出力する信号の流れを示す概念図である。また、適宜、図1を参照して説明した本実施形態の構成例について補足説明を行う。観測信号41は、例えばユーザーの前頭部3箇所の脳波電極22で計測された、2チャンネルであるCh1及びCh2の各々の脳波をデジタル変換したデジタル脳波信号である。図4に示すように、観測信号41は、横軸を時間、縦軸を電圧とした脳波信号の波形として図示している。
Next, the flow of the signal separation process in the present embodiment will be described with reference to FIG. FIG. 4 is a conceptual diagram showing the processing of each part in FIG. 1 and the flow of signals output by each part. In addition, supplementary explanations will be given to the configuration examples of the present embodiment described with reference to FIG. The
刺激提示部11は、例えば、測定者が予め設定した順序に従って、楽曲を自動再生し、スピーカ11aを介してユーザーに、楽曲に対応した提示音を与えて(提示して)、音情報により聴覚に対する刺激を与える。
本実施形態では、人間が聞いた際に与える印象が異なると考えられる、例えば曲調及び曲想の各々が異なる3種類の楽曲を、3回づつ、計9回提示した場合を例として説明する。しかし、すでに述べたように、刺激として音以外、すなわち聴覚以外の刺激を用いる構成としても良い。また、聴覚などの外部刺激のみでなく、ユーザーに対して特定の場面や感情を想起させるよう指示を与えて、内部刺激を与える構成としても良い。
刺激特徴量記憶部14には、ユーザーに対して提示される刺激(本実施形態においては聴覚の刺激)の刺激特徴量45が予め書き込まれて記憶されている。
For example, the
In the present embodiment, an example will be described in which three types of music, which are considered to give different impressions when heard by a human being, for example, have different musical tunes and tunes, are presented three times, for a total of nine times. However, as already described, a configuration in which a stimulus other than sound, that is, a stimulus other than hearing may be used as the stimulus may be used. In addition to external stimuli such as hearing, the user may be instructed to evoke a specific scene or emotion to give an internal stimulus.
In the stimulus feature
図5は、刺激特徴量記憶部14に記憶されている刺激特徴量の構成例を示すテーブルの図である。図5に示すように、本実施形態において、刺激特徴量45は、楽曲の音源番号、快適度、テンポ、音量感、ジャンルであり、9回の刺激提示の時間を示す区間それぞれに対応する値を保持している。音源番号は、3種類の楽曲のいずれであるかを識別するための識別情報である。快適度は、図3に示す質問用紙から得た、ユーザーの刺激に対する主観評価に基づく心理状態を数値化した心理量を示す。テンポは、楽曲の音声信号から物理的に計測された、楽曲において拍をどの程度の速さで打つかを示す物理量である。音量感は、楽曲の音声信号の周波数帯毎の振幅からラウドネス曲線を用いて抽出した心理量である。ジャンルは、例えば、クラッシックを「1」、ポップスを「2」などとした楽曲の分類を示している。
FIG. 5 is a diagram of a table showing a configuration example of the stimulus feature amount stored in the stimulus feature
例えば、テーブルの1行目の区間「1」においては、音源番号が「1」であることを示している。また、音源番号「1」の音源である楽曲は、快適度が「4.2」であり、テンポが「60」であり、音量感が「2.9」であり、ジャンルが「1(クラッシック)」であることを示している。同様に、テーブルの2行目の区間「2」においては、音源番号が「2」であることを示している。また、音源番号「2」の音源である楽曲は、快適度が「2.1」であり、テンポが「80」であり、音量感が「4.1」であり、ジャンルが「1(ポップス)」であることを示している。 For example, in the section "1" in the first row of the table, it is shown that the sound source number is "1". The music that is the sound source of the sound source number "1" has a comfort level of "4.2", a tempo of "60", a volume feeling of "2.9", and a genre of "1 (classical)". ) ”. Similarly, in the section "2" in the second row of the table, it is shown that the sound source number is "2". The music that is the sound source of the sound source number "2" has a comfort level of "2.1", a tempo of "80", a volume feeling of "4.1", and a genre of "1 (pops)". ) ”.
図4に戻り、脳波計測部19は、刺激提示部11がユーザーに刺激を提示した際、ユーザーの頭部に装着された脳波電極22を介して、脳波信号の観測信号41を取得し、信号分離部12に対して出力する。本実施形態においては、脳波計測部19が、脳波電極22の各々に対応したチャンネルCh1及びチャンネルCh2それぞれの観測信号41を取得する。ここで、すでに述べたように、脳波計測部19は、チャンネルCh1及びチャンネルCh2それぞれの観測信号をAD変換器25により所定周期でサンプリングし、アナログ信号をデジタル変換したデジタル信号として観測信号41を信号分離部12へ出力する。また、一般に、頭部における近接した脳波電極から検出される脳波信号には、その特徴(信号波形やパワースペクトルなど)に強い相関が存在する。さらに、脳波電極の数を増加させて、頭部の広い領域をカバーするように脳波信号の観測点を配置することにより、独立成分分析による独立成分がより細かい単位で抽出されるため、得られる情報量も増加することになる。
Returning to FIG. 4, when the
信号分離部12は、脳波特徴量の導出に用いるため、観測信号41を成分信号43に独立成分分析を用いて分離する。本実施形態においては、独立成分分析によって、複数チャンネル(チャンネルCh1、Ch2)の観測信号41を、脳内信号源に対応すると推定される相互に独立した信号である独立成分42に分離する。この独立成分42は、独立成分1(IC1)と独立成分2(IC2)との各々である。一般に、独立成分間の相関は低いため、得られた独立成分1(IC1)と独立成分2(IC2)との各々の相関も低い。
The
図6は、独立成分42における独立成分1(IC1)と独立成分2(IC2)とのデータの一例を示す図である。この図6に示す独立成分1(IC1)のグラフと独立成分2(IC2)のグラフとの各々において、縦軸が独立成分の電圧を示し、横軸が時間を示している。この独立成分1(IC1)と独立成分2(IC2)との各々の信号波形は、ユーザーの脳内信号源の各々が発生した脳波信号の波形、すなわち脳内信号源の各々における脳波信号の特徴を示している。
本実施形態においては、独立成分の信号分離方法として、独立成分分析と周波数フィルタとを用いているが、これらに限るものではなく、独立成分分析の代わりに主成分分析などを用い、周波数フィルタの代わりにフーリエ変換などを用いた構成としてもよい。
FIG. 6 is a diagram showing an example of data of the independent component 1 (IC1) and the independent component 2 (IC2) in the
In the present embodiment, the independent component analysis and the frequency filter are used as the signal separation method of the independent component, but the present invention is not limited to these, and the principal component analysis or the like is used instead of the independent component analysis to obtain the frequency filter. Alternatively, a configuration using a Fourier transform or the like may be used.
図4に戻り、信号分離部12は、求めた独立成分42の独立成分1(IC1)と独立成分2(IC2)との各々を、周波数フィルタ(デジタルフィルタ)により、それぞれアルファ成分、ベータ成分とに分離し、成分信号43を算出する。
Returning to FIG. 4, the
図7は、独立成分42を周波数フィルタにより分離して求めた成分信号43のデータの一例を示す図である。この図6に示すIC1アルファ波及びIC1ベータ波の各々は、独立成分1(IC1)を周波数フィルタにより分離した信号の波形である。また、IC2アルファ波及びIC2ベータ波の各々は、独立成分2(IC2)を周波数フィルタにより分離した信号の波形である。IC1アルファ波、IC1ベータ波、IC2アルファ波及びIC2ベータ波の各々において、縦軸が成分信号の電圧を示し、横軸が時間を示している。上記IC1アルファ波、IC1ベータ波、IC2アルファ波及びIC2ベータ波の各々の信号は、ユーザーの脳内信号源の各々における脳波信号を示す独立成分42の信号に対して、周波数というパラメータを追加して、独立成分42より細分化した特徴としたものである。
FIG. 7 is a diagram showing an example of data of the
図4に戻り、特徴量算出部13は、信号分離部12によって分離された成分信号43におけるIC1アルファ波、IC1ベータ波、IC2アルファ波及びIC2ベータ波の各々から、脳波信号の特徴を示す脳波特徴量44を求める。
Returning to FIG. 4, the feature
図8は、IC1アルファ波、IC1ベータ波、IC2アルファ波及びIC2ベータ波の各々から求めた脳波特徴量44の構成例を示す図である。
図8(a)は、IC1アルファ波とIC1ベータ波との各々の脳波特徴量のグラフを示している。この脳波特徴量のグラフは、横軸が時間(図5に示す時間の区間「1」から区間「9」まで)を示しており、縦軸が時間の区間における信号波形の振幅を平均した平均振幅を示している。
FIG. 8 is a diagram showing a configuration example of an
FIG. 8A shows a graph of each EEG feature amount of the IC1 alpha wave and the IC1 beta wave. In this graph of electroencephalogram features, the horizontal axis indicates time (from the time interval "1" to the interval "9" shown in FIG. 5), and the vertical axis is the average of the amplitudes of the signal waveforms in the time interval. It shows the amplitude.
図8(b)は、IC2アルファ波とIC2ベータ波との各々の脳波特徴量のグラフを示している。この脳波特徴量のグラフは、横軸が時間(図5に示す時間の区間「1」から区間「9」まで)を示しており、縦軸が時間の区間における信号波形の振幅を平均した振幅を示している。 FIG. 8B shows a graph of each electroencephalogram feature amount of the IC2 alpha wave and the IC2 beta wave. In this graph of electroencephalogram features, the horizontal axis indicates time (from the time interval "1" to the interval "9" shown in FIG. 5), and the vertical axis is the amplitude obtained by averaging the amplitudes of the signal waveforms in the time interval. Is shown.
図8(c)は、独立成分間の位相差、すなわちICIアルファ波及びIC2アルファ波の波形における位相差と、ICIベータ波及びIC2ベータ波の波形における位相差との各々を示すグラフである。この脳波特徴量のグラフは、横軸が時間(図5に示す時間の区間「1」から区間「9」まで)を示しており、縦軸が時間の区間における信号波形の位相差を平均した位相差を示している。 FIG. 8C is a graph showing the phase difference between the independent components, that is, the phase difference in the waveforms of the ICI alpha wave and the IC2 alpha wave, and the phase difference in the waveforms of the ICI beta wave and the IC2 beta wave. In this EEG feature graph, the horizontal axis indicates time (from the time interval "1" to the interval "9" shown in FIG. 5), and the vertical axis averages the phase differences of the signal waveforms in the time interval. It shows the phase difference.
図4に戻り、特徴量算出部13は、図8に示すように、時間の区間毎のIC1アルファ波、IC1ベータ波、IC2アルファ波及びIC2ベータ波の各々の振幅と、時間区間毎のICIアルファ波及びIC2アルファ波の位相差と、ICIベータ波及びIC2ベータ波の位相差との合計6個からなる脳波信号の脳波特徴量44を求める。
Returning to FIG. 4, as shown in FIG. 8, the feature
相互情報量算出部15は、特徴量算出部13により求められた脳波信号の脳波特徴量44の各々と、刺激特徴量記憶部14に記憶されている刺激特徴量45との相互情報量46を求める。ここで、相互情報量算出部15は、刺激特徴量45と脳波特徴量44との各々の時間の区間毎を対応させて、相互情報量46の算出を上記(1)式により求める。
本実施形態において、相互情報量算出部15は、例えば、音源番号と快適度とテンポと音量感とジャンルに対応して値が変動する特徴量を脳波特徴量44から特定する場合、5種類の刺激特徴量45と6種類の脳波特徴量44との各々の組合わせについて、脳波特徴量44をXとし、刺激特徴量45をYとして、(1)式により相互情報量46を算出する。
The mutual information
In the present embodiment, when the mutual information
図9は、刺激特徴量45と脳波特徴量44との各々の時間の区間毎を対応させて求めた相互情報量を示す図である。図9(a)に示すテーブルは、音源番号と脳波特徴量44との相互情報量を示している。テーブルの欄101には、音源番号と、脳波特徴量44におけるIC1アルファ波の振幅(図8(a)の左側のグラフ)との相互情報量「0.411」が示されている。テーブルの欄102には、音源番号と、脳波特徴量44におけるIC1ベータ波の振幅(図8(a)の右側のグラフ)との相互情報量「0.351」が示されている。テーブルの欄103には、音源番号と、脳波特徴量44におけるIC2アルファ波の振幅(図8(b)の左側のグラフ)との相互情報量「0.211」が示されている。テーブルの欄104には、音源番号と、脳波特徴量44におけるIC2ベータ波の振幅(図8(b)の右側のグラフ)との相互情報量「0.292」が示されている。テーブルの欄105には、音源番号と、脳波特徴量44におけるIC1アルファ波とIC2アルファ波との位相差(図8(c)の左側のグラフ)との相互情報量「0.177」が示されている。テーブルの欄106には、音源番号と、脳波特徴量44におけるIC1ベータ波とIC2ベータ波との位相差(図8(c)の右側のグラフ)との相互情報量「0.225」が示されている。
FIG. 9 is a diagram showing a mutual information amount obtained by associating each time interval between the
図9(b)に示すテーブルは、快適度と脳波特徴量44との相互情報量を示している。図9(c)に示すテーブルは、テンポと脳波特徴量44との相互情報量を示している。図9(d)に示すテーブルは、音量感と脳波特徴量44との相互情報量を示している。図9(e)に示すテーブルは、ジャンル(ジャンル1)と脳波特徴量44との相互情報量を示している。これら、図9(b)、図9(c)、図9(d)、図9(e)の各々も、図1(a)と同様に、刺激特徴量45の各々と、脳波特徴量44それぞれとの相互情報量がテーブルの各欄に示されている。
The table shown in FIG. 9B shows the amount of mutual information between the comfort level and the
図4に戻り、特徴量選択部16は、図9に示す5種類の刺激特徴量45と6種類の脳波特徴量44との各々の組合わせ毎の相互情報量のなかから、相互情報量が予め設定された判定閾値を超える組合わせを抽出する。この判定閾値は、刺激特徴量と脳波特徴量との関係性の有無を判定する相互情報量の閾値である。実施形態においては、この判定閾値は、例えば、0.5に設定している。このため、特徴量選択部16は、図9に示す5種類の刺激特徴量45と6種類の脳波特徴量44との各々の組合わせ毎の相互情報量のなかから、相互情報量が「0.624」である快適度とIC1アルファ波振幅との組合わせと、相互情報量が「0.656」であるジャンルとIC1アルファ波振幅との組合わせとを、相互情報量が判定閾値を超えているため抽出する。そして、特徴量選択部16は、独立成分におけるIC1アルファ波振幅を、ユーザーの快適度及び楽曲のジャンルの推定に用いる脳波特徴量47として選択する。
Returning to FIG. 4, the feature
特徴量選択部16は、脳波推定データベース20に対して、指標として得た刺激特徴量とこの刺激特徴量に対して関係性がある脳波特徴量とを組として、脳波推定データベース20に対して、書き込んで記憶させる。
ユーザー状態推定部17は、特徴量選択部16が選択した脳波特徴量47、すなわち脳波信号の独立成分におけるIC1アルファ波振幅を、ユーザーの快適度及び楽曲のジャンルの推定の指標として用いる。
The feature
The user
図10は、刺激特徴量と脳波特徴量(IC1アルファ波振幅)との対応関係を示すグラフの図である。図10(a)は、刺激特徴量の快適度と脳波特徴量のアルファ波振幅との対応関係を示している。図10(a)のグラフからは、アルファ波振幅(IC1アルファ波振幅)が大きいほどユーザーが快適な状態であることを推定できることが判る。
同様に、図10(b)は、刺激特徴量のジャンルと脳波特徴量のアルファ波振幅との対応関係を示している。図10(b)のグラフからは、アルファ波振幅(IC1アルファ波振幅)が大きい場合には、ユーザーが視聴している曲のジャンルがジャンル2であると推定できることが判る。この図10の関係から、アルファ波振幅(IC1アルファ波振幅)が大きいほど、ユーザーは快適な状態にあると推定され、アルファ波の振幅が予め設定した推定閾値を超える場合には、ユーザーはジャンル2の楽曲を聴取していると推定することができる。この推定閾値は、例えば、楽曲のジャンルに対する複数のユーザーのIC1アルファ波振幅を取得し、IC1アルファ波振幅の分散などから求める。
FIG. 10 is a graph showing the correspondence between the stimulus feature amount and the electroencephalogram feature amount (IC1 alpha wave amplitude). FIG. 10A shows the correspondence between the comfort level of the stimulus feature and the alpha wave amplitude of the EEG feature. From the graph of FIG. 10A, it can be seen that the larger the alpha wave amplitude (IC1 alpha wave amplitude), the more comfortable the user is.
Similarly, FIG. 10B shows the correspondence between the genre of stimulus features and the alpha wave amplitude of EEG features. From the graph of FIG. 10B, it can be seen that when the alpha wave amplitude (IC1 alpha wave amplitude) is large, it can be estimated that the genre of the song being viewed by the user is
以上から、本実施形態によれば、独立成分のIC1アルファ波振幅が、ユーザーの快適度と聴取している楽曲のジャンルを推定するために有用な脳波特徴量であり、独立成分のIC1アルファ波振幅と、ユーザーの快適度及び聴取している楽曲のジャンルそれぞれとの関係が求められる。
また、本実施形態によれば、上述した脳波推定データベース20を参照することにより、ユーザーが刺激特徴量が未知の楽曲を聴取している場合においても、ユーザー状態推定部17が特徴量算出部13において求められた脳波特徴量から、ユーザーの快適度と聴取している楽曲のジャンルを推定することができる。すなわち、本実施形態においては、脳波推定データベース20を設けて、ユーザーの刺激特徴量と脳波特徴量との関係を蓄積することにより、刺激特徴量記憶部14と相互情報量算出部15と特徴量選択部16との各々の処理を省略し、未知の楽曲に対するユーザーの快適度と聴取している楽曲のジャンルとの推定を行なうことができる。
From the above, according to the present embodiment, the IC1 alpha wave amplitude of the independent component is an electroencephalogram feature amount useful for estimating the comfort level of the user and the genre of the music being listened to, and the IC1 alpha wave of the independent component. The relationship between the amplitude, the comfort level of the user, and the genre of the music being listened to is required.
Further, according to the present embodiment, by referring to the above-mentioned
以上のように、本実施形態のユーザー状態推定装置1(脳波信号処理システム)は、時系列的にユーザーに感覚刺激を提示する刺激提示部11と、頭部に装着された脳波電極22を用いて観測信号41を取得する脳波計測部19と、観測信号41を周波数分割や多変量解析などによって成分信号43に分割する信号分離部12と、成分信号43の脳波特徴量44を求める特徴量算出部13と、脳波特徴量44と刺激特徴量45の相互情報量を求める相互情報量算出部15と、相互情報量に基づいて刺激特徴量と関係性の強い脳波特徴量を特定し選択する特徴量選択部16と、選択された脳波特徴量と刺激特徴量との関係性からユーザーの状態を推定する、ユーザー状態推定部17と、ユーザーの刺激特徴量と脳波特徴量との関係を蓄積する、脳波推定データベース20との各々を備える。
As described above, the user state estimation device 1 (electroencephalogram signal processing system) of the present embodiment uses the
図11は、楽曲と脳波信号の独立成分におけるベータ波の振幅との対応を示すグラフである。図11においては、横軸が時間を示しており、縦軸がベータ波の振幅を示している。
ここで、信号分離部12が4種の楽曲を繰り返して連続聴取した際の脳波信号を19個の脳波電極22で計測し、計測した脳波信号に対して独立成分分析を用いて分離した19個の独立成分の各々のベータ波の振幅変化が図11に示されている。
FIG. 11 is a graph showing the correspondence between the amplitude of the beta wave in the independent component of the music and the electroencephalogram signal. In FIG. 11, the horizontal axis represents time and the vertical axis represents the amplitude of beta waves.
Here, the electroencephalogram signals when the
図12は、図3のアンケート用紙を用いて予め求めた楽曲毎の興奮度を示すテーブルである。楽曲「1」の興奮度が「0.8」であり、楽曲「2」の興奮度が「3.2」であり、楽曲「3」の興奮度が「3.7」であり、楽曲「4」の興奮度が「0.5」であることが示されている。 FIG. 12 is a table showing the degree of excitement for each song obtained in advance using the questionnaire sheet of FIG. The excitement of the song "1" is "0.8", the excitement of the song "2" is "3.2", the excitement of the song "3" is "3.7", and the song "2". It is shown that the degree of excitement of "4" is "0.5".
図13は、図11に示す19個の独立成分の各々のベータ波の振幅と興奮度とから求めた相互情報量のテーブルを示す図である。テーブルの項目は、順位、成分及び相互情報量の各々である。成分は、0から18までの19個の独立成分の各々の区別を示している。相互情報量は、興奮度と各独立成分のベータ波の振幅との相互情報量とを示している。順位は、興奮度と独立成分との相互情報量が大きい方から小さい方への順番を示している。すなわち、最も相互情報量(0.4)が大きい独立成分「9」の順位が「1」であり、最も相互情報量(0.09)が小さい独立成分「3」の順位が「19」である。 FIG. 13 is a diagram showing a table of mutual information obtained from the amplitude and excitement of the beta waves of each of the 19 independent components shown in FIG. The items in the table are each of rank, component and mutual information. The components show the distinction of each of the 19 independent components from 0 to 18. The mutual information indicates the mutual information between the degree of excitement and the amplitude of the beta wave of each independent component. The ranking indicates the order from the one with the largest amount of mutual information between the degree of excitement and the independent component to the one with the smallest amount of mutual information. That is, the rank of the independent component "9" having the largest mutual information (0.4) is "1", and the rank of the independent component "3" having the smallest mutual information (0.09) is "19". is there.
図14は、図11から抽出した、図13において相互情報量が最も大きい独立成分「9」のベータ波の振幅を示す図である。図14においては、横軸が時間を示しており、縦軸がベータ波の振幅を示している。
この図14の結果から、図12に示す取得した楽曲による興奮度が大きい程、ユーザーのベータ波の振幅が大きく、興奮度が小さいほどベータ波の振幅が小さいことが判り、この独立成分「9」であるベータ波の振幅が興奮度の推定に利用できることが判る。
FIG. 14 is a diagram showing the amplitude of the beta wave of the independent component “9” having the largest mutual information in FIG. 13 extracted from FIG. In FIG. 14, the horizontal axis represents time and the vertical axis represents the amplitude of beta waves.
From the results of FIG. 14, it can be seen that the greater the degree of excitement of the acquired music shown in FIG. 12, the greater the amplitude of the user's beta wave, and the smaller the degree of excitement, the smaller the amplitude of the beta wave. It can be seen that the amplitude of the beta wave, which is "", can be used to estimate the degree of excitement.
図15は、図11から抽出した、図13において相互情報量が最も小さい独立成分「3」のベータ波の振幅を示す図である。図15においては、横軸が時間を示しており、縦軸がベータ波の振幅を示している。
この図15の結果から、図12に示す取得した楽曲による興奮度と、ベータ波の振幅との関連性(相関する変化)が観察されないため、この独立成分「3」であるベータ波の振幅が興奮度の推定に利用できないことが判る。
FIG. 15 is a diagram showing the amplitude of the beta wave of the independent component “3” having the smallest mutual information in FIG. 13 extracted from FIG. In FIG. 15, the horizontal axis represents time and the vertical axis represents the amplitude of beta waves.
From the result of FIG. 15, since the relation (correlated change) between the excitement degree by the acquired music shown in FIG. 12 and the amplitude of the beta wave is not observed, the amplitude of the beta wave which is the independent component “3” is It turns out that it cannot be used to estimate the degree of excitement.
なお、本発明における図1の脳波信号処理システムの機能を実現するためのプログラムをコンピュータ読み取り可能な記録媒体に記録して、この記録媒体に記録されたプログラムをコンピュータシステムに読み込ませ、実行することにより脳波信号から受けている刺激を推定するための情報を求める処理を行ってもよい。なお、こでいう「コンピュータシステム」とは、OSや周辺機器等のハードウェアを含むものとする。
また、「コンピュータシステム」は、ホームページ提供環境(あるいは表示環境)を備えたWWWシステムも含むものとする。また、「コンピュータ読み取り可能な記録媒体」とは、フレキシブルディスク、光磁気ディスク、ROM、CD−ROM等の可搬媒体、コンピュータシステムに内蔵されるハードディスク等の記憶装置のことをいう。さらに「コンピュータ読み取り可能な記録媒体」とは、インターネット等のネットワークや電話回線等の通信回線を介してプログラムが送信された場合のサーバやクライアントとなるコンピュータシステム内部の揮発性メモリ(RAM)のように、一定時間プログラムを保持しているものも含むものとする。
The program for realizing the function of the brain wave signal processing system of FIG. 1 in the present invention is recorded on a computer-readable recording medium, and the program recorded on the recording medium is read into the computer system and executed. You may perform the process of obtaining the information for estimating the stimulus received from the brain wave signal. The term "computer system" as used herein includes hardware such as an OS and peripheral devices.
Further, the "computer system" shall also include a WWW system provided with a homepage providing environment (or display environment). Further, the "computer-readable recording medium" refers to a portable medium such as a flexible disk, a magneto-optical disk, a ROM, or a CD-ROM, or a storage device such as a hard disk built in a computer system. Furthermore, a "computer-readable recording medium" is a volatile memory (RAM) inside a computer system that serves as a server or client when a program is transmitted via a network such as the Internet or a communication line such as a telephone line. In addition, it shall include those that hold the program for a certain period of time.
また、上記プログラムは、このプログラムを記憶装置等に格納したコンピュータシステムから、伝送媒体を介して、あるいは、伝送媒体中の伝送波により他のコンピュータシステムに伝送されてもよい。ここで、プログラムを伝送する「伝送媒体」は、インターネット等のネットワーク(通信網)や電話回線等の通信回線(通信線)のように情報を伝送する機能を有する媒体のことをいう。また、上記プログラムは、前述した機能の一部を実現するためのものであっても良い。さらに、前述した機能をコンピュータシステムにすでに記録されているプログラムとの組み合わせで実現できるもの、いわゆる差分ファイル(差分プログラム)であっても良い。 Further, the program may be transmitted from a computer system in which this program is stored in a storage device or the like to another computer system via a transmission medium or by a transmission wave in the transmission medium. Here, the "transmission medium" for transmitting a program refers to a medium having a function of transmitting information, such as a network (communication network) such as the Internet or a communication line (communication line) such as a telephone line. Further, the above program may be for realizing a part of the above-mentioned functions. Further, a so-called difference file (difference program) may be used, which can realize the above-mentioned functions in combination with a program already recorded in the computer system.
ここまで、本発明の実施形態について説明したが、上記実施形態はあくまで一例であり、本発明は上述した実施形態に限定されず、その技術的思想の範囲内において種々異なる形態にて実施されてよいことは言うまでもない。
また、本発明の範囲は、図示され記載された例示的な実施形態に限定されるものではなく、本発明が目的とするものと均等な効果をもたらすすべての実施形態をも含む。さらに、本発明の範囲は、各請求項により画される発明の特徴の組み合わせに限定されるものではなく、すべての開示されたそれぞれの特徴のうち特定の特徴のあらゆる所望する組み合わせを含む。
Although the embodiments of the present invention have been described so far, the above-described embodiments are merely examples, and the present invention is not limited to the above-described embodiments, and is implemented in various different forms within the scope of the technical idea. Needless to say, it's good.
Also, the scope of the present invention is not limited to the exemplary embodiments illustrated and described, but also includes all embodiments that provide an effect equal to that intended by the present invention. Furthermore, the scope of the present invention is not limited to the combination of the features of the invention defined by each claim, but includes any desired combination of the specific features of all the disclosed features.
10…コンピュータ
11…刺激提示部
12…信号分離部
13…特徴量算出部
14…刺激特徴量記憶部
15…相互情報量算出部
16…特徴量選択部
17…ユーザー状態推定部
18…表示部
19…脳波計測部
20…脳波推定データベース
21…接地電極
22…脳波電極
24…増幅器
25…AD変換器
10 ...
Claims (14)
前記刺激に対応した前記脳波信号の複数の前記脳波特徴量を求める特徴量算出部と、
前記刺激の1以上の前記刺激特徴量の各々と、複数の前記脳波特徴量それぞれとの組合わせ毎に、当該組合わせの関連の度合を示す相互関連度を算出する相互関連度算出部と、
前記相互関連度が所定の閾値を超える、前記刺激特徴量と前記脳波特徴量との組合わせを抽出し、前記刺激特徴量に対して関連性の高い前記脳波特徴量を選択する特徴量選択部と、
前記脳波信号を複数系列の脳波の特徴を有する成分信号に分離する信号分離部と
を備え、
前記特徴量算出部が、前記成分信号から前記脳波特徴量を算出する
ことを特徴とする脳波信号処理システム。 It is an electroencephalogram signal processing system that obtains the relationship between the electroencephalogram feature amount, which is the feature amount of the electroencephalogram signal generated in response to a given stimulus, and the stimulus feature amount, which is the feature amount of the stimulus.
A feature amount calculation unit for obtaining a plurality of the brain wave feature amounts of the brain wave signal corresponding to the stimulus, and a feature amount calculation unit.
For each combination of one or more of the stimulus features of the stimulus and each of the plurality of electroencephalogram features, a degree of interrelationship calculation unit for calculating the degree of interrelationship indicating the degree of relevance of the combination, and
A feature amount selection unit that extracts a combination of the stimulus feature amount and the brain wave feature amount whose degree of interrelationship exceeds a predetermined threshold value and selects the brain wave feature amount having a high relevance to the stimulus feature amount. When,
It is provided with a signal separation unit that separates the electroencephalogram signal into component signals having the characteristics of a plurality of series of electroencephalograms.
An electroencephalogram signal processing system characterized in that the feature amount calculation unit calculates the brain wave feature amount from the component signal.
前記刺激の各々に対応して、脳内信号源それぞれが発した前記脳波信号の周波数に基づく複数の前記脳波特徴量を求める特徴量算出部と、
前記刺激の1以上の前記刺激特徴量の各々と、複数の前記脳波特徴量それぞれとの組合わせ毎に、当該組合わせの関連の度合を示す相互関連度を算出する相互関連度算出部と、
前記相互関連度が所定の閾値を超える、前記刺激特徴量と前記脳波特徴量との組合わせを抽出し、前記刺激特徴量に対して関連性の高い前記脳波特徴量を選択する特徴量選択部と
を備えることを特徴とする脳波信号処理システム。 It is an electroencephalogram signal processing system that obtains the relationship between the electroencephalogram feature amount, which is the feature amount of the electroencephalogram signal generated in response to a given stimulus, and the stimulus feature amount, which is the feature amount of the stimulus.
A feature amount calculation unit that obtains a plurality of the electroencephalogram feature amounts based on the frequencies of the electroencephalogram signals emitted by each of the signal sources in the brain in response to each of the stimuli.
For each combination of one or more of the stimulus features of the stimulus and each of the plurality of electroencephalogram features, a degree of interrelationship calculation unit for calculating the degree of interrelationship indicating the degree of relevance of the combination, and
A feature amount selection unit that extracts a combination of the stimulus feature amount and the brain wave feature amount whose degree of interrelationship exceeds a predetermined threshold value and selects the brain wave feature amount having a high relevance to the stimulus feature amount. An electroencephalogram signal processing system characterized by being equipped with.
前記刺激に対応した前記脳波信号の複数の前記脳波特徴量を求める特徴量算出部と、
前記刺激の1以上の前記刺激特徴量の各々と、複数の前記脳波特徴量それぞれとの組合わせ毎に、当該組合わせの関連の度合を示す相互関連度を算出する相互関連度算出部と、
前記相互関連度が所定の閾値を超える、前記刺激特徴量と前記脳波特徴量との組合わせを抽出し、前記刺激特徴量に対して関連性の高い前記脳波特徴量を選択する特徴量選択部と
を備え、
前記複数の前記脳波特徴量のうち、少なくとも2以上の前記脳波特徴量は互いに異なる種類の単位である
ことを特徴とする脳波信号処理システム。 Among the EEG features, which are the features of the EEG signal generated in response to a given stimulus, the EEG that selects the EEG features that are highly relevant to each of the Stimulus features, which are the features of the stimulus. It ’s a signal processing system,
A feature amount calculation unit for obtaining a plurality of the brain wave feature amounts of the brain wave signal corresponding to the stimulus, and a feature amount calculation unit.
For each combination of one or more of the stimulus features of the stimulus and each of the plurality of electroencephalogram features, a degree of interrelationship calculation unit for calculating the degree of interrelationship indicating the degree of relevance of the combination, and
A feature amount selection unit that extracts a combination of the stimulus feature amount and the brain wave feature amount whose degree of interrelationship exceeds a predetermined threshold value and selects the brain wave feature amount having a high relevance to the stimulus feature amount. With and
Of the plurality of electroencephalogram features, at least two or more of the electroencephalogram features are units of different types from each other.
EEG signal processing system comprising a call.
ことを特徴とする請求項1から請求項3のいずれか一項に記載の脳波信号処理システム。 The electroencephalogram signal processing system according to any one of claims 1 to 3, wherein a mutual information amount is used as the degree of interrelationship.
ことを特徴とする請求項1から請求項3のいずれか一項に記載の脳波信号処理システム。 The electroencephalogram signal processing system according to any one of claims 1 to 3, wherein a correlation coefficient is used as the degree of interrelationship.
をさらに有し、
前記刺激提示部が、種類の異なる前記刺激特徴量を有する前記刺激をユーザーに対して提示する
ことを特徴とする請求項1から請求項5のいずれか一項に記載の脳波信号処理システム。 It also has a stimulus presentation unit that presents the stimulus to the user.
The electroencephalogram signal processing system according to any one of claims 1 to 5, wherein the stimulus presenting unit presents the stimulus having a different type of stimulus feature amount to the user.
さらに備え、
前記特徴量選択部が、前記刺激特徴量に対して選択した、当該刺激特徴量に対して関連性の高い前記脳波特徴量を、前記刺激特徴量と組として、前記脳波推定データベースに書き込んで記憶させる
ことを特徴とする請求項1から請求項6のいずれか一項に記載の脳波信号処理システム。 With an additional EEG estimation database
The electroencephalogram feature amount selected by the feature amount selection unit with respect to the stimulus feature amount and having a high relevance to the stimulus feature amount is written and stored in the electroencephalogram estimation database as a set with the stimulus feature amount. The electroencephalogram signal processing system according to any one of claims 1 to 6, wherein the electroencephalogram signal processing system is characterized.
さらに備える
ことを特徴とする請求項7に記載の脳波信号処理システム。 When the brain wave feature amount of the brain wave signal is supplied, the stimulus feature amount corresponding to the brain wave feature amount is extracted from the brain wave estimation database, and a user state estimation unit that estimates the stimulus feature amount of the stimulus is further added. The electroencephalogram signal processing system according to claim 7, further comprising.
与えられた刺激に対応して発生する脳波信号の特徴量である脳波特徴量、前記刺激の特徴量である刺激特徴量の関係性を求める脳波信号処理方法であって、
特徴量算出部が、前記刺激に対応した前記脳波信号の複数の前記脳波特徴量を求める特徴量算出過程と、
相互関連度算出部が、前記刺激の1以上の前記刺激特徴量の各々と、複数の前記脳波特徴量それぞれとの組合わせ毎に、当該組合わせの関連の度合を示す相互関連度を算出する相互関連度算出過程と、
特徴量選択部が、前記相互関連度が所定の閾値を超える、前記刺激特徴量と前記脳波特徴量との組合わせを抽出し、前記刺激特徴量に対して関連性の高い前記脳波特徴量を選択する特徴量選択過程と、
信号分離部が、前記脳波信号を複数系列の脳波の特徴を有する成分信号に分離する信号分離過程と
を含み、
前記特徴量算出部が、前記成分信号から前記脳波特徴量を算出する
ことを特徴とする脳波信号処理方法。 It is an electroencephalogram signal processing method for obtaining the relationship between the electroencephalogram feature amount, which is the feature amount of the electroencephalogram signal generated in response to a given stimulus, and the stimulus feature amount, which is the feature amount of the stimulus.
It is an electroencephalogram signal processing method for obtaining the relationship between the electroencephalogram feature amount, which is the feature amount of the electroencephalogram signal generated in response to a given stimulus, and the stimulus feature amount, which is the feature amount of the stimulus.
A feature amount calculation process in which the feature amount calculation unit obtains a plurality of the brain wave feature amounts of the brain wave signal corresponding to the stimulus, and a feature amount calculation process.
The interrelationship calculation unit calculates the degree of interrelationship indicating the degree of association of the combination for each combination of one or more of the stimulus features of the stimulus and each of the plurality of electroencephalogram features. Interrelationship calculation process and
The feature amount selection unit extracts a combination of the stimulus feature amount and the brain wave feature amount whose degree of interrelationship exceeds a predetermined threshold value, and obtains the brain wave feature amount which is highly related to the stimulus feature amount. Feature selection process to be selected and
The signal separation unit includes a signal separation process of separating the brain wave signal into component signals having the characteristics of a plurality of series of brain waves.
An electroencephalogram signal processing method characterized in that the feature amount calculation unit calculates the brain wave feature amount from the component signal.
特徴量算出部が、前記刺激の各々に対応して、脳内信号源それぞれが発した前記脳波信号の周波数に基づく複数の前記脳波特徴量を求める特徴量算出過程と、
相互関連度算出部が、前記刺激の1以上の前記刺激特徴量の各々と、複数の前記脳波特徴量それぞれとの組合わせ毎に、当該組合わせの関連の度合を示す相互関連度を算出する相互関連度算出過程と、
特徴量選択部が、前記相互関連度が所定の閾値を超える、前記刺激特徴量と前記脳波特徴量との組合わせを抽出し、前記刺激特徴量に対して関連性の高い前記脳波特徴量を選択する特徴量選択過程と
を含むことを特徴とする脳波信号処理方法。 It is an electroencephalogram signal processing method for obtaining the relationship between the electroencephalogram feature amount, which is the feature amount of the electroencephalogram signal generated in response to a given stimulus, and the stimulus feature amount, which is the feature amount of the stimulus.
A feature amount calculation process in which the feature amount calculation unit obtains a plurality of the brain wave feature amounts based on the frequencies of the brain wave signals emitted by each of the signal sources in the brain in response to each of the stimuli.
The interrelationship calculation unit calculates the degree of interrelationship indicating the degree of association of the combination for each combination of one or more of the stimulus features of the stimulus and each of the plurality of electroencephalogram features. Interrelationship calculation process and
The feature amount selection unit extracts a combination of the stimulus feature amount and the brain wave feature amount whose degree of interrelationship exceeds a predetermined threshold, and obtains the brain wave feature amount having a high relevance to the stimulus feature amount. An electroencephalogram signal processing method characterized by including a feature selection process to be selected.
特徴量算出部が、前記刺激に対応した前記脳波信号の複数の前記脳波特徴量を求める特徴量算出過程と、
相互関連度算出部が、前記刺激の1以上の前記刺激特徴量の各々と、複数の前記脳波特徴量それぞれとの組合わせ毎に、当該組合わせの関連の度合を示す相互関連度を算出する相互関連度算出過程と、
特徴量選択部が、前記相互関連度が所定の閾値を超える、前記刺激特徴量と前記脳波特徴量との組合わせを抽出し、前記刺激特徴量に対して関連性の高い前記脳波特徴量を選択する特徴量選択過程と
を含み、
前記複数の前記脳波特徴量のうち、少なくとも2以上の前記脳波特徴量は互いに異なる種類の単位である
ことを特徴とする脳波信号処理方法。 Among the EEG features, which are the features of the EEG signal generated in response to a given stimulus, the EEG that selects the EEG features that are highly relevant to each of the Stimulus features, which are the features of the stimulus. It is a signal processing method
A feature amount calculation process in which the feature amount calculation unit obtains a plurality of the brain wave feature amounts of the brain wave signal corresponding to the stimulus, and a feature amount calculation process.
The interrelationship calculation unit calculates the degree of interrelationship indicating the degree of association of the combination for each combination of one or more of the stimulus features of the stimulus and each of the plurality of electroencephalogram features. Interrelationship calculation process and
The feature amount selection unit extracts a combination of the stimulus feature amount and the brain wave feature amount whose degree of interrelationship exceeds a predetermined threshold value, and obtains the brain wave feature amount which is highly related to the stimulus feature amount. Including the feature selection process to be selected
Of the plurality of electroencephalogram features, at least two or more of the electroencephalogram features are units of different types from each other.
Electroencephalogram signal processing method characterized by and this.
コンピュータを、
前記刺激に対応した前記脳波信号の複数の前記脳波特徴量を求める特徴量算出手段、
前記刺激の1以上の前記刺激特徴量の各々と、複数の前記脳波特徴量それぞれとの組合わせ毎に、当該組合わせの関連の度合を示す相互関連度を算出する相互関連度算出手段、
前記相互関連度が所定の閾値を超える、前記刺激特徴量と前記脳波特徴量との組合わせを抽出し、前記刺激特徴量に対して関連性の高い前記脳波特徴量を選択する特徴量選択手段、
前記脳波信号を複数系列の脳波の特徴を有する成分信号に分離する信号分離手段
として動作させ、
前記信号分離手段が、前記脳波信号を複数系列の脳波の特徴を有する成分信号に分離する
プログラム。 It is a program for finding the relationship between the electroencephalogram feature amount, which is the feature amount of the electroencephalogram signal generated in response to a given stimulus, and the stimulus feature amount, which is the feature amount of the stimulus.
Computer,
A feature amount calculation means for obtaining a plurality of the electroencephalogram feature amounts of the electroencephalogram signal corresponding to the stimulus,
A means for calculating the degree of interrelationship, which calculates the degree of interrelationship indicating the degree of relevance of the combination for each combination of one or more of the stimulus features of the stimulus and each of the plurality of electroencephalogram features.
A feature amount selection means for extracting a combination of the stimulus feature amount and the brain wave feature amount whose degree of interrelationship exceeds a predetermined threshold value and selecting the brain wave feature amount having a high relevance to the stimulus feature amount. ,
It is operated as a signal separation means for separating the electroencephalogram signal into component signals having the characteristics of a plurality of series of electroencephalograms.
A program in which the signal separation means separates the electroencephalogram signal into component signals having the characteristics of a plurality of series of electroencephalograms.
コンピュータを、
前記刺激の各々に対応して、脳内信号源それぞれが発した前記脳波信号の周波数に基づく複数の前記脳波特徴量を求める特徴量算出手段、
前記刺激の1以上の前記刺激特徴量の各々と、複数の前記脳波特徴量それぞれとの組合わせ毎に、当該組合わせの関連の度合を示す相互関連度を算出する相互関連度算出手段、
前記相互関連度が所定の閾値を超える、前記刺激特徴量と前記脳波特徴量との組合わせを抽出し、前記刺激特徴量に対して関連性の高い前記脳波特徴量を選択する特徴量選択手段
として動作させるためのプログラム。 It is a program for finding the relationship between the electroencephalogram feature amount, which is the feature amount of the electroencephalogram signal generated in response to a given stimulus, and the stimulus feature amount, which is the feature amount of the stimulus.
Computer,
A feature amount calculation means for obtaining a plurality of the electroencephalogram feature amounts based on the frequencies of the electroencephalogram signals emitted by each of the intracerebral signal sources in response to each of the stimuli.
A means for calculating the degree of interrelationship, which calculates the degree of interrelationship indicating the degree of relevance of the combination for each combination of one or more of the stimulus features of the stimulus and each of the plurality of electroencephalogram features.
A feature amount selection means for extracting a combination of the stimulus feature amount and the brain wave feature amount whose degree of interrelationship exceeds a predetermined threshold value and selecting the brain wave feature amount having a high relevance to the stimulus feature amount. A program to operate as.
コンピュータを、
前記刺激に対応した前記脳波信号の複数の前記脳波特徴量を求める特徴量算出手段、
前記刺激の1以上の前記刺激特徴量の各々と、複数の前記脳波特徴量それぞれとの組合わせ毎に、当該組合わせの関連の度合を示す相互関連度を算出する相互関連度算出手段、
前記相互関連度が所定の閾値を超える、前記刺激特徴量と前記脳波特徴量との組合わせを抽出し、前記刺激特徴量に対して関連性の高い前記脳波特徴量を選択する特徴量選択手段
として動作させ、
前記複数の前記脳波特徴量のうち、少なくとも2以上の前記脳波特徴量は互いに異なる種類の単位である
プログラム。 To select the electroencephalogram feature that is highly relevant to each of the stimulus features that are the stimulus features, among the electroencephalogram features that are the features of the electroencephalogram signal generated in response to the given stimulus. Program
Computer,
A feature amount calculation means for obtaining a plurality of the electroencephalogram feature amounts of the electroencephalogram signal corresponding to the stimulus,
A means for calculating the degree of interrelationship, which calculates the degree of interrelationship indicating the degree of relevance of the combination for each combination of one or more of the stimulus features of the stimulus and each of the plurality of electroencephalogram features.
A feature amount selection means for extracting a combination of the stimulus feature amount and the brain wave feature amount whose degree of interrelationship exceeds a predetermined threshold value and selecting the brain wave feature amount having a high relevance to the stimulus feature amount. To operate as
A program in which at least two or more of the plurality of electroencephalogram features are units of different types.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2016147559A JP6859622B2 (en) | 2016-07-27 | 2016-07-27 | EEG signal processing system, EEG signal processing method and program |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2016147559A JP6859622B2 (en) | 2016-07-27 | 2016-07-27 | EEG signal processing system, EEG signal processing method and program |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2018015212A JP2018015212A (en) | 2018-02-01 |
| JP6859622B2 true JP6859622B2 (en) | 2021-04-14 |
Family
ID=61076891
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2016147559A Active JP6859622B2 (en) | 2016-07-27 | 2016-07-27 | EEG signal processing system, EEG signal processing method and program |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP6859622B2 (en) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113178195B (en) * | 2021-03-04 | 2022-08-26 | 杭州电子科技大学 | Speaker identification method based on sound-induced electroencephalogram signals |
| CN115517689B (en) * | 2022-09-06 | 2025-04-29 | 上海诺诚电气股份有限公司 | A neuroelectrophysiological monitoring device and method for spinal surgery |
| WO2024100861A1 (en) * | 2022-11-10 | 2024-05-16 | 日本電信電話株式会社 | Presentation device, presentation method, and presentation program |
| WO2025141737A1 (en) * | 2023-12-26 | 2025-07-03 | 日本電信電話株式会社 | Content evaluation device and content evaluation method |
Family Cites Families (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP3044891B2 (en) * | 1991-12-24 | 2000-05-22 | 株式会社デンソー | Vehicle air conditioner |
| JPH09108363A (en) * | 1995-10-18 | 1997-04-28 | Kokusai Electric Co Ltd | Biostimulation alleviation device |
| JP2007329100A (en) * | 2006-06-09 | 2007-12-20 | Toshiba Lighting & Technology Corp | Lighting device |
| US8391966B2 (en) * | 2009-03-16 | 2013-03-05 | Neurosky, Inc. | Sensory-evoked potential (SEP) classification/detection in the time domain |
| KR101130761B1 (en) * | 2010-02-26 | 2012-03-28 | 대구가톨릭대학교산학협력단 | Thematic Apperception Test device Based on BCI. |
| WO2013001836A1 (en) * | 2011-06-30 | 2013-01-03 | パナソニック株式会社 | Discomfort threshold level estimation system, method and program therefor, hearing aid adjustment system, and discomfort threshold level estimation processing circuit |
| US10238310B2 (en) * | 2013-12-16 | 2019-03-26 | Ideal Innovations Incorporated | Knowledge discovery based on brainwave response to external stimulation |
| JP5799351B1 (en) * | 2014-12-09 | 2015-10-21 | 株式会社センタン | Evaluation apparatus and evaluation method |
| US9914508B2 (en) * | 2015-10-16 | 2018-03-13 | Grady-White Boats, Inc. | Stowable folding boat step |
-
2016
- 2016-07-27 JP JP2016147559A patent/JP6859622B2/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| JP2018015212A (en) | 2018-02-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Li et al. | Sensitivity analysis of changes in human physiological indicators observed in soundscapes | |
| JP2593625B2 (en) | Biological information automatic identification device | |
| US8712514B2 (en) | Neurophysiological central auditory processing evaluation system and method | |
| Poling et al. | Characteristics of the 2f1-f2 distortion product otoacoustic emission in a normal hearing population | |
| US20140164095A1 (en) | Neuro-response data synchronization | |
| JP6859622B2 (en) | EEG signal processing system, EEG signal processing method and program | |
| JP2015229040A (en) | Emotion analysis system, emotion analysis method, and emotion analysis program | |
| JP6655242B2 (en) | Music listening experience estimation method, music listening experience estimation device, and music listening experience estimation program | |
| CN111671399A (en) | Method, apparatus and electronic equipment for measuring noise perception intensity | |
| CN104688222A (en) | EEG-based (electroencephalogram-based) tone synthesizer | |
| CN102469961A (en) | Speech sound intelligibility evaluation system, method, and program | |
| Alamir et al. | The effect of age, gender and noise sensitivity on the liking of food in the presence of background noise | |
| Regev et al. | Evidence for linear but not helical automatic representation of pitch in the human auditory system | |
| JP7623716B2 (en) | Information processing device and information processing method | |
| van den Broek et al. | Unobtrusive sensing of emotions (USE) | |
| Manohare et al. | Electroencephalography based classification of emotions associated with road traffic noise using Gradient boosting algorithm | |
| Laufs et al. | Methods and evaluation of physiological measurements with acoustic stimuli—a systematic review | |
| van der Reijden et al. | Frequency-specific objective audiometry: Tone-evoked brainstem responses and steady-state responses to 40 Hz and 90 Hz amplitude modulated stimuli: Audiometría objetiva de frecuencia específica: respuestas de tallo cerebral evocadas por tonos y respuestas de estado estable a 40 y 90 Hz con estímulos de amplitud modulada | |
| Schlittenlacher et al. | A hearing-model-based active-learning test for the determination of dead regions | |
| KR20020015547A (en) | recording medium for sensitivity treatment | |
| Dourou et al. | IoT-enabled analysis of subjective sound quality perception based on out-of-lab physiological measurements | |
| Vitale et al. | Physiological responses related to pleasant and unpleasant sounds | |
| Obeidat et al. | Objective methods to measure vestibular evoked myogenic potential response saccular tuning curves | |
| JP5603281B2 (en) | Hearing examination apparatus, method and computer program, hearing aid adjustment system | |
| Taymanov et al. | Measurement model as a means for studying the process of emotion origination |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A821 Effective date: 20160728 |
|
| RD03 | Notification of appointment of power of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7423 Effective date: 20181102 |
|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20190619 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20200310 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20200311 |
|
| A601 | Written request for extension of time |
Free format text: JAPANESE INTERMEDIATE CODE: A601 Effective date: 20200420 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20200629 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20201020 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20201211 |
|
| TRDD | Decision of grant or rejection written | ||
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20210224 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20210309 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 6859622 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |