(57)【特許請求の範囲】
【請求項１】
快／不快、活性／非活性、および期待感に関連する関心領域から、快／不快、活性／非活性、および期待感の各軸を含む感性多軸モデルの各軸に係る脳生理情報を抽出するステップと、
前記感性多軸モデルの各軸の脳生理情報を用いて感性を評価するステップと、を備え、
期待感に関連する関心領域が頭頂葉、後頭葉、および島皮質を含む領域である、感性評価方法。
【請求項２】
前記感性を評価するステップにおいて、前記感性多軸モデルの各軸の関連性を求め、この相関性と、前記感性多軸モデルの各軸の脳生理情報を用いて感性を評価する、請求項１に記載の感性評価方法。
【請求項３】
快／不快および活性／非活性に関連する関心領域が帯状回を含む領域である、請求項１または請求項２に記載の感性評価方法。
【請求項４】
前記脳生理情報が脳波信号に由来し、
前記脳生理情報を取得するステップは、前記感性多軸モデルの軸ごとに、
被験者の脳波信号を計測するステップと、
前記脳波信号について独立成分分析を行って複数の独立成分を抽出するステップと、
前記複数の独立成分のうち当該軸に関連する独立成分について時間周波数解析を行って時間周波数スペクトラムを算出するステップと、
前記脳生理情報として、前記時間周波数スペクトラムのうち当該軸に関連する関心周波数帯域でのスペクトラム強度から当該軸の脳生理指標値を推定するステップとを含む、請求項１ないし請求項３のいずれかに記載の感性評価方法。
【請求項５】
前記脳生理指標値を推定するステップにおいて、ｆＭＲＩにより計測されるＢＯＬＤ信号を用いて関心領域を特定する、請求項４に記載の感性評価方法。
【請求項６】
快／不快に関連する関心領域の周波数帯域がθ帯域である、請求項４または請求項５に記載の感性評価方法。
【請求項７】
活性／非活性に関連する関心領域の周波数帯域がβ帯域である、請求項４または請求項５に記載の感性評価方法。
【請求項８】
期待感に関連する関心領域の周波数帯域がθ〜α帯域である、請求項４または請求項５に記載の感性評価方法。
【請求項９】
前記脳生理情報がｆＭＲＩにより計測されるＢＯＬＤ信号に由来し、
前記脳生理情報を取得するステップは、前記感性多軸モデルの軸ごとに、
ｆＭＲＩを用いて被験者の全脳のＢＯＬＤ信号を取得するステップと、
前記全脳のＢＯＬＤ信号のうち当該軸に関連するＢＯＬＤ信号を選択するステップと、
前記脳生理情報として、前記選択したＢＯＬＤ信号から当該軸の脳生理指標値を推定するステップとを含む、請求項１ないし請求項３のいずれか記載の感性評価方法。
【請求項１０】
快／不快、活性／非活性、および期待感に関連する関心領域から、快／不快、活性／非活性、および期待感の各軸を含む感性多軸モデルの各軸に係る脳生理情報を抽出するステップと、
前記感性多軸モデルの各軸の脳生理情報から、脳波（EEG）の測定結果に基づく各軸の脳生理指標値（EEG快，EEG活性，EEG期待感）を求めるステップと、
被験者の主観的な統計データから得られ、前記感性多軸モデルの各軸の重み付け係数（a，b，c）を示す主観心理軸を用いて、次の計算式で、感性を評価するステップと、を備
えた感性評価方法。
感性＝[主観心理軸]＊[脳生理指標]＝a＊EEG快＋b＊EEG活性＋c＊EEG期待感
【請求項１１】
被験者の脳波信号を取得するステップと、
取得した脳波信号について独立成分分析を行って、得られた独立成分ごとにダイポールの位置を推定するステップと、
前記独立成分分析で得られた独立成分の主成分分析を行って、独立成分における脳活動データの次元縮約を行うステップと、
次元縮約された独立成分における脳活動データをクラスタリングするステップと、
前記クラスタリングで得られた複数のクラスタから、さまざまな感情あるいは情動のそれぞれを反映する脳活動を示すクラスタを選択するステップと、
選択したクラスタに含まれる成分から各感情あるいは情動の評価値を計算するステップと、
前記各感情あるいは情動について計算された評価値を総合して感性の評価値を計算するステップとを備えた、脳波を用いた感性評価方法。
【請求項１２】
前記感性が、各感情あるいは情動を各軸とする多軸モデルで表され、
前記感性の評価値を計算するステップにおいて、前記感性の評価値を前記各軸の評価値の重み付け加算により計算する、請求項１１に記載の脳波を用いた感性評価方法。
【請求項１３】
前記感性が３軸の多軸モデルで表される、請求項１２に記載の脳波を用いた感性評価方法。
【請求項１４】
前記感性の多軸モデルは、少なくとも快／不快の感情あるいは情動を表す軸を含んでおり、
前記快／不快の感情あるいは情動を反映する脳活動が、後帯状回周辺の脳活動である、請求項１２に記載の脳波を用いた感性評価方法。
【請求項１５】
ｆＭＲＩを用いて被験者の全脳のＢＯＬＤ信号を取得するステップと、
取得した全脳のＢＯＬＤ信号から、さまざまな感情あるいは情動のそれぞれを反映する脳活動を示すボクセルのＢＯＬＤ信号を選択するステップと、
選択したボクセルのＢＯＬＤ信号から各感情あるいは情動の評価値を計算するステップと、
前記各感情あるいは情動について計算された評価値を総合して感性の評価値を計算するステップとを備え、
前記感性の多軸モデルは、少なくとも快予期／不快予期の感情あるいは情動を表す軸を含んでおり、
前記快予期／不快予期の感情あるいは情動を反映する脳活動が、島皮質を含む領域の脳活動である、ｆＭＲＩを用いた感性評価方法。
【請求項１６】
前記感性が各感情あるいは情動を各軸とする多軸モデルで表され、
前記感性の評価値を計算するステップにおいて、前記感性の評価値を前記各軸の評価値の重み付け加算により計算する、請求項１５に記載の脳波を用いた感性評価方法。
【請求項１７】
前記感性が３軸の多軸モデルで表される、請求項１６に記載の脳波を用いた感性評価方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、感性を定量的に評価する方法に関する。
【背景技術】
【０００２】
人が機械やコンピュータなどのモノを操作する場合、手や足などの身体の一部を使ってハンドル、レバー、ボタン、キーボード、マウスなどの補助デバイスを操作したり、発話やジェスチャーによりモノに意思を伝達するのが一般的である。近年、脳と機械とを直接接続して人が思った通りに機械を操作するＢＭＩ（Brain Machine Interface）あるいはＢＣＩ（Brain Computer Interface）と呼ばれる技術が研究開発されている。ＢＭＩあるいはＢＣＩは、人の意思をモノに直接的に伝達できるようになることでモノの使い勝手の向上に期待されるほか、事故や病気によって運動機能や感覚機能などを失った人が自分の意思でモノを操作してモノを通じて他人と意思疎通を図ることができるようになる点で、医療や福祉の分野で期待されている。
【０００３】
上記のように人の意思あるいは能動的な意識といったものでモノを操作すること以外に、人の無意識あるいは潜在意識といったものを読み取ってそれをモノの操作に利用しようとする研究も行われている。例えば、被験者が分類対象情報を視認したときに生じる、当該被験者の主観的知見や感性に関連する脳波データに基づいて、分類対象情報を自動的に分類する技術が知られている（例えば、特許文献１を参照）。また、生体および脳機能計測から推定した操作者個々の内部情報を基に操作者の注意・記憶といった内部情報を推定もしくは予測し、ヒューマンエラーの起こりやすい状態を検出および警告する技術が知られている（例えば、特許文献２を参照）。また、潜在的な意識を示す高度な脳活動の情報を用いて、人が類似する感覚を覚える対象物を検出する技術が知られている（例えば、特許文献３参照）。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２００３−５８２９８号公報
【特許文献２】特開２０１１−１５０４０８号公報
【特許文献３】特開２０１４−１１５９１３号公報
【特許文献４】特開２００６−９５２６６号公報
【特許文献５】特開２０１１−１２０８２４号公報
【特許文献６】特開２００５−５８４４９号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
人の無意識あるいは潜在意識、特に感性といった人の精神活動あるいは心の情報を読み取ることができれば、人に心に優しいモノづくりやサービス提供が可能になる。例えば、対象物に対して人が抱く感性を客観的に検出し、または予測することができれば、そのような感性を発揮させるような対象物を事前にデザインすることができる。さらに、読み取った感性の情報は、人の心のケアや人と人とのコミュニケーションに活かすこともできる。本発明者らは、人の感性を読み取り、感性情報を介して人と人、ヒトとモノを繋ぐＢＥＩ（Brain Emotion Interface）の開発を目指している。
【０００６】
脳波（Electroencephalogram：EEG）などの脳生理情報はＢＥＩを実現するための有力な手段となり得る。そこで、本発明は、脳生理情報を用いて感性を定量的に評価することを課題とする。
【課題を解決するための手段】
【０００７】
本発明の一局面に従った感性評価方法は、快／不快、活性／非活性、および期待感に関連する関心領域から、快／不快、活性／非活性、および期待感の各軸を含む感性多軸モデルの各軸に係る脳生理情報を抽出するステップと、感性多軸モデルの各軸の脳生理情報を用いて感性を評価するステップと、を備えた感性評価方法である。
【０００８】
本発明の別の一局面に従った感性評価方法は、快／不快、活性／非活性、および期待感に関連する関心領域から、快／不快、活性／非活性、および期待感の各軸を含む感性多軸モデルの各軸に係る脳生理情報を抽出するステップと、感性多軸モデルの各軸の脳生理情報から、各軸の脳生理指標値（EEG快，EEG活性，EEG期待感）を求めるステップと、被験者の主観的な統計データから得られ、感性多軸モデルの各軸の重み付け係数（a，b，c）を示す主観心理軸を用いて、次の計算式で、感性を評価するステップと、を備えた感性評価方法である。
感性＝[主観心理軸]＊[脳生理指標]＝a＊EEG快＋b＊EEG活性＋c＊EEG期待感
【０００９】
本発明のさらに別の一局面に従った感性評価方法は、被験者の脳波信号を取得するステップと、取得した脳波信号について独立成分分析を行って、得られた独立成分ごとにダイポールの位置を推定するステップと、独立成分分析で得られた独立成分の主成分分析を行って、独立成分における脳活動データの次元縮約を行うステップと、次元縮約された独立成分における脳活動データをクラスタリングするステップと、クラスタリングで得られた複数のクラスタから、さまざまな感情あるいは情動のそれぞれを反映する脳活動を示すクラスタを選択するステップと、選択したクラスタに含まれる成分から各感情あるいは情動の評価値を計算するステップと、各感情あるいは情動について計算された評価値を総合して感性の評価値を計算するステップとを備えた、脳波を用いた感性評価方法である。
【００１０】
本発明のさらに別の一局面に従った感性評価方法は、ｆＭＲＩを用いて被験者の全脳のＢＯＬＤ信号を取得するステップと、取得した全脳のＢＯＬＤ信号から、さまざまな感情あるいは情動のそれぞれを反映する脳活動を示すボクセルのＢＯＬＤ信号を選択するステップと、選択したボクセルのＢＯＬＤ信号から各感情あるいは情動の評価値を計算するステップと、各感情あるいは情動について計算された評価値を総合して感性の評価値を計算するステップとを備えた、ｆＭＲＩを用いた感性評価方法である。
【発明の効果】
【００１１】
本発明によると、脳生理情報を用いて感性を定量的に評価することができる。これにより、感性情報を用いて、使えば使うほど人の感性に訴えかけて愛着が増すような製品をデザインしたり、感性情報を介して人と人とのコミュニケーションを円滑にすることができる。
【図面の簡単な説明】
【００１２】
【図１】情動、感情、感性の関係を表す模式図
【図２】本発明において採用する感性多軸モデルの模式図
【図３】感性多軸モデルの各軸に関連する関心領域を説明する図
【図４】快反応時のさまざまなｆＭＲＩ画像を示す図
【図５】快反応時のｆＭＲＩ画像およびＥＥＧ信号源をプロットした脳矢状断面を示す図
【図６】関心領域（快反応時の後帯状回）のＥＥＧ信号源の信号を時間周波数解析した結果を示す図
【図７】活性反応性時のｆＭＲＩ画像およびＥＥＧ信号源をプロットした脳矢状断面を示す図
【図８】関心領域（活性反応時の後帯状回）のＥＥＧ信号源の信号を時間周波数解析した結果を示す図
【図９】快／不快の刺激画像呈示実験の概要を説明する図
【図１０】快画像予期時および不快画像予期時の各ｆＭＲＩ画像を示す図
【図１１】快画像予期時と不快画像予期時のＥＥＧ信号の差分における信号源をプロットした脳矢状断面（頭頂葉部分）及び当該部分のＥＥＧ信号の時間周波数分布を示す図
【図１２】快画像予期時と不快画像予期時のＥＥＧ信号の差分における信号源をプロットした脳矢状断面（視覚野）及び当該部分のＥＥＧ信号の時間周波数分布を示す図
【図１３】主観心理軸決定のための自己評価の一例を示す図
【図１４】関心領域の脳波独立成分および周波数帯域を特定するフローチャート
【図１５】脳波信号の独立成分分析で抽出された各独立成分における信号強度分布を表したコンポーネント（脳波トポグラフィ）を示す図
【図１６】独立信号成分の信号源の推定位置をプロットした脳矢状断面図
【図１７】快・不快反応時のｆＭＲＩ画像を示す図
【図１８】ＥＥＧ信号源の信号を時間周波数解析した結果を示す図
【図１９】脳波を用いた感性のリアルタイム評価のフローチャート
【図２０】脳波信号の独立成分分析で抽出された各独立成分における信号強度分布を表したコンポーネント（脳波トポグラフィ）を示す図
【図２１】関心領域に関連する独立成分として特定されたコンポーネントを示す図
【図２２】特定された独立成分についての時間周波数解析の結果を示す図
【図２３】推定された快／不快軸の値を示す模式図
【図２４】推定された活性／非活性軸および期待感軸の各値を示す模式図
【図２５】各クラスタ（１６個）の脳波トポグラフィ
【図２６】各クラスタ（１６個）のダイポール推定位置をプロットした脳内画像
【図２７】快刺激と不快刺激との間に有意差のあったクラスタの脳波トポグラフィおよびダイポール推定位置をプロットした脳内画像
【図２８】（ａ）快画像および（ｂ）不快画像を見たときの主観評価と脳活動との相関を示す図
【図２９】本発明の一実施形態に係る脳波を用いた感性評価方法のフローチャート
【図３０】快画像を予期しているときの脳内画像
【図３１】不快画像を予期しているときの脳内画像
【図３２】本発明の一実施形態に係るｆＭＲＩを用いた感性評価方法のフローチャート
【発明を実施するための形態】
【００１３】
以下、適宜図面を参照しながら、実施の形態を詳細に説明する。ただし、必要以上に詳細な説明は省略する場合がある。例えば、既によく知られた事項の詳細説明や実質的に同一の構成に対する重複説明を省略する場合がある。これは、以下の説明が不必要に冗長になるのを避け、当業者の理解を容易にするためである。
【００１４】
なお、発明者らは、当業者が本発明を十分に理解するために添付図面および以下の説明を提供するのであって、これらによって特許請求の範囲に記載の主題を限定することを意図するものではない。
【００１５】
以下、１．本発明に至った背景、２．関心領域の特定、３．感性の可視化、４．脳波を用いた感性評価方法、５．ｆＭＲＩを用いた感性評価方法の順に説明する。なお、上記４は本願において優先権を主張する特願２０１５-２０４９６３に記載された内容であり、上記５は本願において優先権を主張する特願２０１５−２０４９６９に記載された内容である。
【００１６】
１．本発明に至った背景
人は何かを見たり、聞いたり、あるいは何かに触れたり、触れられたりしたときに、わくわくしたり、うきうきしたり、はらはらしたり、どきどきしたりする。これらは、単なる情動や感情と異なり、運動神経および感覚神経を含む体性神経系を通して脳に入ってくる外受容感覚、交感神経および副交感神経を含む自律神経系、それに基づく内受容感覚、さらには記憶や経験などが深く関与した複雑で高次の脳活動によってもたらされていると考えられる。
【００１７】
本発明では、わくわく、うきうき、はらはら、どきどきなどの感情あるいは情動とは異なる複雑な高次脳機能を広く「感性」として捉える。すなわち、本発明において、感性を、外受容感覚情報（体性神経系）と内受容感覚情報（自律神経系）を統合し、過去の経験、記憶と照らし合わせて生じる情動反応を、より上位のレベルで俯瞰する高次脳機能と定義する。換言すると、感性は、予測（イメージ）と結果（感覚情報）とのギャップを経験・知識と比較することによって直感的に“はっ”と気付く高次脳機能であると言える。
【００１８】
ここで、情動、感情、および感性の３つの概念を整理する。図１は情動、感情、感性の関係を表す模式図である。情動は外界からの刺激などによって引き起こされる無意識的・本能的な脳機能であり、３つの中で最も低次の脳機能である。感情は情動を意識化したより高次の脳機能である。そして、感性は経験・知識も反映したヒト特有の脳機能であり、３つの中で最も高次の脳機能である。
【００１９】
このような高次脳機能である感性の全体像を把握するには、種々の観点あるいは側面から総合的に感性を捉える必要がある。
【００２０】
例えば、人が快い、快適、あるいは心地よいと感じているか、あるいは反対に人が気持ち悪い、不快、あるいは心地よくないと感じているかといった「快／不快」の観点あるいは側面から感性を捉えることができる。
【００２１】
また、例えば、人が覚醒、興奮、あるいは活性状態にあるか、あるいは反対に人がぼんやり、沈静、あるいは非活性状態にあるかといった「活性／非活性」の観点あるいは側面から感性を捉えることができる。
【００２２】
また、例えば、人が何かを期待あるいは予期してわくわくしているか、あるいは期待が外れてがっかりしているかといった「期待感」の観点あるいは側面から感性を捉えることができる。
【００２３】
快／不快および活性／非活性を２軸に表したラッセル（Russell）の円環モデルが知られている。感情はこの円環モデルで表すことができる。しかし、感性は予測（イメージ）と結果（感覚情報）とのギャップを経験・知識と比較する高次脳機能であるので、快／不快および活性／非活性の２軸からなる既存の円環モデルでは十分に表し得ないと本発明者らは考える。そこで、本発明者らは、ラッセルの円環モデルに、時間軸（例えば、期待間）を第３軸として加えた感性多軸モデルを提唱する。
【００２４】
図２は、本発明において採用する感性多軸モデルの模式図である。感性多軸モデルは、例えば、「快／不快」を第１軸、「活性／非活性」を第２軸、「時間（期待感）」を第３軸として表すことができる。感性を多軸モデル化することのメリットは、各軸について評価値を算出し、それらを総合することで、漠然と広い概念の感性を定量的に評価する、すなわち、可視化することができる点にある。
【００２５】
この高次脳機能である感性を正確に評価することができれば、ヒトとモノを繋ぐＢＥＩ技術の確立に繋がる。そして、多様な分野で感性情報を活用して新価値を創造して、新しい価値を生み出すことができる。例えば、使えば使うほどヒトの思いに的確に反応し、喜び、やる気、愛情などの精神的価値が成長する製品・システムの創出を通してＢＥＩ技術の社会実装を実現すると考えられる。
【００２６】
従来、感性について言及した先行文献が多数存在するが、いずれの文献においても感性と感情を明確に区別せずに感性を感情と同意味で使用している。例えば、特許文献４には、感性には感情や意志が含まれるとし、喜怒哀楽などの感性状態を定量的に計測する内容が開示されている。しかし、特許文献４では感性と感情が区別されておらず、時間軸（期待感）の観点から感性を評価することは記載されていない。
【００２７】
特許文献５には、快、不快などに対応する刺激を与えた時の学習用生体情報を主成分分析し、複数の感性を快、不快などとして感性を定量的に評価する内容が開示されている（特に請求項５）。また、特許文献５において先行技術文献としてあげられている特許文献６には、感情データとして２軸モデルや３軸モデルなどの感情モデルを感情パラメータ値による感情を可視化する装置が開示されている。しかし、特許文献６では感性を定量的に評価できないと説明していることから明らかなように、感情と感性を混同しており、時間軸（期待感）についての記載も示唆もない。
【００２８】
なお、特許文献６には、快適と不快とのどちらにより近いかの度合いを示す第１軸と、興奮または緊張とリラックスとのどちらにより近いかの度合いを示す第２軸と、緊張と弛緩とのどちらにより近いかの度合いを示す第３軸との３軸により感情モデルを形成し、当該３軸空間における座標値からなる感情パラメータにより感情の状態を表現することが記載されているが、これはあくまで感情を表すためのモデルであり、ラッセルの円環モデルを複雑にしたに過ぎない。
【００２９】
以上のように、いずれの先行文献も感情分析の域を出ておらず、感性を正確に評価し得ないと考えられる。
【００３０】
２．関心領域の特定
快／不快、活性／非可性、および期待感の各脳反応に伴い、脳のどの部位が活動するかをｆＭＲＩとＥＥＧにより測定した結果について説明する。この測定結果は、感性を可視化、数値化する上での基礎データになり、極めて重要な位置づけにある。
【００３１】
ｆＭＲＩとは、ある心的過程と特定の脳構造を非侵襲的に対応づける脳機能画像法の一つであり、神経活動に伴う局所脳血流の酸素レベルに依存した信号強度を計測するものである。そのためｆＭＲＩはＢＯＬＤ（Blood Oxygen Level Dependent）法とも呼ばれる。
【００３２】
脳の中で神経細胞の活動が生じると多くの酸素が要求されるため、脳血流を通して酸素と結合した酸化ヘモグロビン（oxyhemoglobin）が局所において流れ込んでくる。そのときに神経細胞の酸素摂取を上回る酸素が供給され、結果として酸素を運び終えた還元型ヘモグロビン（deoxyhemoglobin）が局所において相対的に減少することになる。この還元型ヘモグロビンは磁気的性質を持ち、血管周囲の磁場の局所的不均一性を引き起こす。ｆＭＲＩは、このような酸素との結合関係に応じて磁気的性質を変化させるヘモグロビンの特徴を利用して、神経細胞の活動に伴う脳血流の酸素化バランスの局所的変化によって二次的に起こる信号増強を捉えるものである。現在では、局所的な脳血流の変化を全脳にわたり、数ミリメートル程度の空間的解像度で、秒単位で計測することが可能である。
【００３３】
図３は、感性多軸モデルの各軸に関連する関心領域を説明する図であり、各軸に関連する脳反応についてｆＭＲＩとＥＥＧにより測定した結果を示す。図３において、快・不快軸、活性・非活性軸のｆＭＲＩ画像、ＥＥＧ画像は、それぞれ、快反応時と不快反応時、活性反応時と非活性反応時との差分（変化分）を示すものである。また、期待感軸のｆＭＲＩ画像は快画面予期反応時のものであり、ＥＥＧ画像は快画像予期反応時と不快画像予期反応時との差分を示すものである。
【００３４】
図３に示したように、「快・不快」と「活性・非活性」反応時には帯状回が活動していることがｆＭＲＩおよびＥＥＧの測定結果から示され、「期待感」反応時にはｆＭＲＩおよびＥＥＧの測定結果から、頭頂葉、視覚野において脳活動があることが示される。
【００３５】
図３に示した感性多軸モデルの各軸に関連する関心領域は、ｆＭＲＩおよびＥＥＧを用いたさまざまな条件下での脳反応の観測実験を通じて得られた知見である。以下、その観測実験について具体的に説明する。
【００３６】
（１）快／不快時の脳反応について
まず、国際感情画像システム（International Affective Picture System：IAPS）から抽出した快画像（例えば、愛くるしいアザラシの赤ちゃん画像）と不快画像（例えば、危険な産業廃棄物画像）を２７名の実験参加者に提示することにより、実験参加者の快／不快時の脳反応を観察する。
【００３７】
図４は、快反応時のさまざまなｆＭＲＩ画像（脳の矢状断、冠状断、および水平断の各ｆＭＲＩ断面画像）を示す図である。図４において、不快反応時（不快画像を見た場合）と比較して快反応時（快画像を見た場合）に顕著に反応した領域に○印を付している。図４から明らかなように、快反応時には、後帯状回、視野野、線条体、眼窩前頭前野が賦活する。
【００３８】
図５は、快反応時のｆＭＲＩ画像およびＥＥＧ信号源をプロットした脳矢状断面を示す図である。図５において、不快反応時と比較して快反応時に顕著に反応した領域に○印を付している。図５からわかるように、快反応時には後帯状回を含む領域の脳活動がｆＭＲＩとＥＥＧの観測結果で共通している。この結果から、帯状回を含む領域を快／不快反応時の関心領域として特定することができる。
【００３９】
図６は、関心領域（快反応時の後帯状回）のＥＥＧ信号源の信号を時間周波数解析した結果を示す図である。図６左側は、関心領域（快反応時の後帯状回）のＥＥＧ信号源の信号を時間周波数解析した結果を示す図である。図６右側は快反応時と不快反応時の差分を示す。図６右側において色の濃い部分は差分が大きいことを表す。このＥＥＧの測定結果から、快反応時には関心領域のθ帯域の反応が関与していることがわかる。
【００４０】
（２）活性／非活性時の脳反応について
ＩＡＰＳから抽出した活性画像（例えば、美味しそうな寿司の画像）および非活性画像（例えば、静かな田園にたたずむ館の画像）を２７名の実験参加者に提示することにより、実験参加者の活性／非活性時の脳反応を観察する。
【００４１】
図７は、活性反応時のｆＭＲＩ画像およびＥＥＧ信号源をプロットした脳矢状断面を示す図である。図７において、非活性反応時（非活性画像を見た場合）と比較して活性反応時（活性画像を見た場合）に顕著に反応した領域に○印を付している。図７からわかるように、活性反応時には後帯状回を含む領域の脳活動がｆＭＲＩとＥＥＧの観測結果で共通している。この結果から、帯状回を含む領域を活性／非活性反応時の関心領域として特定することができる。
【００４２】
図８は、関心領域（活性反応時の後帯状回）のＥＥＧ信号源の信号を時間周波数解析した結果を示す図である。図８左側は、関心領域（活性反応時の後帯状回）のＥＥＧ信号源の信号を時間周波数解析した結果を示す図である。図８右側は活性反応時と非活性反応時の差分を示す。図８右側において色の濃い部分は差分が大きいことを表す。このＥＥＧの測定結果から、活性反応時には関心領域のβ帯域の反応が関与していることがわかる。
【００４３】
（３）期待時の脳反応について
まず、２７名の実験参加者に対して、情動を喚起する刺激画像を呈示し、画像を視認しているときの実験参加者の感情状態を評定させる実験を行う。刺激画像として、ＩＡＰＳから抽出した情動を喚起するカラー画像８０枚を用いる。そのうち４０枚が快さを喚起する画像（快画像）であり、残りの４０枚が不快を喚起する画像（不快画像）である。
【００４４】
図９は、快／不快の刺激画像呈示実験の概要を説明する図である。刺激画像は、短いトーン音（Ｃｕｅ）を０．２５秒間鳴らして、その３．７５秒後に４秒間だけ呈示する。そして、呈示された画像を快いと感じたか、不快と感じたかを被験者にボタンで回答してもらう。ただし、低いトーン音（５００Ｈｚ）が鳴った後には必ず快画像が呈示される。高いトーン音（４０００Ｈｚ）が鳴った後には必ず不快画像が呈示される。そして、中くらいのトーン音（１５００Ｈｚ）が鳴った後には、５０％の確率で快画像または不快画像が呈示される。
【００４５】
この実験において、いずれかのトーン音がなってから画像が呈示されるまでの４秒間は実験参加者が次に起こるであろうこと（この実験の場合には、快画像または不快画像が呈示されること）を予期している期間であり、この予期時における脳活動を観測した。例えば、低いトーン音が鳴ったとき、実験参加者は快画像が呈示されることを予期する「快画像予期」の状態にあり、高いトーン音が鳴ったとき、不快画像が呈示されることを予期する「不快画像予期」の状態にある。一方、中くらいのトーン音が鳴ったとき、実験参加者は快画像および不快画像のいずれが呈示されるのかがわからない「快・不快予期不可」の状態にある。
【００４６】
図１０は、快画像予期時および不快画像予期時の各ｆＭＲＩ画像（脳の矢状断および水平断の各ｆＭＲＩ断面画像）を示す図である。図１０の○印部分から明らかなように、ｆＭＲＩでは快画像予期時と不快画像予期時には、頭頂葉、視覚野、島皮質を含む脳領域が関与していることがわかる。
【００４７】
図１１は、ＥＥＧによる測定結果を示し、図１１ａは脳の矢状断の断面を示すものであり、不快画像予期時と比較して快画像予期時において顕著に反応した領域に破線○印を付している。また、図１１ｂは、関心領域（快画像予測時の頭頂葉の領域）のＥＥＧ信号源の信号を時間周波数解析した結果を示し、図１１ｃは、関心領域（不快画像予測時の頭頂葉の領域）のＥＥＧ信号源の信号を時間周波数解析した結果を示す。さらに、図１１ｄは、快予測時と不快予測時の差分を示した図であり、図中で丸で囲んだ部分は差分があった領域であり、その他の部分は差分がなかった領域である。このＥＥＧの測定結果から、快画像予測時において頭頂葉のβ帯域の反応が関与していることが理解される。
【００４８】
図１２は、ＥＥＧによる測定結果を示し、図１２ａは脳の矢状断の断面を示すものであり、不快画像予期時と比較して快画像予期時において顕著に反応した領域に破線○印を付している。また、図１２ｂは、関心領域（快画像予測時の視覚野の領域）のＥＥＧ信号源の信号を時間周波数解析した結果を示し、図１２ｃは、関心領域（不快画像予測時の視覚野の領域）のＥＥＧ信号源の信号を時間周波数解析した結果を示す。さらに、図１２ｄは、快予測時と不快予測時の差分を示した図であり、図中で丸で囲んだ部分は差分があった領域であり、その他の部分は差分がなかった領域である。このＥＥＧの測定結果から、快画像予測時視覚野のα帯域の反応が関与していることが理解される。
【００４９】
３．感性の可視化
感性を快／不快の軸、活性／非活性の軸、および期待感（時間）の軸の３軸を含む感性多軸モデルを用いて表すことについては上述したが、次は、具体的に感性をどのように可視化・数値化してＢＥＩ構築に結び付けできるかが課題になる。
【００５０】
本発明者らは、感性を構成する３軸は独立したものでなく相関性があるものであり、各軸の値を実測すると同時に各軸の感性に寄与する関係性を特定する必要があるとの知見に基づき、感性の主観心理軸と感性の脳生理指標を次のように融合させて感性の可視化を図っている。
【００５１】
感性＝［主観心理軸］＊［脳生理指標］＝a＊EEG快＋b＊EEG活性＋c＊EEG期待感…（式１）
ここで、主観心理軸は各軸の重み付け係数（a，b，c）を示し、脳生理指標はＥＥＧの測定結果に基づく各軸の値（EEG快，EEG活性，EEG期待感)を示す。
【００５２】
以下、主観心理軸の決定手順、および脳生理指標の選定手順について順に説明する。
【００５３】
Ａ．主観心理軸の決定
感性の主観心理軸を用いた各軸の寄与率、すなわち重み付けは次の手順で決定することができる。
【００５４】
（１）実験参加者（男女学生：２７名）に対して、上述の快／不快の刺激画像呈示実験を行う。ここでは、各トーン音が鳴ってから画像が呈示されるまでの４秒間（予期時）における能の感性状態を実験参加者の自己評価により評定してもらう。
【００５５】
（２）実験参加者には、３条件（快画像予期時、不快画像予期時、および快・不快予期不可）ごとに、わくわく（感性）度合、快度合（快軸）、活性度合（活性軸）、期待感度合（期待感軸）についてＶＡＳ（Visual Analog Scale）を用いて０から１００までの１０１段階で評定してもらう。図１３は、主観心理軸決定のための自己評価の一例を示す図であり、低いトーン音が鳴ったとき（快画像予期時）の快度合を評定している様子を示す。実験参加者は０から１００の間でカーソルを移動させて評定を行う。評定の結果、例えば、ある実験参加者から、快画像予期に関して、わくわく＝７３、快＝６８、活性＝４５、期待感＝７８といったような主観評定値が得られる。
【００５６】
（３）実験参加者全員から得られた３条件それぞれの主観評定値から、線形回帰により主観心理軸の各係数を算出する。この結果、例えば次式のような主観心理軸における感性評価式が得られる。
【００５７】
感性＝0.38＊主観快＋0.11＊主観活性＋0.51＊主観期待感…（式２）
ただし、主観快、主観活性、主観期待感は、実験参加者が評定した快度合、活性度合、期待感度合の各数値である。
【００５８】
（４）主観心理軸における主観快、主観活性、主観期待感と脳生理指標のEEG快、EEG活性、EEG期待感とはそれぞれ対応関係にある。したがって、主観評定値の線形回帰により算出された主観心理軸の各軸の重み係数は脳生理指標のEEG快、EEG活性、EEG期待感の各重み係数として用いることができる。そこで、式２で得られた各軸の重み係数を式１に適用することで、感性は、時々刻々測定されるEEG快，EEG活性，EEG期待感を用いて次式のように表される。
【００５９】
感性＝0.38＊EEG快＋0.11＊EEG活性＋0.51＊EEG期待感…（式３）
すなわち、式３により感性を数値により可視化することができる。
【００６０】
Ｂ．脳生理指標の選定
脳生理指標は、ＥＥＧの測定結果から計算される感性多軸モデルの各軸の推定値である。しかし、脳活動には個人差があるため、リアルタイムで感性を評価する前にあらかじめ各個人のＥＥＧを計測して各個人の脳波独立成分およびその周波数帯域を特定しておく必要がある。
【００６１】
まず、被験者の快／不快、活性／非活性、期待感の各脳波測定時に用いる周波数帯域を特定する手順について説明する。図１４は、関心領域の脳波独立成分および周波数帯域を特定するフローチャートである。
【００６２】
被験者に例えば快／不快を伴う画像を呈示して視覚刺激を与え、この刺激により誘発されたＥＥＧ脳波信号を計測する（Ｓ１）。なお、計測された脳波信号には、瞬き、目の動き、筋電に伴うノイズ（アーチファクト）が混在しているので、これらノイズを除去する。
【００６３】
計測した脳波信号に対して独立成分分析（ICA：Independent Component Analysis）を行って複数の独立成分（信号源）を抽出する（Ｓ２）。例えば、３２チャンネルで脳波を計測した場合は、チャンネルの数に応じた３２の独立成分が抽出される。計測した脳波の独立成分分析の結果、信号源の位置が特定される（Ｓ３）。
【００６４】
図１５は、ステップＳ２での脳波信号の独立成分分析で抽出された各独立成分における信号強度分布を表したコンポーネント（脳波トポグラフィ）を示す。また、図１６は、独立信号成分の信号源の推定位置をプロットした脳矢状断面を示す。
【００６５】
なお、脳波の測定とは別にｆＭＲＩによる測定を行う。図１７は、快・不快反応時のｆＭＲＩ画像を示す。快・不快の反応時において、図１７中で○印で示すように帯状回が関与していることがわかる。
【００６６】
このように別途行われるｆＭＲＩによる測定により、例えば、「快」の状態では帯状回が関与していることが判明しているため、「快」に関連する独立成分を選定する場合、関心領域の候補として帯状回付近に存在する信号源（独立成分）を選定することができる（Ｓ４）。例えば、３２の独立成分が取捨選択されて１０の独立成分に絞り込まれる。
【００６７】
関心領域の候補となる信号源の信号（例えば１０の独立成分）のそれぞれについて、時間周波数解析を行って、各時間ポイントおよび各周波数ポイントにおけるパワー値を算出する（Ｓ５）。例えば、４０の時間ポイントのそれぞれにおいて２０の周波数のポイントを設定して合計８００ポイントでのパワー値を算出する。
【００６８】
図１８は、ステップＳ５において、ＥＥＧ信号源の信号を時間周波数解析した結果を示す図である。図１８のグラフにおいて縦軸は周波数であり、横軸は時間である。周波数はβ、α、θの順で高い。グラフの色の濃淡は信号強度を表す。実際には時間周波数解析結果のグラフはカラーで表されるが、ここでは便宜上グレースケールで表している。
【００６９】
次に、時間周波数に分解された各独立成分に対して主成分分析（PCA：Principal Component Analysis）を行って、時間および周波帯域の主成分に絞り込みを行う（Ｓ６）。これにより、特徴の数が絞り込まれる。例えば、上記の８００ポイントの特徴から４０の主成分に次元が削減される。
【００７０】
各独立成分において、絞り込まれた時間周波数の主成分を対象に機械学習（SLR：Sparse Logistic Regression）を用いて判別学習を行う（Ｓ７）。これにより、その独立成分（信号源）における軸（例えば快／不快軸）の判別に寄与する主成分（時間周波数）が検出される。例えば、被験者の「快」測定時において、関心領域の信号源ではθ帯域が関係していることが判明する。また、例えば、快または不快の２択での判別精度が７０％であるといったように、その独立成分の周波帯域における判別精度が算出される。
【００７１】
算出された判別精度を元に、有意な判別率をもつ独立成分およびその周波帯域を特定する（Ｓ８）。これにより、関心領域の候補である例えば１０の独立成分の中からトップの独立成分およびその周波帯域が１つ選定される。
【００７２】
上記は快／不快の測定時における手順であるが、活性／非活性および期待感の各測定時においても同様の手順で関心領域の脳波独立成分および周波数帯域の特定を行う。この結果、活性／非活性の場合には関心領域のβ帯域が、期待感の場合には関心領域のθ〜α帯域がそれぞれ関与していることが判明する。
【００７３】
上記手順で得られた結果は、次の感性のリアルタイム評価において空間フィルタとして適用される。
【００７４】
なお、上記のステップＳ３およびＳ４ではすべての独立成分に対する信号源を推定した後にｆＭＲＩ情報に基づいて信号源（独立成分）の絞り込みを行っているが、ｆＭＲＩ情報を用いずにステップＳ３〜Ｓ７を実施し、最後のステップＳ８において有意に判別に寄与する独立成分の中から、ｆＭＲＩ情報を用いて独立成分（信号源）の選定を行い、その中で最も判別に寄与する独立成分を選択してもよい。このようにしても結果は同じになる。
【００７５】
次に、上記手順で特定された独立成分の周波数帯域を用いて、時々刻々変化する被験者の脳活動を推定して感性をリアルタイムで評価する手順について説明する。図１９は、脳波を用いた感性のリアルタイム評価のフローチャートである。
【００７６】
被験者の脳波を計測し、リアルタイムで脳波情報（各チャンネルでの脳活動）を抽出する（Ｓ１１）。なお、計測された各チャンネルの脳波信号には、瞬き、目の動き、筋電に伴うノイズ（アーチファクト）が混在しているので、これらノイズ成分を除去する。
【００７７】
計測した脳波信号に対して独立成分分析を行って複数の独立成分（信号源）を抽出する（Ｓ１２）。例えば、３２チャンネルで脳波を計測した場合は、チャンネルの数に応じた３２の独立成分が抽出される。図２０は、ステップＳ１２での脳波信号の独立成分分析で抽出された各独立成分における信号強度分布を表したコンポーネント（脳波トポグラフィ）を示す。
【００７８】
抽出された３２個の独立成分から、関心領域に関連する独立成分を特定する（Ｓ１３）。ここでは、図１４のフローチャートで示す手順により対象とする独立成分があらかじめ特定されているので、対象のコンポーネントは容易に特定される。図２１は、関心領域に関連する独立成分として特定されたコンポーネントを示す。
【００７９】
次に、特定された独立成分について時間周波数解析を行って、時間周波数スペクトラムを算出する（Ｓ１４）。図２２は、特定された独立成分についての時間周波数解析の結果を示す。
【００８０】
ここで、被験者の「快」測定時においてその独立成分（関心領域の信号源）では、対象となる周波帯域がθ帯域であることが判明しているため、当該帯域でのスペクトラムの信号強度より、ある時点での、快／不快軸の値（脳生理指標値）が推定される（Ｓ１５）。脳生理指標値は、例えば、０〜１００の数値で表される。図２３は、推定された快／不快軸の値を模式的に示す。例えば、図２３に示したように、快／不快軸の値としてEEG快＝６３が推定される。
【００８１】
上記は快／不快の測定時における手順であるが、活性／非活性および期待感の各測定時においても同様の手順で脳生理指標値を推定する。図２４は、推定された活性／非活性軸および期待感軸の各値を模式的に示す。例えば、図２４に示したように、活性／非活性軸の値としてEEG活性＝４２が、期待感軸（時間軸）の値としてEEG期待感＝７２が推定される。
【００８２】
推定した脳生理指標値を式３に代入して感性の評価値を計算する（Ｓ１６）。例えば、EEG快＝６３、EEG活性＝４２、EEG期待感＝７２という推定結果が得られた場合、感性の評価値は６５．２８と計算される。
【００８３】
４．脳波を用いた感性評価方法
次に、脳波を用いた別の感性評価方法について説明する。なお、以下の記載内容は、本願において優先権を主張する特願２０１５−２０４９６３に記載されたものである。
【００８４】
感性多軸モデルにおける各軸は、脳波に基づいて表される脳活動データによって評価することができる。例えば、脳波から得られる快／不快に関連する脳活動データをＥＥＧＶ１，ＥＥＧＶ２，・・・，ＥＥＧＶｉとし、各脳活動データの係数をｖ１，ｖ２，・・・，ｖｉとすると、図２の感性多軸モデルにおける第１軸の評価値Ｖａｌｅｎｃｅは、
Ｖａｌｅｎｃｅ＝ｖ１×ＥＥＧＶ１＋ｖ２×ＥＥＧＶ２＋・・・＋ｖｉ×ＥＥＧＶｉ
と表すことができる。
【００８５】
また、脳波から得られる活性／非活性に関連する脳活動データをＥＥＧＡ１，ＥＥＧＡ２，・・・，ＥＥＧＡｊとし、各脳活動データの係数をａ１，ａ２，・・・，ａｊとすると、図２の感性多軸モデルにおける第２軸の評価値Ａｒｏｕｓａｌは、
Ａｒｏｕｓａｌ＝ａ１×ＥＥＧＡ１＋ａ２×ＥＥＧＡ２＋・・・＋ａｊ×ＥＥＧＡｊ
と表すことができる。
【００８６】
また、脳波から得られる時間に関連する脳活動データをＥＥＧＴ１，ＥＥＧＴ２，・・・，ＥＥＧＴｋとし、各脳活動データの係数をｔ１，ｔ２，・・・，ｔｋとすると、図２の感性多軸モデルにおける第３軸の評価値Ｔｉｍｅは、
Ｔｉｍｅ＝ｔ１×ＥＥＧＴ１＋ｔ２×ＥＥＧＴ２＋・・・＋ｔｋ×ＥＥＧＴｋ
と表すことができる。
【００８７】
そして、感性の評価値Ｅｍｏｔｉｏｎは、各軸の係数をａ，ｂ，ｃとして、
Ｅｍｏｔｉｏｎ＝ａ×Ｖａｌｅｎｃｅ＋ｂ×Ａｒｏｕｓａｌ＋ｃ×Ｔｉｍｅ
と表すことができる。
【００８８】
なお、感性多軸モデルにおける各軸の評価値の計算に使用される係数ｖ，ａ，ｔ、および感性の評価値Ｅｍｏｔｉｏｎの計算に使用される係数ａ，ｂ，ｃはどのような値であってもよく、目的に応じて設定すればよい。例えば、評価値が大きくなるようにしたい場合（増幅）には、係数を１以上に設定すればよく、評価値が小さくなるようにしたい場合（減衰）には係数を０以上１未満に設定すればよい。特に、感性の評価値Ｅｍｏｔｉｏｎの計算に使用される係数ａ，ｂ，ｃは、感性多軸モデルにおける各軸の評価値を重み付け加算できるように、それぞれ０から１の範囲内かつ合計値が１になるように設定してもよい。これにより、感性多軸モデルにおける各軸の重要度あるいは寄与度に応じて各軸の評価値にバイアスをかけて感性を評価することができる。
【００８９】
また、図２では３軸で感性をモデル化しているが、さらに多くの軸を用いて感性をモデル化してもよい。また、上記の快／不快軸や活性／非活性軸などは軸の一例であって、別の軸で感性を多軸モデル化してもよい。
【００９０】
（１）実験例
快／不快と脳波との関連を調べるために実験を行った。実験の概要は図９を参照して説明した通りであるためここでの説明は省略する。
【００９１】
各被験者の頭部に６４個以上の電極を取り付けて、呈示された刺激画像を視認しているときの脳波信号を取得する。刺激画像呈示の試行回数は２００回以上とし、脳波信号の判定を図る。脳波信号のサンプリングレートは１０００Ｈｚ以上とする。
【００９２】
すべての脳波信号に対して１Ｈｚのハイパスフィルタをかける。
【００９３】
電極で測定される脳波信号には、脳（皮質）の活動に関連した電位変動のほかに、まばたきや筋電位などのアーチファクトや、商用電源ノイズなどの外部ノイズが含まれる。脳（皮質）の活動に関連した電位変動は非常に微弱であり、アーチファクトや外部ノイズに比べて小さいためＳＮ比が非常に低い。そのため、測定された脳波信号から純粋な脳反応を反映している可能性が高い信号を抽出するために、可能な限りノイズを除去する。
【００９４】
被験者ごとに、上記のノイズ除去後の脳波信号の独立成分分析を行って、得られた独立成分ごとにダイポールの位置推定を行う。
【００９５】
さらに、独立成分分析で得られたすべての独立成分について主成分分析を行って次元縮約を行う。そして、次元縮約後のすべての独立成分をｋ−ｍｅａｎｓ法により１６のクラスタにクラスタリングする。
【００９６】
図２５は、１６の各クラスタの脳波トポグラフィである。また、図２６は、１６の各クラスタのダイポール推定位置をプロットした脳内画像である。
【００９７】
上記１６のクラスタのうち、快刺激と不快刺激との間に有意差（ｐ＜０．０５）が見られ、かつ、脳活動データと被験者による主観評価値との間に高い相関を示すクラスタが一つあった。図２７は、快刺激と不快刺激との間に有意差のあったクラスタの脳波トポグラフィおよびダイポール推定位置をプロットした脳内画像である。図２８は、（ａ）快画像および（ｂ）不快画像を見たときの主観評価と脳活動との相関を示す。図２７のクラスタは、快画像を見たときおよび不快画像を見たときのいずれにおいても、主観評価と脳活動の高い関連が確認された。このように、今回の実験で、少なくとも後帯状回周辺が人の快／不快といった感情あるいは情動に深く関与していることが確認された。
【００９８】
なお、活性／非活性や時間についても脳内の同じまたは別の特定部分の活動が深く関与することが予想される。
【００９９】
（２）実施形態
今回の実験により得られた知見を基に、本発明に係る脳波を用いた感性評価方法は次のように実施することができる。図２９は、本発明の一実施形態に係る脳波を用いた感性評価方法のフローチャートである。なお、下記の処理フローは、ＰＣなどの汎用のコンピュータを用いて実施することができる。
【０１００】
被験者に何らかの刺激（対象物）を呈示した状態で被験者の脳波信号を取得する（Ｓ２１）。例えば、頭部に６４個の電極を取り付けて６４チャンネルの脳波信号を取得する。
【０１０１】
得られた脳波信号のフィルタリングおよびノイズ除去を行う（Ｓ２２）。フィルタリングは、例えば、１Ｈｚのハイパスフィルタである。ノイズ除去は、アーチファクトの除去、および商用電源ノイズなどの外乱ノイズの除去を含む。
【０１０２】
フィルタリングおよびノイズ除去後の脳波信号について、独立成分分析を行い、得られた独立成分ごとにダイポールの位置推定を行う（Ｓ２３）。すなわち、頭皮上で計測された脳波から脳内の主な活動源をダイポール（電流双極子）として推定する。例えば、６４チャンネルの脳波信号を、６４個の独立成分（ダイポール推定位置）に分解する。
【０１０３】
独立成分分析で得られたすべての独立成分について主成分分析を行って次元縮約を行う（Ｓ２４）。
【０１０４】
次元縮約後のすべての独立成分をｋ−ｍｅａｎｓ法により所定数（例えば、１６）のクラスタにクラスタリングする（Ｓ２５）。
【０１０５】
得られたクラスタのうち、快／不快を反映する脳活動を示すクラスタを選択し、その選択したクラスタに含まれる各成分（ダイポール推定位置を信号源とする脳活動データ）から、感性多軸モデルにおける快／不快軸の評価値Ｖａｌｅｎｃｅを計算する（Ｓ２６ａ）。
【０１０６】
また、得られたクラスタのうち、活性／非活性を反映する脳活動を示すクラスタを選択し、その選択したクラスタに含まれる各成分（ダイポール推定位置を信号源とする脳活動データ）から、感性多軸モデルにおける活性／非活性軸の評価値Ａｒｏｕｓａｌを計算する（Ｓ２６ｂ）。
【０１０７】
また、得られたクラスタのうち、時間を反映する脳活動を示すクラスタを選択し、その選択したクラスタに含まれる各成分（ダイポール推定位置を信号源とする脳活動データ）から、感性多軸モデルにおける時間軸の評価値Ｔｉｍｅを計算する（Ｓ２６ｃ）。
【０１０８】
ステップＳ２６ａ、Ｓ２６ｂ、Ｓ２６ｃでのクラスタ選択の際に感性データベース１００を参照することができる。感性データベース１００には、これまでの試行に係る多数の被験者の脳波信号およびそれら被験者による主観評価値が大量に蓄積されている。感性データベース１００に蓄積された大量のデータを用いてダイポール推定位置およびクラスタリングを行うことで、ステップＳ２６ａ、Ｓ２６ｂ、Ｓ２６ｃで選択すべきクラスタを特定することができる。
【０１０９】
なお、感性データベース１００を拡充するために、ステップＳ２１で取得した脳波信号および呈示刺激に対する被験者の主観評価値も新たに感性データベース１００に記録することが好ましい。
【０１１０】
ステップＳ２６ａ、Ｓ２６ｂ、Ｓ２６ｃで計算された各評価値を総合して、感性の評価値Ｅｍｏｔｉｏｎを計算する（Ｓ２７）。
【０１１１】
５．ｆＭＲＩを用いた感性評価方法
次に、ｆＭＲＩを用いた別の感性評価方法について説明する。なお、以下の記載内容は、本願において優先権を主張する特願２０１５−２０４９６９に記載されたものである。
【０１１２】
感性多軸モデルにおける各軸は、ｆＭＲＩに基づいて表される脳活動データによって評価することができる。例えば、ｆＭＲＩから得られる快／不快に関連する脳活動データをｆＭＲＩＶ１，ｆＭＲＩＶ２，・・・，ｆＭＲＩＶｉとし、各脳活動データの係数をｖ１，ｖ２，・・・，ｖｉとすると、図２の感性多軸モデルにおける第１軸の評価値Ｖａｌｅｎｃｅは、
Ｖａｌｅｎｃｅ＝ｖ１×ｆＭＲＩＶ１＋ｖ２×ｆＭＲＩＶ２＋・・・＋ｖｉ×ｆＭＲＩＶｉ
と表すことができる。
【０１１３】
また、ｆＭＲＩから得られる活性／非活性に関連する脳活動データをｆＭＲＩＡ１，ｆＭＲＩＡ２，・・・，ｆＭＲＩＡｊとし、各脳活動データの係数をａ１，ａ２，・・・，ａｊとすると、図２の感性多軸モデルにおける第２軸の評価値Ａｒｏｕｓａｌは、
Ａｒｏｕｓａｌ＝ａ１×ｆＭＲＩＡ１＋ａ２×ｆＭＲＩＡ２＋・・・＋ａｊ×ｆＭＲＩＡｊ
と表すことができる。
【０１１４】
また、ｆＭＲＩから得られる時間に関連する脳活動データをｆＭＲＩＴ１，ｆＭＲＩＴ２，・・・，ｆＭＲＩＴｋとし、各脳活動データの係数をｔ１，ｔ２，・・・，ｔｋとすると、図２の感性多軸モデルにおける第３軸の評価値Ｔｉｍｅは、
Ｔｉｍｅ＝ｔ１×ｆＭＲＩＴ１＋ｔ２×ｆＭＲＩＴ２＋・・・＋ｔｋ×ｆＭＲＩＴｋ
と表すことができる。
【０１１５】
そして、感性の評価値Ｅｍｏｔｉｏｎは、各軸の係数をａ，ｂ，ｃとして、
Ｅｍｏｔｉｏｎ＝ａ×Ｖａｌｅｎｃｅ＋ｂ×Ａｒｏｕｓａｌ＋ｃ×Ｔｉｍｅ
と表すことができる。
【０１１６】
なお、感性多軸モデルにおける各軸の評価値の計算に使用される係数ｖ，ａ，ｔ、および感性の評価値Ｅｍｏｔｉｏｎの計算に使用される係数ａ，ｂ，ｃはどのような値であってもよく、目的に応じて設定すればよい。例えば、評価値が大きくなるようにしたい場合（増幅）には、係数を１以上に設定すればよく、評価値が小さくなるようにしたい場合（減衰）には係数を０以上１未満に設定すればよい。特に、感性の評価値Ｅｍｏｔｉｏｎの計算に使用される係数ａ，ｂ，ｃは、感性多軸モデルにおける各軸の評価値を重み付け加算できるように、それぞれ０から１の範囲内かつ合計値が１になるように設定してもよい。これにより、感性多軸モデルにおける各軸の重要度あるいは寄与度に応じて各軸の評価値にバイアスをかけて感性を評価することができる。
【０１１７】
また、図２では３軸で感性をモデル化しているが、さらに多くの軸を用いて感性をモデル化してもよい。また、上記の快／不快軸や活性／非活性軸などは軸の一例であって、別の軸で感性を多軸モデル化してもよい。
【０１１８】
（１）実験例
先行研究では、情動を喚起する画像を用いた実験から脳内の島や扁桃体が人の不快予期と性格特性のそれぞれに関与することが報告されている（Simmons et al., 2006；Schuerbeek et al., 2014）。しかし、先行研究では、用いた画像の種類が少なかったり、快画像あるいは不快画像との比較がなかったりするなど、いくつかの制限がある。実験結果を汎化するためにはより多様な状況に即した知見を得ることが重要である。そこで、より多様な刺激画像を用いつつ、さらに各個人の性格特性と、快／不快を予期しているときの脳の特定部分の活動状況（賦活の程度）との関連を調べるために実験を行った。なお、実験の概要は図９を参照して説明した通りであるためここでの説明は省略する。
【０１１９】
被験者の性格特性として「損害回避」の得点を用いた。事前に各被験者に質問紙による性格特性調査の一つである気質性格検査（Temperament and Character Inventory：TCI）を実施して「損害回避」の得点を計算した。損害回避の得点が高いほど、将来起こり得る問題に対する悲観的な心配や、好ましくない事態の予期に対する受動的回避行動（例えば、ある嫌悪刺激に対してそれに近づかないあるいは行動しないことにより回避しようとする行為）などの傾向が高いことを意味する。
【０１２０】
被験者のｆＭＲＩの測定には３ＴＭＲＩ装置を用いた。刺激画像呈示の試行回数は１２０回とし、脳活動由来の信号の安定を図る。
【０１２１】
図３０は、快画像を予期しているとき、すなわち、低いトーン音が鳴った後の脳内画像である。図３１は、不快画像を予期しているとき、すなわち、高いトーンが鳴った後の脳内画像である。快予期時および不快予期時のいずれにおいても、島皮質を含む脳領域の賦活が高まることが確認された。
【０１２２】
さらに、島皮質を含む領域について快予期時と不快予期時の賦活を比較すると、損害回避の得点が高い人ほど不快予期時の賦活が相対的に高くなることが確認された。すなわち、島皮質を含む領域に損害回避の得点と正の相関がある領域があることが確認された。
【０１２３】
なお、快／不快や活性／非活性についても脳内の同じまたは別の特定部分の活動が深く関与することが予想される。
【０１２４】
（２）実施形態
今回の実験により得られた知見を基に、本発明に係るｆＭＲＩを用いた感性評価方法は次のように実施することができる。図３２は、本発明の一実施形態に係るｆＭＲＩを用いた感性評価方法のフローチャートである。なお、下記の処理フローは、ＰＣなどの汎用のコンピュータを用いて実施することができる。
【０１２５】
被験者に何らかの刺激（対象物）を呈示した状態でｆＭＲＩを用いて被験者の全脳のＢＯＬＤ信号を取得する（Ｓ３１）。
【０１２６】
得られた全脳のＢＯＬＤ信号に対して、空間・時間的前処理を行う（Ｓ３２）。
【０１２７】
前処理後の全脳のＢＯＬＤ信号のうち、快／不快を反映する脳活動を示すボクセルのＢＯＬＤ信号、活性／非活性を反映する脳活動を示すボクセルのＢＯＬＤ信号、および時間を反映する脳活動を示すボクセルのＢＯＬＤ信号を抽出する（Ｓ３３）。
【０１２８】
快／不快を反映する脳活動を示すボクセルのＢＯＬＤ信号から、感性多軸モデルにおける快／不快軸の評価値Ｖａｌｅｎｃｅを計算する（Ｓ３４ａ）。
【０１２９】
また、活性／非活性を反映する脳活動を示すボクセルのＢＯＬＤ信号から、感性多軸モデルにおける活性／非活性軸の評価値Ａｒｏｕｓａｌを計算する（Ｓ３４ｂ）。
【０１３０】
また、時間を反映する脳活動を示すボクセルのＢＯＬＤ信号から、感性多軸モデルにおける時間軸の評価値Ｔｉｍｅを計算する（Ｓ３４ｃ）。
【０１３１】
ステップＳ３３でのボクセル選択の際に感性データベース１００を参照することができる。感性データベース１００には、これまでの試行に係る多数の被験者のｆＭＲＩデータおよびそれら被験者による主観評価値が大量に蓄積されている。感性データベース１００に蓄積された大量のデータを用いてステップＳ３３で選択すべきクラスタを特定することができる。
【０１３２】
なお、感性データベース１００を拡充するために、ステップＳ３１で取得した全脳のＢＯＬＤ信号および呈示刺激に対する被験者の主観評価値も新たに感性データベース１００に記録することが好ましい。
【０１３３】
ステップＳ３４ａ、Ｓ３４ｂ、Ｓ３４ｃで計算された各評価値を総合して、感性の評価値Ｅｍｏｔｉｏｎを計算する（Ｓ３５）。
【０１３４】
ステップＳ３４ａ、Ｓ３４ｂ、Ｓ３４ｃ、Ｓ３５における各評価値の計算において、被験者の性格特性を反映することが好ましい。例えば、損害回避傾向の高い人は、ネガティブな事象に対して脳反応が大きいことがわかっている。すなわち、損害回避傾向の高い人は、損害回避傾向の低い人と比べて、感性の評価値Ｅｍｏｔｉｏｎの絶対値が全体的に高くなる可能性がある。そのため、損害回避傾向の低い人と比較して、損害回避傾向の高い人の評価値Ｅｍｏｔｉｏｎをどのように解釈するかが難しくなる。例えば、損害回避傾向の低い人にとって少し嫌だと感じる（例えば評価値Ｅｍｏｔｉｏｎが−２０程度）イベントも非常に嫌だと感じる（例えば評価値Ｅｍｏｔｉｏｎが−１００以下）イベントも、損害回避傾向の高い人は極端な脳反応を示す（例えばいずれのイベントにおいても評価値Ｅｍｏｔｉｏｎが−１００以下となる）ことでその差がわかりにくくなってしまう。そこで、例えばステップＳ３４ａでの快／不快の評価値Ｖａｌｅｎｃｅの計算において、損害回避傾向の低い人の場合には、係数ｖ１，ｖ２，・・・，ｖｉを例えば１以下に設定して評価値Ｖａｌｅｎｃｅが小さな値になるような減衰処理を施すとよい。このような減衰処理などを加えることで、損害回避傾向の高い人の感性の評価値Ｅｍｏｔｉｏｎが、損害回避傾向の低い人の感性の評価値Ｅｍｏｔｉｏｎと比較して解釈しやすくなると考えられる。
【０１３５】
以上のように、本発明における技術の例示として、実施の形態を説明した。そのために、添付図面および詳細な説明を提供した。
【０１３６】
したがって、添付図面および詳細な説明に記載された構成要素の中には、課題解決のために必須な構成要素だけでなく、上記技術を例示するために、課題解決のためには必須でない構成要素も含まれ得る。そのため、それらの必須ではない構成要素が添付図面や詳細な説明に記載されていることをもって、直ちに、それらの必須ではない構成要素が必須であるとの認定をするべきではない。
【０１３７】
また、上述の実施の形態は、本発明における技術を例示するためのものであるから、特許請求の範囲またはその均等の範囲において種々の変更、置き換え、付加、省略などを行うことができる。
【産業上の利用可能性】
【０１３８】
本発明に係る感性評価方法は、脳波などの脳生理情報を用いて感性を定量的に評価することができるため、ヒトとモノを繋ぐＢＥＩを実現するための基礎技術として有用である。
【図１】
0006590411.tif 000002
【図２】
0006590411.tif 000003
【図３】
0006590411.tif 000004
【図４】
0006590411.tif 000005
【図５】
0006590411.tif 000006
【図６】
0006590411.tif 000007
【図７】
0006590411.tif 000008
【図８】
0006590411.tif 000009
【図９】
0006590411.tif 000010
【図１０】
0006590411.tif 000011
【図１１】
0006590411.tif 000012
【図１２】
0006590411.tif 000013
【図１３】
0006590411.tif 000014
【図１４】
0006590411.tif 000015
【図１５】
0006590411.tif 000016
【図１６】
0006590411.tif 000017
【図１７】
0006590411.tif 000018
【図１８】
0006590411.tif 000019
【図１９】
0006590411.tif 000020
【図２０】
0006590411.tif 000021
【図２１】
0006590411.tif 000022
【図２２】
0006590411.tif 000023
【図２３】
0006590411.tif 000024
【図２４】
0006590411.tif 000025
【図２５】
0006590411.tif 000026
【図２６】
0006590411.tif 000027
【図２７】
0006590411.tif 000028
【図２８】
0006590411.tif 000029
【図２９】
0006590411.tif 000030
【図３０】
0006590411.tif 000031
【図３１】
0006590411.tif 000032
【図３２】
0006590411.tif 000033