CN117100277A - 基于脑磁图、近红外脑成像的多模态同步脑成像方法 - Google Patents
基于脑磁图、近红外脑成像的多模态同步脑成像方法 Download PDFInfo
- Publication number
- CN117100277A CN117100277A CN202310789043.1A CN202310789043A CN117100277A CN 117100277 A CN117100277 A CN 117100277A CN 202310789043 A CN202310789043 A CN 202310789043A CN 117100277 A CN117100277 A CN 117100277A
- Authority
- CN
- China
- Prior art keywords
- brain
- activity
- matrix
- magnetoencephalography
- imaging
- 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
- 238000002582 magnetoencephalography Methods 0.000 title claims abstract description 89
- 238000002610 neuroimaging Methods 0.000 title claims abstract description 75
- 238000000034 method Methods 0.000 title claims abstract description 38
- 230000001360 synchronised effect Effects 0.000 title claims abstract description 29
- 238000012360 testing method Methods 0.000 claims abstract description 96
- 210000004556 brain Anatomy 0.000 claims abstract description 93
- 230000000694 effects Effects 0.000 claims abstract description 54
- 230000008035 nerve activity Effects 0.000 claims abstract description 49
- 230000001537 neural effect Effects 0.000 claims abstract description 43
- 230000008878 coupling Effects 0.000 claims abstract description 39
- 238000010168 coupling process Methods 0.000 claims abstract description 39
- 238000005859 coupling reaction Methods 0.000 claims abstract description 39
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 claims abstract description 21
- 239000008280 blood Substances 0.000 claims abstract description 21
- 210000004369 blood Anatomy 0.000 claims abstract description 21
- 229910052760 oxygen Inorganic materials 0.000 claims abstract description 21
- 239000001301 oxygen Substances 0.000 claims abstract description 21
- 238000003333 near-infrared imaging Methods 0.000 claims abstract description 18
- 210000005036 nerve Anatomy 0.000 claims abstract description 17
- 238000003384 imaging method Methods 0.000 claims abstract description 15
- 238000010586 diagram Methods 0.000 claims abstract description 12
- 230000007177 brain activity Effects 0.000 claims abstract description 11
- 230000000007 visual effect Effects 0.000 claims abstract description 7
- 230000003238 somatosensory effect Effects 0.000 claims abstract description 4
- 239000011159 matrix material Substances 0.000 claims description 100
- 230000001054 cortical effect Effects 0.000 claims description 73
- 238000012546 transfer Methods 0.000 claims description 47
- 210000004884 grey matter Anatomy 0.000 claims description 45
- 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 claims description 36
- 108010064719 Oxyhemoglobins Proteins 0.000 claims description 36
- 108010002255 deoxyhemoglobin Proteins 0.000 claims description 36
- 239000013598 vector Substances 0.000 claims description 18
- 108010054147 Hemoglobins Proteins 0.000 claims description 17
- 102000001554 Hemoglobins Human genes 0.000 claims description 17
- 230000002490 cerebral effect Effects 0.000 claims description 14
- 239000000523 sample Substances 0.000 claims description 11
- 230000004424 eye movement Effects 0.000 claims description 7
- 238000001914 filtration Methods 0.000 claims description 7
- 238000012937 correction Methods 0.000 claims description 6
- 238000000926 separation method Methods 0.000 claims description 6
- 230000001419 dependent effect Effects 0.000 claims description 5
- 230000007613 environmental effect Effects 0.000 claims description 5
- 238000004088 simulation Methods 0.000 claims description 5
- 230000000241 respiratory effect Effects 0.000 claims description 4
- 238000013461 design Methods 0.000 claims description 3
- 230000000004 hemodynamic effect Effects 0.000 abstract description 8
- 230000007831 electrophysiology Effects 0.000 abstract description 3
- 238000002001 electrophysiology Methods 0.000 abstract description 3
- 230000001149 cognitive effect Effects 0.000 abstract description 2
- 238000005259 measurement Methods 0.000 abstract description 2
- 210000003128 head Anatomy 0.000 description 25
- 230000003702 neurovascular coupling effect Effects 0.000 description 10
- 238000005516 engineering process Methods 0.000 description 9
- 230000003287 optical effect Effects 0.000 description 7
- 230000007274 generation of a signal involved in cell-cell signaling Effects 0.000 description 6
- 230000001939 inductive effect Effects 0.000 description 6
- 210000004761 scalp Anatomy 0.000 description 6
- 238000011160 research Methods 0.000 description 5
- 208000014644 Brain disease Diseases 0.000 description 2
- 210000004204 blood vessel Anatomy 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000001816 cooling Methods 0.000 description 2
- 230000004886 head movement Effects 0.000 description 2
- 238000013507 mapping Methods 0.000 description 2
- 239000013307 optical fiber Substances 0.000 description 2
- 230000001575 pathological effect Effects 0.000 description 2
- 238000002360 preparation method Methods 0.000 description 2
- 230000011218 segmentation Effects 0.000 description 2
- 230000000638 stimulation Effects 0.000 description 2
- 208000024827 Alzheimer disease Diseases 0.000 description 1
- 108091026922 FnrS RNA Proteins 0.000 description 1
- 206010020772 Hypertension Diseases 0.000 description 1
- 208000006011 Stroke Diseases 0.000 description 1
- 230000032683 aging Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000017531 blood circulation Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 210000001175 cerebrospinal fluid Anatomy 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 230000005684 electric field Effects 0.000 description 1
- 230000004927 fusion Effects 0.000 description 1
- 230000005358 geomagnetic field Effects 0.000 description 1
- 230000002650 habitual effect Effects 0.000 description 1
- 238000010348 incorporation Methods 0.000 description 1
- 238000012880 independent component analysis Methods 0.000 description 1
- 238000007917 intracranial administration Methods 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 230000002503 metabolic effect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003183 myoelectrical effect Effects 0.000 description 1
- 210000002569 neuron Anatomy 0.000 description 1
- 238000002496 oximetry Methods 0.000 description 1
- 230000008092 positive effect Effects 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000029058 respiratory gaseous exchange Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 210000003625 skull Anatomy 0.000 description 1
- 210000001519 tissue Anatomy 0.000 description 1
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
- A61B5/242—Detecting biomagnetic fields, e.g. magnetic fields produced by bioelectric currents
- A61B5/245—Detecting biomagnetic fields, e.g. magnetic fields produced by bioelectric currents specially adapted for magnetoencephalographic [MEG] signals
- A61B5/246—Detecting biomagnetic fields, e.g. magnetic fields produced by bioelectric currents specially adapted for magnetoencephalographic [MEG] signals using evoked responses
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/145—Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue
- A61B5/14542—Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue for measuring blood gases
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/145—Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue
- A61B5/1455—Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue using optical sensors, e.g. spectral photometrical oximeters
- A61B5/14551—Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue using optical sensors, e.g. spectral photometrical oximeters for measuring blood gases
- A61B5/14553—Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue using optical sensors, e.g. spectral photometrical oximeters for measuring blood gases specially adapted for cerebral tissue
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
- A61B5/316—Modalities, i.e. specific diagnostic methods
- A61B5/318—Heart-related electrical modalities, e.g. electrocardiography [ECG]
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
- A61B5/316—Modalities, i.e. specific diagnostic methods
- A61B5/389—Electromyography [EMG]
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
- A61B5/316—Modalities, i.e. specific diagnostic methods
- A61B5/398—Electrooculography [EOG], e.g. detecting nystagmus; Electroretinography [ERG]
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7203—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7225—Details of analog processing, e.g. isolation amplifier, gain or sensitivity adjustment, filtering, baseline or drift compensation
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7235—Details of waveform analysis
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7235—Details of waveform analysis
- A61B5/725—Details of waveform analysis using specific filters therefor, e.g. Kalman or adaptive filters
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Veterinary Medicine (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Signal Processing (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Artificial Intelligence (AREA)
- Physiology (AREA)
- Psychiatry (AREA)
- Optics & Photonics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Neurology (AREA)
- Power Engineering (AREA)
- Ophthalmology & Optometry (AREA)
- Cardiology (AREA)
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
本发明公开了一种基于脑磁图、近红外脑成像的多模态同步脑成像方法。所述方法包括：对测试对象施加一测试任务，并同步采集基于所述测试任务产生的大脑活动信号；其中，所述测试任务由诱发脑神经活动的视觉任务、听觉任务和体感任务中最少一项组成，所述大脑活动信号包括：多通道脑磁信号和多通道近红外成像信号；基于多通道脑磁信号和多通道近红外成像信号，绘制大脑神经活动成像图、大脑血氧活动成像图和大脑神经血管耦合系数图。本发明能够对全脑电生理、血液动力学进行长时间非侵入性测量，实现真实场景下认知神经科学研究。
Description
技术领域
本发明涉及多模态脑成像研究领域，具体涉及一种基于脑磁图、近红外脑成像的多模态同步脑成像方法。
背景技术
多模态神经影像汇集各单一模态的优势，根据不同模态信号特点从多角度描述大脑活动，如结合电生理和血液动力学信号研究神经血管耦合。多模态脑成像要求各模态需要互相兼容、互补，从而形成对基础科学研究和临床应用有效的多模态脑成像平台。多模态脑成像技术在脑科学研究和脑疾病临床实践中占据重要地位，多模态脑成像平台也是脑科学领域大科学装置中不可或缺的一部分。
脑磁图(MEG)是非侵入性的大脑电生理测量技术，能够以毫秒级时间分辨率在头皮上测量大脑神经元电流产生的磁场或电场，从而直接记录大脑电生理活动。然而，传统的MEG的核心部分超导量子干涉装置(SQUID)需要在大型低温冷却系统中工作，这增加了MEG与其他神经成像模态结合的难度。近年来，快速发展的光泵磁强计(OPM)为SQUID-MEG技术提供了一个非常有效的替代方案。OPM可以在室温环境无需低温冷却系统的情况下测量极弱的磁场，大幅降低设备采购和维护成本，并且可以实现更灵活的探测器阵列排布。与SQUID-MEG探测器相比，OPM-MEG探测器可以放置在更靠近头皮的位置，通过减少与颅内皮质源的距离，显著提高信噪比。灵活的探测器排布方式提高了OPM-MEG与其他成像模态的兼容性。具体而言，OPM探测器可以与其他模态传感器在帽子上交错排列或定制针对个体受试者头部形状和传感器排布需求的头盔。
近红外脑成像(fNIRS)是一种非电离、安全、低成本、便携、可穿戴的氧合血液动力学神经影像技术。利用神经血管耦合效应，fNIRS通过测量大脑功能性血液动力学反应(如含氧血红蛋白、无氧血红蛋白浓度的变化)，间接反映大脑神经活动情况。fNIRS使用近红外光照射头皮，将探测器放置在距光源几厘米处，检测头皮出射光强变化，并根据不同的模型计算大脑血红蛋白浓度。目前，MEG已与fNIRS结合，研究神经血管耦合效应，如运动诱发神经血管耦合、神经血管耦合的习惯效应，推动神经血管耦合效应的生理学研究(如在睡眠、正常衰老领域)和病理学研究(如在阿尔茨海默病、高血压和中风领域)。
综上所述，结合脑磁图、近红外脑成像的多模态脑成像技术，可以同步记录大脑电生理和血流动力学变化，拓宽脑科学和脑疾病的研究方向和应用场景。然而，到目前为止，尚无能够实现OPM-MEG和fNIRS同步采集的多模态脑成像系统。
发明内容
针对上述问题，本发明提供了一种一种基于脑磁图、近红外脑成像的多模态同步脑成像方法，可以实现同步记录大脑的电生理、血液动力学活动，对研究神经血管耦合与代谢活动有着积极作用。
本发明采用的技术方案为：
一种基于脑磁图、近红外脑成像的多模态同步脑成像方法，所述方法包括：
对测试对象施加一测试任务，并同步采集基于所述测试任务产生的大脑活动信号；其中，所述测试对象包括：健康受试者和患病的受试者，所述测试任务包括：由诱发脑神经活动的视觉任务、听觉任务和/或体感任务，所述大脑活动信号包括：多通道脑磁信号和多通道近红外成像信号；
基于所述多通道脑磁信号，构建包括皮层灰质在内的头模型，并在所述头模型的基础上构建脑神经活动的源空间后，结合测试对象头部与脑磁图传感器阵列的相对位置计算从源空间到脑磁图阵列的传递矩阵；其中，所述脑磁图传感器用于采集测试对象的多通道脑磁信号；
结合所述测试任务计算脑磁图的事件相关磁场，并根据所述事件相关磁场和所述传递矩阵，计算皮层灰质神经活动强度；
根据所述多通道近红外成像信号和所述测试任务的设计，获取事件相关血氧变化，并将所述事件相关血氧变化投影到所述头模型的皮层灰质上，以得到氧合血红蛋白的活动矩阵和脱氧血红蛋白的活动矩阵；
基于所述皮层灰质神经活动强度、所述氧合血红蛋白的活动矩阵和所述脱氧血红蛋白的活动矩阵，计算神经活动与氧合血红蛋白的耦合强度和神经活动与脱氧血红蛋白的耦合强度；
使用所述皮层灰质神经活动强度绘制大脑神经活动成像图，使用氧合血红蛋白的活动矩阵和脱氧血红蛋白的活动矩阵绘制大脑血氧活动成像图，使用神经活动与氧合血红蛋白的耦合强度和神经活动与脱氧血红蛋白的耦合强度绘制大脑神经血管耦合系数图。
进一步地，在所述头模型的基础上构建脑神经活动的源空间后，结合测试对象头部与脑磁图传感器阵列的相对位置计算从源空间到脑磁图阵列的传递矩阵，包括：
对测试对象的头模型进行离散网格化；
使用有限元仿真方法计算每个皮层灰质的网格点产生偶极子磁场后，将所述偶极子磁场传递到脑磁图传感器阵列处磁场；
将脑磁图传感器阵列处磁场大小除以对应网格点产生的所述偶极子磁场，得到该处皮层灰质网格神经活动产生磁场对传感器阵列的传递向量；
基于所有皮层灰质网格神经活动产生磁场对传感器阵列的传递向量，构成传递矩阵；其中，所述传递矩阵的行代表某个皮层灰质网格神经活动产生磁场的传递向量，所述传递矩阵的列代表不同的皮层灰质网格。
进一步地，根据所述事件相关磁场和所述传递矩阵，计算皮层灰质神经活动强度，包括：
采集无测试对象情况下同步空采光泵磁强计脑磁图、近红外脑成像，并根据同步空采的光泵磁强计脑磁图，计算脑磁图传感器的噪声协方差矩阵NMEG；
结合所述事件相关磁场、所述传递矩阵和所述噪声协方差矩阵NMEG，使用最小范数估计方法进行脑磁图溯源，以得到皮层灰质神经活动强度。
进一步地，基于所述皮层灰质神经活动强度、所述氧合血红蛋白的活动矩阵和所述脱氧血红蛋白的活动矩阵，计算氧合血红蛋白的耦合强度和脱氧血红蛋白的耦合强度，包括：
采集无测试对象情况下同步空采光泵磁强计脑磁图、近红外脑成像，并根据同步空采近红外脑成像，计算近红外脑成像的氧合血红蛋白的噪声协方差矩阵NHbO和脱氧血红蛋白的噪声协方差矩阵NHbT；
通过皮层灰质神经活动强度、氧合血红蛋白的活动矩阵和所述脱氧血红蛋白的活动矩阵，使用最小范数估计获取氧合血红蛋白的耦合强度CHbO和脱氧血红蛋白的耦合强度CHbR；其中，所述最小范数的表达公式为M＝CHbO(NHbO-σ1NHbO)T和M＝CHbR(HHbR-σ2HHbR)T，M表示皮层灰质神经活动强度，HHbO表示氧合血红蛋白的活动矩阵，HHbR表示脱氧血红蛋白的活动矩阵，σ1表示第一权重系数，σ2表示第二权重系数。
进一步地，所述同步采集基于所述测试任务产生的大脑活动信号的同时，还同步采集基于所述测试任务产生的人体电生理信号，所述人体电生理信号包括：眼动信号、心电信号和探头附件肌电信号。
进一步地，结合所述测试任务计算脑磁图的事件相关磁场，并根据所述事件相关磁场和所述传递矩阵，计算皮层灰质神经活动强度，还包括：
对所述多通道脑磁信号进行信号空间分离，以去除脑神经活动之外的环境干扰信号；
和，
时域上去除所述多通道脑磁信号工频的带阻滤波；
和，
去除所述多通道脑磁信号中与所述眼动信号、所述心电信号和探头附件肌电信号相关的成分；
和，
去除所述多通道近红外成像信号中由于头动校正、呼吸伪影和心跳伪影带来的噪声。
一种基于脑磁图、近红外脑成像的多模态同步脑成像系统，其特征在于，所述系统包括：
测试任务呈现装置，用于对测试对象施加一测试任务，所述测试对象包括：健康受试者和患病的受试者，所述测试任务包括：由诱发脑神经活动的视觉任务、听觉任务和/或体感任务；
第一信号生成装置，用于生成基于所述测试任务的多通道脑磁信号；
第二信号生成装置，用于生成基于所述测试任务产生的多通道近红外成像信号；
数据采集装置，用于同步采集多通道脑磁信号和多通道近红外成像信号；
脑成像装置，用于：
基于所述多通道脑磁信号，构建包括皮层灰质在内的头模型，并在所述头模型的基础上构建脑神经活动的源空间后，结合测试对象头部与脑磁图传感器阵列的相对位置计算从源空间到脑磁图阵列的传递矩阵；其中，所述脑磁图传感器用于采集测试对象的多通道脑磁信号；
结合所述测试任务计算脑磁图的事件相关磁场，并根据所述事件相关磁场和所述传递矩阵，计算皮层灰质神经活动强度；
根据所述多通道近红外成像信号，获取事件相关血氧变化，并将所述事件相关血氧变化投影到所述头模型的皮层灰质上，以得到氧合血红蛋白的活动矩阵和脱氧血红蛋白的活动矩阵；
基于所述皮层灰质神经活动强度、所述氧合血红蛋白的活动矩阵和所述脱氧血红蛋白的活动矩阵，计算神经活动与氧合血红蛋白的耦合强度和神经活动与脱氧血红蛋白的耦合强度；
使用所述皮层灰质神经活动强度绘制大脑神经活动成像图，使用氧合血红蛋白的活动矩阵和脱氧血红蛋白的活动矩阵绘制大脑血氧活动成像图，使用神经活动与氧合血红蛋白的耦合强度和神经活动与脱氧血红蛋白的耦合强度绘制大脑神经血管耦合系数图。
进一步地，还包括：磁屏蔽装置，所述磁屏蔽装置用于屏蔽与大脑神经活动无关的环境磁场。
进一步地，所述脑成像装置，还用于：
对测试对象的头模型进行离散网格化；
使用有限元仿真方法计算每个皮层灰质的网格点产生偶极子磁场后，将所述偶极子磁场传递到脑磁图传感器阵列处磁场；
将脑磁图传感器阵列处磁场大小除以对应网格点产生的所述偶极子磁场，得到该处皮层灰质网格神经活动产生磁场对传感器阵列的传递向量；
基于所有皮层灰质网格神经活动产生磁场对传感器阵列的传递向量，构成传递矩阵；其中，所述传递矩阵的行代表某个皮层灰质网格神经活动产生磁场的传递向量，所述传递矩阵的列代表不同的皮层灰质网格。
进一步地，所述脑成像装置，还用于：
采集无测试对象情况下同步空采光泵磁强计脑磁图、近红外脑成像，并根据同步空采的光泵磁强计脑磁图，计算脑磁图传感器的噪声协方差矩阵NMEG；
结合所述事件相关磁场、所述传递矩阵和所述噪声协方差矩阵NMEG，使用最小范数估计方法进行脑磁图溯源，以得到皮层灰质神经活动强度。
进一步地，所述脑成像装置，还用于：
采集无测试对象情况下同步空采光泵磁强计脑磁图、近红外脑成像，并根据同步空采近红外脑成像，计算近红外脑成像的氧合血红蛋白的噪声协方差矩阵NHbO和脱氧血红蛋白的噪声协方差矩阵NHbT；
通过皮层灰质神经活动强度、氧合血红蛋白的活动矩阵和所述脱氧血红蛋白的活动矩阵，使用最小范数估计获取氧合血红蛋白的耦合强度CHbO和脱氧血红蛋白的耦合强度CHbR；其中，所述最小范数的表达公式为M＝CHbO(HHbO-σ1HHbO)T和M＝CHbR(HHbR-σ2NHbR)T，M表示皮层灰质神经活动强度，HHbO表示氧合血红蛋白的活动矩阵，HHbR表示脱氧血红蛋白的活动矩阵，σ1表示第一权重系数，σ2表示第二权重系数。
进一步地，还包括：第三信号生成装置、第四信号生成装置和第五信号生成装置；其中，
所述第三信号生成装置用于生成基于所述测试任务产生的眼电信号；
所述第四信号生成装置用于生成基于所述测试任务产生的心电信号；
所述第五信号生成装置用于生成基于所述测试任务产生的探头附件肌电信号。
进一步地，所述脑成像装置还用于：
对所述多通道脑磁信号进行信号空间分离，以去除脑神经活动之外的环境干扰信号；
和，
时域上去除所述多通道脑磁信号工频的带阻滤波；
和，
去除所述多通道脑磁信号中与所述眼动信号、所述心电信号和探头附件肌电信号相关的成分；
和，
去除所述多通道近红外成像信号中由于头动校正、呼吸伪影和心跳伪影带来的噪声。
相比于现有技术，本发明至少具有以下优点：
1.能够实现对全脑电生理、血液动力学进行长时间非侵入性测量，适用于几乎所有健康人群和临床患者。
2.多模态同步采集和融合分析可以更全面更高效的刻画大脑活动情况，研究大脑在不同生理和病理条件下的特点，为基础科学和临床应用提供有效研究手段。
3.结合可穿戴脑磁图，本系统能够在最大程度还原使用场景的生态效度(ecological validity)，实现真实场景下认知神经科学研究。
附图说明
图1为本发明的硬件系统构成示意图。
具体实施方式
在下文的描述中，通过具体的实施例对本发明的多模态脑成像系统作进一步描述，以便于本领域技术人员对本发明的特征和优点有更为彻底的理解。
本发明的基于脑磁图、近红外脑成像的多模态同步脑成像方法，如图1所示，其硬件系统包括以下模块：测试任务呈现装置、第一信号生成装置、第二信号生成装置、第三信号生成装置、第四信号生成装置、第五信号生成装置、磁屏蔽装置和脑成像装置。
一.测试任务呈现装置。
测试任务呈现装置用于给测试对象呈现诱发脑神经活动的任务信息，任务信息可以通过多种模态呈现(例如视觉、听觉和体感)，其中，测试对象包括：健康受试者和患病的受试者。测试任务呈现装置将任务的触发信号发送给数据采集装置，以实现后者的同步。
二.第一信号生成装置。
第一信号生成装置是一个脑磁信号采集装置，其用于在给测试对象呈现诱发脑神经活动的任务信息(触发信号)时，采集大脑神经活动产生的磁信号的探测器和配套器件(如前/后端放大器等)；脑磁信号采集装置的实现方式为光泵磁强计；脑磁信号采集装置表面配置多个探测器，每个探测器提供一个信号通道，各探测器采集的模拟信号分别发送给数据采集装置。
根据本发明的一个实施例，脑磁信号采集装置使用基于无自旋交换弛豫(SpinExchange Relaxation Free,SERF)原理的光泵磁强计脑磁图，但本系统也适用于其他原理的光泵磁强计脑磁图。
三.第二信号生成装置。
第二信号生成装置是一个近红外脑成像采集装置，其用于在给测试对象呈现诱发脑神经活动的任务信息时，发射近红外光的光源，采集携带大脑功能性血液动力学反应产生的出射光信号的探测器和配套器件(如前/后端放大器等)；近红外脑成像采集装置的实现方式可以通过连续波技术，频域技术，或者时域技术；近红外脑成像采集装置表面配置多个光源和探测器，光源和探测器配对成为一个信号通道，各通道的模拟信号分别发送给数据采集装置。
根据本发明的一个实施例，近红外脑成像采集装置使用连续波技术的近红外脑成像，但本系统也适用于频域技术或者时域技术的近红外脑成像。
四.第三信号生成装置。
第三信号生成装置是一个眼电信号采集装置，其用于在给测试对象呈现诱发脑神经活动的任务信息，采集测试对象当前的眼电信号。
五.第四信号生成装置。
第四信号生成装置是一个心电信号采集装置，其用于在给测试对象呈现诱发脑神经活动的任务信息，采集测试对象当前的心电信号。
六、第五信号生成装置。
第五信号生成装置是一个肌电信号采集装置，其用于在给测试对象呈现诱发脑神经活动的任务信息，采集测试对象当前的探头附件肌电信号。
七.磁屏蔽装置。
磁屏蔽装置用于屏蔽与大脑神经活动无关的环境磁场(例如地磁场，电子设备的干扰磁信号)。
根据本发明的一个实施例，磁屏蔽装置可以用使用屏蔽效果更优的多层磁屏蔽桶或能够提供自由运动空间的磁屏蔽间。对于可穿戴式的光泵磁强计脑磁图，还可以在每个光泵磁强计探测器上附加主动补偿线圈用于抵消由于运动引起的环境剩磁干扰信号。
八.数据采集装置。
数据采集装置可以同步采集多通道脑磁信号、多通道脑电信号、多通道近红外成像信号、触发信号和其他电生理信号(如眼电信号、心电信号)；数据采集装置对采集信号进行预处理和保存。
九.脑成像装置。
脑成像装置是本发明的核心装置，其基于预处理后的脑磁信号、近红外成像信号、触发信号、眼电信号、心电信号、探头附件肌电信号，绘制大脑神经活动成像图、大脑血氧活动成像图和大脑神经血管耦合系数图。
根据本发明的一个实施例，该脑成像装置融合脑成像的步骤如下：
1)采集测试对象的磁共振或CT头部结构图像，使用灰度梯度法、灰度阈值法、区域分割法等头部组织分割算法，构建包括头皮、头骨、脑脊液和皮层灰质在内的头模型。
2)在头模型的基础上构建脑神经活动的源空间，并结合测试对象头部与脑磁图传感器阵列的相对位置计算从源空间到脑磁图阵列的传递矩阵GMEG。具体方法如下：对测试对象的头模型进行离散网格化，使用有限元仿真方法计算每个皮层灰质的网格点产生偶极子磁场后传递到脑磁图传感器阵列的磁场大小。将阵列处磁场大小除以对应网格点产生的磁场，得到该处皮层灰质网格神经活动产生磁场对传感器阵列的传递向量。将所有皮层灰质神经活动的传递向量并在一起构成传递矩阵GMEG，其中行代表某个皮层灰质网格神经活动产生磁场的传递向量，列代表不同的皮层灰质网格。
3)对采集到的脑磁信号进行信号空间分离(Signal-space separation)，用于去除脑神经活动之外的环境干扰信号。对脑磁信号进一步预处理，包括时域上去除工频(50Hz)的带阻滤波，提取所需频率成分的带通滤波，去除眼动、心电和探头附件肌电噪声的独立成分分析，进行头动校正。对与近红外脑成像数据预处理，包括头动校正，呼吸、心跳伪影去除，提取所需频率成分的带通滤波。
4)采集无测试对象情况下同步空采光泵磁强计脑磁图、近红外脑成像，计算脑磁图传感器的噪声协方差矩阵NMEG和近红外脑成像的氧合血红代表和脱氧血红蛋白的噪声协方差矩阵NHbO和NHbT。
5)按照测试任务设计和测试触发信号结合广义线性模型计算脑磁图的事件相关磁场E(E的行代表某个通道事件相关磁场的时间序列，列代表不同的通道)和近红外的事件相关血氧变化。将全脑血氧变化投影到头模型的皮层灰质上，并进行插值，得到皮层灰质血氧活动矩阵HHbO和HHbR，其中H行代表某个皮层灰质网格处血氧活动的时间序列，列代表不同的皮层灰质网格，HHbO代表氧合血红蛋白的活动矩阵，HHbR代表脱氧血红蛋白的活动矩阵。结合事件相关磁场E和传递矩阵GMEG，使用最小范数估计方法进行脑磁图溯源得到皮层灰质神经活动强度M。皮层灰质神经活动强度M的行代表某个皮层灰质网格处神经活动的时间序列，列代表不同的皮层灰质网格。
6)通过皮层灰质神经活动强度M和皮层灰质血氧活动矩阵HHbO和HHbR，使用最小范数估计获取神经血管耦合向量CHbO和CHbR，计算公式如下。神经血管耦合列向量CHbO和CHbR的代表皮层灰质的神经血管耦合强度，每个位置的取值和皮层灰质网格点对应，CHbO代表氧合血红蛋白的耦合强度，CHbR代表脱氧血红蛋白的耦合强度。
M＝CHbO(HHbO-σ1NHbO)T
M＝CHbR(HHbO-σ2NHbR)T
7)使用皮层灰质神经活动强度M绘制大脑神经活动成像图，使用皮层灰质血氧活动矩阵HHbO和HHbR绘制大脑血氧活动成像图，使用神经血管耦合向量CHbO和CHbR绘制大脑神经血管耦合系数图。
综上所述，在一个实现同步采集大脑的电生理、血液动力学活动信号的多模态脑成像系统的应用，可以包括准备阶段和测试阶段两部分。
准备阶段中，根据使用者头部形状定制搭载OPM探测器和fNIRS光纤头的刚性头盔。编写视觉刺激程序，绘制视觉刺激图片。
测试阶段中，测试对象在磁屏蔽环境中佩戴刚性头盔，紧贴头皮安装OPM探测器，拨开测试对象头发安装近红外脑成像光纤头。测试对象准备好后开始测试，播放视觉刺激程序，测试对象在测试过程中注视屏幕，各模态采集装置同步采集视觉刺激诱发的大脑活动信号，以完成大脑神经活动成像图、大脑血氧活动成像图和大脑神经血管耦合系数图的绘制。
上述实施例仅例示性说明本发明的原理及其功效，而非用于限制本发明。任何本领域技术人员均可在不违背本发明的精神及范畴下，对上述实施例进行修饰与改变。因此，本发明的权利保护范围，应如权利要求书所列。
Claims (13)
1.一种基于脑磁图、近红外脑成像的多模态同步脑成像方法，其特征在于，所述方法包括：
对测试对象施加一测试任务，并同步采集基于所述测试任务产生的大脑活动信号；其中，所述测试对象包括：健康受试者和患病的受试者，所述测试任务包括：由诱发脑神经活动的视觉任务、听觉任务和/或体感任务，所述大脑活动信号包括：多通道脑磁信号和多通道近红外成像信号；
基于所述多通道脑磁信号，构建包括皮层灰质在内的头模型，并在所述头模型的基础上构建脑神经活动的源空间后，结合测试对象头部与脑磁图传感器阵列的相对位置计算从源空间到脑磁图阵列的传递矩阵；其中，所述脑磁图传感器用于采集测试对象的多通道脑磁信号；
结合所述测试任务计算脑磁图的事件相关磁场，并根据所述事件相关磁场和所述传递矩阵，计算皮层灰质神经活动强度；
根据所述多通道近红外成像信号和所述测试任务的设计，获取事件相关血氧变化，并将所述事件相关血氧变化投影到所述头模型的皮层灰质上，以得到氧合血红蛋白的活动矩阵和脱氧血红蛋白的活动矩阵；
基于所述皮层灰质神经活动强度、所述氧合血红蛋白的活动矩阵和所述脱氧血红蛋白的活动矩阵，计算神经活动与氧合血红蛋白的耦合强度和神经活动与脱氧血红蛋白的耦合强度；
使用所述皮层灰质神经活动强度绘制大脑神经活动成像图，使用氧合血红蛋白的活动矩阵和脱氧血红蛋白的活动矩阵绘制大脑血氧活动成像图，使用神经活动与氧合血红蛋白的耦合强度和神经活动与脱氧血红蛋白的耦合强度绘制大脑神经血管耦合系数图。
2.如权利要求1所述的方法，其特征在于，在所述头模型的基础上构建脑神经活动的源空间后，结合测试对象头部与脑磁图传感器阵列的相对位置计算从源空间到脑磁图阵列的传递矩阵，包括：
对测试对象的头模型进行离散网格化；
使用有限元仿真方法计算每个皮层灰质的网格点产生偶极子磁场后，将所述偶极子磁场传递到脑磁图传感器阵列处磁场；
将脑磁图传感器阵列处磁场大小除以对应网格点产生的所述偶极子磁场，得到该处皮层灰质网格神经活动产生磁场对传感器阵列的传递向量；
基于所有皮层灰质网格神经活动产生磁场对传感器阵列的传递向量，构成传递矩阵；其中，所述传递矩阵的行代表某个皮层灰质网格神经活动产生磁场的传递向量，所述传递矩阵的列代表不同的皮层灰质网格。
3.如权利要求1所述的方法，其特征在于，根据所述事件相关磁场和所述传递矩阵，计算皮层灰质神经活动强度，包括：
采集无测试对象情况下同步空采光泵磁强计脑磁图、近红外脑成像，并根据同步空采的光泵磁强计脑磁图，计算脑磁图传感器的噪声协方差矩阵NMEG；
结合所述事件相关磁场、所述传递矩阵和所述噪声协方差矩阵NMEG，使用最小范数估计方法进行脑磁图溯源，以得到皮层灰质神经活动强度。
4.如权利要求1所述的方法，其特征在于，基于所述皮层灰质神经活动强度、所述氧合血红蛋白的活动矩阵和所述脱氧血红蛋白的活动矩阵，计算氧合血红蛋白的耦合强度和脱氧血红蛋白的耦合强度，包括：
采集无测试对象情况下同步空采光泵磁强计脑磁图、近红外脑成像，并根据同步空采近红外脑成像，计算近红外脑成像的氧合血红蛋白的噪声协方差矩阵NHbO和脱氧血红蛋白的噪声协方差矩阵NHbT；
通过皮层灰质神经活动强度、氧合血红蛋白的活动矩阵和所述脱氧血红蛋白的活动矩阵，使用最小范数估计获取氧合血红蛋白的耦合强度CHbO和脱氧血红蛋白的耦合强度CHbR；其中，所述最小范数的表达公式为M＝CHbO(NHbO-σ1NHbO)T和M＝CHbR(HHbR-σ2HHbR)T，M表示皮层灰质神经活动强度，HHbO表示氧合血红蛋白的活动矩阵，HHbR表示脱氧血红蛋白的活动矩阵，σ1表示第一权重系数，σ2表示第二权重系数。
5.如权利要求1至4任一项所述的方法，其特征在于，所述同步采集基于所述测试任务产生的大脑活动信号的同时，还同步采集基于所述测试任务产生的人体电生理信号，所述人体电生理信号包括：眼动信号、心电信号和探头附件肌电信号。
6.如权利要求5所述的方法，其特征在于，结合所述测试任务计算脑磁图的事件相关磁场，并根据所述事件相关磁场和所述传递矩阵，计算皮层灰质神经活动强度，还包括：
对所述多通道脑磁信号进行信号空间分离，以去除脑神经活动之外的环境干扰信号；
和，
时域上去除所述多通道脑磁信号工频的带阻滤波；
和，
去除所述多通道脑磁信号中与所述眼动信号、所述心电信号和探头附件肌电信号相关的成分；
和，
去除所述多通道近红外成像信号中由于头动校正、呼吸伪影和心跳伪影带来的噪声。
7.一种基于脑磁图、近红外脑成像的多模态同步脑成像系统，其特征在于，所述系统包括：
测试任务呈现装置，用于对测试对象施加一测试任务，所述测试对象包括：健康受试者和患病的受试者，所述测试任务包括：由诱发脑神经活动的视觉任务、听觉任务和/或体感任务；
第一信号生成装置，用于生成基于所述测试任务的多通道脑磁信号；
第二信号生成装置，用于生成基于所述测试任务产生的多通道近红外成像信号；
数据采集装置，用于同步采集多通道脑磁信号和多通道近红外成像信号；
脑成像装置，用于：
基于所述多通道脑磁信号，构建包括皮层灰质在内的头模型，并在所述头模型的基础上构建脑神经活动的源空间后，结合测试对象头部与脑磁图传感器阵列的相对位置计算从源空间到脑磁图阵列的传递矩阵；其中，所述脑磁图传感器用于采集测试对象的多通道脑磁信号；
结合所述测试任务计算脑磁图的事件相关磁场，并根据所述事件相关磁场和所述传递矩阵，计算皮层灰质神经活动强度；
根据所述多通道近红外成像信号，获取事件相关血氧变化，并将所述事件相关血氧变化投影到所述头模型的皮层灰质上，以得到氧合血红蛋白的活动矩阵和脱氧血红蛋白的活动矩阵；
基于所述皮层灰质神经活动强度、所述氧合血红蛋白的活动矩阵和所述脱氧血红蛋白的活动矩阵，计算神经活动与氧合血红蛋白的耦合强度和神经活动与脱氧血红蛋白的耦合强度；
使用所述皮层灰质神经活动强度绘制大脑神经活动成像图，使用氧合血红蛋白的活动矩阵和脱氧血红蛋白的活动矩阵绘制大脑血氧活动成像图，使用神经活动与氧合血红蛋白的耦合强度和神经活动与脱氧血红蛋白的耦合强度绘制大脑神经血管耦合系数图。
8.如权利要求7所述的系统，其特征在于，还包括：磁屏蔽装置，所述磁屏蔽装置用于屏蔽与大脑神经活动无关的环境磁场。
9.如权利要求7或8所述的系统，其特征在于，所述脑成像装置，还用于：
对测试对象的头模型进行离散网格化；
使用有限元仿真方法计算每个皮层灰质的网格点产生偶极子磁场后，将所述偶极子磁场传递到脑磁图传感器阵列处磁场；
将脑磁图传感器阵列处磁场大小除以对应网格点产生的所述偶极子磁场，得到该处皮层灰质网格神经活动产生磁场对传感器阵列的传递向量；
基于所有皮层灰质网格神经活动产生磁场对传感器阵列的传递向量，构成传递矩阵；其中，所述传递矩阵的行代表某个皮层灰质网格神经活动产生磁场的传递向量，所述传递矩阵的列代表不同的皮层灰质网格。
10.如权利要求7或8所述的系统，其特征在于，所述脑成像装置，还用于：
采集无测试对象情况下同步空采光泵磁强计脑磁图、近红外脑成像，并根据同步空采的光泵磁强计脑磁图，计算脑磁图传感器的噪声协方差矩阵NMEG；
结合所述事件相关磁场、所述传递矩阵和所述噪声协方差矩阵NMEG，使用最小范数估计方法进行脑磁图溯源，以得到皮层灰质神经活动强度。
11.如权利要求7或8所述的系统，其特征在于，所述脑成像装置，还用于：
采集无测试对象情况下同步空采光泵磁强计脑磁图、近红外脑成像，并根据同步空采近红外脑成像，计算近红外脑成像的氧合血红蛋白的噪声协方差矩阵NHbO和脱氧血红蛋白的噪声协方差矩阵NHbT；
通过皮层灰质神经活动强度、氧合血红蛋白的活动矩阵和所述脱氧血红蛋白的活动矩阵，使用最小范数估计获取氧合血红蛋白的耦合强度CHbO和脱氧血红蛋白的耦合强度CHbR；其中，所述最小范数的表达公式为M＝CHbO(HHbO-σ1NHbO)T和M＝CHbR(HHbR-σ2NHbR)T，M表示皮层灰质神经活动强度，HHbO表示氧合血红蛋白的活动矩阵，HHbR表示脱氧血红蛋白的活动矩阵，σ1表示第一权重系数，σ2表示第二权重系数。
12.如权利要求7或8所述的系统，其特征在于，还包括：第三信号生成装置、第四信号生成装置和第五信号生成装置；其中，
所述第三信号生成装置用于生成基于所述测试任务产生的眼电信号；
所述第四信号生成装置用于生成基于所述测试任务产生的心电信号；
所述第五信号生成装置用于生成基于所述测试任务产生的探头附件肌电信号。
13.如权利要求12所述的系统，其特征在于，所述脑成像装置还用于：
对所述多通道脑磁信号进行信号空间分离，以去除脑神经活动之外的环境干扰信号；
和，
时域上去除所述多通道脑磁信号工频的带阻滤波；
和，
去除所述多通道脑磁信号中与所述眼动信号、所述心电信号和探头附件肌电信号相关的成分；
和，
去除所述多通道近红外成像信号中由于头动校正、呼吸伪影和心跳伪影带来的噪声。
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN2022108778773 | 2022-07-25 | ||
| CN202210877877 | 2022-07-25 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN117100277A true CN117100277A (zh) | 2023-11-24 |
| CN117100277B CN117100277B (zh) | 2024-04-30 |
Family
ID=88811655
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310789043.1A Active CN117100277B (zh) | 2022-07-25 | 2023-06-29 | 基于脑磁图、近红外脑成像的多模态同步脑成像方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN117100277B (zh) |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20020058867A1 (en) * | 1999-12-02 | 2002-05-16 | Breiter Hans C. | Method and apparatus for measuring indices of brain activity during motivational and emotional function |
| US20040082862A1 (en) * | 2002-07-10 | 2004-04-29 | Britton Chance | Examination and imaging of brain cognitive functions |
| WO2012162569A2 (en) * | 2011-05-24 | 2012-11-29 | The Regents Of The University Of California | Magnetoencephalography source imaging |
| CN103325119A (zh) * | 2013-06-27 | 2013-09-25 | 中国科学院自动化研究所 | 一种基于模态融合的默认态脑网络中心节点检测方法 |
| CN107550493A (zh) * | 2017-07-14 | 2018-01-09 | 北京大学 | 基于功能磁共振成像的时变约束脑电图或脑磁图溯源方法 |
| CN113918008A (zh) * | 2021-08-30 | 2022-01-11 | 北京大学 | 一种基于源空间脑磁信号解码的脑机接口系统及应用方法 |
-
2023
- 2023-06-29 CN CN202310789043.1A patent/CN117100277B/zh active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20020058867A1 (en) * | 1999-12-02 | 2002-05-16 | Breiter Hans C. | Method and apparatus for measuring indices of brain activity during motivational and emotional function |
| US20040082862A1 (en) * | 2002-07-10 | 2004-04-29 | Britton Chance | Examination and imaging of brain cognitive functions |
| WO2012162569A2 (en) * | 2011-05-24 | 2012-11-29 | The Regents Of The University Of California | Magnetoencephalography source imaging |
| CN103717129A (zh) * | 2011-05-24 | 2014-04-09 | 加利福尼亚大学董事会 | 脑磁图源成像 |
| CN103325119A (zh) * | 2013-06-27 | 2013-09-25 | 中国科学院自动化研究所 | 一种基于模态融合的默认态脑网络中心节点检测方法 |
| CN107550493A (zh) * | 2017-07-14 | 2018-01-09 | 北京大学 | 基于功能磁共振成像的时变约束脑电图或脑磁图溯源方法 |
| CN113918008A (zh) * | 2021-08-30 | 2022-01-11 | 北京大学 | 一种基于源空间脑磁信号解码的脑机接口系统及应用方法 |
Non-Patent Citations (3)
| Title |
|---|
| S.A. HOSSEINI，等: "Hybrid approach in recognition of visual covert selective spatial attention based on MEG signals", 2015 IEEE INTERNATIONAL CONFERENCE ON FUZZY SYSTEMS (FUZZ-IEEE), 5 August 2015 (2015-08-05), pages 1 - 7, XP032818948, DOI: 10.1109/FUZZ-IEEE.2015.7337958 * |
| XINGYU RU，等: "Multimodal neuroimaging with optically pumped magnetometers: A simultaneous MEG-EEG-fNIRS acquisition system", NEUROIMAGE, vol. 259, 31 July 2022 (2022-07-31), pages 1 - 8 * |
| YOUSRA BEKHTI，等: "M/EEG source localization with multi-scale time-frequency dictionaries", 2016 INTERNATIONAL WORKSHOP ON PATTERN, 22 June 2016 (2016-06-22), pages 1 - 3 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN117100277B (zh) | 2024-04-30 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Niedermeyer | Niedermeyer's electroencephalography: basic principles, clinical applications, and related fields | |
| Teplan | Fundamentals of EEG measurement | |
| Yao et al. | Evaluation of different cortical source localization methods using simulated and experimental EEG data | |
| Braeutigam | Magnetoencephalography: fundamentals and established and emerging clinical applications in radiology | |
| Riera et al. | Fusing EEG and fMRI based on a bottom-up model: inferring activation and effective connectivity in neural masses | |
| Velmurugan et al. | Magnetoencephalography recording and analysis | |
| Kristeva-Feige et al. | The bereitschaftspotential paradigm in investigating voluntary movement organization in humans using magnetoencephalography (MEG) | |
| Fabietti et al. | Effectiveness of employing multimodal signals in removing artifacts from neuronal signals: an empirical analysis | |
| Halchenko et al. | Multimodal integration: fMRI, mri, EEG, MEG | |
| Rahman et al. | Fundamental study on brain signal for BCI-FES system development | |
| Mahmood et al. | Introduction to Non-Invasive Biomedical Signals for Healthcare | |
| SITARAM et al. | 18| BCIS THAT USE BRAIN METABOLIC SIGNALS | |
| CN117100277B (zh) | 基于脑磁图、近红外脑成像的多模态同步脑成像方法 | |
| Song et al. | Mobile electroencephalography for studying neural control of human locomotion | |
| Sato et al. | Movement-related cortical activation with voluntary pinch task: simultaneous monitoring of near-infrared spectroscopy signals and movement-related cortical potentials | |
| Wang et al. | Methods for Improving Movement Compatibility of Wearable OPM-MEG: A Review | |
| Ahmad et al. | Optimization and development of concurrent EEG-fMRI data acquisition setup for understanding neural mechanisms of brain | |
| Nakano | EEG measurement as a tool for rehabilitation assessment and treatment | |
| Berestov et al. | Brain–computer interface technologies for monitoring and control of bionic systems | |
| Papadelis | MEG systems for young children and recent developments of pediatric MEG | |
| Hu et al. | Constructing bodily emotion maps based on high-density body surface potentials for psychophysiological computing | |
| Saxena et al. | Understanding Functional Brain Activation using Source Localization of EEG Signals in Motor Imagery Tasks | |
| Lei | Simultaneous eeg-fmri | |
| Suprijanto et al. | Future Trend of EEG Source Localization for Quantifying Repetitive Motor Activity as feedback for Rehabilitation Device | |
| Kose et al. | A Review on Biomedical Signals with Fundamentals of Digital Signal Processing |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |