JP2018046899A - 脳活動状態定量化方法および脳活動状態計測装置 - Google Patents
脳活動状態定量化方法および脳活動状態計測装置 Download PDFInfo
- Publication number
- JP2018046899A JP2018046899A JP2016182662A JP2016182662A JP2018046899A JP 2018046899 A JP2018046899 A JP 2018046899A JP 2016182662 A JP2016182662 A JP 2016182662A JP 2016182662 A JP2016182662 A JP 2016182662A JP 2018046899 A JP2018046899 A JP 2018046899A
- Authority
- JP
- Japan
- Prior art keywords
- brain
- feature amount
- brain activity
- feature
- data
- 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.)
- Granted
Links
- 210000004556 brain Anatomy 0.000 title claims abstract description 58
- 238000000034 method Methods 0.000 title claims abstract description 42
- 238000013139 quantization Methods 0.000 title 1
- 230000004913 activation Effects 0.000 claims abstract description 17
- 230000000694 effects Effects 0.000 claims abstract description 17
- 230000007177 brain activity Effects 0.000 claims description 68
- 238000004364 calculation method Methods 0.000 claims description 25
- 238000004458 analytical method Methods 0.000 claims description 23
- 238000011002 quantification Methods 0.000 claims description 15
- 239000011159 matrix material Substances 0.000 claims description 11
- 238000005259 measurement Methods 0.000 claims description 7
- 238000000513 principal component analysis Methods 0.000 claims description 6
- 238000001228 spectrum Methods 0.000 claims description 6
- 238000013421 nuclear magnetic resonance imaging Methods 0.000 claims description 4
- 230000009467 reduction Effects 0.000 claims description 3
- 230000002490 cerebral effect Effects 0.000 claims 1
- 238000002599 functional magnetic resonance imaging Methods 0.000 description 10
- 238000002595 magnetic resonance imaging Methods 0.000 description 10
- INGWEZCOABYORO-UHFFFAOYSA-N 2-(furan-2-yl)-7-methyl-1h-1,8-naphthyridin-4-one Chemical compound N=1C2=NC(C)=CC=C2C(O)=CC=1C1=CC=CO1 INGWEZCOABYORO-UHFFFAOYSA-N 0.000 description 7
- 108010002255 deoxyhemoglobin Proteins 0.000 description 7
- 108010064719 Oxyhemoglobins Proteins 0.000 description 6
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical group [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 5
- 239000001301 oxygen Substances 0.000 description 5
- 229910052760 oxygen Inorganic materials 0.000 description 5
- 238000012545 processing Methods 0.000 description 5
- 238000012935 Averaging Methods 0.000 description 4
- 239000008280 blood Substances 0.000 description 4
- 210000004369 blood Anatomy 0.000 description 4
- 108010054147 Hemoglobins Proteins 0.000 description 3
- 102000001554 Hemoglobins Human genes 0.000 description 3
- 210000002569 neuron Anatomy 0.000 description 3
- 230000006870 function Effects 0.000 description 2
- 241000282412 Homo Species 0.000 description 1
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 description 1
- 238000004497 NIR spectroscopy Methods 0.000 description 1
- 238000013528 artificial neural network Methods 0.000 description 1
- 230000017531 blood circulation Effects 0.000 description 1
- 210000004027 cell Anatomy 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 239000000470 constituent Substances 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 229940079593 drug Drugs 0.000 description 1
- 239000001257 hydrogen Substances 0.000 description 1
- 229910052739 hydrogen Inorganic materials 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000002372 labelling Methods 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 238000010801 machine learning Methods 0.000 description 1
- 238000006213 oxygenation reaction Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 239000000941 radioactive substance Substances 0.000 description 1
- 238000011946 reduction process Methods 0.000 description 1
- 230000000284 resting effect Effects 0.000 description 1
- 238000012552 review Methods 0.000 description 1
- 230000002269 spontaneous effect Effects 0.000 description 1
- 238000012706 support-vector machine Methods 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Landscapes
- Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
Description
脳の活動状態を、脳の部位間のつながりにおける中心性を示す第1の特徴量と、脳の賦活状態を示す第2の特徴量で表し、
前記平面上の直交座標系の2つの座標軸のうち一方の座標軸を前記第1の特徴量に対応させ、かつ他方の座標軸を前記第2の座標軸に対応させ、脳の活動状態を、前記直交座標系に点として表示することを特徴とする。
(a)前記算出された相関係数をそのまま加え合わせて、任意の部位における第1の特徴量とする。
(b)前記算出された相関係数について、所定の閾値以上の値を示すものについて部位間の結合があると判定し、任意の部位において、相関係数が閾値以上で有る部位の数を算出し、その値を第1の特徴量とする。
(c)全ての部位間の前記相関係数を、同じ部位間の相関係数が対角要素となるように配置した相関係数行列に対して行列の最大固有値に対応する固有ベクトルを求め、その固有ベクトルの各要素の値を任意の部位における第1の特徴量とする。
(a)脳の任意の部位における時系列データを、高速フーリエ変換を用いて周波数スペクトルに変換する。
(b)変換された前記周波数スペクトルの特定の周波数帯の二乗平均平方根を、全周波数帯の二乗平均平方根で除算する。
人間の脳の活動状態を示す時系列データに基づいて、脳の部位間のつながりにおける中心性を示す第1の特徴量を算出する第1特徴量算出部と、
前記時系列データに基づいて、脳の賦活状態を示す第2の特徴量を算出する第2特徴量算出部と、
前記第1特徴量算出部から出力された第1の特徴量のデータおよび前記第2特徴量算出部から出力された第2の特徴量のデータについてそれぞれ線形判別分析、非線形判別分析または主成分分析のいずれかの手法を用いて次元の削減を行うことにより、一次元のデータに変換する分析部と、
前記分析部で一次元のデータに変換された第1および第2の特徴量のデータを、平面上の直交座標系の2つの座標軸に対応させて、脳活動の状態を直交座標系に点として表示する表示部と、を備えたことを特徴とする。
最初に、本発明における脳活動状態の定量化の原理について説明する。本発明における脳活動状態の定量化は、人の様々な脳活動の状態を定量化して平面上の直交座標系に表示するもので、1つの座標軸を、脳の部位間のつながりにおける中心性を示す第1の特徴量に対応させ、もう1つの座標軸を、脳活動の賦活状態を示す第2の特徴量に対応させている。
(a)前記算出された115の相関係数をそのまま加え合わせて任意の部位における第1の特徴量とする。
(b)前記算出された115の相関係数について、所定の閾値以上の値を示すものについて部位間の結合があると判定し、任意の部位において、相関係数が閾値以上で有る部位の数を算出し、その値を第1の特徴量とする。
(c)全ての部位間の前記相関係数を、同じ部位間の相関係数が対角要素となるように配置した相関係数行列に対して行列の最大固有値に対応する固有ベクトルを求め、その固有ベクトルの各要素の値を任意の部位における第1の特徴量とする。
図8に、本発明に係る脳活動状態定量化方法を実現する計測装置の構成を示す。
脳活動状態計測装置(以降、単に「計測装置」という)1は、MRI装置2から出力されたBOLD信号に基づいて前述した第1および第2の特徴量を算出すると共に、その値を平面上の直交座標系にプロットしてディスプレイに表示するものである。
2 MRI装置
11 インターフェース
12 第1特徴量算出部
13 第2特徴量算出部
14 分析部
15 表示部
16 制御部
Claims (9)
- 人間の脳における活動の状態を定量化し、その値を平面上に位置情報として表示する脳活動状態定量化方法であって、
脳の活動状態を、脳の部位間のつながりにおける中心性を示す第1の特徴量と、脳の賦活状態を示す第2の特徴量で表し、
前記平面上の直交座標系の2つの座標軸のうち一方の座標軸を前記第1の特徴量に対応させ、かつ他方の座標軸を前記第2の座標軸に対応させ、脳の活動状態を、前記直交座標系に点として表示することを特徴とする脳活動状態定量化方法。 - 脳の任意の部位における前記第1および第2の特徴量を、fMRI(機能的核磁気共鳴画像法)により得られる時系列データを用いて算出する、請求項1に記載の脳活動状態定量化方法。
- 前記時系列データは、脳の任意のボクセルの時系列データ、または複数のボクセルの時系列データの加算平均により得られる、請求項2に記載の脳活動状態定量化方法。
- 前記第1の特徴量は、脳の任意の部位における前記時系列データを、他の部位の時系列データと比較して相関係数を算出すると共に、当該相関係数の算出を前記他の全ての部位に対して行い、算出された当該相関係数を用いて、下記のいずれかの方法によって求める、請求項3に記載の脳活動状態定量化方法。
(a)前記算出された相関係数をそのまま加え合わせて、任意の部位における第1の特徴量とする。
(b)前記算出された相関係数について、所定の閾値以上の値を示すものについて部位間の結合があると判定し、任意の部位において、相関係数が閾値以上で有る部位の数を算出し、その値を第1の特徴量とする。
(c)全ての部位間の前記相関係数を、同じ部位間の相関係数が対角要素となるように配置した相関係数行列に対して行列の最大固有値に対応する固有ベクトルを求め、その固有ベクトルの各要素の値を任意の部位における第1の特徴量とする。 - 前記第2の特徴量は、以下の手順により求められる、請求項3または4に記載の脳活動状態定量化方法。
(a)脳の任意の部位における時系列データを、高速フーリエ変換を用いて周波数スペクトルに変換する。
(b)変換された前記周波数スペクトルの特定の周波数帯の二乗平均平方根を、全周波数帯の二乗平均平方根で除算する。 - 前記第1の特徴量および第2の特徴量のデータを、線形判別分析、非線形判別分析または主成分分析のいずれかの手法を用いて次元の削減を行うことにより、それぞれ一次元のデータに変換する、請求項4または5に記載の脳活動状態定量化方法。
- 人間の脳の活動状態を示す時系列データに基づいて、脳の部位間のつながりにおける中心性を示す第1の特徴量を算出する第1特徴量算出部と、
前記時系列データに基づいて、脳の賦活状態を示す第2の特徴量を算出する第2特徴量算出部と、
前記第1特徴量算出部から出力された第1の特徴量のデータおよび前記第2特徴量算出部から出力された第2の特徴量のデータについてそれぞれ線形判別分析、非線形判別分析または主成分分析のいずれかの手法を用いて次元の削減を行うことにより、一次元のデータに変換する分析部と、
前記分析部で一次元のデータに変換された第1および第2の特徴量のデータを、平面上の直交座標系の2つの座標軸に対応させて、脳活動の状態を直交座標系に点として表示する表示部と、を備えたことを特徴とする脳活動状態計測装置。 - 脳の任意の部位における前記第1および第2の特徴量を、fMRI(機能的核磁気共鳴画像法)、近赤外分光計測または脳波計測のいずれかにより得られる脳活動の時系列データを用いて算出する、請求項7に記載の脳活動状態計測装置。
- 前記時系列データは、脳の任意のボクセルの時系列データ、または複数のボクセルの時系列データの加算平均により得られる、請求項8に記載の脳活動状態計測装置。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2016182662A JP6821171B2 (ja) | 2016-09-20 | 2016-09-20 | 脳活動状態定量化方法および脳活動状態計測装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2016182662A JP6821171B2 (ja) | 2016-09-20 | 2016-09-20 | 脳活動状態定量化方法および脳活動状態計測装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2018046899A true JP2018046899A (ja) | 2018-03-29 |
JP6821171B2 JP6821171B2 (ja) | 2021-01-27 |
Family
ID=61766756
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2016182662A Active JP6821171B2 (ja) | 2016-09-20 | 2016-09-20 | 脳活動状態定量化方法および脳活動状態計測装置 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP6821171B2 (ja) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2020043951A (ja) * | 2018-09-18 | 2020-03-26 | 株式会社日立製作所 | 多機能神経フィードバックシステム及び多機能神経フィードバック方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2007289469A (ja) * | 2006-04-26 | 2007-11-08 | Hamano Life Science Research Foundation | 脳神経細胞の活動状態の評価方法 |
WO2014178323A1 (ja) * | 2013-05-01 | 2014-11-06 | 株式会社国際電気通信基礎技術研究所 | 脳活動解析装置、脳活動解析方法およびバイオマーカー装置 |
-
2016
- 2016-09-20 JP JP2016182662A patent/JP6821171B2/ja active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2007289469A (ja) * | 2006-04-26 | 2007-11-08 | Hamano Life Science Research Foundation | 脳神経細胞の活動状態の評価方法 |
WO2014178323A1 (ja) * | 2013-05-01 | 2014-11-06 | 株式会社国際電気通信基礎技術研究所 | 脳活動解析装置、脳活動解析方法およびバイオマーカー装置 |
JP2015062817A (ja) * | 2013-05-01 | 2015-04-09 | 株式会社国際電気通信基礎技術研究所 | 脳活動解析装置、脳活動解析方法およびバイオマーカー装置 |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2020043951A (ja) * | 2018-09-18 | 2020-03-26 | 株式会社日立製作所 | 多機能神経フィードバックシステム及び多機能神経フィードバック方法 |
JP7043374B2 (ja) | 2018-09-18 | 2022-03-29 | 株式会社日立製作所 | 多機能神経フィードバックシステム及び多機能神経フィードバック方法 |
Also Published As
Publication number | Publication date |
---|---|
JP6821171B2 (ja) | 2021-01-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Choe et al. | Comparing test-retest reliability of dynamic functional connectivity methods | |
Hou et al. | NIRS-KIT: a MATLAB toolbox for both resting-state and task fNIRS data analysis | |
Satterthwaite et al. | Linked sex differences in cognition and functional connectivity in youth | |
Huster et al. | Methods for simultaneous EEG-fMRI: an introductory review | |
Gramfort et al. | MNE software for processing MEG and EEG data | |
EP2992823B1 (en) | Brain activity training device | |
França et al. | Fractal and multifractal properties of electrographic recordings of human brain activity: toward its use as a signal feature for machine learning in clinical applications | |
Keilholz et al. | Dynamic properties of functional connectivity in the rodent | |
Omidvarnia et al. | Functional bimodality in the brain networks of preterm and term human newborns | |
EP1946701B1 (en) | Brain function analysis method and brain function analysis program | |
JP6884344B2 (ja) | 脳内ネットワークの活動推定システム、脳内ネットワークの活動推定方法、脳内ネットワークの活動推定プログラム、および、学習済み脳活動推定モデル | |
EP2992824A1 (en) | Brain activity analysis device, brain activity analysis method, and biomarker device | |
Pernet et al. | Best practices in data analysis and sharing in neuroimaging using MEEG | |
Storti et al. | Automatic selection of resting-state networks with functional magnetic resonance imaging | |
US20190117164A1 (en) | Methods and systems of de-noising magnetic-field based sensor data of electrophysiological signals | |
Paszkiel et al. | Methods of acquisition, archiving and biomedical data analysis of brain functioning | |
JP6782006B2 (ja) | 脳情報解析装置及び脳健康指標演算装置 | |
JP2008229238A (ja) | 脳活動解析方法および装置 | |
Shi et al. | Realistic models of apparent dynamic changes in resting‐state connectivity in somatosensory cortex | |
Koush et al. | Single voxel proton spectroscopy for neurofeedback at 7 Tesla | |
WO2018155447A1 (ja) | 脳情報解析装置及び脳情報模倣演算装置 | |
Ewald et al. | Wedge MUSIC: A novel approach to examine experimental differences of brain source connectivity patterns from EEG/MEG data | |
Yao et al. | Estimating the number of sources in magnetoencephalography using spiked population eigenvalues | |
JP6821171B2 (ja) | 脳活動状態定量化方法および脳活動状態計測装置 | |
WO2022004841A1 (ja) | 推定システム、推定方法、プログラム、推定モデル、脳活動トレーニング装置、脳活動トレーニング方法、および、脳活動トレーニングプログラム |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20190802 |
|
RD04 | Notification of resignation of power of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7424 Effective date: 20200701 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20201020 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20201110 |
|
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: 20201208 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20201224 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 6821171 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |