大学生文明上网-大学生文明上网
附件
1
华中科技大学博士学位论文创新点概述
论文题目:提高全脑
fNIRS- EEG
双模态联合检测性能的研究
学科(专业)
、研究方向:生物医学光子学
论文中的主要创新点(不够可附页)
(
1
)将小间距去除表层干扰技术应用到光纤
fNIRS
系 统中,并通过在线自适应滤波首次实现了在实
验中实时去除表层干扰。
测 试结果表明该方法可以有效抑制
fNIRS
信号中的表层干扰和全局噪声,
在
实验中可以观察到不受心跳等噪声干扰的血氧信号。
(
2
)通过研制光电一体式探头帽与研究光电打标技术实现 p>
fNIRS
与
EEG
在物理空间和信号层面的联
合。通过混合编程与远程数据共享技术,实现基于
fNIRS-EEG
双模态信号的在线神经反馈,并建立了
闭环的全脑
fNIRS- EEG
双模态联合检测系统。
(
3
)首次使用全脑
fNIRS-EEG
双模态检测技术以及脑连接的分析方法研究工作记忆中脑功能的性别
差异
。发展了基于
fNIRS-EEG
双模态成像技术的脑连接数据处理方法。结果表 明,脑连接技术的应用
在一定程度上解决了工作记忆中性别相关的脑功能偏侧化的争议。
另外,
fNIRS
和
EEG
信号的联合分
析表明男性和女性在工作记忆中使用不同的策略来编码记忆。
博士学位论文主要研究成果的发表或获奖情况
序号
发表论文题目
/
获奖成果名称
刊物名称
/
获奖部门
刊物
/
奖
励级别
发表论文
/
获奖时间
作者署名
名次
Gender differences in brain
networks during verbal
Sternberg tasks: A
1
Journal of
Biophotonics
A
类
2017
年
9
月
1
simultaneous near-infrared
spectroscopy and
electro
‐
encephalography study.
PFC Activity Pattern
During
Verbal WM Task in Healthy
2
Advances in
Experimental
Medicine and
Biology
Journal of
Innovative Optical
Health
Sciences
B
类
2016
年
9
月
1
Male and Female Subjects: A
NIRS Study
Detecting
bilateral
functional connectivity in
3
B
类
2013
年
6
月
3
the prefrontal cortex during
a stroop task by
near-
infrared spectroscopy
The
Role of Phonological
Processing in
Semantic
Advances in
Experimental
Medicine and
Biology
B
类
2016
年
9
月
3
4
Access of Chinese
Characters: A Near-Infrared
Spectroscopy Study
Extracting heartrate from
optical signal of functional
Journal of
Innovative
Optical
Health Sciences,
B
类
2017
年
11
月
3
5
near-infrared spectroscopy
based on mathematical
morphology.
分
类
号
学号
D201277679
学校代码
10487
密级
博士学位论文
提高全脑
fNIRS- EEG
双模态联合检测
性能的研究
p>
学
位
申
请
人
:
高晨阳
学
科
专
业
:
生物医学工程
指
导
教< /p>
师
:
龚辉
教授
答
辩
日
< p>期:
Dissertation Submitted in Partial
Fulfillment of the Requirements
for the
Degree of Doctor of Philosophy in Engineering
Research on the Improvement of Whole
Brain
fNIRS-EEG bimodal Detection
Performance
Ph.D. Candidate: Gao Chenyang
Major
: Biomedical Engineering
Supervisors
: Prof. Gong Hui
Huazhong University of Science and
Technology
Wuhan, Hubei 430074, P. R.
China
March 1, 2018
华
中
科
技
大
学
博
士
学
位
论
文
摘
要
神经成像技术为人们提供了较好地了解大脑如何工作
的手段。功能近红外光谱
术(
functional
Near-infrared
Spectroscopy
, p>
fNIRS
)作为一种无创的光学成像技术,
已越来越多地应
用于脑功能研究中。
fNIRS
通过通道的排布位置可以确定脑功能信号
的来源,
而
fNIRS
技术本身存在时间分辨率不 够高、
无法直接反映神经活动等问题,
这
些
问
题
严
重
影
响
了
fNIRS
用
于
脑
功
能
检
测
的
可
靠
性
。
脑
电
图
(Electr
oencephalo-graphy , EEG)
作为一种较为成熟的脑功能检测技术,具有 较高的时
间分辨率,其缺点是检测信号无法确定来源的脑区。
fNIRS
和
EEG
技术都具有对环
境和被试限制小
的优势,可以在自然和放松的状态下检测被试的脑功能活动。因此
将
fN IRS
与
EEG
技术进行联合,
使用双模态检测技 术可以综合二者空间和时间分辨
率的优势,帮助人们更好地认识在认知心理学任务中脑活
动的神经机制。目前
fNIRS-EEG
双模态检测技术尚处于起步阶段 ,存在信号受噪声干扰严重、系统整合
性较差、联合通道少等问题。提高
fNIRS-EEG
双模态联合检测的性能是将该技术应
用到脑功能研究
的重要内容。
本文围绕
fNIRS
在光 电联合检测中的方法和应用展开研究,主要内容如下:
(
1
)
围绕课题组研制的光纤式
fNIRS
系统进 行提高全脑
fNIRS
检测性能的研究。
首先研究了影响
fNIRS
信号的表层干扰问题,通过理论分析和模型实验建立了小间
< p>距配合自适应滤波的处理方法,然后在体实验验证该方法的有效性,并开发了插件
< br>嵌入到采集软件中实现在线表层干扰的去除。通过理论分析和仿体测试分析影响
f
NIRS
系统稳定性的具体问题,围绕这些问题设计了光源机箱散热模块和光纤探头,
< p>提高了系统在长时间工作的稳定性。针对
fNIRS
采集数据可视化 问题,对系统软件
进行需求分析,使用
LabVIEW
和
Matlab
混合编程实现了高可视化程度的
fNIRS
系
统软件。
最后通过在体
N-back
< p>实验验证了fNIRS
系统可以可靠地检测到工作记忆中
的脑功能活动。
*
*
本研究得到国家 自然科学创新群体项目(批准号:
61421064
)和武汉光电国家实验室主任 基金的资助。
I
华
中
科
技
大
学
博
士
学
位
论
文
(
2
)在物理结构和数字信号层面实现
fNIRS
和
EEG
的系统联合,提高双模态
联合检测的性能。在物理结构联合层面,通过特制的布帽
基底和探头夹具制作光电
联合探头帽将多通道的
fNIRS
探头以及
EEG
电极固定在一起。在信号同步层面,研
究使用光电打标的方式实现高精度的信号同步触发。围绕双模态技术在脑机接口等
领域的
应用对
fNIRS
和
EEG
信号的在线神经反馈功 能开展研究,
通过运动想象实验
证明神经反馈的有效性。
通过远程数据共享技术实现了基于
LabVIEW
平台的闭环双
< br>模态联合采集软件。最后通过一个
Stroop
实验证明全脑双模态检测技 术可以有效地
检测到脑功能血氧活动和神经电活动。
< br>(
3
)为了促进全脑双模态联合检测技术在脑功能研究中的应用,以工作记 忆任
务下的脑功能性别差异为例进行研究。并针对双模态的数据分析方法较为单一的问<
/p>
题,建立了基于双模态信号的脑网络数据分析方法。使用包括功能连接、效应连接
以及小世界网络特性在内的多种数据分析方法来分析
fNIRS
和
EEG
数据,
并分析两
种信号之间的相关
性。
fNIRS
脑激活和脑连接的结果表明女性更加擅长语言类工作记
< p>忆,而男性需要更强的脑激活、更优化的脑网络以及更擅长的策略来达到和女性相
< br>似的表现。该研究进一步证明基于双模态联合检测的脑连接分析方法可以帮助人们
更全面有效地研究脑功能。
本文围绕信号质量、双模态系统整合、软件
可视化以及数据分析方法等多个方
面对提高双模态联合检测的性能进行研究。本文的工作
可以帮助研究人员使用双模
态联合检测技术更加深入地了解脑功能活动的神经机制。
p>
关键词:
功能近红外光谱术
双模态检测
表层干扰
神经反馈
脑网络
工作记忆
II
华
中
科
技
大
学
博
士
学
位
论
文
Abstract
*
Functional neuroimaging
methods provide a new method for learning how the
brain
works.
As
a
noninvasive
optical
imaging
method,
functional
near- infrared
spectroscopy
(fNIRS)
has been widely employed in studies of brain
function.
By the channels’ location
on the
scalp, fNIRS can detect neuro signals’ source from
the cortex. However, fNIRS has
problems
such
as
insufficient
time
resolution
and
can
’t
detect
neuron
activity
directly,
which affect
the effectiveness
of
fNIRS used in
brain
functional detection.
As a mature
brain functional detection technology,
Electroencephalography (EEG) has sufficient time
resolution, while it can not detect the
source of neuroactivity. Both fNIRS and EEG place
few
limita
tions
on
the
subjects
and
environment,
therefore
can
detect
subjects’
brain
activity in a relaxed state. By
combining fNIRS and EEG, bimodal acquisition
technology
can
have
complementary
advantages
of
the
two
methods,
help
people
to
study
neural
mechanism in cognitive psychology
experiments better. The current fNIRS-EEG bimodal
technology
is
still
in
its
infancy,
has
problems
such
as
low
signal-to-noise
ratio,
poor
system integration, and few detection
channels. Improving fNIRS-EEG bimodal detection
performance is very important for
employing this technology in brain function
research.
This
thesis
focused
on
the
improvement
of
the
fNIRS-EEG
bimodal
detection
performance, and the main contents are
as follows:
(1) Studying key issues
affecting whole brain fNIRS detection performance
based on
our homemade optical fiber
fNIRS system. First, study superficial signals
which affect the
fNIRS
signal
quality.
By
theoretical
analysis
and
model
experiments,
this
study
established a method based on short
separation channels with adaptive filtering to
remove
superficial interference and
verified the effectiveness
of the method in
vivo experiment.
Then this study
developed on-line superficial interference
removing function embedded in
*
This study was supported
by Science Fund for Creative Research Group (Grant
No. 61421064) and Director Fund of
WNLO
III
p>
华
中
科
技
大
学
博
士
学
位
论
文 p>
fNIRS software. Second, this
study analyzed key issues affecting fNIRS system
stability
by
theoretical
analysis
and
phantom
tests.
To
solve
the
issues,
this
study
designed
light
sources
heat dissipation system and fiber-optics probe.
Then, focused on the visualization
of
fNIRS acquisition software, this study analyzed
software requirements and developed
fNIRS
acquisition
software
with
high
visualization
by
Matlab- LabVIEW
hybrid
programming.
At
last,
through
an
in
vivo
N-back
experiment
verified
that
the
fNIRS
system can
detect brain functional activity during working
memory reliably.
(2) Combining
fNIRS and EEG system in both physical structure
and digital signal
level,
improving
bimodal
detection
performance.
In
physical
structure
level,
this
study
developed multi-channels
optic-electronic integrated probe cap using the
special material
basement and probe
fixtures to achieve the joint arrangement of fNIRS
probes and EEG
electrodes.
In
signal
synchronization
level,
this
study
developed
photoelectric
marking
method
to
achieve
simultaneous
bimodal
detection
with
high
precision.
Focused
on
bimodal
detection utilized in
brain-computer-interface, this study did
research on online
neural
feedback
based
on
fNIRS
and
EEG
signals,
and
proved
the
efficiency
of
neural
feedback by a motor imagery experiment.
Then this study developed closed-loop bimodal
acquisition software based on LabVIEW
by remote data access method. At last, through
an in
vivo Stroop experiment
verified that the bimodal system
can detect hemodynamic
response and
neuroelectrical activity simultaneously and
effectively.
(3) In order to promote
the application of the whole brain bimodal
detection, studying
gender-related
brain
function
during
working
memory.
Focused
on
few
bimodal
data
analysis
methods,
this
study
established
brain
networks
analysis
methods
based
on
the
bimodal signal. Multiple data analysis
methods including functional connectivity,
effective
connectivity,
and
small
world
network
property
were
employed
to
analyze
fNIRS
and
EEG data, as well as the correlation
between two types of signal. Results of brain
activity
and
connectivity
indicated
that
women
surpassed
men
in
verbal
working
memory.
Men
achieve
similar
performance
to
women
by
enhancing
brain
activation,
optimizing
brain
IV
p>
华
中
科
技
大
学
博
士
学
位
论
文 p>
networks
and
employing
preferred
visuospatial
strategies
to
encode
memory
in
the
memory tasks. This study further proved
that brain networks analysis based on bimodal
detection
can
help
people
to
extract
the
information
of
brain
function
more
comprehensively and effectively.
This study devoted to improving the
performance of fNIRS-EEG bimodal detection
in
various
aspects
including
signal
quality,
bimodal
system
integration,
software
visualization
and
data
analysis
methods.
This
study
can
help
researchers
to
employ
bimodal
detection technology study brain function in more
details.
Keywords:
Functional near-infrared spectroscopy
Superficial interference
Working memory
Bimodal detection
Neuro-feedback
Brain networks
V
华
中
科
技
大
学
博
士
学
位
论
文
本文所用英文缩写名称一览表
fNIRS
HbO
2
Hb
EEG
ERP
MEG
fMRI
functional Near-infrared
Spectroscopy
Oxy-hemoglobin
功能近红外光谱术
氧合血红蛋白
脱氧血红蛋白
脑电图
事件相关电位
脑磁图
功能磁共振成像
血氧水平依赖
正电子发射断层显像
光密度值
差分光程因子
数值孔径
美国国家仪器
发光二极管
激光二极管
通用串行总线
雪崩光电二极管
模数转换器
均方根值
远程数据访问
前额叶
语言类工作记忆
小波相干
影响差值
支持向量机模型
Deoxy-
hemoglobin
Electroencephalography
Event-related Potential
Magnetoencephalography
functional Magnetic Resonance Imaging
Blood Oxygenation Level Dependent
Positron Emission Tomography
Optical Density
BOLD
PET
OD
DPF
NA
NI
LED
LD
USB
APD
ADC
RMS
RDA
PFC
Differential
Path Length Factor
Numerical Aperture
National Instrument
Light Emitting Diode
Laser
Diode
Universal
Serial Bus
Avalanche Photodiode
Analog to Digital Converter
Root Mean Square
Remote Data Access
Prefrontal Cortex
VWM
WTC
DOI
SVM
Verbal
Working Memory
Wavelet Transform Coherence
Differences of Influence
Support Vector Machine
VI
华
中
科
技
大
学
博
士
学
位
论
文
MUA
Mass Univariate Analyses
大量单变量分析
VII
p>
华
中
科
技
大
学
博
士
学
位
论
文 p>
目
录
摘
要
.
................................................. .................................................. ............. I
Abstract
.
................................................. .................................................. ..... III
英文缩写名称一览表
............ .................................................. ..................... VI
1
绪论
1.1
引言
< /p>
.
.................................... .................................................. .....................
1
1.2
脑功能研究概述
............. .................................................. ........................
1
1.3
近红外脑功能研究
............ .................................................. .....................
6
1.4
近红外
-
脑电双模态联合检测研究
................................................. .........
9
1.5
本文的主要内容
............................ .................................................. .......
11
2
提高全脑
fNIRS
检测性能的研究
2.1
引言
.
... .................................................. .................................................. ..
1
3
2.2
使用小间距通道去除表层干扰
...................... .......................................
1
4
2.3
提高
fNIRS
系统稳定性关键问题研究
.
...................................... ..........
2
8
2.4
混合编程提高
fNIRS
系统软件可视化程度
< /p>
.
.................................... ....
3
5
2.5
基于
fNIRS
的
N-back
任务下脑连接研究
...........................................
4
3
2.6
本章小结
.
................................... .................................................. ............
5
0
3
提高全脑双模态联合检测性能的研究
VIII
华
中
科
技
大
学
博
士
学
位
论
文
3.1
引言
.
... .................................................. .................................................. ..
5
1
3.2
双模态系统的空间融合与信号同步
.................... .................................
5
2
3.3
集成神经反馈的闭环双模态联合系统软件
.........................................
6
0
3.4
双模态联合检测的在体
Stroop
实验测试
< /p>
.
.................................... ........
7
1
3.5
本章小结
.
................ .................................................. ...............................
7
4
4
工作记忆中脑功能性别差异的双模态联合检测研究
4.1
引言
.
... .................................................. .................................................. ..
7
6
4.2
双模态联合检测实验方法
........................ .............................................
7
8
4.3
脑功能激活和连接结果
......................... ................................................
8
5
4.4
脑功能性别差异讨论
.......................... .................................................. .
9
1
4.5
本章小结
.
................ .................................................. ...............................
9
5
5
总结与展望
5.1
p>
本文工作总结
................................ .................................................. ........
9
7
5.2
本文工作的创新点
........................... .................................................. ....
9
8
5.3
工作展望
.
................ .................................................. ...............................
9
8
致
谢
.
................................................. .................................................. .........
1
00
参考文献
.................................................. .................................................. ...
1
01
附录
作者攻读博士学位期间发表论文与研究
........................................
117
IX
华
中
科
技
大
学
博
士
学
位
论
文
1
绪论
1.1
引言
大脑是生命诞生以来,生物经过数亿年漫长岁月进
化的结果,是高度智能化的
复杂系统。一直以来人们对于大脑的研究从未停止。早期主要
通过行为研究来了解
个体的大脑活动,随着科技发展,越来越多的神经影像技术可以帮助
我们更加全面
地认识脑功能活动的规律和神经机制。对大脑功能的研究不仅可以帮助我们
探索人
类行与思的本质关系,帮助我们理解社会现象,指导我们对精神类疾病的治疗,更
可以促进生物科学、信息科学、人工智能等多种学科技术的交叉发展。
脑功能研究的重要性正逐渐被全世界认可。随着欧盟和美国宣布开展
“< /p>
人类大脑
计划
”
和
“
脑活动图谱计划
”
,研究脑可以帮助我们有效地治疗各种类型的神经系统 疾
病。目前自闭症、老年痴呆、抑郁症、精神分裂症等神经与精神疾病给个人、家庭
p>
以及社会带来了极大的负担。
为了对这些疾病进行预防和治疗,
需要从分子、
细胞、
网络和脑区各个层次研究和了解大脑。其中,以人
类大脑局部脑区或者大脑整体作
为研究对象的无损脑功能检测技术具有重要的科学研究和
临床应用价值。在众多的
无损脑功能成像技术中,
近红外光谱术
(
functional Near-infrared Spectroscopy, fNIRS
)
作为一种新兴的功能成像技术已被广泛地应用于脑科学研究
中。该技术具有对环境
和被试限制小、使用方便、对人体无任何损伤等优势
[1-3]
,特别适合应用于儿童,老
年人及患者进
行脑功能活动检测。因此本文的
fNIRS-EEG
双模态检测技术研究工作 p>
也是围绕
fNIRS
技术开展。
< p>
绪论部分首先对脑功能研究进行概述,包括脑活动相关的生理响应、脑功能检
测技术等。接着介绍本文的主要研究工具
—
近红外光谱术,以及在此基础上发展 的
fNIRS-EEG
双模态成像技术的研究概况和存在的问题。最后围 绕这些问题介绍本文
的主要研究内容。
1.2
脑功能研究概述
1.2.1
脑功能活动的生理机制
1
p>
华
中
科
技
大
学
博
士
学
位
论
文 p>
大脑是人体高级神经中枢所在地,位于半封闭式的颅腔内,从内到外分别被
软
脑膜、蛛网膜和硬脑膜包裹,并在膜内有脑脊液来缓冲外界震荡。人脑分为左右两
p>
个半球,之间由粗大的神经纤维束
—
胼胝体将其连接。左右两 个半球的皮层又分为
若干不同的区域,
分别负责不同的功能。
< p>大脑在活动时需要血液输送氧来提供能量,
因此,大脑皮层表面和内部分布有复杂的
血管网络。人脑的功能使人类在进化中拥
有了独一无二的高级智慧。
< /p>
人的大脑有千亿数量的神经细胞和胶质细胞,在脑功能活动时,神经元会产生
动作电位。神经元的生物电活动是人脑调控各种生命活动的基础。神经细胞的胞内
和胞外存在大量的带电离子,这种内外的离子浓度梯度形成了膜电位。在没有外界
刺
激的情况下,膜电位维持在一个恒定的水平,一旦受到外部刺激,膜电位就会发
生改变从
而使神经元产生动作电位。不同神经元之间通过突触进行连接,并且通过
突触来进行信息
的传导。
神经胶质细胞包裹在神经细胞表面,
支持、
保护神经细胞 ,
并维持钾离子的浓度。另外,部分神经胶质细胞的一端贴附在毛细血管壁上,为神
p>
经细胞供给氧和能量
[4]
。
当大脑某个区域的神经活动增加时,局部的血流、氧代谢、糖代谢会产生相应
的变化。因此,通过测量局部血流和代谢等神经元活动的间接反应也可以了解该部
位脑功能活动情况,相应为血液动力学参数的变化。大脑皮层的某个部位被激活需
要消耗能量来产生(或恢复)动作电位和细胞膜两侧的离子浓度梯度,因此该区域
的
新
陈
代
谢
会
增 p>
加
,
并
表
现
为
该
区
域
血
管
网
络
中
血
流
量
以
及
氧
合
血
红
蛋
白 p>
(
Oxy-hemoglobin
,
HbO< /p>
2
)浓度的增加
,
而去氧血红蛋白 (
Deoxy-hemoglobin
,
Hb
)< /p>
浓度小幅度降低
[1]
。
大多数研究认为造成这种血氧响应的原因是神经元直接或通过神
经胶质细胞间接地向血管
释放血管活性因子来调节血流
[5]
。
1.2.2
脑功能成像技术
脑功能成像技术可以动态地检测活体脑的生理活动,对当代认知神经科学的发
展产生深刻而巨大的影响。根据所测量的内容,可以把无创的脑功能成像技术分为
两
大类。第一类是非侵入性电生理技术,第二类是通过测量神经元活动引起的次级
2
华
中
科
技
大
学
博
士
学
位
论
文
反应(如局部葡
萄糖代谢和血流、血氧变化等)研究与行为相关的脑局部神经元的
活动情况
[6]
。
(
1 p>
)脑电图(
Electroencephalography
,
EEG
)
大脑皮层上的神经元无时无刻不在进
行着自发、节律的电活动。大量的神经元
活动会产生足够强的电场,在头颅外可以被检测
并记录。
EEG
通过置于头皮表面的
电极记录大脑皮层不
同部位的大群神经元的总体活动
[7]
。
神经电 活动又可以分为自发
脑电活动和皮层诱发电位。自发脑活动是在没有外界特定刺激的条件
下记录到的电
活动。根据
EEG
的频率和振幅,可分为< /p>
α
、
β
、
θ
、
γ
四种基本波形,由这些波形各
自组合成一连串的波称为节律。
比如,正常人在清醒状态和困倦状态下
α
波成分会
有很大
的差异。
在认知神经科学研究中更常用的是基于
EEG
检测的事件相关电位
(
Event-related
< p>Potential
,
ERP
)技术。
ERP
是一种与特定刺激或事件相关联的脑电活动。单次刺
激所产生的
诱发电位往往比较微弱,会湮没在自发脑电的背景下,因此需要通过多
次刺激将诱发电位
通过叠加平均提取出来。叠加平均是提高信号信噪比的常用方法,
比如将
N
次
ERP
信号叠加平均,其信噪比可以提高
倍。在实际应用中,为了得
到一个可靠
的
ERP
信号需要数十次到上百次刺激的叠加。
ERP
< p>波形由一系列大小和
正(
P
)负(
N
)不同的成分组成,这些成分通常根据正负和出现的时间来命名,比
如<
/p>
P300
表示出现在
300 ms
左右的正走向成分。
尽管
EEG< /p>
有很高的时间分辨率,能够直接反映大脑的神经活动,但是
EEG
信
号来源于大脑皮层神经元的集体放电,因此很难得到准确的放电位置。为了推算放
电位置一般的方法是增加电极数量,并根据头皮电极的电位进行溯源分析。尽管如
此,
EEG
技术的低空间分辨率依然是制约该技术应用的关键问题。< /p>
(
2
)脑磁图(
Magn etoencephalography
,
MEG
)
脑功能活动产生了颅内电流,而电流同时产生磁场。但是由于颅内电流密度很
低,所以形成的磁场非常微弱。
MEG
利用超灵敏的超导量子干 涉器记录脑磁信号的
变化,从而反应脑功能产生的电信号
[8]
。与
EEG
相比,
MEG
成像不受颅骨影响,更
3
p>
华
中
科
技
大
学
博
士
学
位
论
文 p>
加清晰,结果更加准确,同时时间分辨率也可以达到毫秒量级。
MEG
的仪
器非常复杂和精密,并且需要建
造一个严密的电磁场屏蔽室,因此造价高昂,只有
为数不多的单位安装有此设备。
(
3
)功能磁共振成像
(
functional Magnetic Resonance Imaging
,
fMRI
)
< br>核磁共振现象是由原子核的磁矩与自旋角动量这两者在外磁场的联合作用下形
成的
。功能核磁共振成像一般是依据血氧含量对比度的原理:神经元的活动会引起
一系列血液
动力学参数的变化,
比如脑血流、
脑血容、
血红蛋白浓度等
[8]
。
其中
HbO
< br>2
是反磁性物质,
Hb
是顺磁性 物质。当
Hb
和
HbO
2
浓度比降低时,磁化率引起的散
相会减少,从而表现为
MR p>
信号的增强。这种现象就是血氧水平依赖效应
(Blood
Oxygenation Level Dependent, BOLD)
fMRI
最大的优势是空间分辨率较高
(
2 ~4 mm
)
,
但是时间分辨率较低,为秒量级。功能磁
共振成像系统体积庞大,需要建立专门的
电磁屏蔽室来放置仪器。在实验中被试需要躺在
狭小的腔体内,扫描过程中不能发
生肢体运动
[9]
,因此不适合老人、小孩以及患有幽闭恐惧症的人群。
(
4
)正电子发射断层显像
(
Positron Emission Tomography
,
PET
)
PET
利用示踪剂测量局部糖代谢、
血流和血氧变化对各种脑活动进行 功能定位。
在成像之前首先注射带有放射性核素标记的代谢生物分子,比如葡萄糖。当带
有放
射性标记的代谢物聚集在相关区域时便开始进行成像
[8]
。
该区域的功能活动时由于代
谢需要就会聚集带
有标记的葡萄糖。具有很短半衰期的放射性核素开始衰变后,质
子转化为一个中子和一个
正电子,而正电子几乎在瞬间与组织内的电子结合发生湮
灭。根据质能守恒,湮灭后产生
一对飞行方向相反的
γ
光子。这一对光子被环绕在
生物体
外的探测器探测到,根据探测到光子的位置、数量、飞行时间,可以得到这
些代谢生物分
子在体内的分布情况并绘制分布图像。
PET
用于脑功 能检测时,由于人脑在活动时需要消耗能量,活动区域的葡萄糖
分子浓度增加,表现为该
区域图像信号增强。
PET
不仅可以测出大脑活动的区域,
并且还可以定量地测出激活的强度。
PET
利用示踪剂检测组织糖代谢和血液动 力学
参数变化的时间分辨率和空间分辨率都不如
fMRI
。
由于
PET
具有极高的化学选择性
和灵
敏度,
能够分辨出
10
-9
<
/p>
mol
浓度的化合物,
目前主要用于神经递质等物质的测量 。
4
华
中
科
技
大
学
博
士
学
位
论
文
p>
与
fMRI
类似,
PET
对于 环境和被试的限制也非常大。
(
5
)功能近红外光谱术
(
functional Near-infrared Spectroscopy
,
fNIRS
)
p>
1977
年,
Jobsis
发现了生物组织对 于近红外波段(
780-1100 nm
)的光具有低吸
收、高散射的光学窗特性
[10]
。人体组织中水的含量占了很
大比例,而水对近红外波
段的光吸收非常少。因此,近红外光照射人脑时,可以穿透头皮
和颅骨进入大脑皮
层,并被脑组织散射出来。散射出来的光相比于入射光的衰减程度和脑
组织中相关
吸收色团的浓度相关。具体关系可以通过修正的比尔朗博定律来表示。入射到
大脑
的近红外光在具有不同光学属性的多层组织中传播,光子的分布近似
“
香蕉形
”
。在
近红外光学窗波段,
p>
HbO
2
和
Hb
的吸 收光谱曲线有明显差异(图
1.1
)
,因此可以使
用该波段的近红外光对深层组织的功能活动进行检测。与
fMRI
类似,
fNIRS
通过检
测生物组织中
HbO< /p>
2
和
Hb
浓度的相对变化,从而间 接反映脑功能活动
[11, 12]
。
图
1.1
水,
HbO
2
和
Hb
< p>在近红外波段的吸收光谱曲线
[2,
13]
。
Fig. 1.1 The
absorption spectra of water, HbO and Hb in the
near-infrared spectral range
[2,
13]
.
相比于其他的脑功能检测技术,
fN IRS
同时兼顾了空间分辨率和时间分辨率
[3, 9]
,
能够检测多种血氧参数,
并且具有对环境和被试限制小, p>
制作和运行成本低等优势。
5
华
中
科
技
大
学
博
士
学
位
论
文
但
fNIRS p>
技术的局限在于近红外光的穿透深度有限,只能探测大脑灰质部分,而无
法检
测深层脑组织活动
[14]
。因此该技术主要应用于认知心理学
的研究中。
1.3
脑功能的近红外研究
1.3.1
近红外脑功能成像仪器
近红外脑功能成像技术是根据大
脑组织在光学窗波段低吸收高散射的特性来检
测血氧动力学参数,
从而反 映大脑的功能活动。
根据使用不同的近红外光检测方法,
用于脑功能成像
的仪器有三类:时域仪器、频域仪器和连续光仪器。三种仪器的工
作原理如图
所示。
图
1.2
三种近红外技术检测原理示意图
。
Fig. 1.2 Illustration of the three
near-infrared spectroscopy
techniques
[15]
.
[15]
美国宾夕法尼亚大学的
Britton
Chance
于
1988
年使用百皮秒量级的超短 光脉冲
照射大脑,获取出射光表征的包含了组织光学特性信息的时间点扩散函数,再通过
逆向求解时间点扩散函数获取被检测组织的约化散射系数
μ
s
’
和吸收系数
μ
a
。这就是
最早的时域仪器技术
[16,
17]
。在
1990
年,
Britton
Chance
又发明了对生物组织进行测
量的频域方法:使用射频调制的光源照射组织,测量出射光信号的强度和相位,再
将这
些信息带入到由玻尔兹曼辐射传输方程简化得到的标准扩散近似方程中
[18]
,能
够得到
μ
s
’
和
μ
a
[16]
p>
。由于时域仪器和频域仪器的复杂程度较高,大部分仍处于实验样
机
阶段,商业化程度低,因此较少应用于脑功能的研究中。
连续光仪器使
用直流稳态或者低频调制的光源照射大脑,通过探测大脑出射光
信号的光强变化来反映血
液动力学参数的相对变化,
从而间接地评估大脑功能活动。
连续光仪器采
集到的光强数据通过修正的比尔朗伯定律转化为血氧参数,而且只能
6
华
中
科
技
大
学
博
士
学
位
论
文
对血氧参数的相
对变化进行测量
[19]
。一般情况下脑功能研究关注的是大脑
活动也就
是变化过程,因此连续光技术适用于脑功能活动研究。目前,包括哈佛大学、华
中
科技大学和中科院自动化所在内的多家研究单位都在研究连续光近红外技术和相关
p>
应用
[11, 12, 20]
。日本的日
立和岛津以及美国的
NIRx
等公司也于
2003
年开始陆续推出
多款商用的多通道连续光近红外脑功能成像仪器。
连续光 近红外仪器复杂程度较低,
制作和使用成本较低,可以进行便携式设计,因此被广泛应用
于脑科学研究和临床
诊断中。
图
1.3
连续光
fN IRS
检测技术从
1992
年的单一通道仪器到多通道检测系统的 发展历程
。
Fig. 1.3 Sketch of
the development of continuous-wave fNIRS
instrumentation from single channel
instrument in 1992 up to the multi-
channel systems
[16]
.
[16]
早期的连续光
fNIRS
系统使用灯泡加滤波片作为光源,将光源和探测器直接与
头
皮接触,因此只能探测没有头发覆盖的额头区域,而且通道数量受到限制
[21]
。随
着科技的发展,连续光
fNIRS
仪器由单一通道发展为高密度多通道,从而具有全脑
成像的能力(图
1. 3
)
。为了对头发覆盖的脑区(比如运动皮层和视觉皮层)进行近
7
华
中
科
技
大
学
博
士
学
位
论
文
红外脑功能测量,光源和探测器需要通过光纤器件作为光信号入射和接收的设计方
案。
基于此方案近几年发展了多通道的光纤式
fNIRS
脑功能检测技术。光纤器件具
有较小的接触端面,可以透过头发的间隙与头皮良好接触。另外光纤系统可以将雪
崩光电二极管模块用于光电转换,进一步提高微弱光信号的检测灵敏度。目前的光
纤式
fNIRS
系统发展方向为多通道数目、高时间分辨率和高检测灵 敏度,从而可以
应用于光快神经信号以及脑网络连接的研究中。华中科技大学的骆清铭小
组研发了
一套多通道高灵敏度的光纤式
fNIRS
成像系 统,系统各项性能均达到世界领先水平
并可以提取光快神经信号
[22]
。本文的研究工作也是基于光纤式连续光近红外脑功能
检测技术开展的。
1.3.2
近红外脑功能研究的发展现状
从
20< /p>
世纪
70
年代开始,
近红外光谱术就被应用于生理检 测和监护中。
作为
fNIRS
的开端,在
20
世纪
90
年代初,日本、德国和美国的研究小组纷纷开展了使 用近红
外光谱术检测脑功能活动的实验,并发表相关研究
[23
-25]
。
1996
年骆清铭在
Brit ton
Chance
实验室成功研制出了具有
16 p>
个探测通道的近红外脑功能成像仪器
[26]
,
p>
并利用
该仪器对视觉和运动皮层的脑功能活动进行检测
[4, 26-28]
。
该仪器在对指运动范式下测
< br>量到运动皮层的响应结果和
fMRI
测量结果一致,
被 p>
Science
期刊以
The brain lights up p>
为标题进行了报道
[29]
。
在脑功能成像领域,
fNIRS
技术主要应用于认知 心理学的研究以及临床精神类疾
病的研究与诊断中
[12, <
/p>
30-32]
。对于认知心理学的研究主要关注前额叶大脑皮层。比如
p>
人们在进行记忆、情感、判断以及逻辑思维类的任务时,前额叶对应的大脑皮层会
出现血氧激活现象,即
HbO
2
浓
度的显著上升,任务结束后
HbO
2
浓度会逐渐
下降恢
复至基线状态
[33]
。由于前
额叶区域的颅骨较薄,探测难度较低,因此
fNIRS
在该领
域的应用非常广泛。而大脑在认知任务中的表现和许多神经和精神类疾病(比如精
神
分裂、阅读障碍、癫痫等)密切相关
[34,
35]
< br>。通过研究正常人和患者在认知心理学
任务中大脑活动的差异能够帮助研究人员深
入了解这类疾病的神经机制,从而指导
这类疾病的预防、诊断和治疗。随着光纤式
fNIRS
检测技术的发展,目前在临床上
8
华
中
科
技
大
学
博
士
学
位
论
文
fNIRS
p>
技术还可以应用于躯体运动皮层神经活动的研究中,
指导运动功能受损的患者
进行康复训练
[36]
。近几年,
fNIRS
在躯体运动与感觉皮层的脑功能检测开始应用于
脑机接口技
术的研究中,作为辅助定位的手段来提高脑机接口的识别准确率
[37,
38]
。
另外,随着
fNIRS
系统通道数目的增加和脑连接分析技术的发展,全脑的
fNIRS
检
< p>测技术开始应用于神经科学的脑网络研究中,帮助人们深入了解大脑的功能分化和
< br>功能整合
[39-41]
。
< p>
尽管
fNIRS
技术已经广泛应用于脑功能的研究中,目前仍存在制 约
fNIRS
技术
发展的三个主要问题:
(
1
)
由于近红外光需要穿透头皮和颅骨,
fNI RS
检测到的血氧
信号受表层血流影响严重,头皮、颅骨等浅表的信号会
对大脑皮层信号产生干扰,
降低脑功能信号的可靠程度
[42-
44]
;
(
2
)由于光纤式
fNIRS
成像技术发展较晚,系统
的长时间稳定性较差,
(
3
)
fNI RS
数据分析方法较为
单一。
目前的
fN IRS
研究一般只分析
HbO
2
在任务段的激活情况,
而激活情况只能反
映大脑的功能分化情况而不能反
映大脑的功能整合情况
[45]
。研究并解决以上问题是
提高
fNIRS
技术在脑功能研究中可靠性,进一步发掘
fNIRS
应用潜力的关键,也是
fNIRS-
EEG
双模态联合检测技术的基础。
1.4
近红外
-
脑电双模态联合检测的研究
< p>
fNIRS
技术是通过检测神经活动引发的氧代谢来反映脑功能活动,
对神经活动的
响应速度较慢,因此不能确定处于毫秒量级的神经活动的时间特征。大脑
是一个极
其复杂的动态系统,大脑的活动涉及到从分子到大脑皮层整体的多个层次的变化
,
因此需要更加多元化的脑功能检测技术。随着各种脑功能检测技术的不断发展,近
p>
几年发展出了将多种技术联合在一起检测脑功能的多模态成像技术。将
fNI RS
与其
他脑功能成像技术进行联合检测,可以弥补
fN IRS
技术的缺陷,实现优势互补,帮
助研究人员更加系统、深入地了解
认知任务中的神经加工过程和大脑功能机制。
目前已有一些国外的课题
组将
fNIRS
与
EEG
联合,
f NIRS
与
fMRI
联合以及
fNIRS
与
PET
联合进行认知科学和临床方面的研究,
并 取得了一系列有价值的成果
[46,
47]
< br>。其中
EEG
作为一种发展了多年的脑功能检测技术,商业化仪器和信号处 理方法
9
华
中
科
技
大
学
博
士
学
位
论
文
p>
都较为成熟,具有很高的时间分辨率,直接记录大脑皮层的神经电活动可以在毫秒
量级区分认知任务不同阶段的神经活动
[48]
。
fNIRS
和
EEG
都具有对被试限制小、探 p>
头佩戴方便的优势,因此应用对象可以是儿童或者老年人。它们还都具有良好的环
境兼容性可以用于床边检测(
fMRI
、
PET p>
、
CT
等只能在特定房间内使用)
,能够提
< p>供更加贴近生活场景的检测环境,因此适合应用于认知心理学和社会心理学的研究
< br>中
[49]
。
另外,
fN IRS
检测的内源光信号与
EEG
检测的电生理信号之间不会产生 串扰,
因此这两种技术非常适合联合在一起检测脑功能活动。将这两种技术联合在一起对
脑功能进行同步检测可以综合
fNIRS
和
EEG
的优势,
同时具备很高的时间分辨率以
及可以接
受的空间分辨率。另外,这两种技术采集到的脑功能信号可以互相补充、
互相印证,在神
经科学和临床研究中具有非常大的应用潜力
[50]
。
21
世纪初开始一些国外的研究小组使用
fNIRS-E EG
双模态检测技术来研究认知
心理学任务中的大脑活动,比如语言逻辑
、工作记忆、静息与觉醒状态等
[51-53]
。早
期研究目的主要是通过
EEG
信号来佐证
fNIRS
的结果有效性,
随着人们对
fNIRS
技
术的认可,研究人员开始通过分析这两种信号之间的相关性来解释大脑在认知任务
中的神经活动机制。这些研究一般使用商用的
fNIRS
和
EEG
系统,将二者的探头同
时摆放在被试头部实现联合检测。然而目
前光电联合的研究很少有
fNIRS
探头通道
和
EEG
电极同时覆盖全脑。大多是情况下是二者分开排步,即测试不同的区域,而
且
fNIRS
和
EEG
的探测通道 数都比较少(一般在
10
通道以内)
。这不仅影响了数据
有效性,而且局限了
fNIRS
和
EEG
数据分析方法。比如使用二者进行共定位研究或
者溯源分析时,就需要电极和探头
排布在相同的区域;而分析脑连接的时候需要较
多的通道做相关分析。另外,大部分双模
态检测研究中
fNIRS
采样率都在
10
Hz< /p>
左
右,
与
EEG
采样率差距 较大因此无法实现两种信号时域和频域的相关分析。
在国内,
双模态检测
技术起步较晚。
华中科技大学在
2009
年开始使用自主研发的< /p>
fNIRS
系统
和商用
EEG
系统开展双模态检测的研究
[11, 53]
。然而,这种双
模态检测技术也只是将
两个仪器的探测部分联合在一起,
在实际应用中存 在实验过程复杂、
信号质量较差、
数据分析方法单一等问题。主要原因是
自主研发的光纤式
fNIRS
系统仍处于实验机
阶段,系
统稳定性较低,信号质量较差,而且缺乏从物理层面到信号层面的双模态
10
华
中
科
技
大
学
博
士
学
位
论
文
联合手段。一直
以来,国内外各研究小组使用不同的
fNIRS
和
EEG
系统,不同的信
号同步方法以及不同的探头排布位置来进行双模态检测的研究,因
此缺乏规范的系
统联合技术以及双模态信号融合方法,导致双模态联合检测技术难以在认
知神经科
学研究中得到广泛应用。
最近几年
fNIRS-EEG
双模态成像技术开始转向脑 机接口的研究领域
[54-56]
。传统
的脑机接口只用到了
EEG
信号,
双模态的脑机接口技术通过加 入
fNIRS
信号作为参
考来辅助定位,可以提高脑机接
口的识别准确率。另外,双模态的脑机接口不仅可
以提高模式识别的准确率,还能够提供
更多的分类信号,进一步扩展脑机接口的应
用。比如
2016
< p>年韩国的研究小组成功使用fNIRS-EEG
混合脑机接口成功实现了输< /p>
出八种指令来控制飞行器
[57]
。
p>
然而该技术目前并没有实现真正的光电联合脑机接口,
只是将
fNIRS
和眼电信号提取出来进行分类,而且缺乏两种信号的反馈环节,降低
< br>了被试的训练效果从而影响识别准确率。
总体来说,目前的
fNIRS-EEG
双模态检测技术仍 处于起步阶段,尚无成熟的系
统联合技术以及配套的采集软件,缺少可靠的双模态物理和
信号层面联合手段,存
在采集通道数量少、实验实施困难以及系统非闭环等问题。在双模
态检测技术中
fNIRS
检测技术由于发展时间较晚,
还 额外存在系统稳定性差、
信号受表层干扰严重
等问题,需要开展研究工作
来补齐该技术短板。另外双模态检测技术应用于认知心
理学研究时缺乏合适的数据分析方
法,需要开展更多的工作进行方法建立和验证。
这些问题影响了双模态联合检测的可靠性
,限制了
fNIRS-EEG
双模态检测技术在认
知神经
科学和脑机接口领域中的应用。本文围绕以上影响双模态联合检测性能的因
素展开研究,
并将全脑双模态联合检测技术应用于大脑功能活动的研究中。
1.5
本文的主要内容
如前所述,
fNIRS-EEG
双模态成像技术对于脑功 能的研究有着非常重要的意义。
本文主要基于自主研发的
fNIRS p>
系统研究光电联合脑功能检测的方法,
在信号质量、
系统稳定
性、联合方式以及数据分析方法等方面对影响双模态检测性能的问题展开
研究。
11
华
中
科
技
大
学
博
士
学
位
论
文
p>
第一章,绪论。首先说明了脑功能研究的意义,并介绍了常用的几种脑功能检
测技术。基于这些检测技术引出本文的主要研究对象:
fNIRS-EEG
双模态联合成像
技术。介绍了双模态联合成像技术的意义和目前存在的一些问题,然后
引出本文的
主要研究内容。
第二章,提高全脑
fNIRS
检测性能的研究。首先研究了影响
fNIRS
信号的表层
干扰问题,使用小间距配合自适应滤波的方式进行去除。围绕影响系统
稳定性的问
题研究并解决了光源发热问题、重新设计光纤探头。在光纤式
fNIRS
系统基础上开
发了一套功能完善、使用方便的采集软件以及配
套的可视化人机交互界面。最后通
过一个
N-back
实 验验证
fNIRS
技术的在体检测有效性。
第三章,
提高全脑双模态联合检测性能的研究。
首先研制了光电一体式探头帽 ,
通过特制的布帽基底和探头夹具将多通道的
fNIRS
探头以及
EEG
电极固定在一起。
使用光电打标的方式实
现高精度信号同步。围绕双模态脑机接口的应用研究了基于
fNIRS
和
EEG
信号的在线反馈功能并将这两种信号的采集与显示融合在一个软件平
台中。最后通过一个
Stroop
实验验证双模态联合技术的在 体检测有效性。
第四章,工作记忆中脑功能性别差异的双模态成像研究
。本章使用全脑
fNIRS-EEG
双模态成像技术研究 了语言类工作记忆任务下脑功能的性别差异。然后
使用包括功能连接、
效 应连接和小世界网络特性等方法来分析数据。
最后结合
fNIRS
和
EEG
的结果对于性别差异的神经机制进行了讨论。
第五章,总结与展望。对本文工作进行归纳总结,介绍本文的研究内容、研究
成果和主要创新点,并对接下来的研究工作进行展望。
12
华
中
科
技
大
学
博
士
学
位
论
文
2
提高全脑
fNIRS
检测性能的研究
2.1
引言
fNIRS
技术具有良好的定位能力、
无损检测以及对被试和环境限制低等优势,
已
经被广泛应用于认知心理学的研究以及临床精神类疾病的诊断中。然而
该技术仍存
在无法反映精细的神经活动、信号受表层干扰严重等局限。将
fNIRS
与
EEG
联合进
行同步检测,
可以实现优势互补。
fNIRS-EEG
双模态检测技术通过不同结果的互相 p>
验证,在一定程度上避免由单一检测技术的干扰信号造成的误导性结果,从而提高
脑功能检测的可靠性。在双模态联合检测中由于
fNIRS
技术发展 较晚,相比于
EEG
系统在检测技术和信号的降噪处理等方面都不够成熟
,因此为了实现全脑双模态联
合检测首先需要对
fNIRS
系统中存在的影响检测信号有效性的问题进行研究。
本章的研究基于
课题组自主研发的多通道高探测灵敏度的光纤式
fNIRS
成像系
统
*
。该系统有
8
个光 源和
32
个探测器,最多可以组合出
48
个通道< /p>
[22]
。较多的探测
资源能够覆盖更多
脑区,满足多样化的脑功能检测需求,而且可以使用更多的信号
处理方法(比如脑网络连
接)来分析数据。本课题组曾使用该系统进行了一些认知
脑功能的研究,包括工作记忆、
儿童阅读障碍等
[11, 12]
,并在
201 3
年成功地检测到了
视觉皮层光快神经信号
[5
8]
。
然而目前的光纤式系统仍存在一些问题
影响着
fNIRS
探测脑功能活动的能力。
首先是
fNIRS
的信号质量。
在
fNIRS-EEG
双模态技术的应用中为了对
fNIRS
和
EEG
< p>信号做相关分析,需要两种信号都有比较高的信噪比。
EEG
系统中 的信号调理和噪
声抑制技术较为成熟,比如去除眼电、双极参考、工频陷波、小波变换滤
波等
[59]
,
已经具有了较高的信噪
比。而
fNIRS
系统采集到的信号除了包含脑功能相关的目标
< br>信号外,
还包含仪器噪声、
生理噪声,
生理噪声包括呼吸、
心跳、
以及表层干扰
[15, 60]
< br>等。近红外光是穿透头皮和颅骨到达大脑皮层,因此
fNIRS
探测的血氧 信号中包含
头皮和颅骨中的血流信号。这种干扰信号来源于浅表组织,基于组织反射式数
据采
*
该成像系统是在导师指导下,由张垒、张中兴、高晨阳等人共同完成
13
华
中
科
技
大
学
博
士
学
位
论
文
p>
集模式的
fNIRS
对于该信号更加敏感,因此由认知任务引 起的血氧变化信号有可能
被湮没在表层干扰中,从而严重影响
fNIRS
检测脑功能的准确性
[42-44]
。仪器噪声
以
及呼吸、心跳等生理噪声一般可以通过锁相放大以及滤波的方式加以去除
[22,
61]
,而
表层干扰作为
一种自发的血氧信号,混杂在大脑皮层神经活动引发的血液动力学参
数中,去除难度较大
。
其次在系统的长时间稳定性方面:光纤式
fNIRS
系统使用的
LD
光源发热、探
头固定不牢
固等,都会在实验中引入干扰噪声,影响目标信号的提取
[62]
。最后在采
集软件方面:已有的系统软件功能单一,可视化程度低,显示的信号对于研
究人员
而言不直观,从而影响研究人员的判断导致产生无效的实验数据。这些问题不仅影
响
fNIRS
检测脑功能活动的性能,同时影响双模态联 合检测技术的可靠性。因此本
章针对以上问题对于提高
fNIRS
检测技术的有效性开展研究,通过小间距通道配合
自适应滤波去除表层干扰;研
究并解决影响系统稳定性的若干问题;开发基于
LabVIEW
和
Matlab
混合编程的可视化系统软件;最后通过在体
N-back< /p>
实验来对
fNIRS
检测前额叶脑功能活动的有效性进行验 证。
2.2
使用小间距通道去除表层干扰
本章的研究是基于课题组
自主研发的连续光
fNIRS
系统,
整体结构如图
2.1
所示。
该系统的光源和探测器探头部分使用光纤设计,可以使探头
末端避开头发直接接触
到头皮。系统使用的探测器为雪崩光电二极管(
A valanche photodiode
,
APD
)
< p>,其探
测灵敏度为
0.5
A/W
(
800nm
)
,电流
/
电压放大 器增益为
1.5×
10
6
p>
V/A
,能够检测到
非常微弱的光信号。光源使用激光二极管
(
Laser diode
,
LD
)
,并通过捆绑合束将两
个波长耦合在一起。
相比于
LE D
,
LD
的功耗更小,
发出光信号的频谱半高宽更 窄
(约
为
1 nm
)
,使 血液动力学参数的测量更加准确。另外
LD
出射光的准直性非常好,不
< p>仅易于与光纤进行高效的耦合,而且在保证人体安全光功率的前提下可以提高在体
< br>测量的有效信号强度。
LD
光源采用不同频率的正弦波调制,配合上位机的 数字锁定
技术可以在很大程度上抑制信号传输中的噪声。采集卡使用
NI
公司的
PXIe-6368
,
一块采集卡
有
16
个模拟输入的通道,并行采样率可达
2 M
,即采样率不会随着通道
14
华
中
科
技
大
学
博
士
学
位
论
文
数目增加而下降。
整个光纤系统经过大点数锁定后的时间分辨率最高可以达到
5 ms
。
图
2.1
近红外脑功能仪器的整体结构示意图。
Fig.
2.1 Block diagram of the fNIRS instrument.
目前已有一些
fNIRS
课题组围绕
fNIRS
表层干扰的去除方法开展研究。这些方
法通常可以分为纯软件的方法和引入参考信
号的方法。纯软件的方法是通过离线的
数据处理方法,包括低通滤波、主成分分析和小波
滤波等去除包含表层干扰在内的
生理噪声,并将任务相关的血氧信号提取出来
[60,
63-65]
。其中常用的滤波方法是根据
p>
干扰信号的频谱特性将噪声滤除。然而由于表层干扰的成分复杂,其频率范围和大
脑皮层中血氧信号(目标信号)的频率范围很大一部分是混叠在一起的
[6
6]
,因此只
通过软件滤波的方法很难在保留目标信号的前提下
去除表层干扰。
考虑到干扰信号与目标信号的来源不同,引入参考信号
的方法更加有效。该方
法首先需要加入额外的探测通道采集干扰信号,然后将干扰信号作
为参考信号通过
软件算法进行滤除。根据表层干扰信号来源于浅表的头皮和颅骨,可以仅
提取表层
组织中不包含大脑功能活动的血氧信号,然后将该信号作为参考来去除
fNIRS
信号
中的表层干扰成分
[67
]
。根据课题组之前使用蒙特卡洛模拟光子在大脑皮层中的传输
分布,
fNIRS
探测深度和光源与探测器的距离有关,当间距只有
1 cm
时该通道的探
测深度为
~0.5 cm
,因此只能接收到来自头皮和颅骨的血氧信号
[68]
。这种比正常间距
15
p>
华
中
科
技
大
学
博
士
学
位
论
文 p>
短的通道称为小间距探测通道。小间距探测通道和正常间距探测通道的探测
原理如
图
2.2
所示。
其中黄色线条代表 头皮,
红色香蕉形曲线代表不同间距通道探测的光子
分布。通过放置小间
距和正常间距探测器可以分别记录到表层血氧信号和包含表层
血氧信号和脑功能活动的混
合信号。由于小间距通道和正常间距通道一样可以探测
人体组织的血氧浓度变化,只是无
法探测大脑活动的信号。因此小间距通道不仅可
以探测到表层的干扰,还可以接收到来自
系统或者人体的全局噪声,比如光源的温
度漂移、脉搏和呼吸波动等。这类噪声同样可以
通过将小间距通道作为参考的方式
进行滤除。
图
2.2
小间距通道去除表层干扰原理图。
Fig. 2.2
Schematic diagram of short separation channels
removing superficial interface.
2.2.1
构建小间距探测通道
p>
相比于单独使用滤波等信号处理方法,使用小间距通道引入参考信号的方法能
够在不影响目标信号的前提下对表层干扰进行去除或抑制。
本节通过设计新的光纤、 p>
研究信号处理算法、改进采集软件等工作将该方法引入到光纤式
fNIRS< /p>
系统中,并
对其去除表层干扰的效果进行测试。小间距去除表层干扰的方法
首先需要一个或多
个短间距通道将该区域的表层干扰提取出来,然后通过构建数学模型,
将表层信号
作为参考来去除正常通道信号中的表层干扰。
为了小间距去除表层干扰方法应用于光纤式
fNIRS
系统中,首先需 要在系统中
构建小间距探测通道。当光源和探测器间距只有
1 cm p>
时光子在组织中的平均光程长
度将大大缩短,
因此组织的散射 和吸收对入射光造成衰减相比于正常通道大大减小。
由于光纤式
fNIR S
仪器使用了高灵敏度的
APD
进行光电转换,
如 果使用现有的光源
和探测器来构成小间距通道,
探测器检测到的光强经过
APD
转换后有可能达到饱和。
在组织光学仿体上测试<
/p>
3
cm
和
1
cm
间距探测器接收到的光强,该仿体是使用凝胶
16
华
中
科
技
大
学
博
士
学
位
论
文
作为基底材料制
作而成的用来替代生物组织的实验用品,并在制作时掺杂一定比例
的墨水和白色颗粒来模
拟生物组织对光学的吸收和散射
[69]
。
结果 表明
1 cm
间距下探
测器接收到的光强约为
3 cm
< p>间距探测器接受光强的120
倍。为了让小间距探测器工
作在正常的动态范围内,
有两种方案可供选择:
(
1
)
为小间距通道加入额外的光源,
该光源的出射光功率低于正常光
源。
(
2
)设计制作新的探测光纤,新探测光纤的数
值孔径或者探测面积小于正常探测光纤。这两种方案的最终实现效果是相同,但由
于光纤式
fNIRS
系统使用正弦调制的光源,第一种方案会增加系统复杂程 度,而且
临近光源可能会引起周围通道的串扰,因此本节采取了第二种方案来构建小间距
探
测通道。
根据本课题组之前使用蒙特卡洛模
拟光在组织中传输的研究,探测器接受的有
效信号强度和光纤的探测面积成正比,和数值
孔径的平方成正比。因此可以通过减
少小间距探测光纤的探测面积来降低它接收到的有效
信号强度。应用于小间距的探
测光纤探头部分如图
2.3
所示,
光纤材料选用玻璃束,
数值孔径与正常光纤相同,
为
0.55
,芯径设置为
0.2
mm
< p>。根据探测面积比例,在理论上小间距光纤接收光强幅值
为正常探测光纤(芯径
p>
2.5 mm
)接收光强的
156
分之一,可以满足探 测光强在
APD
动态范围内的要求。
SMA 905
90
°
18mm
Φ
1mm
图
2.3
小间距探头示意图。
Fig. 2.3 Sketch map of short separation
detection probe.
小间距光纤末端设计为
SMA90 5
接头,与正常间距探测光纤相同,因此可以根
据实际需求灵活地替换机
箱上的正常探测光纤。制作好的小间距光纤首先要测试它
是否能够工作在正常的动态范围
,测试方法是将光源探头和正常探测探头固定在组
织光学仿体上,间距设置为
。将光源调到正常工作状态,然后在距离光源
1 cm
17
华
中
科
技
大
学
博
士
学
位
论
文
p>
处放置小间距光纤,并观察它和正常通道检测到的光信号。测试结果如图
2. 4
所示,
在相同光源下,两种不同间距通道探测到的信号幅值基本一致,
说明制作的小间距
探测光纤可以满足在
1 cm
源探间距 下
APD
输出电压在合适的动态范围内。
图
2.4
以组织光纤仿体为探测对象的
(a)
小间距通道和
(b)
正常通道的信号波形图。
Fig. 2.4 Waveforms detected from
phantom of short (a) separation channel and (b)
normal channel.
为了能在近红外脑功能实验中应用去除表层干扰的方
法,需要将小间距探头固
定在以低温热塑板为基底探头帽中。由于该探头直径只有
1 mm
,之前使用的探头夹
具无法将它固定。解决方案是在距
离光源位点
1 cm
处打一个直径
1 mm
的小孔 ,将
小间距探头从孔中插入。为了让小间距探头固定得更加稳固,在小孔上下分别粘贴<
/p>
一小块
2 mm
厚的黑胶棉。
由于小间距光 纤探头非常细,
能够直接接触到头皮,
因此
不需要拨开头
发,但需要在小间距探头末端包裹乳胶管以降低其对头皮的刺痛感。
方法确立之后测试小
间距测量血氧信号的效果,将探头佩戴在头上进行测试,可以
清晰地看到心跳。由于表层
干扰信号具有局部一致性,即在一块较小的头皮区域内
的表层血氧参数波动情况是一致的
,因此一个小间距通道可以为附近的所有正常通
道提供参考。为了节省探测器资源,本节
设计的探头帽在每个光源附近安装一个小
间距探头,该通道的信号用于该光源四周通道的
表层干扰去除。在
fNIRS
采集软件
中,
小间距通道和正常间距通道的数据使用相同的方法采集并保存在同一个文件中。
18
华
中
科
技
大
学
博
士
学
位
论
文
2.2.2
表层干扰的去除方法研究和模型实验
小间距通道设计完
成后需要研究使用参考信号去除表层干扰的算法。常用的方
法有静态估计、卡尔曼估计和
自适应滤波。静态估计的方法是从线性最小均方估计
方法中发展而来。
在 这种模型中,
把表层干扰作为残差项从正常通道的信号中去除
[70]<
/p>
。
也就是认为短间距信号成正比地包含在正常间距信号中,通过线
性回归的方法将这
种比例关系计算出来。尽管静态估计的方法容易实现,然而大脑是个复
杂的动态系
统,实验中短间距所包含的表层信号比例并非一直不变的,而是会随着时间变
化。
为了适应这种变化的趋势,在静态估计的基础上发展出了卡尔曼滤波估计方法。卡<
/p>
尔曼滤波是利用线性系统状态方程,通过系统输入输出观测数据,对系统状态进行
最优估计的算法。该方法使用静态估计的结果作为首次估计值。然后对信号运行两
次卡尔曼滤波得出最终估计值
[67,
71]
。该方法可以根据参考信号的变化来调整滤波系
数。
作为一种智能的信号处理手段,自适应滤波同样适用于研究不确定系统或信息
过程。将原始信号和误差信号输入自适应滤波器,滤波参数受误差信号控制,输出
信
号在原始信号基础上将误差信号去除。在表层干扰去除的应用中一般使用线性自
适应滤波
器,首先将光强信号转换为血氧信号,对于短间距和正常间距通道的信号
分别根据它们的
标准差做归一化处理,然后将短间距信号作为参考信号来对正常间
距信号进行自适应滤波
的处理,以去除正常间距信号中的短间距信号成分
[72]
。自
适
应滤波方法首先将正常间距通道和小间距通道的信号根据它们的标准差进行归一化,<
/p>
然后进行逐点计算。根据公式
2.1
来计算去除表层干扰后 的信号
e
。
[n]
,
[n-k]
2.1
其中
代表正常间距通道的数据,
代表小间距通道的数据。
N
< br>a
是滤波器的阶
数,即参考之前多少个小间距数据点。<
/p>
,
是自适应滤波的相关性系数,它的值根
据公式
2.2
< p>来计算。
,
=
,
2u
e[n-1]
[n-1]
2.2
19
华
中
科
技
大
学
博
士
学
位
论
文
初始值
,
设置为
[1
0
0
…]
T
,数组大小与滤波器的阶数有关。
u
为增益系数,
一般可
以设置为
0.0001
。对于
n
阶自适应滤波,每 次计算需要用到前
n
个时间点的
小间距数据和自适应滤波
后的数据。根据公式
2.1
和
2.2
,最后输出的
e
值即去除表
层干扰后的结果。
图
2.4
加入静息态噪声后正常通道的仿真信号。
Fig.
2.4 Simulation signal of normal channels with
noise of resting state.
通过构建
了模拟的小间距通道信号和正常通道信号来测试以上方法对于表层干
扰的去除效果。小间
距通道信号中只包含表层干扰,正常通道信号中包含任务相关
的血氧响应和表层干扰。为
了尽可能地模拟在体实验的情况,任务相关的血氧响应
使用之前实验中叠加平均出来的任
务激活波形,
在测试中只模拟了
HbO
2
的浓度变化,
在任务开始时
HbO
2<
/p>
的浓度上升,
任务结束后浓度下降。
使用在被试静息态下小 间距
通道采集的一段数据作为表层干扰信号,静息态下的信号只包含呼吸、心跳等各种<
/p>
生理噪声,没有任何脑功能活动引发的血氧变化。模拟正常通道中的信号时,将典
型的血氧激活波形(任务开始后上升,任务结束后下降)复制
15
段,加入到静息态
数据中,来模拟一个有
15
个
Block
的实验中的血氧信号。正常通道的模拟信号如图
2.4
所示,任务响应基本上被湮没在噪声中。
由于卡尔曼估
计的方法是在静态估计的基础上发展而来,因此本次测试只使用
卡尔曼估计和自适应滤波
两种方法,将模拟的表层干扰作为参考信号来处理正常通
20
华
中
科
技
大
学
博
士
学
位
论
文
道的仿真信号。为了能够清楚地比
较这两种方法的效果,在使用这两种处理方法之
前都没有用到其他处理方法,比如低通或
带通滤波,而是直接对每个
Block
的信号
进行叠加平
均。
处理后的信号结果对照如图
2.5
所示,
原始 信号作为参照即模拟任务
中的标准血氧响应信号,不包含任何干扰噪声。相比于没有任何
处理的结果,使用
小间距做滤波的信号更为平滑,噪声成分更小。卡尔曼估计的结果原始
信号强度有
一定的削减,相比较而言自适应滤波的结果更加接近原始信号。而且卡尔曼估
计的
程序较为复杂,需要采集完整的数据段进行离线处理,并且只能处理
Block
设计的
实验。相比于卡尔曼估计的方法,自适应滤波不仅可以
应用到
Block
设计的认知实
验任务,还可以应用到比
如静息状态、觉醒状态等非
Block
设计任务中的脑功能研
究中。而且该方法只需要参考前
n
个时间点的数据,因此可以对数据进行逐点 地处
理,具有将小间距去除表层干扰的方法应用到在线处理的潜力。因此,本节选择了<
/p>
自适应滤波作为小间距通道去除表层干扰的处理算法。
图
2.5
不同方法处理后的结果对比图。
Fig. 2.5
Results of different processing methods.
自适应滤波方法在理论上可以去除与参考信号来源相同的信号成分,不仅可以
滤除信
号中的其他全局干扰,比如心跳、呼吸等生理噪声,而且可以抑制由于光源
本身的波动而
引入的噪声。在光学仿体上测试它对光源稳定性的提升效果。使用一
21
华
中
科
技
大
学
博
士
学
位
论
文
个双波长光源(
785
nm
和
850
nm
)
< p>,一个正常间距探测器和一个小间距探测器同时
检测,
测试时间为< /p>
40
分钟。
测试完成后使用自适应滤波离线处理数据。
结果如图
2.6
所示,滤波后的信号(绿色曲线)稳定性有明显提升。
一般使用变异系数作为表征
信号稳定性的指标。变异系数通过标准差除以均值计算得到,
因此可以比较测量尺
度不同的两组数据的波动水平。分别计算两个光源测试信号滤波前后
的变异系数,
具体计算方法为整段数据的标准差除以均值。滤波前两个波长的变异系数为
2.09%
和
1.77%
,滤波后为 p>
0.49%
和
0.35%
。结果表明小间距作为参考信 号可以有效抑制全
局噪声,降低光源波动造成的干扰。
图
2.6
双波长光源稳定性自适应滤波处理结果图。
Fig.
2.6 Double wavelengths light source stability
figures of adaptive filter processing.
2.2.3
小间距去除表层干扰在体测试
根据上节的研究,加入小
间距通道并配合自适应滤波的方法在仿真信号和仿体
测试中有较好的抑制干扰的效果。<
/p>
然而模拟数据中的干扰情况较为简单:
每个
Block
中任务相关的血氧激活都是一致的,小间距通道和正常通道中的模拟干扰信号也是
相同的。
而且在实际实验测量中的情况中,
脑组织结构比光学仿体的复杂度高 很多,
干扰信号成分也更复杂和不可预测,为了进一步验证
fNIRS< /p>
光纤系统的小间距通道
在实际应用中去除表层干扰并还原目标信号的效果,
本节通过一个典型的认知实验
进行在体测试。
实验采用经典的工作记忆范式
2-back
。
实验任务为
Block
设计,
每个
Block
包含 p>
22
华
中
科
技
大
学
博
士
学
位
论
文
20
个字母刺激,任务期持续
40 s
,休息期
30 s
。整个实验包含
6
< p>个Block
。被试需要
在呈现的字母序列中判断当
前出现的字母和之前隔一个出现的字母是否一致,如一
致按键盘上的左箭头,不一致则按
右箭头。作为典型的工作记忆任务,负责信息编
码和存储的前额叶区域在任务中会有显著
的血氧激活现象
[12, 73]
。
为了方便实 验进行,
本次实验使用较为
简易的探头设计。探头帽包 含
2
个光源,
2
个与之对应的小间距
通道和
8
个正常间距通道。
探头排布如图
2.7
所示。
红色圆形表示光源,
蓝色方框代
< br>表正常探头,绿色三角代表小间距探头。摆放时根据脑电帽的电极位点进行定位,
两个光源分别在
FP1
和
FP2
(黑色圆形)的正 下方。实验开始前首先给被试讲明实
验要求和注意事项,然后佩戴探头帽,观察
fNIRS
信号保证每个通道都可以看到清
晰的心跳波动。由于小
间距探头探测面积非常小,本身体积也非常小,在实验中可
能产生晃动而引入额外的噪声
,因此在实验准备中使用柔性胶带来辅助固定。
FP1
FP2
图
2.7
小间距测试实验探头排布图 。圆形代表光源,正方形代表正常探测器,三角形代表小间
距探测器。
E EG
电极
FP1
和
FP2
用于辅助 定位。
Fig. 2.7 Configuration of
fNIRS channels for experiment. The circulars
represents the light sources, the
squares represents the normal detectors
and the triangles represents the short sepration
detectors. EEG
electodes FP1 and FP2
were used for probes positioning.
23
华
中
科
技
大
学
博
士
学
位
论
文
图
2.8
通道
7
和通道
8
整段
HbO
2
信号自适应滤波前后对比图
Fig. 2.8 Hemoglobin
signal time series for before and after adaptive
filter.
将采集完的
10
个通道
(
8
个正常间距通道,
2
个小间距通道)
的数据导入
Matlab
软件进行离线处理,对每通道的数据进行
滤波、血氧转换。然后再使用小间距通道
的信号作为参考,用自适应滤波算法去除每个正
常通道血氧信号中的表层干扰。选
取了两个典型通道滤波前后的整段
Hb O
2
信号作对比
(图
2.8 p>
)
,
图中可以看出自适应
滤波处理后的信号(
Filtered Ch7 and Filtered Ch8
)更加平滑,基 线更加平稳,干扰
信号受到了有效的抑制。
图
2.9
滤波前后血氧参数叠加平均时间序列对比图。
子图中左右两个黑色竖线分别代表任务
开始和结束时间。
Fig. 2.9 Averaged hemoglobin signal
time series for before and after adaptive filter.
Left and right
24
华
中
科
技
大
学
博
士
学
位
论
文
black vertical
lines in each subgraph represent begin and end of
task.
之后,分别对去除表层干扰前后的血氧信号进行叠加平均观察任务中的血氧
响
应(图
2.9
)
。对比后可以看出滤波 后的信号(实线)噪声波动和伪迹更小,任务中
的血氧响应也更加显著。
正常情况下
HbO2
响应曲线在任务结束阶段应该回到基线附
近,而滤波后的曲线在任务结束后更接近基线,因此表明滤波是有效果的。滤波前
信
号在任务结束后没有回到基线可能是因为血氧的全局波动导致的基线漂移。在体
实验的结
果证明本节的小间距通道自适应滤波方法具有较好的噪声抑制效果,能够
还原出任务中的
血氧响应,可以应用到
fNIRS
脑功能成像的研究中。
2.2.4
在线自适应滤波去除表层干扰
目前已有的表层干扰去除
方法均是在实验结束后将
fNIRS
数据全部导出离线处
理完成。然而在很多应用中,需要在实验过程中就可以看到去除表层干扰的血氧信
号。去
除表层干扰后的血氧信号对于实验的指导意义更加可靠。尤其在基于
fNIRS
< br>的神经反馈应用中,反馈信号长度相对较短,一般在
2~4 s
,在较短的 时间内表层干
扰和生理噪声比如心跳对于反馈信号的影响将变得非常大,而且无法通过叠
加平均
的方法来提高信噪比,从而可能产生虚假的反馈信号干扰被试的判断。因此,在神
经反馈应用中实时表层干扰去除更加必要。反馈信号需要实时地呈现给被试,没有
办法进行比如伪迹剔除、叠加平均等离线处理。使用小间距通道配合自适应滤波的
方法具有在线处理的能力,因此本节将该方法融合到光纤
fNIRS
系 统的上位机采集
软件中,在人机交互界面实时显示去除表层干扰后的信号。
y
1
(1)
,< /p>
y
3
(1)
Initia
lization
y
1
(1),w(1),e(
1)
y
1
(n-1),w(n-1),
e(n-1)
Adaptive filter
y
1
(n),y
3
(n)
p>
y
1
(n),w(n),e(n)
图
2.10
实时自适应滤波程序框图。
Fig. 2.10
Block diagram of real-time adaptive filter.
25
华
中
科
技
大
学
博
士
学
位
论
文
p>
实时的表层干扰去除第一步是要在程序中嵌入自适应滤波器。出于简化处理过
程,提高程序运行速度的目的,使用一阶的自适应滤波器。根据公式
2.1
和
2.2
,只
需要当前时间点采集的数据和上一个时间
点采集的数据即可进行得出去除表层干扰
后的数据。由于自适应滤波在执行中需要设置初
始值,因此在程序中加入一个采集
循环运行次数的计数器,将计数器的值作为后续自适应
滤波的执行条件。实时自适
应滤波的程序框图如图
2.10
所示。在第一次采集循环对数据进行初始化。之后的每
一次循环调用上个循环的小间距
通道数据
,滤波系数
w p>
以及去除表层干扰后的数
据
e
值,自适应滤波 处理完后更新这个三个值并保存到一个数组中供下一个时间点
的自适应滤波程序来调用。
自适应滤波运算涉及较多矩阵乘法运算,如果使用
LabVIEW
自带的 乘法和
for
循环来实现会拖慢程序的运行速度,从而降低表层干扰
去除的实时性,因此在线的自适应滤波部分使用
Matlab
脚本实现 。每个时间点的滤
波处理之后,将数组
e
按照通道次序分 别送到每个通道的波形图中显示。
使用图
2.7
排布的探头帽进行测试。
测试时让被试保持静息状态。
在静息态中主 p>
要的可见生理噪声为脉搏波
[74]
,因此
该测试主要观察实时自适应滤波对于心跳波动
的去除效果。测试实验开始前先将探头帽固
定在被试的额头部位,然后依次放入长
间距光纤、小间距光纤和光源。静息态下可以观察
到通道中显著的脉搏波动,因此
在测试中主要观察自适应滤波对心跳的去除效果。首次实
验中的实时效果如图
2.11
所示。上面一排波形图显示四个正常通道信
号,最右边为小间距通道信号,下面一
排为去除后的信号。每个波形图中的红色曲线代表
850 nm
波长光信号,蓝色曲线代
表
785 nm
< p>波长光信号。可以看出正常通道有明显的心跳波动,而自适应滤波后的信
号中心跳波
动收到抑制。
26
华
中
科
技
大
学
博
士
学
位
论
文
图
2.11
使用在线自适应滤波的表层干扰去除结果。 粗线为
HbO
2
信号,细线为
H b
信号。
Fig. 2.11
Superficial removing results using online adaptive
filter. The thick line for
HbO
2
and the thin
line for Hb.
然而本次测试中只有在通道
1
和通道
3
表现出心跳波动滤除的效果,通道
< p>2和
通道
4
滤波后的心跳信号仍比较 显著,而且滤波后信号的直流量也有很大程度的降
低。
观察信号的幅值可 以看出通道
1
、
3
和小间距通道的数量级更加接近 ,
而通道
2
,
4
幅值与小 间距通道幅值差距较大,
因此影响了自适应滤波的去除效果。
为了让信号
更加真实地反映血氧变化,在不改变信号幅值的前提下滤除心跳波动,在进行自适
应滤波之前需要对参考信号进行矫正。首先将小间距信号的心跳波动水平和正常间
距信号统一起来,具体实现方法为:将小间距通道信号对应每个正常间距通道信号
乘一
个增益系数作为该正常通道信号的参考信号。该增益系数是通过被参考信号的
标准差除以
参考信号的标准差来得到。为了不影响正常间距通道的幅值,小间距通
道信号在自适应滤
波计算中需要矫正到基线值,即去除该信号的直流成分。由于该
软件为实时采集,个体的
血流、血氧等生理参数会随着时间变化,与之对应的增益
系数和基线也会随着时间改变,
因此在程序中需要对这两个参数进行调整。程序中
的自适应滤波实时执行,只能根据之前
时间段的数据对这两个参数进行估算,用以
矫正接下来的数据。
为了增加 软件的运行效率,
程序设定为开始采集后先取
6 s
信号
计算增益系数和基线,之后每隔
12 s
取当前
12 s
的信号计算增益系数和基线,然后
把计算出的结果实时更新到一个数组中,随后的自适应滤波调用该数组,通过公式
2.3
将小间距通道信号通过公式
3.3
矫正后作为参考信号 去除表层干扰。
[n]
[n]-
)
2.3
其
中
[n]
表
示
小
间
距
原
始
信
号
,
[n]
表
示
< p>矫正
后
用
作
参
考
的
信
号
。
和
分别代表当前时间段的 基线值和增益系数,
这两个参数在程序运行时根据当
前数据段进行更新。
27
p>
华
中
科
技
大
学
博
士
学
位
论
文 p>
图
2.12
使 用矫正参考信号的在线自适应滤波的表层干扰去除结果。粗线为
HbO
2
信号,细线为
Hb
信号。
Fig. 2.12 Superficial removing results
using online adaptive filter with modified
reference signal. The
thick line for
HbO
2
and the thin line for
Hb.
使用加入参考信号矫正的在线滤波程序程序测试静息态被试的血氧变化,实验<
/p>
中的实时效果如图
2.12
所示,上面四个通道为原始信号 ,下面四个通道为滤波后的
信号。可以看出,滤波后每个通道的心跳波动均有明显抑制,
而且滤波后信号的幅
值相比于原始信号没有太大改变。
2.3
提高
fNIRS
系统稳定性关键问题研究
作为生理信号检测系统,为了提取可靠的目标生理信号,需要
fNIRS
系统能够
在长时间工作在稳定的状态
[75]
。而系统在实际工作中由于外部信号干扰、系统发热
以及探测位点接触不
良等原因造成检测信号的无规则波动,这种波动会叠加到任务
相关的血氧响应信号中从而
影响反映脑功能活动的目标信号提取。为了得到有效的
检测信号,首先对这些影响系统稳
定性的问题进行研究。其中外部信号干扰的问题
在课题组研制的光纤式
f NIRS
系统中通过加装电磁屏蔽盒以及设计稳压电路的方式
进行了解决
。因此本节主要围绕系统发热和探测位点接触不良的问题进行研究与解
决。
2.3.1
光源温度监测与控制
fNIRS
仪器光源所使用的
LD
器件在 正常工作时都会产生热量,
热量无法及时散
去就会造成
L D
表面温度上升,不仅影响
LD
光功率的输出并造成探测信号的波 动,
28
华
中
科
技
大
学
博
士
学
位
论
文
p>
而且
LD
长时间工作在高温下其使用寿命会缩短。对于整个光 源机箱而言,
LD
是主
要的发热器件,在工作中还会导致
整个光源机箱温度上升,从而影响机箱内其他电
子器件的性能,并加速器件的老化。目前
的商业
fNIRS
系统均没有解决光源发热以
及控制系统
温度的相关措施。
为了研究并解决该问题,
本节首先在 系统中加入了温度传感器
LM35D
对光源机
箱的温度进
行实时监测。该传感器根据环境温度输出线性的电压信号,探测灵敏度
为
10 mV/
℃,误差为
±
0.5
℃。把温度传 感器固定在光源机箱内,通过导线与采集卡
相连,在实验中将采集到的电压信号由采集卡
发送到上位机,然后转换为温度值在
人机交互界面中实时显示。在实验中如果机箱温度过
高或者有异常发热现象可以提
醒实验人员让系统暂时休息一段时间或对异常情况进行排查
。采集软件可以将实验
中的温度信号和生理信号保存在同一个文件中,可供离线做误差分
析,或者将温度
变化用作信号校准的参考。经测试,在系统正常工作时机箱内的温度逐渐
上升,可
以达到
43
℃。在同一时间点下,
< p>LD表面的温度比机箱温度高
3~5
℃。
为了研究
LD
的工作温度对于
LD
输出光功率的影响,在光学仿体上测试
LD
输
出
光功率与温度之间的关系。
首先将温度传感器与测试所用的
LD
固 定在一起,
二者
之间通过导热硅胶粘合在一起。然后开机将光源功率调至
正常工作范围,然后同时
记录长时间的光功率和
LD
表明 温度,测试时间为
8000
s
,测试时的环境温度为
27.5
℃。
LD
温度和输出光强的测试结果曲线如图
2.13
所示。
29
p>
华
中
科
技
大
学
博
士
学
位
论
文 p>
图
2.13 LD
输出功率和表面温度的长时间曲线。
Fig.
2.13 Long time curves of LD power and surface
temperature.
测试结果表明随着
LD
表面 温度的升高,
LD
输出功率不断下降,变化速度会随
着时
间减慢。二者之间有很强的负相关性,相关系数为
-0.989
(
p<0.0001
)
,说明
LD
温度的上
升直接导致了输出功率的下降。
造成这种现象的原因主要有两个:
(
1
)
. LD
的阈值电流会随着温度升高而变大,从而
导致用于发光的电流降低。
(
2
)
. LD
的电
-
光转换效率(外微分量子效率)会随着温度升高而下降。
即
LD
温度越高,发光效率
越低。
图 p>
3.1
中可以看出,
在开机大约
1
小时 之后温度和输出功率仍在缓慢变化这种
变化一直持续到测试结束。
LD< /p>
输出功率的变化会影响实验的结果。
为了增加光源在工作
中的稳定性,需要对光源部分进行散热处理。在仪器中进
行散热的方法可以分为两种:被
动散热和主动散热。被动散热是指加装金属散热片
或者是风扇,保证仪器产生的热量可以
更快地消散到环境中。被动散热的方式实现
较为简单,但是不能将温度降低至室温以下而
且没有办法控制温度。主动散热是指
使用水冷、小型空气压缩机或者半导体制冷器(
p>
Thermo Electric Cooler
,
TEC
< p>)等降
低仪器或者特定空间内的温度。主动散热的方式可以使仪器工作时的温度降低
至室
温以下,并且通过温度反馈模块可以将温度控制在所设定的范围内。水冷和小型空<
/p>
30
华
中
科
技
大
学
博
士
学
位
论
文
p>
气压缩机的方式散热效率高,但实现起来比较复杂,需要增加额外的管道、控制电
路以及发动机等,适合应用于大型仪器的降温制冷。
TEC
的制冷原 理是半导体材料
的帕尔贴效应:当直流电流通过两种半导体材料组成的电偶时,其一端吸
热,一端
放热。其吸热的一段朝向机箱可以进行制冷,放热的一面一般需要加装风扇对它
进
行散热,整个主动散热模块(
TEC
和风扇)正常工作 时只需要外部提供一路
12 V
的
直流供电。由于
fNIRS
仪器的光源机箱本身体积较小,考虑到降低系统复杂性的原
< br>则,本节使用
TEC
作为主动散热方式。
图
2.14
散热片对于光源长时间稳定性影响的测试结果。
Fig. 2.14 The test result of the
influence of the heat sink on the long time
stability of the light source
对
L D
进行降温处理,首先需要将
LD
器件温度散到环境中,然后再通 过
TEC
控制光源机箱整体的温度,即结合被动散热和主动散热两种方式
。尽管光源机箱顶
部有散热板,
但是
LD
本身面积很小,
难以将热量散发到机箱空间中。
对于小型器件
的被动散热一般选择使用散热片,
首先为每个
LD
加装一块金 属散热片,
将金属散热
片卡在
LD
四周,
增加
LD
的散热面积。
在光学仿体上进行测试散热 片对于光源稳定
性的提高效果:首先将光源和探测器固定在仿体上,源探间距为
3
cm
,然后分别在
加散热片以及不加散热片的
情况下测试光源的长时间稳定性。为了避免无关因素的
影响,
测试使用同 一个
LD
和探测器,
两次测试前均保证光源处于室温状态,
测试时
间为
2000
s
。测试 结果如图
2.14
所示。可以看出加散热片的光源工作较为稳定,而
不加散热片的光源需要较长的时间才能进入稳定期。用变异系数衡量两种情况下的
31
华
中
科
技
大
学
博
士
学
位
论
文
光源稳定性,分别为
3.6%
(不加散热片)和
1.4%
(加散热片)
。测试结果说明加散
热片可以有效提高光源稳定性
。从图中可以看出散热片对于光源稳定性的提升主要
是减少光源预热的时间,即加装散热
片后光源可以更快地进入稳定的工作状态。
图
2.15
加入温控后的
LD
光强和表明温度变化曲线。黑色竖线代表开启温控的时间点。
Fig. 2.15 Long time curve of LD power
and surface temperature with temperature control.
The black vertical
line indicate the
time point of opening temperature controller.
加装散热片只能帮助
LD
快速地将表面热量散到环境中,
而在长时间多光源同时
工作的情况下光源机箱的整体温度仍然会有显著的上升。为
了让光源工作在稳定的
温度范围内,
使用
TEC
和风扇对光源机箱进行主动散热,
风扇的作用是将
TEC
发 热
面的热量散到外部。同时使用温度传感器和上位机构成一个闭环的温度反馈模块来
p>
控制
TEC
和风扇的工作情况,使实验中光源机箱的温度稳定 在提前设定好的温度范
围内。温度反馈模块可以接受上位机发送的设定温度的指令,并通
过电磁继电器控
制
TEC
和风扇供电电源的通断。整个主 动散热模块,包括
TEC
,风扇以及温度反馈
模块,只需
要系统自带的
12
V
电源进行供电,不需要加入额外的控制或电 源电路。
测试加入温控后的光强和温度随时间变化曲线。测试时先关闭温控模块让系统正
常
运行,
然后在大约
1
小时后打开温控,
并将温度设定为
32
℃,
记录光强和温度信号< /p>
的变化。二者的长时间变化曲线如图
2.15
所示,打开温 控后
LD
温度开始下降,输
出光强信号幅值随之上升,最
终温度温度稳定在
32
℃左右。打开温控到
LD
温度稳
32
华
中
科
技
大
学
博
士
学
位
论
文
定需要一段时间
,因此在实验中也需要先打开
LD
和温控,等待
5-10
分钟再进行正
式地实验数据采集。
综上
,本节将被动散热与主动降温结合来控制系统光源部分的温度,提高了光
源工作的稳定性
,由于温度和发光功率线性相关,同时也提高了光源的工作效率。
加入温度显示还可以在
实验中实时观测仪器的运行状态,避免系统过热运行。
2.3.2
探测光纤设计
< p>
光纤式
fNIRS
系统的光源和探测器通过光纤与被试的头皮接触, 可以探测头发
覆盖的位置。作为信号获取的关键部件,光纤的传光特性和制作工艺直接影
响到
fNIRS
检测信号的质量
[76]
。制作光纤的材料一般可以选择塑料、玻璃或者石英。塑
料光纤尽管制
作工艺简单,但是传光效率较低,不适用于微弱生理信号探测。综合
考虑传光效率、
p>
柔韧性、
制作成本等因素,
fNIRS
系统的光源光纤 材质选择石英光纤,
探测器选择玻璃光纤束。由于
fNIRS
< p>系统光纤需要传输两个波长的光信号,且传输
距离相对较短
(
1.5~5 m
)
,
因此光源光纤和探测光纤均选择 p>
62.5/125
的多模光纤,
数
值孔径为<
/p>
0.6
。目前使用的光纤探头在使用时固定不方便,容易晃动,而且长时间使
用时包层容易磨损,影响
fNIRS
探测信号的质量。光纤的结 构和材料直接影响到被
试在实验中的体验,本节对已有的光纤进行了重新设计,设计的基
本原则是舒适、
轻便和耐用。在实验中光源光纤的使用数量较少,因此本节的工作主要围
绕探测光
纤开展。
探测用的玻璃光纤束由数千根直径
50
μm
< p>的玻璃丝组成。光纤束的一端通过
FC/APC
接头与
APD
模块耦合在一起,另一端熔融为镜面接收光信号。细微的光纤
束
封装在护套中,使其能够弯曲而不至于断裂。护套一般选择黑色材料实现遮光的
目的。护
套材料的选择在很大程度上影响光纤的柔韧性。为了让被试更加舒适地进
行实验,希望光
纤在保证传光效率的基础上尽可能柔软。原光纤的包裹材料使用的
是热缩管,然而在使用
中发现这种材料导致光纤过于僵硬,对于固定在被试头皮的
光纤有很大的拉扯力。这种拉
扯不仅会降低信号质量,缩短光纤的使用寿命,还会
降低被试在实验中的舒适程度,从而
影响实验结果。因此,本节更换了更加柔软而
33
华
中
科
技
大
学
博
士
学
位
论
文
有韧性的橡胶管来包裹光纤。
探测光纤的末端
有一个
90
°的弯角,弯折后伸出一段接触到头皮。原光纤的伸
< br>出部分较长,不仅增加光纤的重量,使被试佩戴不舒服,而且使探头在实验中容易
发生倾斜。新版光纤通过缩减弯折的长度和倒角的体积以减轻重量。原光纤与被试
头皮接
触的部分为圆形平面,边缘比较尖锐,压迫被试头皮时会产生刺痛,从而影
响被试的配合
程度和脑功能信号。为了解决尖锐边缘产生的压迫感,新版光纤中在
末端用玻璃材料做了
一个弧状倒角处理。另外,探测光纤的末端需要长时间与头皮
或毛发接触,人体分泌的汗
液会对其产生腐蚀效果。为了提高探测光纤的耐用性,
探头末端部分使用重量更轻、硬度
更低、并且抗酸碱腐蚀的聚四氟乙烯作为包裹材
料。
< br>热缩管包层
橡胶管包层
10mm
7mm
普通塑料
22mm
13mm
p>
聚四氟乙烯
原探测光纤
新探测光纤
图
2.16
两版光纤对比图。
Fig. 2.16
Compared graph of two probe fibers.
新旧两
版光纤对比如图
2.16
所示,左侧为原版光纤,右侧为新版光纤。底部的
黄色曲线代表与探头接触的头皮部分。新版光纤更加轻便、耐用、贴附性更好,而
且佩戴更加舒适。被试在舒适的实验环境下能够更好地配合实验,从而能够得到更
可靠的实验数据。
为了测试新版光纤的性能,
在光 学仿体上分别固定光源探头和两版探测器探头,
源探间距均为
3
cm
,然后将光源打开并设定为正常工作状态,观察两个通道检测到
的解调前的波形。结果表明在相同的光源和源探间距下,新版光纤探测到的信号与
原
光纤相比在数量级上保持一致,而且波形稳定。在长时间的实验中新版光纤能够
34
华
中
科
技
大
学
博
士
学
位
论
文
更加稳定地贴合
在头皮表面。说明本节设计的新版探测光纤可以替代原探测光纤用
于血氧参数的检测。<
/p>
在实验中需要用到多根探测光纤和光源光纤,光纤数量一般是左右侧脑区
对称
分布。这些两米长的光纤需要搭在支撑架上,光纤之间会相互挤压,在下面的光纤<
/p>
玻璃束有可能会折断。
而且无规则的放置也为实验中安放对应编号的光纤带 来麻烦。
为了方便光纤的安放以及保护柔软的探测光纤束不被折断,使用黑色套管包裹多
根
光纤。每一个套管中包含左侧或者右侧脑区用到的所有探测光纤和光源光纤,然后
p>
将两个包含光纤的套管搭放在仪器左右两侧的支架上。
以上的研究工作能够良好地抑制系统发热以及探头与头皮接触不良带来的系统
不稳定
噪声,但是在实际实验中仍有可能因为市电电压的跳动或者外部触碰光纤而
导致探测信号
发生波动,这种波动由于偶然性和不可预测性而难以通过特定的方式
来避免。由于这些波
动一般都具有全局性,即对于每个探测通道的影响是一致的,
针对这类系统噪声本节给出
的解决方案是使用一个额外的探测通道来专门记录这类
异常波动。具体做法是将一路探测
光纤进行遮光处理,在暗环境下记录信号。正常
情况下此通道信号应该稳定在
V
附近,当有全局的异常波动产生时该通道可以将
波动情况记录下来供后续离线处理,进而在离线数据分析中可以对这路信号进行噪
声分
析或者将它作为参考信号来进行噪声滤除。比较常规的做法是使用自适应滤波
的方法来去
除探测通道中与参考信号相关的噪声成分。
将探头固定在光学仿体上测
量系统的长时间稳定性,分别使用系统稳定改进前
的方案和系统稳定性改进后的方案进行
测量,并对比结果。在测试时使用遮光布包
裹光学仿体和探头,测试前均对系统进行预热
。使用改进后的方案时加入额外的探
测通道来记录异常波动,并在离线处理时将该通道作
为参考信号使用自适应滤波的
方式去除异常噪声。测试时间为
30
分钟,测试时随机选取光源和探测器配成探测通
道,为了保证结果的鲁棒性两种
实验条件各测试
10
组并取平均。使用变异系数表征
系统
的长时间稳定性,测试结果在表
2.1
中,平均后的结果分别为
0 .78%
(系统改进
前)和
0.11%
( 系统改进后,
p<0.01
)
。结果表明本章的研究工作显著提高 了
fNIRS
系统的长时间稳定性。
35
华
中
科
技
大
学
博
士
学
位
论
文
表
2.1
系统长时间稳定性测试结果。
Table 2.1 Results of system long-term
stability.
测试组
原系统(
%
)
1
2
3
4
5
6
7
8
9
10
0.654
0.544
0.922
0.854
1.027
0.774
0.953
0.497
0.816
0.715
改进后系统
(
%
)
0.072
0.089
0.122
0.204
0.131
0.049
0.062
0.117
0.143
0.081
2.4
混合编程提高
fNIRS
系统软件可视化程度
2.4.1
系统软件需求分析
光纤式
fNIRS
系统的上位机软件是系统中非常重要的组成部 分。原有的系统采
集软件功能单一,采集参数设置困难,一般只能由软件开发者操作,未
经培训的研
究人员难以使用,从而限制了
fNIRS
技术 的应用。另外在实验过程中,研究人员需
要实时地观测被试的脑功能活动并判断他们所处
的状态以及是否配合实验,并根据
在人机交互界面中观察到的信息对被试进行相应的指导
或者对仪器的设置进行调整,
从而避免产生无效的实验数据。因此,提高
fNIRS
系统软件的可视化程度可以提升
fNIRS
的 脑功能检测性能
[77]
。为了使
fNIRS< /p>
系统使用更加方便,功能更加完善,检
测的信号更加有效,本章开发了一套
配合光纤式
fNIRS
系统使用的上位机软件。在
软件中
实现数字锁定、采集卡参数设置、信号显示与保存等功能,使该软件能够在
实验过程中提
供反映被试脑功能活动的可靠信息。
fNIRS
采集软件的整体设计原则是
可视化程度高、使用灵活方便。
本课题组的光纤系统
使用
NI
采集卡来采集光电转换后的信号并送到上位机,
因
此上位机软件使用了和
NI
采集卡可以无缝联接的
LabVIEW
软件编写。
LabVIEW
作
< br>为广泛应用于虚拟仪器领域的编程工具,相比于其他开发工具在仪器控制软件编写
上更加高效,
而且该软件可以方便地实现并行处理。
同时在
Lab VIEW
中使用
Matlab
脚本进行混合编程,进一步
提高程序的执行效率和可扩展性。
Matlab
作为一种数学
软件被广泛应用于
fNIRS
信号处理中,有大量成熟的信号处理范例可供借 鉴
[78]
。
Matlab
LabVIEW
中的子函数,脚本使用
Matl ab
语言编写,通过脚本
的输入输出接口与
LabVIE W
主程序互联,脚本中调用
Matlab
软件对输入数据进行 p>
36
华
中
科
技
大
学
博
士
学
位
论
文
p>
处理。使用混合编程可以直接移植编写好的
Matlab
信号 处理程序,弥补
LabVIEW
难以进行复杂数据操作(比如矩阵乘法、
二维插值等)的缺陷,使软件使用更加灵
活,功能更加完善。
光纤式
fNIRS
系统的人机交互软件作为一种生理信号检测软件 ,除了应具备基
本的数据采集、显示与保存之外,还需要满足以下五项需求:
1
)软件中可以实时
采集并记录每通道的数据以及
任务起止信号,保存路径可以使用默认或者人为设置;
(
2
)
软件中可以实时显示采集到的光强信号以及经过修正比尔朗博定律转换后的血
氧信号,二者可以根据需求实时切换;
(
3
)系统运行的 基本参数,包括采样率、光
源选择、调制频率等可以在界面中进行配置;
(
4
)软件中可以实现一些基本的数据
处理功能,
比如带通滤波和叠加平均,
能够在软件中看到任务中的血氧激活情况;
(
5
)
为了满足光电联合检测的应用需求,软件需要有同步
触发以及与其他上位机通信的
功能。另外,该软件还需要具有较强的可移植性,对于一些
关键性的功能提供二次
开发的接口。
2.4.2
可视化
fNIRS
系统软件设计与实现
信号显示
初始参数设置
信号采集
数据保存
图
2.17 fNIRS
系统软件框架。
Fig. 2.17
Software framework of fNIRS system.
p>
根据上一节的软件设计,
fNIRS
系统软件的基本框架如图
2.17
所示,主要包含
参数设置、信号采集、信号显示
与数据保存四项功能。
(
1
)
初始参数设置
上位机软件在运行之前需要设置一些初始
的运行参数。包括采样率、使用的光
源、被试信息、同步信号等。由于光纤系统使用调制
光源,调制频率在
1 kHz
到
100
kHz
之间,在程序中需要对采集到的原始信号进行数字锁定运算。为了保证数字锁
37