波谱学杂志, 2022, 39(3): 337-344 doi: 10.11938/cjmr20222975

研究论文

基于原子磁力计的穿戴式脑磁图动态测量研究

陈春巧1,2,3, 张欣2,3, 郭清乾2,3, 徐佳玉1,2, 冯晓宇2,3, 常严2,3, 胡涛,2,3, 杨晓冬,1,2,3

1. 长春理工大学 电子信息工程学院,吉林 长春 130022

2. 中国科学院苏州生物医学工程技术研究所,江苏 苏州 215163

3. 季华实验室,广东 佛山 528200

Moving Wearable Magnetoencephalography Measurement Study Based on Optically-pumped Magnetometer

CHEN Chun-qiao1,2,3, ZHANG Xin2,3, GUO Qing-qian2,3, XU Jia-yu1,2, FENG Xiao-yu2,3, CHANG Yan2,3, HU Tao,2,3, YANG Xiao-dong,1,2,3

1. School of Electronic and Information Engineering, Changchun University of Science and Technology, Changchun 130022, China

2. Suzhou Institute of Biomedical Engineering and Technology, Chinese Academy of Sciences, Suzhou 215163, China

3. Jihua Laboratory, Foshan 528200, China

通讯作者: 胡涛, Tel: 17625326300, E-mail:hutao@sibet.ac.cn杨晓冬, Tel: 18900616030, E-mail:xiaodong.yang@sibet.ac.cn

收稿日期: 2022-02-15  

基金资助: 苏州市基础研究试点项目.  SJC2021024
季华实验室项目.  X190131TD190
江苏省自然科学基金青年项目.  BK20200215

Received: 2022-02-15  

摘要

脑磁图作为一种无创的脑功能成像技术,依靠超高的时间及空间溯源分辨率,在脑科学研究和临床应用领域中有着极其重要的价值.本文介绍了自主搭建的基于原子磁力计的穿戴式脑磁图系统,通过设计匀场补偿线圈组并结合参考传感器阵列,实现被试头部运动区域内剩磁在±1 nT以内,保证动态测量过程中传感器输出维持在动态范围以内;同时提出了一种虚拟合成梯度去噪方法,显著抑制了环境共模噪声;最终在被试者头部自然运动状态下,成功检测到高信噪比的α节律信号与听觉诱发磁场信号,证实了该系统的有效性,为穿戴式脑磁图应用推广提供更多的可能性.

关键词: 穿戴式脑磁图 ; 原子磁力计 ; 动态测量 ; 匀场补偿 ; 噪声抑制

Abstract

Magnetoencephalography is a non-invasive technology for brain function imaging, which is of enormous value to brain science research and clinical application due to its ultra-high temporal and spatial trace resolution. In this paper, we introduce a self-built and atomic magnetometer based wearable magnetoencephalography system. By designing bi-planar coils system and combining with reference sensor array, the residual magnetic field in the subject's head movement area is controlled to be within ±1 nT, which ensures the sensors are maintained within their dynamic range during the moving measurement. At the same time, a virtual gradiometer-based noise reduction method is proposed to suppress the common-mode magnetic-field noise. Finally, the alpha rhythm and auditory evoked magnetic field signals with high signal-to-noise ratio are successfully detected under the subject's natural head movement and the effectiveness of the system is confirmed. This study could provide more possibilities for the application and promotion of moving wearable magnetoencephalography.

Keywords: wearable magnetoencephalography ; optically-pumped magnetometer ; moving measurement ; field nulling ; noise reduction

PDF (867KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

陈春巧, 张欣, 郭清乾, 徐佳玉, 冯晓宇, 常严, 胡涛, 杨晓冬. 基于原子磁力计的穿戴式脑磁图动态测量研究. 波谱学杂志[J], 2022, 39(3): 337-344 doi:10.11938/cjmr20222975

CHEN Chun-qiao. Moving Wearable Magnetoencephalography Measurement Study Based on Optically-pumped Magnetometer. Chinese Journal of Magnetic Resonance[J], 2022, 39(3): 337-344 doi:10.11938/cjmr20222975

引言

随着现代物理的迅速发展,近年来脑功能成像技术取得了巨大进步,功能磁共振成像(functional magnetic resonance imaging,fMRI)、脑电图(electroencephalography,EEG)、近红外光谱(near infrared spectrum,NIRS)、导航经颅磁刺激(navigated transcranial magnetic stimulation,nTMS)和脑磁图(magnetoencephalography,MEG)等无创脑功能成像技术相继诞生[1-3],极大地促进了脑科学领域深入化研究进程.其中大多数脑功能成像技术都能够提供毫秒范围的时间分辨率,而MEG技术同时也兼具良好的空间分辨率优势.作为一种高时空分辨的探测大脑神经元电活动所产生磁场的神经影像学技术,MEG在脑科学研究和临床诊断方面具有重要的研究价值和广泛的应用前景[4-6].当前商用MEG采用超导量子干涉仪(superconductor quantum interference devices,SQUID),具有1 fT/$\sqrt {{\rm{Hz}}} $的灵敏度[7, 8].但是其正常工作需要在4 K超低温环境下,以维持其超导状态,因此传感器需固定浸泡在液氦杜瓦中[9, 10],从而导致SQUID-MEG设备体积庞大且要求被试者在测量时保持静止状态.同时由于杜瓦的存在,传感器离被试头皮较远(~2 cm),脑磁信噪比受到限制,这些问题制约了SQUID-MEG的广泛应用[11, 12].

近年来,一种新型的量子弱磁传感器——无自旋交换弛豫(spin-exchange relaxation-free,SERF)原子磁力计(optically-pumped magnetometer,OPM)得到了快速发展[13-15],其利用碱金属原子光泵浦自旋极化和零场磁共振检测原理,通过抑制自旋交换弛豫,它无需工作在低温环境就能实现和SQUID相当的探测灵敏度,并且具有轻便、可小型化、探测距离近等优势[16].通过阵列形式将小型原子磁力计排布在脑磁测量帽插槽内可实现穿戴式脑磁图——OPM-MEG,并可在被试者自然运动状态下进行测量,大大提升了系统的灵活性和便利性[17, 18].2018年,英国诺丁汉大学团队首次利用OPM-MEG实现了被试者运动过程中β波(13~30 Hz)脑磁信号测量[19].2019年,Barry等[20]利用OPM-MEG实现了过程中海马体θ波(4~8 Hz)信号测量动态.2021年,Seymour等[21]利用OPM-MEG开展了动态过程中听觉诱发信号的检测.轻量化穿戴式脑磁图动态测量的实现在提升被试者体验感的同时扩大了其临床应用范围,尤其是对癫痫、多动症、帕金森病等无法控制自身活动的患者,弥补了传统脑磁图技术的不足.

本文详细介绍了自主搭建的基于原子磁力计的穿戴式脑磁图系统,为实现动态脑磁信号测量,设计了匀场补偿线圈组并结合参考传感器阵列,用于被试头部运动区域内剩磁补偿,保证动态测量过程中传感器输出维持在动态范围内.此外,研究了一种虚拟合成梯度去噪方法,用于抑制环境共模噪声.最后,设计了两组动态脑磁信号测量实验,用于验证系统的可靠性和有效性.

1 穿戴式脑磁图系统

1.1 系统整体介绍

本文自主搭建的基于原子磁力计的穿戴式脑磁图系统基本结构如图 1(a)所示,其中OPM传感器均使用美国QuSpin公司二代原子磁力计(QZFM Gen-2.0),该OPM最多支持双通道(z轴与y轴)脑磁信号采集,其测量方向如图 1(b)所示.系统主体在磁屏蔽室中,磁屏蔽室内部尺寸为1.75 m×1.95 m×2.23 m,可将地磁场屏蔽至10 nT以内,使传感器处于正常工作状态,同时利用左右两侧匀场补偿线圈组进一步补偿屏蔽房中心区域剩磁.被试头部佩戴柔性脑磁测量帽,表面共排布32个插槽用于放置OPM探测传感器阵列,测量帽可根据任意头型调整大小,并使传感器紧贴头皮. 4个OPMs作为参考传感器阵列固定在左右两侧支架上,用于测量屏蔽房内部剩磁及环境噪声. OPM参考传感器阵列和探测传感器阵列信号经屏蔽房外的OPM电子学系统,输出到控制采集系统,进行数据处理.被试接受脑磁诱发刺激模块产生的听觉、视觉或触觉等相关刺激[22],并由脑磁诱发刺激控制系统所控制.控制采集系统采用C++语言编写,以1 024 Hz采样率记录原始数据并与刺激序列保持同步.

图1

图1   (a) 自主搭建的基于原子磁力计(OPM)的穿戴式脑磁图仪系统整机图;(b) OPM传感器双轴测量方向

Fig.1   (a) Photograph of the self-built wearable magnetoencephalography system based on OPM; (b) Dual measurement axes of OPM sensor


1.2 匀场补偿

由于磁屏蔽室内部残余剩磁的大小约为10 nT,且在空间中存在较大的梯度场,在OPM-MEG动态实验过程中,头部的轻微晃动都有可能超出OPM传感器自身的动态范围(±1.5 nT),使其无法正确探测脑磁信号.因此,如何使用匀场补偿线圈组进一步补偿磁屏蔽室内部剩磁及梯度场至关重要[19, 23]. 本文基于平面型载流面[24]自主设计了7组对称补偿线圈组(3组匀场线圈BxByBz,4组一阶梯度线圈dBz/dx、dBz/dy、dBz/dz和dBy/dx),补偿中心区域为40 cm×40 cm×40 cm,7组补偿线圈张叠加粘贴在左右两侧固定支架上,线圈尺寸均为165 cm×140 cm,左右两侧线圈间距120 cm.7组补偿线圈系数实测值如表 1所示.

表1   补偿线圈系数实测值

Table 1  Measured coefficients of field nulling coils

匀场线圈一阶梯度线圈
BxByBzdBz/dxdBz/dydBz/dzdBy/dx
补偿线圈系数1.37 nT/mA1.84 nT/mA9.74 nT/mA11.62 nT/m/mA13.17 nT/m/mA52.24 nT/m/mA5.22 nT/m/mA

新窗口打开| 下载CSV


采用4个OPMs作为参考传感器阵列摆放在被试头部左右两侧,两侧间距为40 cm.具体位置如图 2(a)所示,4个OPMs参考传感器围成一个30 cm×30 cm×40 cm区域,可涵盖被试头部自然运动范围.OPMs参考传感器均工作在三轴直流测量(field-zero)模式,实时探测所在位置处三轴方向剩磁,参考传感器位置摆放包含各个方向梯度磁场(例如1号与2号参考传感器作差可得到x方向的梯度场),便于最优化补偿.匀场补偿需建立优化目标函数,目标函数F定义为:

$F = {\left\| {{\mathit{\boldsymbol{B}}_0} - \mathit{\boldsymbol{A}} \cdot \mathit{\boldsymbol{I}}} \right\|_2}$

图2

图2   (a) OPM参考传感器位置摆放示意图;(b)被试自然头动状态下,探测传感器所测磁场变化(细虚线:未施加匀场补偿;粗实线:施加匀场补偿)

Fig.2   (a) Configuration of the OPM reference sensors; (b) Variation of the measured magnetic field under the subject's natural head movement (thin dashed line: before field nulling; thick solid line: after field nulling)


其中,B0为4个OPMs参考传感器分别测量到所在位置处的三轴背景磁场所组成的12×1的矩阵;A为12×7的补偿线圈系统矩阵,其通过计算分别给7组补偿线圈施加特定的电流,记录每个电流值下4个OPMs参考传感器所测得12个磁场值,再利用线性拟合得出电流和磁场间比例系数,组成12×7的线圈系统矩阵;I为施加给补偿线圈的7×1电流值矩阵;${\left\| \sim \right\|_2}$表示矩阵的二阶范数.匀场补偿通过施加补偿电流,使目标函数达到最小.根据目标函数定义,补偿电流值计算可采用最小二乘法,即

$\mathit{\boldsymbol{I}} = {\mathit{\boldsymbol{A}}^ + } \cdot {\mathit{\boldsymbol{B}}_0}$

其中,A+为线圈系统矩阵的伪逆.实际补偿过程中,将参考传感器实际测得磁场值输入补偿优化模型,计算得出7组匀场线圈的补偿电流值,通过串口通信将电流值发送给多通道电流源,并进行多次迭代优化,直到补偿效果达到理想状态.

为测试匀场补偿效果,我们将一个OPM探测传感器插入脑磁测量帽中,传感器工作在单轴弱磁测量模式.被试头部在30 cm×30 cm×40 cm区域内各个方向自由晃动,测量补偿前后剩磁变化,结果如图 2(b)所示.可以看出,未施加匀场补偿情况下,被试头动极易导致被测磁场(细虚线)超出传感器动态范围(±1 500 pT)达到饱和,无法正常进行脑磁测量;而在完成匀场补偿后(粗实线),被试头部任意晃动,被测磁场仍在±1 000 pT(即±1 nT)以内,传感器可正常工作.

1.3 虚拟合成梯度去噪

环境噪声的抑制将有利于获取高信噪比的脑磁信号,本文提出一种虚拟合成梯度去噪方法用于环境共模噪声抑制.该系统配置由脑磁测量帽上的OPMs探测传感器阵列和远离头皮的参考传感器阵列组成,参考传感器阵列摆放位置与图 2(a)保持一致.在完成匀场补偿后,4个参考传感器工作在单轴弱磁测量模式,虚拟合成梯度去噪具体方法如下:

${\mathit{\boldsymbol{D}}_r} = \mathit{\boldsymbol{D}} - \mathit{\boldsymbol{b}} \cdot \mathit{\boldsymbol{R}}$

$\mathit{\boldsymbol{D}} = \left\{ {{D_1}, {D_2}, {D_3} \cdots {D_i} \cdots {D_N}} \right\}$

$\mathit{\boldsymbol{R}} = \left\{ {{R_1}, {R_2}, {R_3}, {R_4}} \right\}$

其中,Dr为去除环境噪声后的OPMs探测传感器阵列信号矩阵. D为去噪前探测传感器阵列信号N×T的矩阵,N为探测传感器通道数;T为时间点数,也是一维矩阵DiRi的长度. ${{D_1}, {D_2}, {D_3} \cdots {D_i} \cdots {D_N}}$分别表示1~N通道探测传感器数据.R为参考传感器阵列所测环境噪声4×T的矩阵.R1, R2, R3R4分别表示各参考传感器所测环境噪声分量.bN×4的权重系数矩阵,取决于探测传感器所在位置及所测磁场方向.权重系数矩阵b计算同样采用最小二乘法,表示如下:

$\mathit{\boldsymbol{b}} = \mathit{\boldsymbol{D}} \cdot {\mathit{\boldsymbol{R}}^ + }$

其中,R+为参考传感器信号矩阵的伪逆.

为验证虚拟合成梯度去噪效果,我们将脑磁测量帽固定在头模上,如图 3(a)所示,并使用3通道OPMs探测传感器任意插入测量帽插槽中,在空屏蔽房中进行环境噪声测量.头模高度及位置与脑磁实验被试头部接近.实验中探测传感器阵列与4个参考传感器阵列均工作在单轴测量模式(z轴),探测传感器测量轴垂直头皮方向,测量结果如图 3(b)所示.从图中可以看出,3通道OPMs探测传感器所测原始本底噪声,即粉红色功率谱密度(power spectral density,PSD)曲线,维持在0.03~0.04 pT/$\sqrt {{\rm{Hz}}} $;经过去噪修正后,本底噪声(蓝色PSD曲线)可降低至0.01~0.02 pT/$\sqrt {{\rm{Hz}}} $,接近传感器自身灵敏度极限.而50 Hz工频噪声干扰也从约50 pT/$\sqrt {{\rm{Hz}}} $大幅降低至0.8 pT/$\sqrt {{\rm{Hz}}} $以下,降噪比达到62以上.此外,在25 Hz附近也具有一定的降噪效果.然而,对于接近75~80 Hz区间的OPM传感器固有电子学噪声干扰,去噪结果并不明显,因为传感器固有电子学噪声干扰为非共模噪声.

图3

图3   (a) 头模照片;(b)空屏蔽房中,3通道传感器的虚拟合成梯度去噪结果

Fig.3   (a) The head model; (b) Virtual gradiometer-based noise reduction results of 3-channel OPMs in the empty shielded room


2 脑磁图动态测量

为检验穿戴式脑磁图系统的可靠性和实用性,我们开展了两组被试者自然头动状态下的脑磁测量实验,分别是α节律光刺激和听觉诱发刺激,每组实验被试者4人,实验设计均相同,α节律实验重复三次,听觉实验重复两次,本文展示其中一名被试实验结果[其余三名被试者实验结果见图S1和图S2(扫描文章首页二维码或在论文网页版查看)].在两组实验中,被试头部在30 cm×30 cm×40 cm区域内自由摆动.脑磁图动态测量实验开始前,参考OPMs传感器阵列在三轴直流测量模式下工作用于匀场补偿.补偿完毕后,补偿线圈组电流保持恒定;此时,参考传感器阵列切换到单轴弱磁测量模式(z轴)测量环境噪声,用于梯度去噪,梯度去噪覆盖整个动态脑磁测量过程并与探测传感器保持同步.

对于α节律光刺激实验,5通道OPMs探测传感器配置在被试者的枕骨区域,因为α节律(8~13 Hz范围)信号源主要位于枕骨区[4, 25, 26].光刺激任务设计如图 4所示,被试者通过空气耳机指示在整个实验中每5 s交替睁眼和闭眼.传感器测量得到的全频段时域信号经过虚拟合成梯度去噪及8~13 Hz带通滤波后,结果如图 5(a)所示,灰色和黄色区域代表睁眼和闭眼的状态.从图中可以看出5通道时域数据均有明显的强弱交替的α节律信号变化,当被试者睁眼时,信号幅度明显减弱.各通道间信号幅度强弱则取决于OPM传感器与颅内信号源相对距离,越靠近信号源则幅度越强.

图4

图4   α节律光刺激实验设计

Fig.4   Experimental design of α-rhythm light stimulation


图5

图5   (a) 被试者自然头动状态下5通道α节律信号时域图;(b)被试自然头动状态下10通道听觉诱发磁场信号时域图

Fig.5   (a) Time-domain of the 5-channel α-rhythm signals under the subject's natural head movement; (b) Time-domain of the 10-channel auditory evoked magnetic field signals under the subject's natural head movement


听觉信号源主要位于大脑皮层颞叶,因此听觉诱发磁场信号测量的实验中将10通道OPMs探测传感器放置在被试者头部的左侧和右侧颞叶区域,两侧各5通道.被试者佩戴空气耳机听取来自刺激装置产生的声音刺激,实验总共施加400次听觉刺激,持续时间为0.3 s,刺激时间间隔为1.7 s.同时,听觉刺激频率随机呈现1 000 Hz或1 200 Hz,以保证被试不易产生听觉疲劳,其中1 000 Hz刺激的数量占刺激总数的80%.听觉诱发磁场信号所处波段为2~40 Hz,我们同样采用虚拟合成梯度去噪方法对原始信号进行噪声抑制并使用带通滤波截取2~40 Hz范围内信号,对400个听觉刺激信号做叠加取平均,得到如图 5(b)所示的听觉诱发磁场时域信号,正负信号反映两侧脑磁信号方向相对传感器测量方向相反.从图中可以看出,在刺激产生大约100 ms后,10通道OPM传感器均出现明显的磁场强度峰值,该峰值即为M100信号.M100是一个典型的听觉刺激诱发反应尖峰信号,产生于听觉刺激发生100 ms以后,该反应是一个瞬态过程,后续不再有明显反应[9, 27].

两组实验均表明该系统可实现高信噪比的脑磁图动态测量,证实了该系统的有效性和实用性.

3 结论

本文详细描述了自主搭建的基于原子磁力计的穿戴式脑磁图动态测量系统,通过配置参考传感器阵列和匀场补偿线圈组,结合相应的补偿算法,可实现30 cm×30 cm×40 cm区域内剩磁在±1 nT以内,保证动态测量过程中传感器输出维持在动态范围以内.同时本文提出了一种虚拟合成梯度去噪方法,实现了环境共模噪声的有效抑制,使信号本底噪声可降低至传感器灵敏度附近.最后设计了两组实验用于验证穿戴式脑磁图系统动态测量性能,结果表明被试在头部自然运动状态下,可以成功检测到高信噪比的α节律信号与听觉诱发磁场信号,证实了该系统的有效性.该研究为穿戴式脑磁图在癫痫、多动症、帕金森病等患者动态测量的应用推广提供了更多的可能性.

利益冲突


参考文献

YANG L Q , LIN F C , LEI H .

Resting state functional connectivity in brain studied by fMRI approach

[J]. Chinese J Magn Reson, 2010, 27 (3): 326- 340.

DOI:10.3969/j.issn.1000-4556.2010.03.005      [本文引用: 1]

杨丽琴, 林富春, 雷皓.

静息状态下脑功能连接的磁共振成像研究

[J]. 波谱学杂志, 2010, 27 (3): 326- 340.

DOI:10.3969/j.issn.1000-4556.2010.03.005      [本文引用: 1]

CHENG L W , WANG L L , ZHONG K .

Application of fMRI in transcranial direct current stimulation researches

[J]. Chinese J Magn Reson, 2020, 37 (4): 533- 546.

URL    

程力维, 王璐璐, 钟凯.

fMRI在经颅直流电刺激研究中的应用进展

[J]. 波谱学杂志, 2020, 37 (4): 533- 546.

URL    

VARONE G, HUSSAIN Z, SHEIKH Z, et al. Real-time artifacts reduction during TMS-EEG co-registration: A comprehensive review on technologies and procedures[J]. Sensors (Basel, Switzerland), 21(2): 637.

[本文引用: 1]

COHEN D .

Magnetoencephalography: detection of the brain's electrical activity with a superconducting magnetometer

[J]. Science, 1972, 175 (4022): 664- 666.

DOI:10.1126/science.175.4022.664      [本文引用: 2]

COHEN .

Magnetoencephalography: evidence of magnetic fields produced by alpha-rhythm currents

[J]. Science, 1968, 161 (3843): 784- 786.

DOI:10.1126/science.161.3843.784     

GRATTA C D , PIZZELLA V , TECCHIO F , et al.

Magnetoencephalography - a noninvasive brain imaging method with 1 ms time resolution

[J]. Rep Prog Phys, 2001, 64 (12): 1759- 1814.

DOI:10.1088/0034-4885/64/12/204      [本文引用: 1]

VRBA J . Multichannel SQUID biomagnetic systems. applications of superconductivity[M]. Springer Netherlands, 2000, 61- 138.

[本文引用: 1]

CHEN M J , LIAO S X , YANG H C , et al.

Nuclear magnetic resonance and imaging of hyperpolarized 3He using high-Tc superconducting quantum interference device in microtesla magnetic fields

[J]. Chinese J Magn Reson, 2010, 27 (3): 386- 395.

DOI:10.3969/j.issn.1000-4556.2010.03.010      [本文引用: 1]

陈名杰, 廖书贤, 杨鸿昌, .

采用超导量子干涉组件在微特斯拉磁场下获取超极化3He的核磁共振波谱和影像

[J]. 波谱学杂志, 2010, 27 (3): 386- 395.

DOI:10.3969/j.issn.1000-4556.2010.03.010      [本文引用: 1]

HÄMÄLÄINEN M , HARI R , ILMONIEMI R J , et al.

Magnetoencephalography—theory, instrumentation, and applications to noninvasive studies of the working human brain

[J]. Rev Mod Phys, 1993, 65 (2): 413- 497.

DOI:10.1103/RevModPhys.65.413      [本文引用: 2]

WEINSTOCK H .

SQUID sensors: fundamentals, fabrication and applications

[J]. Springer Netherlands, 1996,

[本文引用: 1]

VRBA J .

Magnetoencephalography: The art of finding a needle in a haystack

[J]. J Psychophysiol, 2003, 17 (4): 237- 237.

[本文引用: 1]

MUKAMEL R , GELBARD H , ARIELI A , et al.

Coupling between neuronal firing, field potentials, and FMRI in human auditory cortex

[J]. Science, 2005, 309 (5736): 951- 954.

DOI:10.1126/science.1110913      [本文引用: 1]

WANG X F , SUN X P , ZHAO X C , et al.

Progress in biomagnetic signal measurements with ultra-sensitive atomic magnetometers

[J]. Chinese Journal of Lasers, 2018, 45 (2): 0207012.

URL     [本文引用: 1]

王晓飞, 孙献平, 赵修超, .

超灵敏原子磁力计在生物磁应用中的研究进展

[J]. 中国激光, 2018, 45 (2): 0207012.

URL     [本文引用: 1]

MA H F , WU Y T , ZHAO W , et al.

Research progress of magnetoencephalography in the functional mechanism of bilingual brain

[J]. Chinese Journal of Biomedical Engineering, 2021, 40 (4): 477- 484.

DOI:10.3969/j.issn.0258-8021.2021.04.12     

马恒芬, 吴云涛, 赵文, .

双语脑功能机制的脑磁图研究进展

[J]. 中国生物医学工程学报, 2021, 40 (4): 477- 484.

DOI:10.3969/j.issn.0258-8021.2021.04.12     

XU W J , JIANG M , PENG X H .

Study on ultra-low-field nuclear magnetic resonance spectroscopy based on high-sensitivity atomic magnetometer

[J]. Journal of University of Science and Technology of China, 2020, 50 (8): 1138- 1143.

URL     [本文引用: 1]

徐文杰, 江敏, 彭新华.

基于高灵敏度原子磁力计的超低场核磁共振谱学研究

[J]. 中国科学技术大学学报, 2020, 50 (8): 1138- 1143.

URL     [本文引用: 1]

ZHANG S L , CAO N .

A synthetic optically pumped gradiometer for magnetocardiography measurements

[J]. Chinese Physics B, 2020, 29 (4)

[本文引用: 1]

BOTO E , MEYER S S , SHAH V , et al.

A new generation of magnetoencephalography: Room temperature measurements using optically-pumped magnetometers

[J]. Neuroimage, 2017, 149, 404- 414.

DOI:10.1016/j.neuroimage.2017.01.034      [本文引用: 1]

ZHANG X , CHEN C Q , ZHANG M K , et al.

Detection and analysis of MEG signals in occipital region with double-channel OPM sensors

[J]. J Neurosci Methods, 2020, 346, 108948.

DOI:10.1016/j.jneumeth.2020.108948      [本文引用: 1]

BOTO E , HOLMES N , LEGGETT J , et al.

Moving magnetoencephalography towards real-world applications with a wearable system

[J]. Nature, 2018, 555 (7698): 657- 661.

DOI:10.1038/nature26147      [本文引用: 2]

BARRY D N , TIERNEY T M , HOLMES N , et al.

Imaging the human hippocampus with optically-pumped magnetoencephalography

[J]. Neuroimage, 2019, 203, 116192.

DOI:10.1016/j.neuroimage.2019.116192      [本文引用: 1]

SEYMOUR R A , ALEXANDER N , MELLOR S , et al.

Using OPMs to measure neural activity in standing, mobile participants

[J]. Neuroimage, 2021, 244, 118604.

DOI:10.1016/j.neuroimage.2021.118604      [本文引用: 1]

SUN W , WANG H , ZHANG Y , et al.

Optimal design for quantification of gas concentration based olfactory stimulator

[J]. Chinese J Magn Reson, 2021, 38 (1): 12- 21.

URL     [本文引用: 1]

孙韦, 王慧, 张寅, .

基于气体浓度定量的嗅觉刺激器优化设计

[J]. 波谱学杂志, 2021, 38 (1): 12- 21.

URL     [本文引用: 1]

IIVANAINENJ, ZETTERR, GRONM, .

On-scalp MEG system utilizing an actively shielded array of optically-pumped magnetometers

[J]. Neuroimage, 2019, 194, 244- 258.

URL     [本文引用: 1]

HOLMES N , LEGGETT J , BOTO E , et al.

A bi-planar coil system for nulling background magnetic fields in scalp mounted magnetoencephalography

[J]. Neuroimage, 2018, 181, 760- 774.

DOI:10.1016/j.neuroimage.2018.07.028      [本文引用: 1]

KAWABATA N .

Nonstationary power spectrum analysis of the photic alpha blocking

[J]. Kybernetik, 1972, 12 (1): 40- 44.

DOI:10.1007/BF00289235      [本文引用: 1]

BRENNER D , WILLIAMSON S J , KAUFMAN L .

Visually evoked magnetic fields of the human brain

[J]. Science, 1975, 190 (4213): 480- 482.

DOI:10.1126/science.170683      [本文引用: 1]

JOHNSON C N , SCHWINDT P D , WEISEND M .

Multi-sensor magnetoencephalography with atomic magnetometers

[J]. Phys Med Biol, 2013, 58 (17): 6065- 6077.

DOI:10.1088/0031-9155/58/17/6065      [本文引用: 1]

/