JP6093422B1 - Brain age presentation device - Google Patents

Brain age presentation device Download PDF

Info

Publication number
JP6093422B1
JP6093422B1 JP2015203358A JP2015203358A JP6093422B1 JP 6093422 B1 JP6093422 B1 JP 6093422B1 JP 2015203358 A JP2015203358 A JP 2015203358A JP 2015203358 A JP2015203358 A JP 2015203358A JP 6093422 B1 JP6093422 B1 JP 6093422B1
Authority
JP
Japan
Prior art keywords
brain
data
skin temperature
component
facial skin
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
Application number
JP2015203358A
Other languages
Japanese (ja)
Other versions
JP2017074215A (en
Inventor
潤一郎 新井
潤一郎 新井
後藤 尚志
尚志 後藤
誠 岩亀
誠 岩亀
健一 日野
健一 日野
平野 智也
智也 平野
泰則 小谷
泰則 小谷
淑美 大上
淑美 大上
太郎 戸松
太郎 戸松
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Daikin Industries Ltd
Tokyo Institute of Technology NUC
Original Assignee
Daikin Industries Ltd
Tokyo Institute of Technology NUC
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Daikin Industries Ltd, Tokyo Institute of Technology NUC filed Critical Daikin Industries Ltd
Priority to JP2015203358A priority Critical patent/JP6093422B1/en
Priority to PCT/JP2016/080711 priority patent/WO2017065313A1/en
Priority to US15/767,067 priority patent/US10709338B2/en
Priority to EP16855564.7A priority patent/EP3363352B1/en
Priority to CN201680060244.6A priority patent/CN108135498B/en
Application granted granted Critical
Publication of JP6093422B1 publication Critical patent/JP6093422B1/en
Publication of JP2017074215A publication Critical patent/JP2017074215A/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Abstract

【課題】簡便に脳年齢を提示し得る脳年齢提示装置を提供する。【解決手段】脳年齢提示装置301は、顔面皮膚温度取得手段20及び/又は画像データ取得手段120と、脳活動推定手段30,130と、状態可視化手段200とを備える。顔面皮膚温度取得手段20は、顔面皮膚温度データを時系列で取得する。画像データ取得手段120は、顔面の撮像画像データを時系列で取得する。脳活動推定手段30,130は、顔面皮膚温度データ及び/又は顔面血行量データを、特異値分解、主成分分析或いは独立成分分析により分解することで得られた複数の成分に基づき、対象者の脳活動を推定する。状態可視化手段200は、脳活動推定手段30により推定された対象者の脳活動に基づき、対象者の脳年齢を提示する。【選択図】図23A brain age presentation device capable of simply presenting brain age is provided. A brain age presentation device 301 includes facial skin temperature acquisition means 20 and / or image data acquisition means 120, brain activity estimation means 30, 130, and state visualization means 200. The facial skin temperature acquisition means 20 acquires facial skin temperature data in time series. The image data acquisition unit 120 acquires captured image data of the face in time series. Based on a plurality of components obtained by decomposing facial skin temperature data and / or facial blood flow data by singular value decomposition, principal component analysis, or independent component analysis, the brain activity estimating means 30, 130 Estimate brain activity. The state visualization unit 200 presents the brain age of the subject based on the brain activity of the subject estimated by the brain activity estimation unit 30. [Selection] Figure 23

Description

本発明は、脳年齢提示装置に関する。   The present invention relates to a brain age presentation device.

近年、脳年齢を把握したいというニーズが高まっている。脳年齢を把握するために、特許文献1(特開2013−176406号公報)に開示されているような脳波計測法(EEG)、磁気共鳴画像法(fMRI: functional Magnetic Resonance Imaging)、又は近赤外線分光法(NIRS)などが用いられることがある。   In recent years, there has been an increasing need to ascertain brain age. In order to grasp the brain age, an electroencephalogram measurement method (EEG), magnetic resonance imaging (fMRI), or near infrared ray as disclosed in Patent Document 1 (Japanese Patent Laid-Open No. 2013-176406) Spectroscopy (NIRS) may be used.

しかしながら、上述した各方法による脳波計測に基づいて脳年齢を推定する場合、対象者に対して前処理の必要な電極を装着する必要がある等、準備段階での作業が繁雑であったり、所定の計測室でなければ計測できなかったり、多大なコストを要したりするといった問題がある。   However, when estimating the brain age based on the electroencephalogram measurement by each method described above, it is necessary to wear an electrode that requires pretreatment for the subject, and the work in the preparation stage is complicated, There is a problem that it cannot be measured unless it is a measurement room, and it requires a great deal of cost.

本発明の課題は、簡便に脳年齢を提示し得る脳年齢提示装置を提供することにある。   The subject of this invention is providing the brain age presentation apparatus which can present brain age simply.

本発明の第1観点に係る脳年齢提示装置は、顔面皮膚温度取得手段及び/又は顔面血行量取得手段と、脳活動推定手段と、脳年齢提示手段と、を備える。顔面皮膚温度取得手段は、顔面皮膚温度データを時系列で取得する。顔面皮膚温度データは、対象者の顔面から検出した皮膚温度データ及びその検出部位の位置データを含む。顔面血行量取得手段は、顔面血行量データを時系列で取得する。顔面血行量データは、対象者の顔面を時系列で撮影した撮影画像データに対してRGB処理を行い、RGB処理により得られた撮影画像データのRGBデータに基づき取得される。脳活動推定手段は、顔面皮膚温度データ及び/又は顔面血行量データを、特異値分解、主成分分析或いは独立成分分析により分解することで得られた複数の成分に基づき、対象者の脳活動を推定する。脳年齢提示手段は、脳活動推定手段により推定された対象者の脳活動に基づき、対象者の脳年齢を提示する。   The brain age presentation device according to the first aspect of the present invention comprises facial skin temperature acquisition means and / or facial blood flow acquisition means, brain activity estimation means, and brain age presentation means. The facial skin temperature acquisition means acquires facial skin temperature data in time series. The facial skin temperature data includes skin temperature data detected from the subject's face and position data of the detected part. The facial blood flow acquisition means acquires facial blood flow data in time series. The facial blood flow data is obtained based on the RGB data of the captured image data obtained by the RGB processing by performing RGB processing on the captured image data obtained by capturing the face of the subject in time series. The brain activity estimating means determines the brain activity of the subject based on a plurality of components obtained by decomposing facial skin temperature data and / or facial blood flow data by singular value decomposition, principal component analysis, or independent component analysis. presume. The brain age presentation means presents the brain age of the subject based on the brain activity of the subject estimated by the brain activity estimation means.

本発明の第1観点の係る脳年齢提示装置では、顔面皮膚温度取得手段及び/又は顔面血行量取得手段によって取得された時系列の顔面皮膚温度データ及び/又は顔面血行量データに基づき対象者の脳活動が推定される。このため、この脳年齢提示装置では、脳波電極等の装着前に処理が必要なセンサを装着しなくても、対象者の脳活動を推定することができる。結果として、簡便に脳年齢を提示することができる。   In the brain age presentation device according to the first aspect of the present invention, based on the time-series facial skin temperature data and / or facial blood flow data acquired by the facial skin temperature acquisition means and / or facial blood flow volume acquisition means, Brain activity is estimated. For this reason, in this brain age presentation device, it is possible to estimate the brain activity of the subject even without wearing a sensor that requires processing before wearing an electroencephalogram electrode or the like. As a result, it is possible to simply present the brain age.

本発明の第2観点に係る脳年齢提示装置は、第1観点の脳年齢提示装置において、脳活動推定手段は、複数の成分のうち、その成分波形の振幅が脳の賦活時及び脳の非賦活時の変化と相関関係にある成分を判定用成分として抽出する。そして、脳活動推定手段は、判定用成分に基づき対象者の脳活動を推定する。この脳年齢提示装置では、複数の成分のうちの脳の賦活/非賦活と相関関係にある成分が対象者の脳活動を推定する判定用成分として抽出される。このため、対象者の脳活動と関連性が高いと推測される成分から脳活動を推定することができる。結果として、脳年齢の推定精度を高めることができる。   The brain age presentation device according to a second aspect of the present invention is the brain age presentation device according to the first aspect, wherein the brain activity estimation means includes a plurality of components whose amplitude of the component waveform is not activated when the brain is activated. A component having a correlation with a change at the time of activation is extracted as a determination component. Then, the brain activity estimating means estimates the brain activity of the subject based on the determination component. In this brain age presentation device, a component having a correlation with brain activation / deactivation among a plurality of components is extracted as a determination component for estimating the brain activity of the subject. For this reason, it is possible to estimate the brain activity from a component that is presumed to be highly related to the brain activity of the subject. As a result, it is possible to improve the estimation accuracy of brain age.

本発明の第3観点に係る脳年齢提示装置は、第1観点又は第2観点の脳年齢提示装置において、顔面皮膚温度取得手段及び/又は顔面血行量取得手段は、対象者の顔面の副鼻腔周辺及び/又は前額部の顔面皮膚温度データ及び/又は顔面血行量データを少なくとも取得する。   The brain age presentation apparatus according to a third aspect of the present invention is the brain age presentation apparatus according to the first aspect or the second aspect, wherein the facial skin temperature acquisition means and / or the facial blood flow acquisition means is the sinus of the face of the subject. At least facial skin temperature data and / or facial blood flow data of the peripheral and / or forehead are obtained.

ここで、脳には、選択的脳冷却機構(Selective Brain Cooling System)という体温とは独立して脳を冷却する仕組みがある。選択的脳冷却機構としては、脳活動によって生じた熱を、副鼻腔及び前額部周辺を用いて排熱していることが知られている。また、顔面の皮膚温度は、顔面の血行量に比例すると考えられる。   Here, the brain has a mechanism called a selective brain cooling system that cools the brain independently of the body temperature. As a selective brain cooling mechanism, it is known that heat generated by brain activity is exhausted using the sinuses and the periphery of the forehead. Further, the skin temperature of the face is considered to be proportional to the blood flow volume of the face.

本発明の第3観点に係る脳年齢提示装置では、脳活動を反映していると推測される副鼻腔及び/又は前額部における時系列の顔面皮膚温度データ及び/又は顔面血行量データに基づき、対象者の脳活動が推定される。このため、この脳年齢提示装置では、対象者の脳活動を精度よく推定することができる。結果として、脳年齢の推定精度を高めることができる。   In the brain age presentation device according to the third aspect of the present invention, based on time series facial skin temperature data and / or facial blood flow data in the paranasal sinuses and / or the forehead portion, which is presumed to reflect brain activity. The brain activity of the subject is estimated. For this reason, in this brain age presentation apparatus, a subject's brain activity can be estimated accurately. As a result, it is possible to improve the estimation accuracy of brain age.

本発明の第4観点に係る脳年齢提示装置は、第1観点から第3観点のいずれかの脳年齢提示装置であって、脳年齢記憶部をさらに備える。脳年齢記憶部は、脳活動の推定結果と脳年齢を示す脳年齢データとを関連付けて記憶する。また、脳活動推定手段は、脳活動の推定結果を算出する。また、脳年齢提示手段は、推定結果に基づいて脳年齢記憶部から対象者の脳年齢データを抽出して提示する。   A brain age presentation device according to a fourth aspect of the present invention is the brain age presentation device according to any one of the first to third aspects, and further includes a brain age storage unit. The brain age storage unit stores an estimation result of brain activity and brain age data indicating the brain age in association with each other. The brain activity estimating means calculates a brain activity estimation result. The brain age presenting means extracts and presents the brain age data of the subject from the brain age storage unit based on the estimation result.

本発明の第4観点に係る脳年齢提示装置では、脳年齢記憶部が、脳活動の推定結果と脳年齢を示す脳年齢データとを関連付けて記憶するので、脳活動と脳年齢との相関関係を逐次最適化することができる。   In the brain age presentation device according to the fourth aspect of the present invention, the brain age storage unit stores the brain activity estimation result and the brain age data indicating the brain age in association with each other, so that the correlation between the brain activity and the brain age is stored. Can be optimized sequentially.

本発明の第5観点に係る脳年齢提示装置は、第1観点から第4観点のいずれかの脳年齢提示装置であって、顔面皮膚温度取得手段及び/又は顔面血行量取得手段と、脳年齢提示手段とがスマートデバイスに格納されているものである。   A brain age presentation device according to a fifth aspect of the present invention is the brain age presentation device according to any one of the first to fourth aspects, comprising facial skin temperature acquisition means and / or facial blood flow acquisition means, and brain age. The presenting means is stored in the smart device.

本発明の第5観点に係る脳年齢提示装置では、スマートデバイスにより情報の入出力がなされるので、脳年齢を簡便に提示することができる。また、脳活動推定手段をサーバに設けることで、演算負荷を軽減することができる。   In the brain age presentation apparatus according to the fifth aspect of the present invention, information is input / output by the smart device, so that the brain age can be presented easily. In addition, by providing the brain activity estimation means in the server, the calculation load can be reduced.

本発明の第6観点に係る脳年齢提示装置は、第1観点から第5観点のいずれかの脳年齢提示装置であって、顔面皮膚温度取得手段及び/又は顔面血行量取得手段と、脳活動推定手段と、脳年齢提示手段とがスマートデバイスに格納されているものである。   A brain age presentation device according to a sixth aspect of the present invention is the brain age presentation device according to any one of the first to fifth aspects, comprising facial skin temperature acquisition means and / or facial blood flow acquisition means, and brain activity. The estimation means and the brain age presentation means are stored in the smart device.

本発明の第6観点に係る脳年齢提示装置では、スマートデバイスにより情報の入出力及び演算がなされるので、脳年齢を簡便に提示することができる。   In the brain age presentation apparatus according to the sixth aspect of the present invention, information input / output and computation are performed by the smart device, so that the brain age can be presented easily.

本発明の第7観点に係る脳年齢提示装置は、第1観点から第6観点のいずれかの脳年齢提示装置であって、顔面皮膚温度取得手段及び/又は顔面血行量取得手段が、ゲーム機に搭載されるものである。また、顔面皮膚温度取得手段及び/又は顔面血行量取得手段は、ゲーム機によるゲーム開始後の所定期間に、顔面皮膚温度データ及び/又は顔面血行量データを時系列で取得する。   A brain age presentation device according to a seventh aspect of the present invention is the brain age presentation device according to any one of the first to sixth aspects, wherein the facial skin temperature acquisition means and / or the facial blood flow acquisition means is a game machine. It is to be mounted on. Further, the facial skin temperature acquisition means and / or the facial blood flow volume acquisition means acquires facial skin temperature data and / or facial blood flow volume data in time series during a predetermined period after the game is started by the game machine.

本発明の第7観点に係る脳年齢提示装置では、上記構成を具備することにより、ゲーム実行中に脳年齢を提示することができる。また、ゲーム機によるゲーム自体が脳機能賦活課題を構成することで、当該ゲームの実行によって脳年齢を提示できる。   In the brain age presentation device according to the seventh aspect of the present invention, the brain age can be presented during the game execution by having the above-described configuration. Moreover, since the game itself by the game machine constitutes the brain function activation task, the brain age can be presented by executing the game.

本発明の第1観点に係る脳年齢提示装置では、簡便に脳年齢を提示することができる。   In the brain age presentation device according to the first aspect of the present invention, the brain age can be presented simply.

本発明の第2,3観点に係る脳年齢提示装置では、脳年齢の推定精度を高めることができる。   In the brain age presentation device according to the second and third aspects of the present invention, it is possible to improve the estimation accuracy of the brain age.

本発明の第4観点に係る脳年齢提示装置では、脳活動と脳年齢との相関関係を逐次最適化できる。   In the brain age presentation device according to the fourth aspect of the present invention, the correlation between brain activity and brain age can be sequentially optimized.

本発明の第5,6観点に係る脳年齢提示装置では、スマートデバイスを利用して脳年齢を簡便に提示できる。   In the brain age presentation apparatus according to the fifth and sixth aspects of the present invention, the brain age can be simply presented using a smart device.

本発明の第7観点に係る脳年齢提示装置では、ゲーム実行中に脳年齢を提示することができる。   In the brain age presentation device according to the seventh aspect of the present invention, the brain age can be presented during game execution.

撮影画像データ及びこれを解析した結果の一例を示す図。The figure which shows an example of picked-up image data and the result of having analyzed this. 顔面皮膚温度データを解析した結果の一部を示す図。The figure which shows a part of result of having analyzed facial skin temperature data. 顔面皮膚温度データを解析した結果の一部を示す図。The figure which shows a part of result of having analyzed facial skin temperature data. 成分2の成分波形の振幅と、測定された脳波のうちのβ波の振幅とを示す図。The figure which shows the amplitude of the component waveform of the component 2, and the amplitude of (beta) wave among the measured brain waves. 成分3の成分波形の振幅と、測定された脳波のうちのβ波の振幅とを示す図。The figure which shows the amplitude of the component waveform of the component 3, and the amplitude of (beta) wave among the measured brain waves. 対照実験で得られた顔面皮膚温度データを解析した結果の一部を示す図。The figure which shows a part of result of having analyzed the facial skin temperature data obtained by the control experiment. 顔面の撮影画像データに基づく成分波形と、測定された脳波のうちのβ波の振幅とを示す図。The figure which shows the component waveform based on the picked-up image data of a face, and the amplitude of (beta) wave among the measured brain waves. 顔面皮膚温度データに基づく成分波形と、測定された脳波のうちのβ波の振幅とを示す図。The figure which shows the component waveform based on facial skin temperature data, and the amplitude of (beta) wave among the measured brain waves. 顔面の撮影画像データに基づく成分波形と、測定された脳波のうちのβ波の振幅とを示す図。The figure which shows the component waveform based on the picked-up image data of a face, and the amplitude of (beta) wave among the measured brain waves. 顔面皮膚温度データに基づく成分波形と、測定された脳波のうちのβ波の振幅とを示す図。The figure which shows the component waveform based on facial skin temperature data, and the amplitude of (beta) wave among the measured brain waves. 顔面の撮影画像データに基づく成分波形と、測定された脳波のうちのβ波の振幅とを示す図。The figure which shows the component waveform based on the picked-up image data of a face, and the amplitude of (beta) wave among the measured brain waves. 顔面皮膚温度データに基づく成分波形と、測定された脳波のうちのβ波の振幅とを示す図。The figure which shows the component waveform based on facial skin temperature data, and the amplitude of (beta) wave among the measured brain waves. 顔面の撮影画像データに基づく成分波形と、測定された脳波のうちのβ波の振幅とを示す図。The figure which shows the component waveform based on the picked-up image data of a face, and the amplitude of (beta) wave among the measured brain waves. 顔面皮膚温度データに基づく成分波形と、測定された脳波のうちのβ波の振幅とを示す図。The figure which shows the component waveform based on facial skin temperature data, and the amplitude of (beta) wave among the measured brain waves. 顔面の撮影画像データに基づく成分波形と、測定された脳波のうちのβ波の振幅とを示す図。The figure which shows the component waveform based on the picked-up image data of a face, and the amplitude of (beta) wave among the measured brain waves. 顔面皮膚温度データに基づく成分波形と、測定された脳波のうちのβ波の振幅とを示す図。The figure which shows the component waveform based on facial skin temperature data, and the amplitude of (beta) wave among the measured brain waves. 顔面の撮影画像データに基づく成分波形と、測定された脳波のうちのβ波の振幅とを示す図。The figure which shows the component waveform based on the picked-up image data of a face, and the amplitude of (beta) wave among the measured brain waves. 顔面皮膚温度データに基づく成分波形と、測定された脳波のうちのβ波の振幅とを示す図。The figure which shows the component waveform based on facial skin temperature data, and the amplitude of (beta) wave among the measured brain waves. 本発明の一実施形態に係る脳活動可視化装置の概略図。1 is a schematic diagram of a brain activity visualization apparatus according to an embodiment of the present invention. 脳活動可視化装置において脳機能を反映した皮膚温度の変化を示す成分を同定する際の処理の流れの一例を示すフローチャート。The flowchart which shows an example of the flow of a process at the time of identifying the component which shows the change of the skin temperature which reflected the brain function in the brain activity visualization apparatus. 本発明の一実施形態に係る脳活動可視化装置の概略図。1 is a schematic diagram of a brain activity visualization apparatus according to an embodiment of the present invention. 脳活動可視化装置において脳機能を反映した顔面のRGB変化を示す成分を同定する際の処理の流れの一例を示すフローチャート。The flowchart which shows an example of the flow of a process at the time of identifying the component which shows the RGB change of the face reflecting a brain function in a brain activity visualization apparatus. 本発明の一実施形態に係る脳年齢提示装置の概略図。The schematic diagram of the brain age presentation device concerning one embodiment of the present invention.

本発明の実施形態を説明する前に、まず、本発明者らが本発明を為すにあたって重要な基礎となった、本発明者らによる知見について説明する。   Before describing the embodiments of the present invention, first, the knowledge obtained by the present inventors, which is an important basis for the inventors of the present invention, will be described.

(1)本発明者らによる知見の要点
人間の脳活動には、人間の知的活動(認知活動等)及び情動活動(快/不快等の活動)が反映されていることが知られている。そして、従来、人間の脳活動を推定する試みがされているが、この場合、脳波計測法、磁気共鳴画像法及び近赤外線分光法のいずれかの方法によって検出されたデータが利用されることが多い。
(1) Key points of knowledge by the present inventors It is known that human brain activities reflect human intellectual activities (cognitive activities, etc.) and emotional activities (activity such as pleasure / discomfort). . Conventionally, attempts have been made to estimate human brain activity. In this case, data detected by any of the electroencephalography, magnetic resonance imaging, and near-infrared spectroscopy may be used. Many.

ここで、検出方法として、例えば、脳波計測法が採用される場合には、被験者に対して脳波電極を装着する必要がある。そして、脳波電極を装着する際には皮膚と電極との間の抵抗を小さくする必要があるため、皮膚を研磨する処理を行ったり電極にペーストを塗布したりする等の作業が必要になる。また、磁気共鳴画像法が採用される場合には、MRI室以外での測定が不可能であるとともに、測定室内に金属を持ち込むことができない等の測定条件に制約がある。さらに、近赤外線分光法が採用される場合には、被験者に対してプローブを装着する必要があるが、プローブを長時間装着することで被験者が痛みを感じたり、被験者の髪とプローブとの接触具合によっては正確に検出できなかったりすることがある。このように、人間の脳活動を測定するために従来の検出方法が採用される場合、脳波電極やプローブ等を装着する際の前処理が必要であったり、測定条件が限定されたりする等、被験者に与える負担が大きくなる。   Here, for example, when an electroencephalogram measurement method is employed as the detection method, it is necessary to attach an electroencephalogram electrode to the subject. And when attaching an electroencephalogram electrode, since it is necessary to make resistance between skin and an electrode small, the process of grind | polishing the skin or apply | coating a paste to an electrode is needed. In addition, when magnetic resonance imaging is adopted, there are restrictions on measurement conditions such as that measurement outside the MRI room is impossible and metal cannot be brought into the measurement room. Furthermore, when near-infrared spectroscopy is used, it is necessary to wear a probe on the subject. However, if the probe is worn for a long time, the subject feels pain or the subject's hair contacts the probe. Depending on the condition, it may not be detected accurately. In this way, when a conventional detection method is adopted to measure human brain activity, pre-processing when wearing an electroencephalogram electrode or a probe is necessary, measurement conditions are limited, etc. The burden on the subject increases.

したがって、被験者の負担を軽減し、かつ簡便に人間の脳活動を推定できる手段の開発が求められている。   Therefore, there is a need for development of means that can reduce the burden on the subject and can easily estimate human brain activity.

そして、本発明者らは、人間の顔面の皮膚温度又は顔面の皮膚温度に比例すると考えられている顔面の血行状態に基づき人間の脳活動を推定することができるのではないか、と考えた。人間の顔面の皮膚温度であればサーモグラフィ等の測定装置を用いることで取得することができ、顔面の血行状態すなわち顔面の血行量であれば撮影装置を利用して得られる顔面の撮影画像のRGBデータから推定することができる。このように、顔面の皮膚温度や顔面の撮影画像であれば、脳波電極やプローブ等の装着前に処理が必要なセンサを装着することなく取得することができる。   And the present inventors thought that human brain activity could be estimated based on the blood circulation state of the face that is considered to be proportional to the skin temperature of the human face or the skin temperature of the face. . If it is the skin temperature of the human face, it can be obtained by using a measuring device such as a thermography, and if the blood circulation state of the face, that is, the blood circulation amount of the face, RGB of the captured image of the face obtained using the imaging device It can be estimated from the data. As described above, a facial skin temperature or a photographed image of a face can be acquired without wearing a sensor that requires processing before wearing an electroencephalogram electrode or a probe.

一方で、人間の顔面の皮膚温度は、外気温度及び/又は自律神経の活動等の様々な要因の影響を受けて変化することが知られている。このため、顔面の皮膚温度に基づいて又は顔面の皮膚温度に比例すると考えられる顔面の血行量に基づいて脳活動を推定しようとすると、取得したデータが脳活動のみを反映しているかどうかを判断することは、非常に困難であると考えられる。   On the other hand, it is known that the skin temperature of a human face changes under the influence of various factors such as the outside air temperature and / or autonomic nerve activity. Therefore, when trying to estimate brain activity based on facial skin temperature or based on facial blood flow, which is considered to be proportional to facial skin temperature, it is determined whether the acquired data reflects only brain activity. It seems to be very difficult to do.

本発明者らは、鋭意検討した結果、顔面の皮膚温度を検出し、検出した温度データ及び検出部位の位置データ(座標データ)を含む時系列の顔面皮膚温度データを、或いは、時系列の顔面の撮影画像データから得られるRGBデータに基づき算出された時系列の顔面の血行量データを、特異値分解法、主成分分析法若しくは独立成分分析法を用いて複数の成分に分解し、分解した複数の成分について解析を行うことで、脳活動を反映した顔面の皮膚温度の変化或いは顔面の血行量の変化を示す成分を同定することができることを見いだした。そして、本発明者らは、対象者の脳活動を推定し、これを解析することで、推定した脳活動に基づき対象者の生理状態を可視化することのできる本発明に到達した。   As a result of intensive studies, the present inventors have detected the skin temperature of the face, time-series facial skin temperature data including detected temperature data and position data (coordinate data) of the detected part, or time-series face Time-series facial blood flow data calculated based on RGB data obtained from the captured image data was decomposed into multiple components using the singular value decomposition method, principal component analysis method, or independent component analysis method, and decomposed. By analyzing a plurality of components, it was found that components showing changes in facial skin temperature or facial blood flow reflecting brain activity can be identified. The inventors of the present invention have arrived at the present invention that can visualize the subject's physiological state based on the estimated brain activity by estimating and analyzing the subject's brain activity.

(2)顔面の各種データの取得方法、及び取得した各種データの解析方法
(2−1)顔面皮膚温度データの取得方法、及び顔面皮膚温度データの解析方法
次に、本発明者らが上記の知見を得るに際して用いた顔面皮膚温度データの取得方法、及び顔面皮膚温度データの解析方法について説明する。
(2) Method for acquiring various facial data, and method for analyzing various acquired data (2-1) Method for acquiring facial skin temperature data, and method for analyzing facial skin temperature data Next, the inventors described above A method of acquiring facial skin temperature data and a method of analyzing facial skin temperature data used for obtaining knowledge will be described.

この試験では、6名の被験者から顔面皮膚温度データを取得した。具体的には、室温25℃を維持した人工気象室内に設置した椅子に被験者を座らせて、赤外線サーモグラフィ装置を用いて、被験者の顔面全体から顔面皮膚温度データを取得した。赤外線サーモグラフィ装置は、対象物から出ている赤外線放射エネルギーを赤外線カメラで検出し、検出した赤外線放射エネルギーを対象物表面の温度(ここでは、摂氏での温度)に変換して、その温度分布を顔面皮膚温度データ(例えば、温度分布を表した画像データ)として表示、蓄積することが可能な装置である。なお、この試験では、赤外線サーモグラフィ装置として、NEC Avio 赤外線テクノロジー株式会社製のR300を使用した。また、赤外線カメラは、被験者の正面であって、被験者から1.5m離れた地点に設置した。そして、顔面皮膚温度データは、30分間取得した。   In this test, facial skin temperature data was obtained from six subjects. Specifically, the subject was seated on a chair installed in an artificial weather room maintained at a room temperature of 25 ° C., and facial skin temperature data was obtained from the entire face of the subject using an infrared thermography device. The infrared thermography device detects the infrared radiation energy emitted from the object with an infrared camera, converts the detected infrared radiation energy into the temperature of the object surface (in this case, the temperature in degrees Celsius), and converts the temperature distribution. It is a device that can display and accumulate as facial skin temperature data (for example, image data representing a temperature distribution). In this test, NEC Avio Infrared Technology Co., Ltd. R300 was used as an infrared thermography apparatus. The infrared camera was installed in front of the subject at a point 1.5 m away from the subject. The facial skin temperature data was acquired for 30 minutes.

また、この試験では、顔面皮膚温度データを取得している間に、被験者に対して脳機能賦活課題を与えた。これにより、脳の非賦活時の顔面皮膚温度データ、及び脳の賦活時の顔面皮膚温度データを取得した。脳機能賦活課題としては、被験者が表示装置等に表示された映像に基づいて、計算、又は、数値、形状及び色の認知、或いは、記号、文字ないし言語の記憶などの心理的作業が挙げられる。この試験では、脳機能賦活課題として「かけ算の暗算」を採用し、被験者に、表示装置に筆算形式で表示される数字を計算させ、その回答をキーボードに入力させる作業を課した。なお、この試験では、顔面皮膚温度データの取得開始から5分経過後から10分間継続して、被験者に対して脳機能賦活課題を与えた。   In this test, a brain function activation task was given to the subject while acquiring facial skin temperature data. As a result, facial skin temperature data when the brain was not activated and facial skin temperature data when the brain was activated were acquired. Examples of the brain function activation task include calculation or recognition of numerical values, shapes and colors, or psychological work such as memory of symbols, characters or languages based on the image displayed on the display device or the like by the subject. . In this test, “Mental calculation of multiplication” was adopted as a brain function activation task, and subjects were asked to calculate numbers displayed in a writing format on the display device and to input the answers to the keyboard. In this test, a brain function activation task was given to the subject continuously for 10 minutes after the start of acquisition of facial skin temperature data for 5 minutes.

顔面皮膚温度データの解析としては、取得した顔面皮膚温度データを対象として、MATLAB(登録商標)のSVD(Singular Value Decomposition)を分析ツールとして用いて特異値分解を行った。特異値分解では、時系列で取得した全ての顔面皮膚温度データ(30分間のデータ)を対象とし、要因を30秒毎の時間データ(30分間で60 time point)とし、測度をその期間(30秒間)における顔面皮膚温度データ(240×320 pixels)とした。そして、特異値分解により、顔面皮膚温度データXを、複数の成分に分解し、それぞれの成分の時間分布Vと、空間分布Uと、各成分の大きさを示す特異値Sとを算出した。なお、これらの関係は、以下の式で表される。また、V’は、Vの行と列とを入れ替えた行列である。   For analysis of facial skin temperature data, singular value decomposition was performed on the acquired facial skin temperature data using SVLAB (Singular Value Decomposition) of MATLAB (registered trademark) as an analysis tool. In singular value decomposition, all facial skin temperature data (data for 30 minutes) acquired in time series are targeted, the factor is time data every 30 seconds (60 time points in 30 minutes), and the measure is the period (30 Second) facial skin temperature data (240 × 320 pixels). Then, the facial skin temperature data X was decomposed into a plurality of components by singular value decomposition, and a temporal distribution V of each component, a spatial distribution U, and a singular value S indicating the size of each component were calculated. These relationships are expressed by the following equations. V ′ is a matrix in which the rows and columns of V are exchanged.

そして、特異値分解によって求められた各成分の時間分布V及び空間分布Uをグラフにプロットし、各成分の成分波形図と温度分布図とを作成した。   And the time distribution V and the spatial distribution U of each component calculated | required by the singular value decomposition were plotted on the graph, and the component waveform diagram and temperature distribution diagram of each component were created.

さらに、作成した各成分の成分波形図及び温度分布図について、脳活動を反映した皮膚温度の変化を示す成分を同定するための解析を行った。   Furthermore, the analysis for identifying the component which shows the change of the skin temperature reflecting brain activity was performed about the created component waveform diagram and temperature distribution diagram of each component.

各成分の成分波形図については、その成分波形の振幅と、脳の非賦活時及び脳の賦活時との相関関係の有無について解析した。具体的には、各成分の成分波形図に示された振幅と、脳の非賦活期間/脳の賦活期間との間に相関関係があるか否かを評価した。この試験では、顔面皮膚温度データを取得している期間のうち、被験者に対して脳機能賦活課題が与えられていない期間であるデータ取得開始時点から5分が経過した時点までの5分間の期間、及びデータ取得開始時から15分が経過した時点からデータ取得終了時点までの15分間の期間を脳の非賦活時とし、被験者に対して脳機能賦活課題が与えられている期間であるデータ取得開始時から5分が経過した時点から10分が経過した時点までの10分間の期間を脳の賦活時とした。そして、各成分の成分波形図に示された振幅と、脳の非賦活時及び脳の賦活時との相関関係の有無について評価した。なお、相関関係の有無については、統計的相関分析を行い、有意水準(α)が0.05以下の場合に相関があると判断した。   About the component waveform figure of each component, it analyzed about the presence or absence of the correlation of the amplitude of the component waveform, and the time of brain non-activation and the time of brain activation. Specifically, it was evaluated whether there was a correlation between the amplitude shown in the component waveform diagram of each component and the brain non-activation period / brain activation period. In this test, among the period of acquiring facial skin temperature data, a period of 5 minutes from the start of data acquisition, which is a period during which no brain function activation task is given to the subject, Data acquisition is a period in which a brain function activation task is given to a subject, with a period of 15 minutes from the time when 15 minutes have elapsed from the start of data acquisition to the end of data acquisition being the non-activation time of the brain The period of 10 minutes from the time when 5 minutes passed from the start to the time when 10 minutes passed was defined as the time of brain activation. Then, the presence or absence of a correlation between the amplitude shown in the component waveform diagram of each component and the time of brain non-activation and the time of brain activation was evaluated. In addition, about the presence or absence of correlation, statistical correlation analysis was conducted and it was judged that there was a correlation when a significance level ((alpha)) is 0.05 or less.

各成分の温度分布図については、顔面の所定部位における温度変化の有無について解析した。ここで、脳には、選択的脳冷却機構(Selective Brain Cooling System)という体温とは独立して脳を冷却する仕組みがある。選択的脳冷却機構としては、脳活動によって生じた熱を前額部及び副鼻腔周辺(眉間及び鼻部周辺を含む)を用いて排熱していることが知られている。そこで、この試験では、各成分の温度分布図において、副鼻腔周辺及び前額部における温度変化があるか否かを評価した。なお、温度分布図における副鼻腔周辺及び前額部の温度変化の有無については、目視(visual inspection)による温度変化の有無、もしくは副鼻腔周辺及び前額部の温度が測定データ全体の平均温度から1標準偏差(SD)以上異なるか否かを温度変化の有無の基準とした。   Regarding the temperature distribution chart of each component, the presence or absence of temperature change at a predetermined part of the face was analyzed. Here, the brain has a mechanism called a selective brain cooling system that cools the brain independently of the body temperature. As a selective brain cooling mechanism, it is known that the heat generated by the brain activity is exhausted using the forehead and the sinuses (including the eyebrows and the nose). Therefore, in this test, it was evaluated whether or not there was a temperature change around the sinuses and the forehead in the temperature distribution chart of each component. In addition, regarding the presence or absence of temperature changes around the sinuses and the forehead in the temperature distribution chart, the presence or absence of temperature changes by visual inspection, or the temperature around the sinuses and the forehead is determined from the average temperature of the entire measurement data Whether there was a difference of 1 standard deviation (SD) or more was used as a reference for the presence or absence of temperature change.

なお、空間分布U、特異値S及び時間分布Vの値の関係で、顔面皮膚温度データXの極性(プラスマイナス)が決定するため、各成分の成分波形図及び温度分布図において極性が反転して現れることがある。このため、成分波形図及び温度分布図の評価に関して、極性については評価対象としないこととした。   Since the polarity (plus or minus) of the facial skin temperature data X is determined by the relationship between the spatial distribution U, the singular value S, and the time distribution V, the polarity is inverted in the component waveform diagram and temperature distribution diagram of each component. May appear. For this reason, regarding the evaluation of the component waveform diagram and the temperature distribution diagram, the polarity is not considered as an evaluation target.

ここで、この赤外線サーモグラフィ装置では、上述しているように、対象物から検出された赤外線放射エネルギーを温度に変換して、その温度分布を顔面皮膚温度データとしている。ところで、人間を対象として赤外線サーモグラフィ装置を用いて顔面の皮膚温度を取得する場合、顔面の動き及び/又は自律神経の活動等の様々な脳活動とは関連しない温度変化(いわゆるノイズ)についても顔面皮膚温度データとして取得してしまう。そこで、このような脳活動とは関連しない温度変化を検出するために、30秒毎の顔面皮膚温度データに含まれる温度データの全平均値を「0」とした相対的な顔面皮膚温度データを作成し、作成した顔面皮膚温度データについても、MATLAB(登録商標)のSVDを分析ツールとして用いて特異値分解を行い、特異値Sに応じた各成分の成分波形図と温度分布図とを作成し、脳活動を反映した皮膚温度の変化を示す成分を同定するための解析を行った。   Here, in this infrared thermography apparatus, as described above, infrared radiation energy detected from an object is converted into temperature, and the temperature distribution is used as facial skin temperature data. By the way, when the skin temperature of the face is acquired by using an infrared thermography apparatus for human subjects, temperature changes (so-called noise) that are not related to various brain activities such as facial movement and / or autonomic nerve activity are also faced. It is acquired as skin temperature data. Therefore, in order to detect such a temperature change not related to the brain activity, relative facial skin temperature data in which the total average value of the temperature data included in the facial skin temperature data every 30 seconds is “0” is calculated. The facial skin temperature data created is also subjected to singular value decomposition using the MATLAB® SVD as an analysis tool, and a component waveform diagram and temperature distribution diagram of each component according to the singular value S are created. Then, an analysis was performed to identify a component showing a change in skin temperature reflecting brain activity.

なお、以下より、説明の便宜上、赤外線サーモグラフィ装置で取得した顔面皮膚温度データを「温度換算データに応じた顔面皮膚温度データ」といい、所定時間毎(この試験では30秒毎)の温度換算データに応じた顔面皮膚温度データに含まれる温度データの全平均値を「0」とした相対的な顔面皮膚温度データを「相対温度換算データに応じた顔面皮膚温度データ」という。   In the following, for convenience of explanation, the facial skin temperature data acquired by the infrared thermography device is referred to as “facial skin temperature data according to the temperature converted data”, and the temperature converted data every predetermined time (in this test every 30 seconds). The relative facial skin temperature data in which the total average value of the temperature data included in the facial skin temperature data corresponding to “0” is “0” is referred to as “facial skin temperature data according to relative temperature conversion data”.

また、6名の被験者のうちの1名に対しては、赤外線サーモグラフィ装置による顔面皮膚温度の検出の他に、被験者の頭皮上に電極を接続して脳波を測定し、覚醒時や意識が緊張した時に現れる波形として知られているβ波(14〜30Hzの周波数の脳波)の振幅と、成分波形図の振幅との間の相関関係についても評価した。なお、脳波測定では、国際式10−20法に基づき、6つの部位(F3、F4,C3、C4、Cz、Pz)に電極を配置した。   For one of the six subjects, in addition to detecting the facial skin temperature using an infrared thermography device, electrodes were connected to the subject's scalp to measure brain waves, and during wakefulness and consciousness became tense. The correlation between the amplitude of the β wave (the brain wave having a frequency of 14 to 30 Hz), which is known as the waveform appearing at the time, and the amplitude of the component waveform diagram was also evaluated. In the electroencephalogram measurement, electrodes were arranged at six parts (F3, F4, C3, C4, Cz, Pz) based on the international method 10-20.

ところで、被験者に脳機能賦活課題が与えられている間、被験者の頭が上下に動くことが考えられる。そうすると、赤外線カメラに対する被験者の顔面の位置が変化することになる。この顔面の位置の変化が皮膚温度の変化に影響しているか否かを検証するために、被験者1名に対して対照試験を行った。顔面皮膚温度データを取得する際の被験者の動きの影響を検証するための対照試験では、上記試験と同様に赤外線サーモグラフィ装置を用いて被験者の顔面皮膚温度データを取得するが、脳機能賦活課題が与えられていない間(すなわち、脳の非賦活時)についてもランダムなタイミングでキーボードを押す作業を被験者に課した。この対照実験によって得られた温度換算データに応じた顔面皮膚温度データ及び相対温度換算データに応じた顔面皮膚温度データについても、MATLAB(登録商標)のSVDを分析ツールとして用いて特異値分解を行い、特異値Sに応じた各成分の成分波形図と温度分布図とを作成し、脳活動を反映した皮膚温度の変化を示す成分を同定するための解析を行った。   By the way, it is conceivable that the subject's head moves up and down while the subject is given a brain function activation task. Then, the position of the subject's face with respect to the infrared camera changes. In order to verify whether or not the change in the position of the face affects the change in the skin temperature, a control test was performed on one subject. In the control test to verify the influence of the subject's movement when acquiring the facial skin temperature data, the facial skin temperature data of the subject is acquired using the infrared thermography device as in the above test, but the brain function activation task is The subject was forced to press the keyboard at random times even when it was not given (ie, when the brain was not activated). The facial skin temperature data according to the temperature conversion data obtained by the control experiment and the facial skin temperature data according to the relative temperature conversion data are also subjected to singular value decomposition using the MATLAB® SVD as an analysis tool. Then, a component waveform diagram and a temperature distribution diagram of each component according to the singular value S were created, and an analysis for identifying a component showing a change in skin temperature reflecting brain activity was performed.

(2−2)顔面撮影画像データの取得方法、及び顔面撮影画像データの解析方法
図1(a)は、撮影装置にて撮影した被験者の顔面の副鼻腔周辺の撮影画像データの一例を示す図である。図1(b)は、血行量分布図(画像マップ)の一例を示す図である。
(2-2) Method for Acquiring Face Shooting Image Data and Method for Analyzing Face Shooting Image Data FIG. 1A is a diagram showing an example of shooting image data around the paranasal sinuses of a subject's face taken by a photographing device. It is. FIG. 1B is a diagram illustrating an example of a blood flow distribution map (image map).

次に、本発明者らが上記の知見を得るに際して用いた顔面撮影画像データの取得方法、及び顔面撮影画像データの解析方法について説明する。   Next, a method for acquiring face-captured image data and a method for analyzing face-captured image data used by the inventors to obtain the above knowledge will be described.

この試験では、6名の被験者から顔面の撮影画像データを取得した。具体的には、室温25℃を維持した人工気象室内に設置した椅子に被験者を座らせて、時系列で画像を取得可能な撮影装置を用いて、被験者の顔面全体の副鼻腔周辺の撮影画像データを時系列で取得した。   In this test, face image data was obtained from six subjects. Specifically, a photographed image around the sinus of the entire face of the subject using a photographing device that allows the subject to sit in a chair installed in an artificial weather room maintained at room temperature of 25 ° C. and acquire images in time series Data was acquired in time series.

また、上述した選択的脳冷却機構に基づくと、脳活動に伴う顔面皮膚温度に比例すると考えられる顔面の血行量の変化は、前額部及び/又は副鼻腔周辺に出現すると考えられる。このことから、本発明者らは、少なくとも前額部及び/又は副鼻腔周辺の顔面の血行量の変化を捉えることができれば、精度良く脳活動を推定することができる、と考えた。そして、この試験では、被験者の顔面の副鼻腔周辺の撮影画像データを時系列で取得した。   Further, based on the selective brain cooling mechanism described above, a change in facial blood flow that is considered to be proportional to the facial skin temperature accompanying brain activity is considered to appear around the forehead and / or the sinuses. From this, the present inventors considered that brain activity can be estimated with high accuracy if at least changes in blood flow in the face around the forehead and / or sinuses can be captured. In this test, captured image data around the paranasal sinuses of the subject's face was acquired in time series.

また、この試験では、撮影装置として、Apple社製のiPad Air(登録商標)の備える液晶画面側の撮影装置を使用し、時系列の撮影画像データとしてカラーの動画データを取得した。また、撮影装置を、被験者の正面側であって、被験者から1.0m離れた地点に設置した。そして、撮影装置によって、30フレーム/秒の撮影周期で時間軸に沿って30分間の撮影画像データを連続撮影することで、顔面の動画データを得た。   In this test, a photographing device on the liquid crystal screen side of the iPad Air (registered trademark) manufactured by Apple Inc. was used as a photographing device, and color moving image data was acquired as time-series photographed image data. In addition, the photographing apparatus was installed on the front side of the subject and at a point 1.0 m away from the subject. Then, the moving image data of the face was obtained by continuously capturing the captured image data for 30 minutes along the time axis at a capturing period of 30 frames / second by the capturing device.

さらに、この試験では、顔面の動画データを取得している間に、被験者に対して脳機能賦活課題を与えた。これにより、脳の非賦活時の顔面の動画データ、及び脳の賦活時の顔面の動画データを取得した。この試験では、上記試験と同様に、脳機能賦活課題として「かけ算の暗算」を採用し、被験者に、表示装置に筆算形式で表示される数字を計算させ、その回答をキーボードに入力させる作業を課した。なお、この試験では、顔面の動画データの取得開始から5分経過後から10分間継続して、被験者に対して脳機能賦活課題を与えた。   Furthermore, in this test, a brain function activation task was given to the subject while acquiring facial moving image data. Thereby, the moving image data of the face when the brain was not activated and the moving image data of the face when the brain was activated were acquired. In this test, as in the case of the above test, the task of activating the brain function is `` multiplying mental arithmetic '', allowing the subject to calculate the numbers displayed in the writing format on the display device, and inputting the answers to the keyboard. Imposed. In this test, a brain function activation task was given to the subject continuously for 10 minutes after 5 minutes from the start of facial moving image data acquisition.

顔面の動画データの解析としては、撮影した顔面の動画データより得られたRGBデータに基づき血行量データを算出し、算出した時系列の血行量データを対象として、MATLAB(登録商標)のSVDを分析ツールとして用いて特異値分解を行った。ここでは、CIE−L***表色系に従って、画像のRGBデータより演算される肌の赤みやヘモグロビン量と相関のある紅斑指数「a*」を求め、これを血行量データとした。また、特異値分解では、時系列で取得した全ての動画データ(30分間のデータ)から得られたRGBデータに基づく血行量データ(ここでは、紅斑指数)を対象とし、要因を30秒毎の時間データ(30分間で60 time point)とし、測度をその期間(30秒毎)におけるRGBデータから演算した紅斑指数(30秒毎に1秒間のフレームデータを取り出し、該フレームデータから得られるRGB値の平均値から演算した紅斑指数;240×320 pixels)とした。そして、特異値分解により、顔面の動画データより得られたRGBデータに基づく時系列の血行量データを、複数の成分に分解し、それぞれの成分の時間分布Vと、空間分布Uと、各成分の大きさを示す特異値Sとを算出した。なお、これらの関係は、上記式(数1)と同様の式で表される。 For analysis of facial moving image data, blood flow data is calculated based on RGB data obtained from the captured moving image data of the face, and MATLAB (registered trademark) SVD is used for the calculated time-series blood flow data. Singular value decomposition was performed as an analytical tool. Here, according to the CIE-L * a * b * color system, an erythema index “a * ” correlated with the redness and hemoglobin amount calculated from the RGB data of the image is obtained, and this is used as blood flow data. . In the singular value decomposition, blood flow volume data (here, erythema index) based on RGB data obtained from all moving image data (30-minute data) acquired in time series is targeted, and the factor is set every 30 seconds. Erythema index calculated from time data (60 time points in 30 minutes) and RGB data in that period (every 30 seconds) (taken out frame data for 1 second every 30 seconds, and RGB values obtained from the frame data) Erythema index calculated from the average value of 240 × 320 pixels). Then, by singular value decomposition, time-series blood flow data based on RGB data obtained from facial moving image data is decomposed into a plurality of components, and each component's time distribution V, spatial distribution U, and each component And a singular value S indicating the magnitude of. In addition, these relationships are represented by the same formula as the above formula (Formula 1).

そして、特異値分解によって求められた各成分の時間分布V及び空間分布Uをグラフにプロットし、各成分の成分波形図と血行量分布図とを作成した。   Then, the time distribution V and the spatial distribution U of each component obtained by singular value decomposition were plotted on a graph, and a component waveform diagram and a blood flow distribution diagram of each component were created.

さらに、作成した各成分の成分波形図及び血行量分布図について、脳活動を反映した顔面の血行量の変化すなわち顔面のRGB変化を示す成分を同定するための解析を行った。   Further, analysis was performed on the component waveform diagram and the blood flow distribution diagram of each component that was created to identify a component indicating a change in facial blood flow that reflects brain activity, that is, a component that indicates a change in facial RGB.

各成分の成分波形図については、その成分波形の振幅と、脳の非賦活時及び脳の賦活時との相関関係の有無について解析した。具体的には、各成分の成分波形図に示された振幅と、脳の非賦活期間/脳の賦活期間との間に相関関係があるか否かを評価した。この試験では、顔面の撮影画像データを取得している期間のうち、被験者に対して脳機能賦活課題が与えられていない期間であるデータ取得開始時点から5分が経過した時点までの5分間の期間、及びデータ取得開始時から15分が経過した時点からデータ取得終了時点までの15分間の期間を脳の非賦活時とし、被験者に対して脳機能賦活課題が与えられている期間であるデータ取得開始時から5分が経過した時点から10分が経過した時点までの10分間の期間を脳の賦活時とした。そして、各成分の成分波形図に示された振幅と、脳の非賦活時及び脳の賦活時との相関関係の有無について評価した。なお、相関関係の有無については、統計的相関分析を行い、有意水準(α)が0.01以下の場合に相関があると判断した。   About the component waveform figure of each component, it analyzed about the presence or absence of the correlation of the amplitude of the component waveform, and the time of brain non-activation and the time of brain activation. Specifically, it was evaluated whether there was a correlation between the amplitude shown in the component waveform diagram of each component and the brain non-activation period / brain activation period. In this test, during the period when the captured image data of the face is acquired, the subject is not given a brain function activation task, and the 5-minute period from the start of data acquisition to the time when 5 minutes have elapsed. Data that is a period in which a brain function activation task is given to a subject, with a period of 15 minutes from the start of data acquisition and the period of 15 minutes from the start of data acquisition to the end of data acquisition being the non-activated state of the brain A period of 10 minutes from the time when 5 minutes passed from the start of acquisition to the time when 10 minutes passed was defined as the time of brain activation. Then, the presence or absence of a correlation between the amplitude shown in the component waveform diagram of each component and the time of brain non-activation and the time of brain activation was evaluated. In addition, regarding the presence or absence of the correlation, statistical correlation analysis was performed, and it was determined that there was a correlation when the significance level (α) was 0.01 or less.

各成分の血行量分布図については、顔面の所定部位における血行量変化の有無について解析した。血行量分布図は、ピクセル毎に算出された空間分布Uを各ピクセルの位置に並べることで作成される。このように作成された各成分の血行量分布図において、副鼻腔周辺及び前額部における血行量の変化があるか否かを評価した。なお、血行量分布図における副鼻腔周辺及び前額部の血行量変化の有無については、目視(visual inspection)による血行量変化の有無、もしくは図1(b)に示す副鼻腔周辺及び前額部の血行量の値が「0.000」でないことを血行量変化の有無の基準とした。   Regarding the blood flow distribution chart of each component, the presence or absence of a change in blood flow at a predetermined part of the face was analyzed. The blood circulation distribution map is created by arranging the spatial distribution U calculated for each pixel at the position of each pixel. In the blood flow distribution chart of each component thus created, it was evaluated whether there was a change in blood flow around the sinuses and in the forehead. In addition, regarding the presence or absence of changes in blood flow around the sinuses and the forehead in the blood flow distribution diagram, the presence or absence of changes in blood flow by visual inspection, or around the sinuses and the forehead shown in FIG. The reference value for the presence or absence of a change in blood flow was that the value of the blood flow was not “0.000”.

なお、空間分布U、特異値S及び時間分布Vの値の関係で、血行量データXの極性(プラスマイナス)が決定するため、各成分の成分波形図及び血行量分布図において極性が反転して現れることがある。このため、成分波形図及び血行量分布図の評価に関して、極性については評価対象としないこととした。   Note that the polarity (plus or minus) of the blood flow data X is determined by the relationship between the spatial distribution U, the singular value S, and the time distribution V. Therefore, the polarity is reversed in the component waveform diagram and the blood flow distribution diagram of each component. May appear. For this reason, regarding the evaluation of the component waveform diagram and the blood flow distribution diagram, the polarity is not considered as an evaluation target.

さらに、顔面の皮膚温度と顔面の血行量との相関関係を検証するために、6名の被験者から顔面の撮影画像データを時系列で取得している間、赤外線サーモグラフィ装置により顔面皮膚温度データも時系列で取得し、取得した顔面皮膚温度データについてもMATLAB(登録商標)のSVDを分析ツールとして用いて特異値分解を行い、特異値Sに応じた各成分の成分波形図を作成し、その成分波形の振幅と、脳の非賦活時及び脳の賦活時との相関関係の有無について解析した。なお、この試験では、赤外線サーモグラフィ装置として、上記試験と同様の装置を用いた。また、赤外線カメラは、被験者の正面であって、被験者から1.5m離れた地点に設置した。   Furthermore, in order to verify the correlation between facial skin temperature and facial blood flow, facial skin temperature data was also obtained by infrared thermography equipment while acquiring facial image data in time series from six subjects. Time-series acquisition, and the acquired facial skin temperature data is also subjected to singular value decomposition using the MATLAB (registered trademark) SVD as an analysis tool, and a component waveform diagram of each component according to the singular value S is created. The presence or absence of a correlation between the amplitude of the component waveform and the brain non-activation and brain activation was analyzed. In this test, an apparatus similar to the above test was used as the infrared thermography apparatus. The infrared camera was installed in front of the subject at a point 1.5 m away from the subject.

ところで、撮影装置を用いて顔面の撮影画像データを取得する場合、撮影中に太陽の光等が顔に当たることで光が顔で反射し、この反射光が撮影装置のレンズに入り込んでしまうことがある。そうすると、撮影された顔面の撮影画像データにはこの反射光が記録されてしまうことになる。ここで、撮影画像データから得られるRGBデータにおいて、顔面の血行量に基づく明度の変化は反射光に基づく明度の変化よりも小さいため、反射光の記録された撮影画像データから得られるRGBデータに基づいて算出された血行量が解析されると、脳活動とは関連しない顔面のRGB変化(いわゆるノイズ)が混入してしまう可能性があると考えられた。そこで、このような脳活動とは関連しない顔面のRGB変化の混入を防ぐために、30秒毎のRGBデータの全平均値を「0」とした相対的なRGBデータから相対的な血行量データを作成し、作成した血行量データについても、MATLAB(登録商標)のSVDを分析ツールとして用いて特異値分解を行い、特異値Sに応じた各成分の成分波形図と血行量分布図とを作成し、脳活動を反映した顔面のRGB変化を示す成分を同定するための解析を行った。   By the way, when acquiring captured image data of a face using a photographing device, light may be reflected by the face when the sunlight hits the face during photographing, and this reflected light may enter the lens of the photographing device. is there. Then, the reflected light is recorded in the photographed image data of the photographed face. Here, in the RGB data obtained from the photographed image data, the change in brightness based on the blood flow amount on the face is smaller than the change in brightness based on the reflected light. Therefore, the RGB data obtained from the photographed image data in which the reflected light is recorded. When the blood flow volume calculated based on the analysis was analyzed, it was considered that there was a possibility that a face RGB change (so-called noise) that is not related to the brain activity might be mixed. Therefore, in order to prevent the mixing of RGB changes in the face that are not related to such brain activity, relative blood flow data is obtained from the relative RGB data with the total average value of RGB data every 30 seconds being “0”. The created blood flow data is also subjected to singular value decomposition using the MATLAB (registered trademark) SVD as an analysis tool, and a component waveform diagram and blood flow distribution diagram of each component according to the singular value S are created. Then, an analysis was performed to identify the component showing the RGB change of the face reflecting the brain activity.

なお、以下より、説明の便宜上、所定時間毎(この試験では30秒毎)のRGBデータの全平均値を「0」とした相対的なRGBデータに基づく相対的な血行量データを「相対換算血行量データ」といい、相対的なRGBデータに換算する前のRGBデータに基づく血行量データを単に「血行量データ」という。   In the following, for convenience of explanation, relative blood flow data based on relative RGB data with a total average value of RGB data every predetermined time (every 30 seconds in this test) being “0” is referred to as “relative conversion. The blood flow data is referred to as “blood flow data”, and the blood flow data based on the RGB data before conversion into relative RGB data is simply referred to as “blood flow data”.

また、6名の被験者に対して撮影装置によって顔面の時系列の撮影画像データを取得している間、各被験者の頭皮上に電極を接続して脳波を測定し、覚醒時等の脳細胞が活動している時に現れる波形として知られているβ波(13〜30Hzの周波数の脳波)の振幅と、成分波形図の振幅との間の相関関係についても評価した。なお、脳波測定では、国際式10−20法に基づき、頭皮上19の部位(Fp1、Fp2、F3、F4、C3、C4、P3、P4、O1、O2、F7、F8、T3、T4、T5、T6、Fz、Cz及びPz)に電極を配置した。   In addition, while acquiring time-series captured image data of the face with the imaging device for 6 subjects, electrodes are connected to the scalp of each subject to measure brain waves, and brain cells at the time of waking etc. The correlation between the amplitude of the β wave (the brain wave having a frequency of 13 to 30 Hz), which is known as a waveform appearing when active, and the amplitude of the component waveform diagram were also evaluated. In the electroencephalogram measurement, the region of the scalp 19 (Fp1, Fp2, F3, F4, C3, C4, P3, P4, O1, O2, F7, F8, T3, T4, T5, based on the international method 10-20) , T6, Fz, Cz and Pz).

さらに、被験者に脳機能賦活課題が与えられている間、被験者の頭が上下に動くことが考えられる。そうすると、撮影装置に対する被験者の顔面の位置が変化することになる。この顔面の位置の変化が顔面のRGB変化に影響しているか否かを検証するために、被験者1名に対して対照試験を行った。対照試験では、上記試験と同様に撮影装置を用いて被験者の顔面の時系列の撮影画像データを取得するが、脳機能賦活課題が与えられていない間(すなわち、脳の非賦活時)についてもランダムなタイミングでキーボードを押す作業を被験者に対して課した。この対照実験によって撮影された顔面の時系列の撮影画像データから得られたRGBデータに基づく時系列の血行量データについても、MATLAB(登録商標)のSVDを分析ツールとして用いて特異値分解を行い、特異値Sに応じた各成分の成分波形図を作成し、その成分波形の振幅と、脳の非賦活時及び脳の賦活時との相関関係の有無について解析した。また、各成分波形の振幅と、実際の顔面の動きとの相関関係の有無について解析した。実際の顔面の動きは、撮影画像データから顔の同一箇所の2次元座標を取得し、対照実験開始時の撮影画像データを基準として撮影時における30秒毎の顔面の移動距離を算出することで評価した。さらに、各成分波形の振幅と、撮影中のキーボードの入力数との相関関係の有無についても解析した。撮影中のキーボードの入力数は、時系列の撮影画像データにおける30秒毎の単純移動平均を算出することで評価した。   Furthermore, it is conceivable that the subject's head moves up and down while the subject is given a brain function activation task. Then, the position of the subject's face with respect to the imaging device changes. In order to verify whether or not the change in the position of the face affects the RGB change in the face, a control test was performed on one subject. In the control test, time-series captured image data of the subject's face is acquired using the imaging device in the same manner as in the above test, but also during the time when the brain function activation task is not given (that is, when the brain is not activated) The subject was forced to press the keyboard at random timing. For chronological blood flow data based on RGB data obtained from time-series photographed image data of the face photographed in this control experiment, singular value decomposition was performed using the SVLAB of MATLAB (registered trademark) as an analysis tool. Then, a component waveform diagram of each component corresponding to the singular value S was created, and the presence or absence of a correlation between the amplitude of the component waveform and the brain non-activation and brain activation was analyzed. In addition, the presence or absence of a correlation between the amplitude of each component waveform and the actual facial movement was analyzed. The actual facial movement is obtained by obtaining the two-dimensional coordinates of the same part of the face from the photographed image data and calculating the movement distance of the face every 30 seconds during photographing with reference to the photographed image data at the start of the control experiment. evaluated. Furthermore, the presence / absence of a correlation between the amplitude of each component waveform and the number of keyboard inputs during shooting was also analyzed. The number of keyboard inputs during shooting was evaluated by calculating a simple moving average every 30 seconds in time-series shot image data.

(3)解析結果
以下、特異値の大きいものから順次「成分1,2,・・・」としている。各成分は全体の分散に占める割合を示しており、変化の大きい影響が成分1に反映される。具体的には、被験者の動作や外界ノイズの影響が成分1に反映される傾向がある。結果として、成分1より成分2の成分波形に脳活動の影響が反映されることが少なくない。
(3) Analysis results Hereinafter, “components 1, 2,... Each component indicates the ratio of the total dispersion, and the effect of large change is reflected in component 1. Specifically, the effects of the subject's movement and external noise tend to be reflected in the component 1. As a result, the influence of brain activity is often reflected in the component waveforms of component 1 and component 2.

(3−1)顔面皮膚温度データの解析結果
図2は、温度換算データに応じた顔面皮膚温度データを解析した結果の一部を示す図である。図2(a)は、被験者1の成分2の成分波形図を示している。図2(b)は、被験者1の成分2の温度分布図を示している。図3(a)は、被験者1の成分3の成分波形図を示している。図3(b)は、被験者1の成分3の温度分布図を示している。図4及び図5は、成分波形の振幅と、脳波との関係を示す図である。図4は、被験者1の成分2の成分波形の振幅と、測定された脳波のうちのβ波の振幅とを示す図である。図5は、被験者1の成分3の成分波形の振幅と、測定された脳波のうちのβ波の振幅とを示す図である。図6は、対照実験で得られた顔面皮膚温度データを解析した結果の一部を示す図である。図6(a)は、成分3の成分波形図を示している。図6(b)は、成分3の温度分布図を示している。
(3-1) Analysis result of facial skin temperature data FIG. 2 is a diagram showing a part of the result of analyzing facial skin temperature data according to temperature conversion data. FIG. 2A shows a component waveform diagram of the component 2 of the subject 1. FIG. 2B shows a temperature distribution diagram of the component 2 of the subject 1. FIG. 3A shows a component waveform diagram of the component 3 of the subject 1. FIG. 3B shows a temperature distribution diagram of the component 3 of the subject 1. 4 and 5 are diagrams showing the relationship between the amplitude of the component waveform and the electroencephalogram. FIG. 4 is a diagram illustrating the amplitude of the component waveform of the component 2 of the subject 1 and the amplitude of the β wave among the measured brain waves. FIG. 5 is a diagram illustrating the amplitude of the component waveform of the component 3 of the subject 1 and the amplitude of the β wave among the measured brain waves. FIG. 6 is a diagram showing a part of the result of analyzing the facial skin temperature data obtained in the control experiment. FIG. 6A shows a component waveform diagram of component 3. FIG. FIG. 6B shows a temperature distribution diagram of the component 3.

表1は、各被験者に対する顔面皮膚温度データの解析結果を示したものである。   Table 1 shows the analysis results of facial skin temperature data for each subject.

上記の顔面皮膚温度データの解析によって得られた結果から、時系列の顔面皮膚温度データを特異値分解により分解して得られた複数の成分のうち、成分2及び/又は成分3と、人間の脳活動との間に有意な相関があることが確認された。   Of the plurality of components obtained by decomposing time-series facial skin temperature data by singular value decomposition from the results obtained by the above analysis of facial skin temperature data, component 2 and / or component 3 and human It was confirmed that there was a significant correlation with brain activity.

また、図4及び図5に示すように、脳波解析の結果から、成分2及び成分3の各成分波形の振幅と、脳波のβ波の振幅との間に有意な相関が確認された。   Also, as shown in FIGS. 4 and 5, a significant correlation was confirmed between the amplitude of each component waveform of component 2 and component 3 and the amplitude of the β wave of the electroencephalogram from the results of the electroencephalogram analysis.

さらに、対照実験では、脳機能賦活課題が与えられていない間に被験者に動きがある状態であっても、成分3と人間の脳活動との間に有意な相関があった(図6参照)。このことから、複数の成分のうち、成分3については、顔面皮膚温度データを取得する際の被験者の動きが影響していないことが認められた。   Furthermore, in the control experiment, there was a significant correlation between component 3 and human brain activity even when the subject was in motion while no brain function activation task was given (see FIG. 6). . From this, it was recognized that among the plurality of components, the movement of the subject at the time of acquiring facial skin temperature data was not affected for component 3.

これらの結果から、本発明者らは、以下の知見を得た。   From these results, the present inventors obtained the following knowledge.

被験者から取得した時系列の顔面皮膚温度データを特異値分解により複数の成分に分解し、分解した各成分について解析した結果、複数の成分のうちの成分3が脳活動に関連する成分であると認められた。すなわち、時系列の顔面皮膚温度データを特異値分解により複数の成分に分解し、分解した複数の成分から脳の賦活/非賦活と相関のある成分を抽出し、抽出した成分について選択的脳冷却機構を利用した解析を行うことで、複数の成分から脳活動を反映した皮膚温度の変化を示す成分を同定することができることが判明した。このことから、本発明者らは、人間の顔面の皮膚温度に基づいて、脳活動を推定することができる、という知見を得た。   The time-series facial skin temperature data obtained from the test subject is decomposed into a plurality of components by singular value decomposition, and as a result of analyzing each decomposed component, component 3 of the plurality of components is a component related to brain activity. Admitted. In other words, time-series facial skin temperature data is decomposed into multiple components by singular value decomposition, components that correlate with brain activation / deactivation are extracted from the decomposed multiple components, and selective brain cooling is performed on the extracted components By analyzing using the mechanism, it was found that a component showing a change in skin temperature reflecting brain activity can be identified from a plurality of components. From this, the present inventors have obtained the knowledge that brain activity can be estimated based on the skin temperature of the human face.

(3−2)顔面の撮影画像データの解析結果
図7〜図18は、顔面の撮影画像データ(血行量データ)又は顔面皮膚温度データに基づく成分波形図を解析した結果の一部を示す図である。図7は、被験者1の撮影画像データに基づく成分2の成分波形の振幅と、測定された被験者1の脳波のうちのβ波の振幅とを示す図である。図8は、被験者1の顔面皮膚温度データに基づく成分2の成分波形の振幅と、測定された被験者1の脳波のうちのβ波の振幅とを示す図である。図9は、被験者2の撮影画像データに基づく成分2の成分波形の振幅と、測定された被験者2の脳波のうちのβ波の振幅とを示す図である。図10は、被験者2の顔面皮膚温度データに基づく成分2の成分波形の振幅と、測定された被験者2の脳波のうちのβ波の振幅とを示す図である。図11は、被験者3の撮影画像データに基づく成分4の成分波形の振幅と、測定された被験者3の脳波のうちのβ波の振幅とを示す図である。図12は、被験者3の顔面皮膚温度データに基づく成分3の成分波形の振幅と、測定された被験者3の脳波のうちのβ波の振幅とを示す図である。図13は、被験者4の撮影画像データに基づく成分3の成分波形の振幅と、測定された被験者4の脳波のうちのβ波の振幅とを示す図である。図14は、被験者4の顔面皮膚温度データに基づく成分2の成分波形の振幅と、測定された被験者4の脳波のうちのβ波の振幅とを示す図である。図15は、被験者5の撮影画像データに基づく成分2の成分波形の振幅と、測定された被験者5の脳波のうちのβ波の振幅とを示す図である。図16は、被験者5の顔面皮膚温度データに基づく成分2の成分波形の振幅と、測定された被験者5の脳波のうちのβ波の振幅とを示す図である。図17は、被験者6の撮影画像データに基づく成分4の成分波形の振幅と、測定された被験者6の脳波のうちのβ波の振幅とを示す図である。図18は、被験者6の顔面皮膚温度データに基づく成分3の成分波形の振幅と、測定された被験者6の脳波のうちのβ波の振幅とを示す図である。
(3-2) Analysis Results of Photographed Image Data of Face FIGS. 7 to 18 are diagrams showing a part of the result of analyzing a component waveform diagram based on photographed image data (blood flow data) or facial skin temperature data of the face. It is. FIG. 7 is a diagram illustrating the amplitude of the component 2 component waveform based on the captured image data of the subject 1 and the measured β wave amplitude of the brain wave of the subject 1. FIG. 8 is a diagram showing the amplitude of the component 2 component waveform based on the facial skin temperature data of the subject 1 and the measured β wave amplitude of the brain wave of the subject 1. FIG. 9 is a diagram showing the amplitude of the component 2 component waveform based on the captured image data of the subject 2 and the measured amplitude of the β wave in the brain wave of the subject 2. FIG. 10 is a diagram illustrating the amplitude of the component 2 component waveform based on the facial skin temperature data of the subject 2 and the measured β wave amplitude of the brain wave of the subject 2. FIG. 11 is a diagram illustrating the amplitude of the component waveform of component 4 based on the captured image data of the subject 3 and the amplitude of the β wave among the measured brain waves of the subject 3. FIG. 12 is a diagram illustrating the amplitude of the component 3 component waveform based on the facial skin temperature data of the subject 3 and the measured β wave amplitude of the brain wave of the subject 3. FIG. 13 is a diagram illustrating the amplitude of the component 3 component waveform based on the captured image data of the subject 4 and the measured β-wave amplitude of the brain wave of the subject 4. FIG. 14 is a diagram showing the amplitude of the component 2 component waveform based on the facial skin temperature data of the subject 4 and the measured β wave amplitude of the brain wave of the subject 4. FIG. 15 is a diagram illustrating the amplitude of the component 2 component waveform based on the captured image data of the subject 5 and the measured β wave amplitude of the brain wave of the subject 5. FIG. 16 is a diagram showing the amplitude of the component 2 component waveform based on the facial skin temperature data of the subject 5 and the measured β wave amplitude of the brain wave of the subject 5. FIG. 17 is a diagram illustrating the amplitude of the component waveform of component 4 based on the captured image data of the subject 6 and the amplitude of the β wave among the measured brain waves of the subject 6. FIG. 18 is a diagram illustrating the amplitude of the component 3 component waveform based on the facial skin temperature data of the subject 6 and the measured β wave amplitude of the brain wave of the subject 6.

図7〜図18に示すように、各成分波形と脳波解析との結果から、顔面の皮膚温度と顔面の血行量とが相関関係にあることが確認された。なお、顔面の皮膚温度データ及び顔面の血行量データのいずれのデータに基づく解析においても、各成分波形の振幅と、頭頂部又は後頭部に装着した電極が測定した脳波のβ波の振幅との間に有意な相関が確認された。   As shown in FIG. 7 to FIG. 18, it was confirmed from the results of the component waveforms and the electroencephalogram analysis that the facial skin temperature and the facial blood flow volume are in a correlation. In the analysis based on any of the skin temperature data of the face and the blood flow data of the face, the amplitude of each component waveform and the amplitude of the β wave of the electroencephalogram measured by the electrode mounted on the top or back of the head A significant correlation was confirmed.

以下に示す表2は、各被験者に対する顔面の撮影画像データの解析結果を示したものである。   Table 2 shown below shows the analysis results of the captured image data of the face for each subject.

表2に示すように、上記の顔面の撮影画像データの解析によって得られた結果から、顔面の撮影画像データに基づく時系列の血行量データを特異値分解により分解して得られた複数の成分のうち、成分1,2,3,4,5と人間の脳活動との間に有意な相関があることが確認された。なお、ここでは、血行量データに基づく相関において有意な相関が見られかつ相対換算血行量データに基づく相関において有意な相関が見られた成分だけでなく、血行量データに基づく相関においては有意な相関が見られなかったが相対換算血行量データに基づく相関において有意な相関が見られた成分も、人間の脳活動と有意な相関があると認めるようにした。   As shown in Table 2, a plurality of components obtained by decomposing time-series blood flow data based on the photographed image data of the face by singular value decomposition based on the results obtained by analyzing the photographed image data of the face. Of these, it was confirmed that there was a significant correlation between components 1, 2, 3, 4, 5 and human brain activity. In addition, here, not only the component in which the correlation based on the blood flow data is significant and the correlation based on the relative blood flow data is significant, but also the correlation based on the blood flow data is significant. Components that were not correlated but showed significant correlation in the correlation based on relative blood flow data were also recognized as having significant correlation with human brain activity.

また、以下に示す表3は、対照実験の結果を示したものである。   Table 3 shown below shows the results of the control experiment.

表3に示すように、対照実験では、脳機能賦活課題が与えられていない間に被験者に動きがある場合、成分1,2の成分波形の振幅と脳の非賦活時及び脳の賦活時との間に有意な相関があった。さらに、成分2については、脳機能賦活課題が与えられていない間に被験者に動きがある場合でも、移動距離及びキーボード入力数それぞれとの間に有意な相関が認められなかった。このことから、顔面の撮影画像データから取得したRGBデータに基づく血行量データを特異値分解することで得られる複数の成分において、脳活動との間に有意な相関がある成分については、顔面の時系列の撮影画像データを取得する際の被験者の動きによる影響を受けたとしても、その影響は脳の脳活動による影響(脳の賦活や非賦活による影響)よりも遙かに小さいことが確認された。   As shown in Table 3, in the control experiment, when the subject moves while the brain function activation task is not given, the amplitude of the component waveforms of components 1 and 2 and the brain activation time and the brain activation time There was a significant correlation between Furthermore, for Component 2, no significant correlation was observed between the movement distance and the number of keyboard inputs even when the subject moved while the brain function activation task was not given. From this, among components obtained by singular value decomposition of blood flow data based on RGB data acquired from captured image data of a face, components having a significant correlation with brain activity are Even if it is affected by the movement of the subject when acquiring time-series captured image data, it is confirmed that the effect is much smaller than the effect of brain activity (effect of brain activation or non-activation). It was done.

これらの結果から、本発明者らは、以下の知見を得た。   From these results, the present inventors obtained the following knowledge.

被験者から取得した時系列の顔面の撮影画像データに基づく顔面のRGBデータから得られる血行量データを特異値分解により複数の成分に分解し、分解した各成分について解析した結果、複数の成分のうちの成分1,2,3,4,5が脳活動に関連する成分であると認められた。すなわち、時系列の顔面の撮影画像データに基づく顔面のRGBデータから得られる血行量データを特異値分解により複数の成分に分解し、分解した複数の成分から脳の賦活/非賦活と相関のある成分を抽出し、抽出した成分について解析することで、複数の成分から脳活動を反映した顔面のRGB変化を示す成分を同定することができることが判明した。このことから、本発明者らは、人間の顔面の時系列の撮影画像データに基づいて、脳活動を推定することができる、という知見を得た。   Blood flow volume data obtained from facial RGB data based on time-series facial captured image data obtained from the subject was decomposed into a plurality of components by singular value decomposition, and each decomposed component was analyzed. Components 1, 2, 3, 4 and 5 were found to be components related to brain activity. That is, blood flow data obtained from facial RGB data based on time-series facial image data is decomposed into a plurality of components by singular value decomposition, and there is a correlation between activation / inactivation of the brain from the decomposed components. It was found that by extracting the components and analyzing the extracted components, it is possible to identify a component showing the RGB change of the face reflecting brain activity from a plurality of components. From this, the present inventors have found that brain activity can be estimated based on time-series captured image data of the human face.

(4)脳活動可視化装置
次に、上記に説明した知見に基づいて、本発明者らが完成するに至った本発明の一実施形態に係る脳活動可視化装置10,110について説明する。なお、本発明に係る脳活動可視化装置は、以下の実施形態に限定されるものではなく、要旨を逸脱しない範囲で適宜変更可能である。
(4) Brain Activity Visualization Device Next, brain activity visualization devices 10 and 110 according to an embodiment of the present invention that have been completed by the present inventors based on the knowledge described above will be described. In addition, the brain activity visualization apparatus according to the present invention is not limited to the following embodiments, and can be appropriately changed without departing from the gist.

本発明の一実施形態に係る脳活動可視化装置10,110は、顔面皮膚温度データに基づき脳活動を推定する脳活動推定手段30、及び/又は顔面の撮影画像データに基づき脳活動を推定する脳活動推定手段130を備えている。以下では、本発明の一実施形態に係る脳活動可視化装置10,110を説明する前に、各脳活動推定手段30,130について説明する。   The brain activity visualization devices 10 and 110 according to an embodiment of the present invention include a brain activity estimation unit 30 that estimates brain activity based on facial skin temperature data, and / or a brain that estimates brain activity based on face image data. The activity estimation means 130 is provided. Below, before demonstrating the brain activity visualization apparatuses 10 and 110 which concern on one Embodiment of this invention, each brain activity estimation means 30 and 130 are demonstrated.

(4−1)顔面皮膚温度データに基づき脳活動を推定する脳活動推定手段30
図19は、本発明の一実施形態に係る脳活動可視化装置10の概略図である。図20は、脳活動可視化装置10において脳機能を反映した皮膚温度の変化を示す成分を同定する際の処理の流れを示すフローチャートである。
(4-1) Brain activity estimating means 30 for estimating brain activity based on facial skin temperature data
FIG. 19 is a schematic diagram of the brain activity visualization apparatus 10 according to an embodiment of the present invention. FIG. 20 is a flowchart showing the flow of processing when the brain activity visualization apparatus 10 identifies a component indicating a change in skin temperature reflecting the brain function.

脳活動可視化装置10の備える脳活動推定手段30は、個人(被験者)の顔面の皮膚温度から、個人の脳活動を推定する。脳活動可視化装置10は、図19に示すように、顔面皮膚温度取得手段20と、脳活動推定手段30と、状態可視化手段200と、を備える。   The brain activity estimating means 30 provided in the brain activity visualization device 10 estimates the brain activity of the individual from the skin temperature of the face of the individual (subject). As shown in FIG. 19, the brain activity visualization device 10 includes a facial skin temperature acquisition unit 20, a brain activity estimation unit 30, and a state visualization unit 200.

顔面皮膚温度取得手段20は、個人の顔面の少なくとも一部の皮膚温度を検出し、検出した温度データ及びその検出部位の位置データを含む顔面皮膚温度データを時系列で取得する(ステップS1)。なお、ここでは、顔面皮膚温度取得手段20は、赤外線サーモグラフィ装置であり、図19に示すように、赤外線カメラ21と、処理部22と、を有する。赤外線カメラ21は、個人の顔面から出ている赤外線放射エネルギーを検出するためのものである。そして、ここでは、赤外線カメラ21は、個人の顔面全体から赤外線放射エネルギーを検出しているものとする。処理部22は、赤外線カメラ21によって検出した赤外線放射エネルギーを温度に変換して温度データとし、赤外線放射エネルギーの検出された部位を位置データ(座標データ)とした顔面全体における顔面皮膚温度の温度分布図を作成し、作成した温度分布図を温度換算データに応じた顔面皮膚温度データとして処理する。温度換算データに応じた顔面皮膚温度データは、処理部22の有する記憶部(図示せず)に蓄積される。   The facial skin temperature acquisition means 20 detects the skin temperature of at least a part of the individual's face, and acquires facial skin temperature data including the detected temperature data and the position data of the detected part in time series (step S1). Here, the facial skin temperature acquisition means 20 is an infrared thermography apparatus, and includes an infrared camera 21 and a processing unit 22 as shown in FIG. The infrared camera 21 is for detecting infrared radiant energy emitted from an individual's face. Here, it is assumed that the infrared camera 21 detects infrared radiation energy from the entire face of the individual. The processing unit 22 converts the infrared radiant energy detected by the infrared camera 21 into temperature to obtain temperature data, and the temperature distribution of the facial skin temperature over the entire face using the position where the infrared radiant energy is detected as position data (coordinate data). A diagram is created, and the created temperature distribution diagram is processed as facial skin temperature data corresponding to the temperature conversion data. The facial skin temperature data corresponding to the temperature conversion data is accumulated in a storage unit (not shown) included in the processing unit 22.

ここでは、処理部22において、顔面全体における顔面皮膚温度の温度分布図が作成されているが、これに限定されず、少なくとも副鼻腔周辺及び/又は前額部を含む顔面皮膚温度の温度分布図が作成され、これが温度換算データに応じた顔面皮膚温度データとされてもよい。   Here, the temperature distribution chart of the facial skin temperature in the entire face is created in the processing unit 22, but the present invention is not limited to this, and the temperature distribution chart of the facial skin temperature including at least the periphery of the sinuses and / or the forehead portion. May be created, and this may be the facial skin temperature data corresponding to the temperature conversion data.

また、ここでは、顔面皮膚温度取得手段20により温度換算データに応じた顔面皮膚温度データが取得されている間に、個人に対して脳機能賦活課題が一定期間与えられる。すなわち、顔面皮膚温度取得手段20により取得される温度換算データに応じた顔面皮膚温度データには、個人に対して脳機能賦活課題が与えられている期間のデータが含まれていることになる。なお、個人に対して与えられる脳機能賦活課題としては、脳が賦活状態になると推定されるものであれば特に限定されるものではなく、例えば、脳活動可視化装置10の利用目的に応じてその内容が適宜決定されるよう構成されていてもよい。   Also, here, while the facial skin temperature data corresponding to the temperature conversion data is acquired by the facial skin temperature acquisition means 20, a brain function activation task is given to the individual for a certain period. That is, the facial skin temperature data corresponding to the temperature conversion data acquired by the facial skin temperature acquisition means 20 includes data for a period during which a brain function activation task is given to an individual. The brain function activation task given to an individual is not particularly limited as long as it is estimated that the brain is activated, for example, depending on the purpose of use of the brain activity visualization device 10 You may be comprised so that the content may be determined suitably.

脳活動推定手段30は、顔面皮膚温度取得手段20により取得された温度換算データに応じた顔面皮膚温度データに基づき、人間の脳活動を推定する。具体的には、脳活動推定手段30は、図19に示すように、換算部31と、解析部32と、推定部33と、を有する。   The brain activity estimation unit 30 estimates human brain activity based on the facial skin temperature data corresponding to the temperature conversion data acquired by the facial skin temperature acquisition unit 20. Specifically, as shown in FIG. 19, the brain activity estimation means 30 includes a conversion unit 31, an analysis unit 32, and an estimation unit 33.

換算部31は、温度換算データに応じた顔面皮膚温度データに含まれる温度データを相対的な温度データに換算し、換算した相対的な温度データに基づく顔面皮膚温度データすなわち相対温度換算データに応じた顔面皮膚温度データを作成する(ステップS2)。具体的には、換算部31は、所定時間毎(例えば、30秒)の温度換算データに応じた顔面皮膚温度データに含まれる温度データの平均値を基準値として、該温度データを相対的な温度データに換算する。そして、換算部31は、換算した相対的な温度データ及び位置データを利用して、相対温度換算データに応じた顔面皮膚温度データを作成する。   The conversion unit 31 converts the temperature data included in the facial skin temperature data corresponding to the temperature conversion data into relative temperature data, and responds to the facial skin temperature data based on the converted relative temperature data, that is, relative temperature conversion data. Facial skin temperature data is created (step S2). Specifically, the conversion unit 31 uses the average value of the temperature data included in the facial skin temperature data corresponding to the temperature conversion data every predetermined time (for example, 30 seconds) as a reference value, and converts the temperature data into a relative value. Convert to temperature data. And the conversion part 31 produces the face skin temperature data according to relative temperature conversion data using the converted relative temperature data and position data.

解析部32は、時系列の温度換算データに応じた顔面皮膚温度データ及び相対温度換算データに応じた顔面皮膚温度データのそれぞれを、特異値分解、主成分分析或いは独立成分分析により複数の成分に分解する(ステップS3)。ここでは、解析部32は、取得した温度換算データに応じた顔面皮膚温度データ及び換算した相対温度換算データに応じた顔面皮膚温度データのそれぞれを対象として、MATLAB(登録商標)のSVDを分析ツールとして用いて、特異値分解を行う。特異値分解は、時系列で取得した温度換算データに応じた顔面皮膚温度データ及び相対温度換算データに応じた顔面皮膚温度データについて、要因を所定期間(例えば、30秒)毎の時間データとし、測度をその期間における温度換算データに応じた顔面皮膚温度データ及び相対温度換算データに応じた顔面皮膚温度データとして行われる。そして、特異値分解により、温度換算データに応じた顔面皮膚温度データ及び相対温度換算データに応じた顔面皮膚温度データのそれぞれを複数の成分に分解し、時間分布と、空間分布と、各成分の大きさを示す特異値とを算出する。   The analysis unit 32 converts the facial skin temperature data according to the time-series temperature conversion data and the facial skin temperature data according to the relative temperature conversion data into a plurality of components by singular value decomposition, principal component analysis, or independent component analysis. Decompose (step S3). Here, the analysis unit 32 analyzes the SVLAB of MATLAB (registered trademark) for each of the facial skin temperature data according to the acquired temperature conversion data and the facial skin temperature data according to the converted relative temperature conversion data. To perform singular value decomposition. The singular value decomposition is performed with the factor of time data for a predetermined period (for example, 30 seconds) for the facial skin temperature data according to the temperature conversion data acquired in time series and the facial skin temperature data according to the relative temperature conversion data, The measure is performed as facial skin temperature data corresponding to the temperature conversion data in that period and facial skin temperature data corresponding to the relative temperature conversion data. Then, by singular value decomposition, each of the facial skin temperature data according to the temperature conversion data and the facial skin temperature data according to the relative temperature conversion data is decomposed into a plurality of components, and the time distribution, spatial distribution, and A singular value indicating the magnitude is calculated.

また、解析部32は、特異値分解によって分解した複数の成分から脳活動を反映した皮膚温度の変化を示す成分を同定するために、各成分が第1条件及び第2条件を満たすか否かを判定する(ステップS4a、ステップS4b、ステップS5a、及びステップS5b)。なお、ここでは、解析部32において、まず、温度換算データに応じた顔面皮膚温度データに基づく各成分について第1条件が満たされているか否かが判定され(ステップS4a)、ステップS4aにおいて第1条件が満たされていると判定された温度換算データに応じた顔面皮膚温度データに基づく成分について第2条件が満たされているか否かが判定される(ステップS4b)。そして、相対温度換算データに応じた顔面皮膚温度データに基づく各成分のうちステップS4a及びステップS4bにおいて第1条件及び第2条件を満たすと判定された成分と一致する成分についてのみ第1条件が満たされているか否かが判定され(ステップS5a)、その後、ステップS5aにおいて第1条件が満たされていると判定された相対温度換算データに応じた顔面皮膚温度データに基づく成分について第2条件が満たされているか否かが判定される(ステップS5b)。しかしながら、解析部32における該判定の順序はこれに限定されるものではなく、例えば、温度換算データに応じた顔面皮膚温度データに基づく各成分と、相対温度換算データに応じた顔面皮膚温度データに基づく各成分とが、第1条件及び第2条件を満たすか否かがそれぞれ判定され、判定結果の一致する成分が最終的に抽出されてもよい。   In addition, the analysis unit 32 determines whether or not each component satisfies the first condition and the second condition in order to identify a component indicating a change in skin temperature reflecting brain activity from a plurality of components decomposed by singular value decomposition. (Step S4a, Step S4b, Step S5a, and Step S5b). Here, the analysis unit 32 first determines whether or not the first condition is satisfied for each component based on the facial skin temperature data corresponding to the temperature conversion data (step S4a), and the first in step S4a. It is determined whether or not the second condition is satisfied for the component based on the facial skin temperature data corresponding to the temperature conversion data determined to satisfy the condition (step S4b). Of the components based on the facial skin temperature data corresponding to the relative temperature conversion data, the first condition is satisfied only for components that match the components determined to satisfy the first condition and the second condition in step S4a and step S4b. (Step S5a), and then the second condition is satisfied for the component based on the facial skin temperature data corresponding to the relative temperature conversion data determined to satisfy the first condition in Step S5a. It is determined whether or not it has been performed (step S5b). However, the order of the determination in the analysis unit 32 is not limited to this. For example, each component based on the facial skin temperature data corresponding to the temperature conversion data and the facial skin temperature data corresponding to the relative temperature conversion data are used. It may be determined whether or not each component based on the first condition and the second condition is satisfied, and a component having a matching determination result may be finally extracted.

第1条件とは、特異値分解によって分解した成分の成分波形の振幅が、脳の非賦活時及び脳の賦活時の変化と相関関係にある、という条件である。解析部32は、複数の成分のうち、第1条件を満たす成分を、判定用成分として抽出する。なお、ここでは、温度換算データに応じた顔面皮膚温度データを取得している間に、個人に対して脳機能賦活課題が与えられている期間が一定期間ある。解析部32は、個人に対して脳機能賦活課題が与えられていない期間を脳の非賦活時とし、個人に対して脳機能賦活課題が与えられている期間を脳の賦活時として、脳機能賦活課題が与えられている期間及び与えられていない期間と、各成分の成分波形とを比較解析する。解析部32は、成分波形データに基づく比較解析結果を利用して、各成分の成分波形と脳の非賦活時及び脳の賦活時とが相関関係にあるか否かを評価し、複数の成分のうち相関関係にあると評価した成分を、第1条件を満たす判定用成分として抽出する。一方、解析部32は、複数の成分のうち相関関係にないと評価した成分を、第1条件を満たさず人間の脳活動を反映した温度変化を示す成分ではないと判定する(ステップS6)。   The first condition is a condition that the amplitude of the component waveform of the component decomposed by the singular value decomposition is correlated with a change at the time of non-activation of the brain and activation of the brain. The analysis unit 32 extracts a component that satisfies the first condition from the plurality of components as a determination component. Here, there is a certain period during which the brain function activation task is given to an individual while acquiring facial skin temperature data corresponding to the temperature conversion data. The analysis unit 32 uses the period when the brain function activation task is not given to the individual as the non-activation time of the brain, and sets the period when the brain function activation task is given to the individual as the brain activation time. The period in which the activation task is given and the period in which it is not given and the component waveform of each component are compared and analyzed. The analysis unit 32 uses the comparison analysis result based on the component waveform data to evaluate whether the component waveform of each component and the brain non-activation time and the brain activation time are correlated, and a plurality of components Among these, the component evaluated as having a correlation is extracted as a determination component that satisfies the first condition. On the other hand, the analysis unit 32 determines that a component evaluated as having no correlation among the plurality of components is not a component that does not satisfy the first condition and shows a temperature change reflecting human brain activity (step S6).

ここでは、温度換算データに応じた顔面皮膚温度データの取得時に個人に対して脳機能賦活課題が一定期間与えられており、これに基づき解析部32は判定用成分を抽出しているが、第1条件の内容、すなわち解析部32における判定用成分の抽出手段はこれに限定されない。例えば、予め実験等がされていることで複数の成分のうち脳の非賦活時及び脳の賦活時と相関関係にある成分波形を示す成分が特定されている場合には、解析部32は、複数の成分から特定されている該成分を判定用成分として抽出する。また、本脳活動可視化装置において眼球運動又はまたたき等の脳の賦活/非賦活に関連することが知られている人間の動作が検出される場合には、解析部32が、この検出結果と各成分の成分波形とを比較解析及び評価することで、複数の成分から判定用成分を抽出してもよい。なお、解析部32による第1条件を満たすか否かの判定の基準は、脳活動可視化装置10の利用目的等に応じて、シミュレーションや実験、机上計算等によって適宜決定される。   Here, a brain function activation task is given to an individual for a certain period when acquiring facial skin temperature data corresponding to temperature conversion data, and based on this, the analysis unit 32 extracts a determination component. The contents of one condition, that is, the determination component extraction means in the analysis unit 32 is not limited to this. For example, when a component showing a component waveform correlated with the non-activation time and the activation time of the brain among a plurality of components by performing an experiment or the like in advance, the analysis unit 32 The component specified from the plurality of components is extracted as a determination component. Further, in the case where the human activity known to be related to activation / deactivation of the brain such as eye movement or tapping is detected in the present brain activity visualization device, the analysis unit 32 analyzes the detection result and each A determination component may be extracted from a plurality of components by comparative analysis and evaluation of the component waveform of the component. Note that the criterion for determining whether or not the first condition is satisfied by the analysis unit 32 is appropriately determined by simulation, experiment, desktop calculation, or the like according to the purpose of use of the brain activity visualization device 10 or the like.

第2条件は、抽出した判定用成分において、人間の顔面の所定部位における温度変化がある、という条件である。解析部32は、判定用成分のうち、第2条件を満たす成分を、人間の脳活動に関連している可能性の高い成分と判定し、候補成分として抽出する。すなわち、解析部32は、人間の顔面の所定部位における温度変化の有無に基づき、判定用成分が人間の脳活動に関連しているか否かを判定する。具体的には、解析部32は、抽出した判定用成分の温度分布データに基づき、副鼻腔周辺及び/又は前額部において温度変化が生じているか否かを判定し、温度変化が生じている場合には該判定用成分が第2条件を満たす人間の脳活動に関連する可能性の高い成分であると判定し、候補成分として抽出する。一方で、解析部32は、副鼻腔周辺及び/又は前額部において温度変化が生じていない場合には、該判定用成分は第2条件を満たさず脳活動を反映した皮膚温度の変化を示す成分ではない、と判定する(ステップS6)。なお、解析部32による第2条件を満たすか否かの判定の基準は、脳活動可視化装置10の利用目的等に応じて、シミュレーションや実験、机上計算等によって適宜決定される。   The second condition is a condition that there is a temperature change in a predetermined part of the human face in the extracted determination component. The analysis unit 32 determines a component that satisfies the second condition among the determination components as a component that is highly likely to be related to human brain activity, and extracts it as a candidate component. That is, the analysis unit 32 determines whether or not the determination component is related to human brain activity based on the presence or absence of a temperature change at a predetermined part of the human face. Specifically, the analysis unit 32 determines whether a temperature change has occurred around the sinuses and / or the forehead based on the extracted temperature distribution data of the determination component, and the temperature change has occurred. In this case, it is determined that the determination component is a component that is highly likely to be related to human brain activity that satisfies the second condition, and is extracted as a candidate component. On the other hand, when no temperature change occurs around the sinuses and / or the forehead part, the analysis unit 32 does not satisfy the second condition and shows a change in skin temperature reflecting brain activity. It is determined that it is not a component (step S6). Note that the criterion for determining whether or not the second condition is satisfied by the analysis unit 32 is appropriately determined by simulation, experiment, desktop calculation, or the like according to the purpose of use of the brain activity visualization device 10 or the like.

そして、解析部32は、ステップS5bにおいて第2条件を満たすと判定した成分を、脳活動を反映した皮膚温度の変化を示す成分として同定する(ステップS7)。すなわち、ステップS7において脳活動を反映した皮膚温度の変化を示す成分として同定される成分は、温度換算データに応じた顔面皮膚温度データを特異値分解により分解し解析することで抽出された候補成分と、相対温度換算データに応じた顔面皮膚温度データを特異値分解により分解し解析することで抽出された候補成分と、の間で一致している成分ということになる。なお、両解析で一致していない候補成分については、ステップS6において脳活動を反映した皮膚温度の変化を示す成分ではない、と判定されている。   And the analysis part 32 identifies the component determined to satisfy | fill 2nd conditions in step S5b as a component which shows the change of the skin temperature reflecting brain activity (step S7). That is, the component identified as the component indicating the change in skin temperature reflecting the brain activity in step S7 is the candidate component extracted by decomposing and analyzing the facial skin temperature data corresponding to the temperature conversion data by singular value decomposition And the candidate skin component extracted by analyzing the facial skin temperature data corresponding to the relative temperature conversion data by decomposing and analyzing by singular value decomposition. Note that candidate components that do not match in both analyzes are determined not to be components indicating changes in skin temperature reflecting brain activity in step S6.

推定部33は、解析部32において人間の脳活動を反映した皮膚温度の変化を示す成分として同定された成分に基づいて、人間の脳活動を推定する。具体的には、推定部33は、解析部32において同定された成分の成分波形データに基づいて、顔面皮膚温度データの取得時における脳活動量を推定する。   The estimation unit 33 estimates the human brain activity based on the component identified by the analysis unit 32 as the component indicating the change in skin temperature reflecting the human brain activity. Specifically, the estimation unit 33 estimates the amount of brain activity when acquiring facial skin temperature data based on the component waveform data of the component identified by the analysis unit 32.

(4−1−1)変形例1A
上記脳活動推定手段30は換算部31を有しており、換算部31によって相対温度換算データに応じた顔面皮膚温度データが作成されている。そして、解析部32が、顔面皮膚温度取得手段20により取得された温度換算データに応じた顔面皮膚温度データだけでなく、相対的な温度データに換算された温度データに基づく相対温度データに応じた顔面皮膚温度データについても、特異値分解により複数の成分に分解し、各成分についての解析を行っている。
(4-1-1) Modification 1A
The brain activity estimating means 30 has a conversion unit 31, and facial skin temperature data corresponding to the relative temperature conversion data is created by the conversion unit 31. And the analysis part 32 respond | corresponded not only the face skin temperature data according to the temperature conversion data acquired by the face skin temperature acquisition means 20, but the relative temperature data based on the temperature data converted into relative temperature data. The facial skin temperature data is also decomposed into a plurality of components by singular value decomposition, and each component is analyzed.

これに代えて、脳活動推定手段30が換算部31を有していなくてもよい。この場合、相対温度換算データに応じた顔面皮膚温度データを作成したり、相対温度換算データに応じた顔面皮膚温度データに基づくデータの解析を行ったりする処理を省くことができる。   Instead of this, the brain activity estimating means 30 may not have the conversion unit 31. In this case, it is possible to omit processing for creating facial skin temperature data corresponding to the relative temperature conversion data or analyzing data based on the facial skin temperature data corresponding to the relative temperature conversion data.

ただし、人間の脳活動に関連する成分を精度よく同定するためには、上記実施形態のように脳活動推定手段30が換算部31を有しており、解析部32によって、顔面皮膚温度取得手段20により取得された温度換算データに応じた顔面皮膚温度データだけでなく、相対的な温度データに換算された温度データに基づく相対温度データに応じた顔面皮膚温度データについても、特異値分解により複数の成分に分解され、各成分についての解析が行われるほうが望ましい。   However, in order to accurately identify a component related to human brain activity, the brain activity estimation means 30 has a conversion unit 31 as in the above embodiment, and the analysis unit 32 performs facial skin temperature acquisition means. A plurality of facial skin temperature data corresponding to the relative temperature data based on the temperature data converted to the relative temperature data as well as the facial skin temperature data corresponding to the temperature converted data acquired by 20 are obtained by singular value decomposition. It is desirable that each component is analyzed after being decomposed into these components.

(4−1−2)変形例1B
また、上記顔面皮膚温度取得手段20は、対象物と非接触の状態で温度データを取得することができる赤外線サーモグラフィ装置である。
(4-1-2) Modification 1B
The facial skin temperature acquisition means 20 is an infrared thermography device that can acquire temperature data in a non-contact state with an object.

しかしながら、個人の顔面の少なくとも一部の皮膚温度を検出し、検出した温度データ及びその検出部位の位置データを含む顔面皮膚温度データを時系列で取得することができれば、顔面皮膚温度取得手段は赤外線サーモグラフィ装置に限定されない。   However, if the skin temperature of at least a part of an individual's face is detected and the facial skin temperature data including the detected temperature data and the position data of the detected part can be acquired in time series, the facial skin temperature acquisition means is infrared. It is not limited to a thermography device.

例えば、顔面皮膚温度取得手段が温度センサを含む装置であってもよい。具体的には、個人の顔面の所定部位に温度センサを装着し、温度センサによって検出される温度データと、温度センサを装着した部位の位置データとに基づいて、時系列の顔面皮膚温度データが取得されてもよい。このように、温度センサにより対象となる個人に接触した状態で顔面皮膚温度データが取得される場合であっても、温度センサは脳波電極等のように装着前の処理が必要ではないため、脳波計測法、磁気共鳴画像法、及び近赤外線分光法等の従来の検出方法と比較して、簡便にデータを取得することができる。これにより、簡便に人間の脳活動を推定することができる。   For example, the facial skin temperature acquisition means may be a device including a temperature sensor. Specifically, a temperature sensor is attached to a predetermined part of an individual's face, and time-series facial skin temperature data is based on temperature data detected by the temperature sensor and position data of the part where the temperature sensor is attached. May be acquired. As described above, even when facial skin temperature data is acquired in a state in which the temperature sensor is in contact with the target individual, the temperature sensor does not require pre-wearing processing like an electroencephalogram electrode. Compared with conventional detection methods such as measurement, magnetic resonance imaging, and near infrared spectroscopy, data can be easily acquired. Thereby, human brain activity can be estimated simply.

(4−2)顔面の撮影画像データに基づき脳活動を推定する脳活動推定手段130
図21は、本発明の実施形態に係る脳活動可視化装置110の概略図である。図22は、脳活動可視化装置110において脳機能を反映した顔面のRGB変化を示す成分を同定する際の処理の流れの一例を示すフローチャートである。
(4-2) Brain activity estimating means 130 for estimating brain activity based on face image data.
FIG. 21 is a schematic diagram of the brain activity visualization apparatus 110 according to the embodiment of the present invention. FIG. 22 is a flowchart illustrating an example of a processing flow when the component indicating the RGB change of the face reflecting the brain function is identified in the brain activity visualization apparatus 110.

脳活動可視化装置110の備える脳活動推定手段130は、個人(被験者)の顔面の撮影画像データから、個人の脳活動を推定するための装置である。脳活動可視化装置110は、図21に示すように、画像データ取得手段120と、脳活動推定手段130と、状態可視化手段200と、を備える。   The brain activity estimating means 130 included in the brain activity visualization device 110 is a device for estimating the brain activity of the individual from the captured image data of the face of the individual (subject). As shown in FIG. 21, the brain activity visualization apparatus 110 includes image data acquisition means 120, brain activity estimation means 130, and state visualization means 200.

画像データ取得手段120は、個人の顔面の少なくとも一部の撮影画像データを時系列で取得する(ステップS101)。なお、画像データ取得手段120は、少なくとも撮影装置を有するものであれば特に限定されるものではなく、例えば、スマートフォンやタブレット(例えば、iPad:登録商標)等の撮影装置内蔵型ポータブル端末等が挙げられる。ここでは、画像データ取得手段120は、図21に示すように、撮影装置としてのカメラ121と、記憶部122とを有する。カメラ121は、個人の顔面の撮影画像データを時系列で取得するためのものである。ここでは、カメラ121は、個人の顔面全体の動画を撮影し、撮影した動画データを取得する。記憶部122は、撮影装置により撮影された時系列の撮影画像データを蓄積する。ここでは、記憶部122は、カメラ121によって取得された動画データを蓄積する。   The image data acquisition unit 120 acquires captured image data of at least a part of an individual's face in time series (step S101). The image data acquisition unit 120 is not particularly limited as long as it has at least a photographing device, and examples thereof include a portable terminal with a photographing device such as a smartphone or a tablet (for example, iPad: registered trademark). It is done. Here, the image data acquisition unit 120 includes a camera 121 as a photographing apparatus and a storage unit 122, as shown in FIG. The camera 121 is for acquiring captured image data of an individual's face in time series. Here, the camera 121 captures a moving image of the entire face of the individual and acquires the captured moving image data. The storage unit 122 accumulates time-series photographed image data photographed by the photographing apparatus. Here, the storage unit 122 accumulates moving image data acquired by the camera 121.

なお、ここでは、カメラ121によって顔面全体の動画が撮影されているが、これに限定されず、顔面の少なくとも前額部及び/又は副鼻腔周辺の画像を含む動画が撮影されていればよい。   Here, the moving image of the entire face is captured by the camera 121, but the present invention is not limited to this, and it is only necessary to capture a moving image including an image of at least the forehead portion and / or the sinuses around the face.

また、ここでは、画像データ取得手段120により顔面の時系列の撮影画像データが取得されている間に、個人に対して脳機能賦活課題が一定期間与えられる。すなわち、画像データ取得手段120により取得される撮影画像データには、個人に対して脳機能賦活課題が与えられている期間のデータが含まれていることになる。なお、個人に対して与えられる脳機能賦活課題としては、脳が賦活状態になると推定されるものであれば特に限定されるものではなく、例えば、脳活動可視化装置110の利用目的に応じてその内容が適宜決定されるよう構成されていてもよい。   Also, here, while the time-series captured image data of the face is acquired by the image data acquisition means 120, a brain function activation task is given to the individual for a certain period. That is, the captured image data acquired by the image data acquisition unit 120 includes data for a period in which a brain function activation task is given to an individual. The brain function activation task given to an individual is not particularly limited as long as it is estimated that the brain is activated, for example, depending on the purpose of use of the brain activity visualization device 110 You may be comprised so that the content may be determined suitably.

脳活動推定手段130は、画像データ取得手段120により取得された顔面の時系列の撮影画像データに基づき、人間の脳活動を推定する。具体的には、脳活動推定手段130は、図21に示すように、RGB処理部131と、換算部132と、血行量算出部133と、解析部134と、推定部135と、を有する。なお、図21では、脳活動推定手段130が、RGB処理部131、換算部132、血行量算出部133、解析部134及び推定部135を有する1つの装置として存在している態様が示されているが、本発明はこれに限定されるものではなく、RGB処理部131、換算部132、血行量算出部133、解析部134及び推定部135の一部或いはそれぞれが独立した装置として存在していてもよい。また、ここでは、画像データ取得手段120、RGB処理部131、換算部132、及び血行量算出部133により顔面血行量取得手段が構成されている。   The brain activity estimation unit 130 estimates human brain activity based on time-series captured image data of the face acquired by the image data acquisition unit 120. Specifically, as shown in FIG. 21, the brain activity estimation unit 130 includes an RGB processing unit 131, a conversion unit 132, a blood flow rate calculation unit 133, an analysis unit 134, and an estimation unit 135. FIG. 21 shows an aspect in which the brain activity estimation means 130 exists as one device having the RGB processing unit 131, the conversion unit 132, the blood flow calculation unit 133, the analysis unit 134, and the estimation unit 135. However, the present invention is not limited to this, and part or each of the RGB processing unit 131, the conversion unit 132, the blood flow volume calculation unit 133, the analysis unit 134, and the estimation unit 135 exists as an independent device. May be. In addition, here, the facial blood flow volume acquisition unit is configured by the image data acquisition unit 120, the RGB processing unit 131, the conversion unit 132, and the blood flow volume calculation unit 133.

RGB処理部131は、画像データ取得手段120により取得された撮影画像データに対して、R成分、G成分及びB成分の3つの色成分に分解するRGB処理を行う(ステップS102)。ここで、顔面全体の撮影画像データに対してRGB処理を行ってもよいが、ここでは、演算処理量及びノイズを減らすために、撮影画像データから前額部及び/又は副鼻腔周辺のデータを抽出し、抽出したデータについてのみRGB処理を行うものとする。   The RGB processing unit 131 performs RGB processing that decomposes the captured image data acquired by the image data acquisition unit 120 into three color components, an R component, a G component, and a B component (step S102). Here, the RGB processing may be performed on the captured image data of the entire face, but here, in order to reduce the calculation processing amount and noise, the data on the forehead part and / or the sinuses around the captured image data is reduced. It is assumed that RGB processing is performed only on the extracted data.

換算部132は、RGB処理により得られた撮影画像データのRGBデータを相対的なRGBデータに換算する(ステップS103)。具体的には、換算部132は、取得された所定時間毎(例えば、30秒)の撮影画像データから得られるRGBデータの平均値を基準値として、該RGBデータを相対的なRGBデータに換算する。   The conversion unit 132 converts the RGB data of the captured image data obtained by the RGB processing into relative RGB data (step S103). Specifically, the conversion unit 132 converts the RGB data into relative RGB data using the average value of the RGB data obtained from the acquired captured image data every predetermined time (for example, 30 seconds) as a reference value. To do.

血行量算出部133は、RGB処理により得られた撮影画像データのRGBデータに基づき、顔面の時系列の血行量データを算出する(ステップS104)。   The blood flow volume calculation unit 133 calculates time-series blood flow volume data of the face based on the RGB data of the captured image data obtained by the RGB processing (step S104).

解析部134は、時系列の相対換算血行量データを、特異値分解、主成分分析或いは独立成分分析により複数の成分に分解する(ステップS105)。ここでは、解析部134は、相対換算血行量データに対して、MATLAB(登録商標)のSVDを分析ツールとして用いて、特異値分解を行う。具体的には、特異値分解は、時系列の相対換算血行量データを対象として、要因を所定期間(例えば、30秒)毎の時間データとし、測度をその期間毎における相対的なRGBデータから演算したピクセル毎の相対換算血行量データとして行われる。そして、特異値分解により、時系列の相対換算血行量データを複数の成分に分解し、時間分布と、空間分布と、各成分の大きさを示す特異値とを算出する。   The analysis unit 134 decomposes the time-series relative converted blood flow data into a plurality of components by singular value decomposition, principal component analysis, or independent component analysis (step S105). Here, the analysis unit 134 performs singular value decomposition on the relative converted blood flow data using the SVLAB of MATLAB (registered trademark) as an analysis tool. Specifically, the singular value decomposition is performed on time-series relative converted blood flow data, with the factor being time data for each predetermined period (for example, 30 seconds), and the measure is calculated from relative RGB data for each period. This is performed as the calculated relative blood flow volume data for each pixel. Then, by singular value decomposition, the time-series relative converted blood flow data is decomposed into a plurality of components, and a time distribution, a spatial distribution, and a singular value indicating the size of each component are calculated.

また、解析部134は、特異値分解によって分解した複数の成分から脳活動を反映した顔面のRGB変化を示す成分を同定するために、各成分が所定条件を満たすか否かを判定する(ステップS106)。ここで、所定条件としては、例えば、特異値分解によって分解した成分の成分波形の振幅が、脳の非賦活時及び脳の賦活時の変化と相関関係にあるという条件(以下、第1条件という)や、特異値分解によって分解した成分において人間の顔面の所定部位に血行量変化があるという条件(以下、第2条件という)等が含まれる。解析部134において判定される所定条件としては、1又は複数の条件が設定されていればよく、ここでは、所定条件として第1条件が設定されているものとする。   Further, the analysis unit 134 determines whether or not each component satisfies a predetermined condition in order to identify a component indicating the RGB change of the face reflecting the brain activity from the plurality of components decomposed by the singular value decomposition (step) S106). Here, as the predetermined condition, for example, the condition that the amplitude of the component waveform of the component decomposed by the singular value decomposition is correlated with the change at the time of non-activation of the brain and the activation of the brain (hereinafter referred to as the first condition). ), A condition that there is a change in blood flow in a predetermined part of the human face in a component decomposed by singular value decomposition (hereinafter referred to as a second condition), and the like. As the predetermined condition determined by the analysis unit 134, one or a plurality of conditions may be set. Here, it is assumed that the first condition is set as the predetermined condition.

そして、解析部134は、複数の成分のうち所定条件を満たす成分を、判定用成分として抽出する。さらに、解析部134は、抽出した判定用成分のうち所定条件に含まれる全ての条件を満たす成分を、脳活動を反映した顔面のRGB変化を示す成分として同定する(ステップS107)。一方、解析部134は、複数の成分のうち所定条件に含まれる少なくとも1つの条件を満たさないと判定した成分を、脳活動を反映した顔面のRGB変化を示す成分ではないと判定する(ステップS108)。   And the analysis part 134 extracts the component which satisfy | fills predetermined conditions among several components as a component for determination. Further, the analysis unit 134 identifies a component satisfying all the conditions included in the predetermined condition among the extracted determination components as a component indicating the RGB change of the face reflecting the brain activity (step S107). On the other hand, the analysis unit 134 determines that the component determined not to satisfy at least one condition included in the predetermined condition among the plurality of components is not a component indicating the RGB change of the face reflecting the brain activity (step S108). ).

ここでは、上述のように所定条件として1つの条件(第1条件)のみが設定されており、顔面の時系列の撮影画像データを取得している間に、個人に対して脳機能賦活課題が与えられている期間が一定期間ある。このため、解析部134は、個人に対して脳機能賦活課題が与えられていない期間を脳の非賦活時とし、個人に対して脳機能賦活課題が与えられている期間を脳の賦活時として、脳機能賦活課題が与えられている期間及び与えられていない期間と、各成分の成分波形とを比較解析する。そして、解析部134は、成分波形データに基づく比較解析結果を利用して、各成分の成分波形と脳の非賦活時及び脳の賦活時とが相関関係にあるか否かを評価し、複数の成分のうち相関関係にあると評価した成分を、所定条件を満たす判定用成分として抽出すると共に、脳活動を反映した顔面のRGB変化を示す成分として同定する。一方、解析部134は、複数の成分のうち相関関係にないと評価した成分を、所定条件を満たさず人間の脳活動を反映した顔面のRGB変化を示す成分ではないと判定する。   Here, only one condition (first condition) is set as the predetermined condition as described above, and the brain function activation task is given to the individual while acquiring time-series captured image data of the face. There is a given period of time. For this reason, the analysis unit 134 sets the period when the brain function activation task is not given to the individual as the non-activation time of the brain, and sets the period when the brain function activation task is given to the individual as the brain activation time The period in which the brain function activation task is given and the period in which the brain function activation task is not given, and the component waveform of each component are compared and analyzed. And the analysis part 134 evaluates whether the component waveform of each component and the time of non-activation of a brain and the time of activation of a brain have a correlation using the comparative analysis result based on component waveform data, Among these components, a component evaluated as having a correlation is extracted as a determination component that satisfies a predetermined condition, and is identified as a component indicating a change in RGB of the face reflecting brain activity. On the other hand, the analysis unit 134 determines that a component evaluated as having no correlation among the plurality of components is not a component that does not satisfy the predetermined condition and indicates a change in RGB of the face that reflects human brain activity.

ここでは、顔面の時系列の撮影画像データが取得される際に個人に対して脳機能賦活課題が一定期間与えられており、これに基づき解析部134が判定用成分を抽出しているが、第1条件の内容、すなわち解析部134における判定用成分の抽出手段はこれに限定されない。例えば、予め実験等がされていることで複数の成分のうち脳の非賦活時及び脳の賦活時と相関関係にある成分波形を示す成分が特定されている場合には、解析部134は、複数の成分から特定されている該成分を判定用成分として抽出する。また、脳活動可視化装置110において眼球運動又はまたたき等の脳の賦活/非賦活に関連することが知られている人間の動作についても検出される場合には、解析部134が、この検出結果と各成分の成分波形とを比較解析及び評価することで、複数の成分から判定用成分を抽出してもよい。なお、解析部134による第1条件を満たすか否かの判定の基準は、脳活動可視化装置110の利用目的等に応じて、シミュレーションや実験、机上計算等によって適宜決定される。   Here, a brain function activation task is given to an individual for a certain period when time-series captured image data of the face is acquired, and the analysis unit 134 extracts a determination component based on this, The content of the first condition, that is, the determination component extraction means in the analysis unit 134 is not limited to this. For example, when a component indicating a component waveform that is correlated with the non-activation time and the activation time of the brain among a plurality of components is identified by performing an experiment or the like in advance, the analysis unit 134 The component specified from the plurality of components is extracted as a determination component. In addition, when the human activity known to be related to activation / deactivation of the brain, such as eye movement or tapping, is also detected in the brain activity visualization device 110, the analysis unit 134 A determination component may be extracted from a plurality of components by comparative analysis and evaluation of the component waveform of each component. Note that the criterion for determining whether or not the first condition is satisfied by the analysis unit 134 is appropriately determined by simulation, experiment, desktop calculation, or the like according to the purpose of use of the brain activity visualization device 110 or the like.

また、所定条件として第2条件が設定されている場合には、解析部134は、人間の顔面の所定部位における顔面の血行量変化の有無に基づき、判定用成分を抽出する。具体的には、解析部134は、特異値分解によって分解された複数の成分に応じた血行量分布図に基づき、副鼻腔周辺及び/又は前額部において血行量の変化が生じているか否かを判定し、血行量の変化が生じている場合には該成分が第2条件を満たしていると判定する。一方で、副鼻腔周辺及び/又は前額部において血行量の変化が生じていない場合には、解析部134は、該成分が第2条件を満たしていないと判定する。なお、解析部134による第2条件を満たすか否かの判定の基準は、脳活動可視化装置110の利用目的等に応じて、シミュレーションや実験、机上計算等によって適宜決定されるものとする。   When the second condition is set as the predetermined condition, the analysis unit 134 extracts a determination component based on whether or not there is a change in the blood flow rate of the face at a predetermined part of the human face. Specifically, the analysis unit 134 determines whether or not there is a change in blood flow around the sinuses and / or the forehead based on the blood flow distribution map corresponding to a plurality of components decomposed by the singular value decomposition. If a change in blood flow occurs, it is determined that the component satisfies the second condition. On the other hand, when there is no change in the blood flow volume around the sinuses and / or the forehead, the analysis unit 134 determines that the component does not satisfy the second condition. Note that the criterion for determining whether or not the second condition is satisfied by the analysis unit 134 is appropriately determined by simulation, experiment, desktop calculation, or the like according to the purpose of use of the brain activity visualization device 110 or the like.

さらに、血行量算出部133によって相対的なRGBデータに換算される前のRGBデータに基づく時系列の血行量データが算出される場合には、解析部134によって、該血行量データを特異値分解等することで得られた複数の成分についても、上記第1条件及び/又は第2条件が満たされるか否かが判定され、判定用成分が抽出されてもよい。   Further, when time-series blood flow data based on RGB data before being converted into relative RGB data is calculated by the blood flow calculation unit 133, the blood flow data is subjected to singular value decomposition by the analysis unit 134. With respect to a plurality of components obtained by performing the same, it may be determined whether or not the first condition and / or the second condition is satisfied, and a determination component may be extracted.

推定部135は、解析部134において人間の脳活動を反映した顔面のRGB変化を示す成分として同定された成分に基づいて、人間の脳活動を推定する。具体的には、推定部135は、解析部134において同定された成分の成分波形データに基づいて、顔面の撮影画像データの取得時における脳活動量を推定する。   The estimation unit 135 estimates the human brain activity based on the component identified by the analysis unit 134 as the component indicating the RGB change of the face reflecting the human brain activity. Specifically, the estimation unit 135 estimates the amount of brain activity when acquiring facial image data based on the component waveform data of the component identified by the analysis unit 134.

(4−2−1)変形例2A
なお、上記血行量算出部133において、RGBデータに含まれる各画素のうちの主にR成分を用いて顔面の血行量データが算出されてもよい。また、RGBデータに基づき血行量データを算出できるのであれば、血行量データは必ずしも紅斑指数に限定されるものではない。
(4-2-1) Modification 2A
The blood flow rate calculation unit 133 may calculate the blood flow rate data of the face mainly using the R component of each pixel included in the RGB data. In addition, blood flow volume data is not necessarily limited to the erythema index as long as blood flow volume data can be calculated based on RGB data.

(4−2−2)変形例2B
上記血行量算出部133は、換算部132によって換算された相対的なRGBデータに基づき相対換算血行量データを算出するが、これに代えて或いはこれに加えて、相対的なRGBデータに換算される前のRGBデータに基づき血行量データが算出されてもよい。ここで、相対的なRGBデータに換算される前のRGBデータに基づき算出された血行量データには、脳活動と相関する成分が出やすい(検定力が高い)ため、例えば、相対的なRGBデータに換算される前のRGBデータに基づき算出された血行量データを、相対的なRGBデータに基づき算出された相対換算血行量データよりも先行して解析してもよい。また、例えば、まず、血行量データを解析して有意な相関のある成分を抽出し、相対換算血行量データに関しては、前記抽出した成分に対応する成分のみを解析することで、演算処理量を減らすことができる。
(4-2-2) Modification 2B
The blood flow volume calculation unit 133 calculates relative converted blood flow volume data based on the relative RGB data converted by the conversion unit 132, but instead of or in addition to this, the blood flow volume calculation unit 133 converts the blood flow volume data into relative RGB data. The blood flow volume data may be calculated based on the RGB data before the check. Here, since the blood flow volume data calculated based on the RGB data before being converted into the relative RGB data easily includes a component correlated with the brain activity (the test power is high), for example, the relative RGB data Blood flow volume data calculated based on RGB data before conversion into data may be analyzed prior to relative converted blood flow volume data calculated based on relative RGB data. In addition, for example, first, the blood flow data is analyzed to extract a component having a significant correlation, and the relative conversion blood flow data is calculated by analyzing only the component corresponding to the extracted component. Can be reduced.

(4−3)状態可視化手段200
状態可視化手段200は、脳活動推定手段30及び/又は脳活動推定手段130により推定された対象者の脳活動に基づき、対象者の生理状態を表示することにより可視化する。例えば、状態可視化手段200が、対象者の脳活動量の変化を解析することで、対象者の生理状態を解析する解析部201を有していてもよい。具体的には、解析部201が、対象者に対して与えられた刺激(視覚刺激、聴覚刺激、触覚刺激、臭覚刺激或いは味覚刺激等)に対する脳活動量の変化を解析することで、対象者の生理状態を判定する。なお、生理状態の種類やレベルについては、脳活動量の上昇度合い及び/又は持続時間に基づき、脳活動可視化装置10,110の用途に応じて適宜設置可能になっていてもよい。そして、解析部201により解析された対象者の生理状態を状態可視化手段200の表示部202から管理者へと出力されることで、管理者は対象者の生理状態を知ることができる。表示部202としては、画像やメッセージを表示する表示デバイス等、解析した対象者の生理状態に関する情報を管理者に対して可視化できるものであればどのようなものであっても採用することができる。
(4-3) State visualization means 200
The state visualizing means 200 is visualized by displaying the physiological state of the subject based on the brain activity of the subject estimated by the brain activity estimating means 30 and / or the brain activity estimating means 130. For example, the state visualization unit 200 may include an analysis unit 201 that analyzes the physiological state of the subject by analyzing changes in the amount of brain activity of the subject. Specifically, the analysis unit 201 analyzes the change in the amount of brain activity with respect to the stimulus (visual stimulus, auditory stimulus, tactile stimulus, olfactory stimulus, taste stimulus, etc.) given to the subject, so that the subject Physiological state is determined. In addition, about the kind and level of a physiological state, based on the raise degree and / or duration of brain activity amount, according to the use of the brain activity visualization apparatuses 10 and 110, you may be able to install suitably. Then, the physiological state of the subject analyzed by the analysis unit 201 is output from the display unit 202 of the state visualization means 200 to the administrator, so that the administrator can know the physiological state of the subject. As the display unit 202, any display device that displays an image or a message can be adopted as long as it can visualize information on the analyzed physiological state of the subject person to the administrator. .

また、解析部32,134において脳活動を反映する成分が同定された後に、さらに顔面皮膚温度取得手段20及び/又は画像データ取得手段120により時系列の各種データが取得される場合には、脳活動可視化装置10,110において、さらに取得された各種データが特異値分解により複数の成分に分解され、同定された成分のみが解析されることで、対象者の生理状態をリアルタイムで知ることができる。   In addition, when various components of time series are acquired by the facial skin temperature acquisition unit 20 and / or the image data acquisition unit 120 after the components reflecting the brain activity are identified in the analysis units 32 and 134, In the activity visualization apparatuses 10 and 110, the acquired various data are further decomposed into a plurality of components by singular value decomposition, and only the identified components are analyzed, so that the physiological state of the subject can be known in real time. .

さらに、被験者の顔面の皮膚温度や撮影した画像から被験者の心拍情報や生体情報等を取得する技術が従来よりあるが、顔面皮膚温度取得手段20及び/又は画像データ取得手段120から得られた各種データが特異値分解等されることで得られる成分に対して従来の技術を採用することで、心拍情報や生体情報を精度良く取得することができる。したがって、特異値分解した複数の成分を解析して心拍情報や生体情報を取得する機能を、解析部32及び/又は解析部134に持たせ、取得した心拍情報や生体情報に基づき交換神経/副交感神経の働きを推定する機能を上記実施形態の推定部33,135に持たせてもよい。   Furthermore, there is a conventional technique for acquiring the subject's facial skin temperature or captured image from heartbeat information, biological information, or the like. However, various techniques obtained from the facial skin temperature acquisition means 20 and / or the image data acquisition means 120 are available. By adopting a conventional technique for a component obtained by performing singular value decomposition or the like on the data, heart rate information and biological information can be obtained with high accuracy. Therefore, the analysis unit 32 and / or the analysis unit 134 have a function of analyzing a plurality of components subjected to singular value decomposition to acquire heart rate information and biological information, and based on the acquired heart rate information and biological information, the exchange nerve / parasympathy You may give the estimation part 33,135 of the said embodiment the function which estimates the function of a nerve.

(5)特徴
(5−1)
本実施形態では、顔面皮膚温度取得手段20及び/又は画像データ取得手段120によって取得された時系列の顔面皮膚温度データ及び/又は顔面血行量データに基づき人間の脳活動が推定される。このため、脳波電極等の装着前に処理が必要なセンサを装着しなくても、人間の脳活動を推定することができる。したがって、簡便に人間の脳活動を推定し、推定した脳活動に基づき対象者の生理状態を可視化することができている。
(5) Features (5-1)
In the present embodiment, human brain activity is estimated based on time-series facial skin temperature data and / or facial blood flow data acquired by the facial skin temperature acquisition means 20 and / or the image data acquisition means 120. For this reason, human brain activity can be estimated without mounting a sensor that requires processing before mounting an electroencephalogram electrode or the like. Therefore, it is possible to easily estimate human brain activity and visualize the physiological state of the subject based on the estimated brain activity.

(5−2)
ここで、時系列の顔面の皮膚温度データ及び/又は画像データが取得される際に、人間に対して実際に脳機能賦活課題が与えられたり与えられなかったりすることにより、人間の脳が賦活化したり賦活化しなかったりする状況が作られている場合、各成分の成分波形と脳の賦活時及び非賦活時との間に相関関係のある成分は、脳活動を反映した皮膚温度及び/又は血行量の変化を示す成分である可能性が高い成分であるといえる。
(5-2)
Here, when time-series facial skin temperature data and / or image data is acquired, a human brain function is activated by whether or not a human brain function activation task is actually given. If the situation is activated or not activated, the component having a correlation between the component waveform of each component and the activation and non-activation of the brain is the skin temperature reflecting the brain activity and / or It can be said that it is a component with high possibility of being a component which shows the change of blood flow volume.

本実施形態では、顔面皮膚温度取得手段20及び/又は画像データ取得手段120により時系列の顔面の皮膚温度データ及び/又は画像データが取得されている間に、個人に対して脳機能賦活課題が一定期間与えられている。すなわち、本実施形態では、個人に対して実際に脳機能賦活課題を与えたり与えなかったりすることにより、人間の脳が賦活化したり賦活化しなかったりする状況が作られている。そして、このように取得された時系列の各種データが特異値分解により複数の成分に分解され、各成分についてその成分波形と脳の賦活時及び非賦活時との相関関係が評価され、相関関係にある成分が判定用成分として複数の成分から抽出される。このため、例えば、予め実験等により特定された所定の成分が抽出用成分として複数の成分から抽出される場合と比較して、人間の脳活動と関連性の低い成分が抽出用成分として複数の成分から抽出されるおそれを低減することができている。   In this embodiment, while the facial skin temperature acquisition means 20 and / or the image data acquisition means 120 are acquiring time-series facial skin temperature data and / or image data, the brain function activation task is performed on the individual. Given for a certain period. That is, in this embodiment, a situation is created in which the human brain is activated or not activated by actually giving or not giving a brain function activation task to an individual. The time series data acquired in this way is decomposed into a plurality of components by singular value decomposition, and the correlation between the component waveform and activation and non-activation of the brain is evaluated for each component. Are extracted from a plurality of components as determination components. For this reason, for example, compared with a case where a predetermined component specified in advance by experiments or the like is extracted from a plurality of components as an extraction component, a plurality of components that are less relevant to human brain activity are extracted as a plurality of extraction components. The possibility of being extracted from the components can be reduced.

(5−3)
ここで、脳には、選択的脳冷却機構という体温とは独立して脳を冷却する仕組みがある。選択的脳冷却機構としては、脳活動によって生じた熱を前額部及び副鼻腔周辺を用いて排熱していることが知られている。そうすると、脳活動に伴う顔面皮膚温度や顔面皮膚温度に相関する顔面の血行量の変化は、前額部及び/又は副鼻腔周辺に出現することになる。
(5-3)
Here, the brain has a mechanism that cools the brain independently of the body temperature, which is a selective brain cooling mechanism. As a selective brain cooling mechanism, it is known that heat generated by brain activity is exhausted using the forehead portion and the periphery of the sinuses. If it does so, the change of the blood flow volume of the face correlating with the facial skin temperature and facial skin temperature accompanying brain activity will appear in a forehead part and / or a sinus periphery.

本実施形態では、前額部及び/又は副鼻腔周辺の各種データが解析されて、判定用成分が抽出されている。このため、人間の脳活動に関連する成分を精度よく抽出することができている。   In the present embodiment, various data around the forehead and / or the sinuses are analyzed, and the determination component is extracted. For this reason, it is possible to accurately extract components related to human brain activity.

(6)脳活動可視化装置の用途例
次に、本発明に係る脳活動可視化装置の用途例について説明する。
(6) Application example of brain activity visualization apparatus Next, an application example of the brain activity visualization apparatus according to the present invention will be described.

(6−1)脳年齢の推定に用いる場合
上記実施形態又は上記変形例の脳活動可視化装置10,110を、対象者の脳年齢の推定に用いる場合の一例を説明する。例えば、上述した脳活動可視化装置10,110を脳年齢提示装置301に組み込むことで、検査対象の感情の判定が可能になる。
(6-1) When Used for Estimating Brain Age An example of using the brain activity visualization devices 10 and 110 of the embodiment or the modified example for estimating the brain age of the subject will be described. For example, by incorporating the brain activity visualization devices 10 and 110 described above into the brain age presentation device 301, it is possible to determine the emotion to be examined.

詳しくは、図23に示すように、上記状態可視化手段200に脳年齢記憶部204を設ける。この脳年齢記憶部204は、脳活動の推定結果と脳年齢を示す「脳年齢データ」とを関連付けて記憶するメモリである。ここでは、例えば、複数の対象者に同一の脳機能賦活課題を一定期間与えたときの脳活動量のデータを収集する。そして、対象者の年齢に応じて脳活動の推定結果を分類し、分類した脳活動の推定結果と脳年齢データとを関連付けて感情記憶部204に格納する。なお、脳活動の推定結果と脳年齢を示す脳年齢データとの相関関係は逐次最適化される。   Specifically, as shown in FIG. 23, a brain age storage unit 204 is provided in the state visualization means 200. The brain age storage unit 204 is a memory that stores an estimation result of brain activity and “brain age data” indicating the brain age in association with each other. Here, for example, data on the amount of brain activity when the same brain function activation task is given to a plurality of subjects for a certain period is collected. Then, the brain activity estimation results are classified according to the age of the subject, and the classified brain activity estimation results and brain age data are associated with each other and stored in the emotion storage unit 204. The correlation between the estimation result of the brain activity and the brain age data indicating the brain age is sequentially optimized.

そして、状態可視化手段200が、脳活動推定手段30,130により算出された脳活動の推定結果から脳年齢データを脳年齢記憶部204から抽出する。これにより、脳年齢提示装置301は、対象者の脳年齢を表示部202に提示することが可能となる。   Then, the state visualization unit 200 extracts brain age data from the brain age storage unit 204 from the estimation result of the brain activity calculated by the brain activity estimation units 30 and 130. Thereby, the brain age presentation device 301 can present the brain age of the subject on the display unit 202.

要するに、脳年齢提示装置301は、顔面皮膚温度取得手段20及び/又は画像データ取得手段120と、脳活動推定手段30,130と、状態可視化手段(脳年齢提示手段)200とを備える。顔面皮膚温度取得手段20は、顔面皮膚温度データを時系列で取得する。画像データ取得手段120は、顔面の撮像画像データを時系列で取得する。脳活動推定手段30,130は、顔面皮膚温度データ及び/又は顔面血行量データを、特異値分解、主成分分析或いは独立成分分析により分解することで得られた複数の成分に基づき、対象者の脳活動を推定する。状態可視化手段200は、脳活動推定手段30により推定された対象者の脳活動に基づき、対象者の脳年齢を提示する。したがって、脳年齢提示装置301は、脳波電極等の装着前に処理が必要なセンサを装着しなくても、簡便に脳年齢を提示することができる。   In short, the brain age presentation device 301 includes facial skin temperature acquisition means 20 and / or image data acquisition means 120, brain activity estimation means 30 and 130, and state visualization means (brain age presentation means) 200. The facial skin temperature acquisition means 20 acquires facial skin temperature data in time series. The image data acquisition unit 120 acquires captured image data of the face in time series. Based on a plurality of components obtained by decomposing facial skin temperature data and / or facial blood flow data by singular value decomposition, principal component analysis, or independent component analysis, the brain activity estimating means 30, 130 Estimate brain activity. The state visualization unit 200 presents the brain age of the subject based on the brain activity of the subject estimated by the brain activity estimation unit 30. Therefore, the brain age presentation device 301 can easily present the brain age without wearing a sensor that needs to be processed before wearing an electroencephalogram electrode or the like.

(6−2)ゲーム機を利用する場合
上記脳年齢提示装置301は、ゲーム機に格納されるものであってもよい。具体的には、顔面皮膚温度取得手段20及び/又は画像データ取得手段120が、ゲーム機に搭載される。ここでは、顔面皮膚温度取得手段20及び/又は画像データ取得手段120は、ゲーム機によるゲーム開始後の所定期間に、顔面皮膚温度データ及び/又は顔面の撮像画像データを時系列で取得する。このような構成により、脳年齢提示装置301は、ゲーム実行中に脳年齢を提示することができる。
(6-2) When using a game machine The brain age presentation device 301 may be stored in a game machine. Specifically, the facial skin temperature acquisition means 20 and / or the image data acquisition means 120 are mounted on a game machine. Here, the facial skin temperature acquisition means 20 and / or the image data acquisition means 120 acquire facial skin temperature data and / or captured image data of the face in time series for a predetermined period after the game is started by the game machine. With such a configuration, the brain age presentation device 301 can present the brain age during the execution of the game.

また、ゲーム機によるゲーム自体が脳機能賦活課題を構成するものであってもよい。この場合、ゲームの内容に応じて、複数の被験者から収集したデータに基づいて、年齢と脳活動昂進の推定式を事前に作成しておく。そして、ゲームが実行されるたびに、ゲーム毎に作成された推定式に基づいて脳年齢が推定される。これにより、ゲームが実行されるたびに「脳年齢x歳」等の提示がなされることになる。   Further, the game itself by the game machine may constitute a brain function activation task. In this case, an estimation formula for age and brain activity progression is created in advance based on data collected from a plurality of subjects according to the contents of the game. Each time the game is executed, the brain age is estimated based on the estimation formula created for each game. Thus, every time the game is executed, “brain age x year old” or the like is presented.

なお、脳年齢提示装置301は、脳年齢の提示とともにアンチエイジングの診断情報等を提供するものであってもよい。   The brain age presentation device 301 may provide anti-aging diagnostic information and the like together with the presentation of the brain age.

(7)スマートデバイスの利用
本実施形態に係る脳活動可視化装置10,110は、スマートデバイスに格納されるものであってもよい。具体的には、顔面皮膚温度取得手段20及び/又は画像データ取得手段120と、脳活動推定手段30,130と、状態可視化手段200とが多機能機能携帯電話等のスマートデバイスに格納されていてもよい。さらに、脳年齢提示装置として利用する脳活動可視化装置10,110をスマートデバイスに格納する場合、ユーザは携帯しているスマートデバイスにより情報の入出力及び演算が可能となるので、脳年齢を簡便に取得できるようになる。
(7) Use of Smart Device The brain activity visualization apparatuses 10 and 110 according to the present embodiment may be stored in a smart device. Specifically, facial skin temperature acquisition means 20 and / or image data acquisition means 120, brain activity estimation means 30, 130, and state visualization means 200 are stored in a smart device such as a multi-function mobile phone. Also good. Furthermore, when the brain activity visualization devices 10 and 110 used as the brain age presentation device are stored in the smart device, the user can input and output information and perform calculations with the smart device that is carried, so that the brain age can be easily set. You can get it.

なお、上述の形態に代えて、顔面皮膚温度取得手段20及び/又は画像データ取得手段120と状態可視化手段200とを多機能機能携帯電話等のスマートデバイスに格納し、脳活動推定手段30,130をサーバ装置に格納する構成を採用してもよい。この場合、スマートデバイス本体の演算負荷を軽減することができる。また、サーバ装置上でデータ管理することで、多数のデータを収集し最適化していくことで、出力結果の精度を高めていくことができる。例えば脳活動の推定結果と検査対象の脳年齢データとの関連を適時最適化していくことで、脳年齢判定の精度を高めることができる。   Instead of the above-described embodiment, the facial skin temperature acquisition means 20 and / or the image data acquisition means 120 and the state visualization means 200 are stored in a smart device such as a multi-function mobile phone, and the brain activity estimation means 30, 130 are stored. May be employed in the server device. In this case, the calculation load on the smart device body can be reduced. Further, by managing data on the server device, it is possible to improve the accuracy of the output result by collecting and optimizing a large number of data. For example, by accurately optimizing the relationship between the estimation result of brain activity and the brain age data to be examined, the accuracy of the brain age determination can be improved.

<付記>
なお、本発明は、上記実施形態そのままに限定されるものではない。本発明は、実施段階ではその要旨を逸脱しない範囲で構成要素を変形して具体化できる。また、本発明は、上記実施形態に開示されている複数の構成要素の適宜な組み合わせにより種々の発明を形成できるものである。例えば、実施形態に示される全構成要素から幾つかの構成要素は削除してもよいものである。さらに、異なる実施形態に構成要素を適宜組み合わせてもよいものである。
<Appendix>
In addition, this invention is not limited to the said embodiment as it is. The present invention can be embodied by modifying the components without departing from the scope of the invention in the implementation stage. In addition, the present invention can form various inventions by appropriately combining a plurality of constituent elements disclosed in the above embodiments. For example, some components may be deleted from all the components shown in the embodiment. Furthermore, constituent elements may be appropriately combined in different embodiments.

本発明によれば、簡便に脳活動を推定することができるため、脳活動に基づき対象者の生理状態を可視化する脳活動可視化装置への適用が有効である。   According to the present invention, since brain activity can be easily estimated, application to a brain activity visualization device that visualizes the physiological state of a subject based on brain activity is effective.

20 顔面皮膚温度取得手段
200 状態可視化手段
201 解析部
204 脳年齢記憶部
301 脳年齢提示装置
10,110 評価装置
30,130 脳活動推定手段
DESCRIPTION OF SYMBOLS 20 Facial skin temperature acquisition means 200 State visualization means 201 Analysis part 204 Brain age memory | storage part 301 Brain age presentation apparatus 10,110 Evaluation apparatus 30,130 Brain activity estimation means

特開2013−176406号公報JP 2013-176406 A

Claims (7)

対象者の顔面から検出した皮膚温度データ及びその検出部位の位置データを含む顔面皮膚温度データを時系列で取得する顔面皮膚温度取得手段(20)及び/又は対象者の顔面を時系列で撮影した撮影画像データに対して行われるRGB処理により得られた前記撮影画像データのRGBデータに基づき顔面血行量データを時系列で取得する顔面血行量取得手段と、
前記顔面皮膚温度取得手段により取得された前記顔面皮膚温度データ及び/又は前記顔面血行量取得手段により取得された前記顔面血行量データを、特異値分解、主成分分析或いは独立成分分析により分解することで得られた複数の成分に基づき、前記対象者の脳活動を推定する脳活動推定手段(30,130)と、
前記脳活動推定手段により推定された前記対象者の脳活動に基づき、前記対象者の脳年齢を提示する脳年齢提示手段(200)と、
を備える脳年齢提示装置(301)。
Facial skin temperature acquisition means (20) for acquiring the skin temperature data detected from the face of the subject and the face skin temperature data including the position data of the detected part in time series and / or the face of the subject was photographed in time series. Facial blood flow volume acquisition means for acquiring facial blood flow volume data in time series based on RGB data of the captured image data obtained by RGB processing performed on the captured image data;
Decomposing the facial skin temperature data acquired by the facial skin temperature acquisition means and / or the facial blood flow volume data acquired by the facial blood flow volume acquisition means by singular value decomposition, principal component analysis, or independent component analysis. Brain activity estimating means (30, 130) for estimating the brain activity of the subject based on a plurality of components obtained in
A brain age presenting means (200) for presenting the brain age of the subject based on the brain activity of the subject estimated by the brain activity estimating means;
A brain age presentation device (301) comprising:
前記脳活動推定手段は、
前記複数の成分のうち、その成分波形の振幅が脳の安静時及び脳の賦活時の変化と相関関係にある成分を判定用成分として抽出し、
前記判定用成分に基づき、前記対象者の脳活動を推定する、
請求項1に記載の脳年齢提示装置。
The brain activity estimating means includes
Among the plurality of components, the component waveform amplitude is correlated with a change at rest and activation of the brain as a component for determination,
Based on the determination component, the brain activity of the subject is estimated.
The brain age presentation apparatus according to claim 1.
前記顔面皮膚温度取得手段及び/又は前記顔面血行量取得手段は、前記対象者の顔面の副鼻腔周辺及び/又は前額部の前記顔面皮膚温度データ及び/又は前記顔面血行量データを少なくとも取得する、
請求項1又は2に記載の脳年齢提示装置。
The facial skin temperature acquisition means and / or the facial blood flow acquisition means acquire at least the facial skin temperature data and / or the facial blood flow data in the vicinity of the sinuses and / or the forehead of the subject's face. ,
The brain age presentation device according to claim 1 or 2.
脳活動の推定結果と脳年齢を示す脳年齢データとを関連付けて記憶する脳年齢記憶部(204)をさらに備え、
前記脳活動推定手段は、脳活動の推定結果を算出し、
前記脳年齢提示手段は、前記推定結果に基づいて前記脳年齢記憶部から前記対象者の脳年齢データを抽出して提示する、
請求項1から3のいずれか1項に記載の脳年齢提示装置。
A brain age storage unit (204) for storing the brain activity estimation result and the brain age data indicating the brain age in association with each other;
The brain activity estimating means calculates a brain activity estimation result,
The brain age presentation means extracts and presents the brain age data of the subject from the brain age storage unit based on the estimation result,
The brain age presentation apparatus according to any one of claims 1 to 3.
前記顔面皮膚温度取得手段及び/又は前記顔面血行量取得手段と、前記脳年齢提示手段とがスマートデバイスに格納されている、
請求項1から4のいずれか1項に記載の脳年齢提示装置。
The facial skin temperature acquisition means and / or the facial blood flow acquisition means and the brain age presentation means are stored in a smart device.
The brain age presentation apparatus according to any one of claims 1 to 4.
前記顔面皮膚温度取得手段及び/又は前記顔面血行量取得手段と、前記脳活動推定手段と、前記脳年齢提示手段とがスマートデバイスに格納されている、
請求項1から5のいずれか1項に記載の脳年齢提示装置。
The facial skin temperature acquisition means and / or the facial blood flow acquisition means, the brain activity estimation means, and the brain age presentation means are stored in a smart device.
The brain age presentation apparatus according to any one of claims 1 to 5.
前記顔面皮膚温度取得手段及び/又は前記顔面血行量取得手段は、ゲーム機に搭載されるものであり、前記ゲーム機によるゲーム開始後の所定期間に、前記顔面皮膚温度データ及び/又は前記顔面血行量データを時系列で取得する、
請求項1から6のいずれか1項に記載の脳年齢提示装置。
The facial skin temperature acquisition means and / or the facial blood flow amount acquisition means are mounted on a game machine, and the facial skin temperature data and / or the facial blood circulation during a predetermined period after the game is started by the game machine. Get quantity data in time series,
The brain age presentation apparatus according to any one of claims 1 to 6.
JP2015203358A 2015-10-15 2015-10-15 Brain age presentation device Active JP6093422B1 (en)

Priority Applications (5)

Application Number Priority Date Filing Date Title
JP2015203358A JP6093422B1 (en) 2015-10-15 2015-10-15 Brain age presentation device
PCT/JP2016/080711 WO2017065313A1 (en) 2015-10-15 2016-10-17 Useful information presentation device
US15/767,067 US10709338B2 (en) 2015-10-15 2016-10-17 Useful information presentation device
EP16855564.7A EP3363352B1 (en) 2015-10-15 2016-10-17 Useful information presentation device
CN201680060244.6A CN108135498B (en) 2015-10-15 2016-10-17 Useful information presentation device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2015203358A JP6093422B1 (en) 2015-10-15 2015-10-15 Brain age presentation device

Publications (2)

Publication Number Publication Date
JP6093422B1 true JP6093422B1 (en) 2017-03-08
JP2017074215A JP2017074215A (en) 2017-04-20

Family

ID=58261955

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2015203358A Active JP6093422B1 (en) 2015-10-15 2015-10-15 Brain age presentation device

Country Status (1)

Country Link
JP (1) JP6093422B1 (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113057594A (en) * 2021-03-19 2021-07-02 疲劳科学研究(广州)有限公司 Judgment processing system for fatigue degree evaluation and fatigue judgment method thereof

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR102607398B1 (en) * 2022-01-19 2023-11-29 이화여자대학교 산학협력단 Method and apparatus for monitoring of age-associated cognitive decline based on brain age

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006271815A (en) * 2005-03-30 2006-10-12 Kao Corp Method and system for measuring oxygen saturation of skin
JP2012239661A (en) * 2011-05-20 2012-12-10 Fujitsu Ltd Heart rate/respiration rate detection apparatus, method and program
JP2013198576A (en) * 2012-03-23 2013-10-03 Kao Corp Method for measuring skin pigment concentration

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006271815A (en) * 2005-03-30 2006-10-12 Kao Corp Method and system for measuring oxygen saturation of skin
JP2012239661A (en) * 2011-05-20 2012-12-10 Fujitsu Ltd Heart rate/respiration rate detection apparatus, method and program
JP2013198576A (en) * 2012-03-23 2013-10-03 Kao Corp Method for measuring skin pigment concentration

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113057594A (en) * 2021-03-19 2021-07-02 疲劳科学研究(广州)有限公司 Judgment processing system for fatigue degree evaluation and fatigue judgment method thereof

Also Published As

Publication number Publication date
JP2017074215A (en) 2017-04-20

Similar Documents

Publication Publication Date Title
US10667738B2 (en) Brain activity estimation device
JP6196402B2 (en) Determination result output device, determination result providing device, and determination result output system
US11786165B2 (en) Physiological state determination device
JP6503327B2 (en) Physiological condition determination apparatus and physiological condition determination method
WO2017065313A1 (en) Useful information presentation device
JP7111450B2 (en) Brain activity estimation device
JP2017200632A (en) Operator state determination device and operator state determination method
JP6463392B2 (en) Fatigue state determination device and fatigue state determination method
JP6093422B1 (en) Brain age presentation device
JP6158887B2 (en) Useful information presentation device
JP6240289B2 (en) Evaluation device, market research device, and learning evaluation device
JP6096857B1 (en) Emotion judgment device
JP2016182504A (en) Brain activity estimation device
JP2017086992A (en) Feeling determination device

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20161026

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: 20170117

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20170210

R150 Certificate of patent or registration of utility model

Ref document number: 6093422

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

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250