关键词不能为空

当前您在: 大学查询网 > 大学 >

休斯顿大学怎么样华中科技大学博士论文

作者:高考题库网
来源:https://bjmy2z.cn/daxue
2020-11-27 23:08
tags:博士论文

大学生文明上网-大学生文明上网

2020年11月27日发(作者:伍晋南)



附件

1


华中科技大学博士学位论文创新点概述


论文题目:提高全脑

fNIRS- EEG

双模态联合检测性能的研究


学科(专业)

、研究方向:生物医学光子学


论文中的主要创新点(不够可附页)


1

)将小间距去除表层干扰技术应用到光纤

fNIRS

系 统中,并通过在线自适应滤波首次实现了在实


验中实时去除表层干扰。

测 试结果表明该方法可以有效抑制

fNIRS

信号中的表层干扰和全局噪声,


实验中可以观察到不受心跳等噪声干扰的血氧信号。


2

)通过研制光电一体式探头帽与研究光电打标技术实现

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>期






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

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>似的表现。该研究进一步证明基于双模态联合检测的脑连接分析方法可以帮助人们


更全面有效地研究脑功能。


本文围绕信号质量、双模态系统整合、软件 可视化以及数据分析方法等多个方


面对提高双模态联合检测的性能进行研究。本文的工作 可以帮助研究人员使用双模


态联合检测技术更加深入地了解脑功能活动的神经机制。


关键词:

功能近红外光谱术

双模态检测

表层干扰

神经反馈

脑网络

工作记忆



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





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





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






.

................................................. .................................................. ............. 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

本文工作总结

................................ .................................................. ........

9

7


5.2

本文工作的创新点

........................... .................................................. ....

9

8


5.3

工作展望

.

................ .................................................. ...............................

9

8


.

................................................. .................................................. .........

1

00


参考文献

.................................................. .................................................. ...

1

01


附录

作者攻读博士学位期间发表论文与研究

........................................

117




IX





1

绪论


1.1

引言


大脑是生命诞生以来,生物经过数亿年漫长岁月进 化的结果,是高度智能化的


复杂系统。一直以来人们对于大脑的研究从未停止。早期主要 通过行为研究来了解


个体的大脑活动,随着科技发展,越来越多的神经影像技术可以帮助 我们更加全面


地认识脑功能活动的规律和神经机制。对大脑功能的研究不仅可以帮助我们 探索人


类行与思的本质关系,帮助我们理解社会现象,指导我们对精神类疾病的治疗,更


可以促进生物科学、信息科学、人工智能等多种学科技术的交叉发展。


脑功能研究的重要性正逐渐被全世界认可。随着欧盟和美国宣布开展

“< /p>

人类大脑


计划

脑活动图谱计划

,研究脑可以帮助我们有效地治疗各种类型的神经系统 疾


病。目前自闭症、老年痴呆、抑郁症、精神分裂症等神经与精神疾病给个人、家庭


以及社会带来了极大的负担。

为了对这些疾病进行预防和治疗,

需要从分子、

细胞、


网络和脑区各个层次研究和了解大脑。其中,以人 类大脑局部脑区或者大脑整体作


为研究对象的无损脑功能检测技术具有重要的科学研究和 临床应用价值。在众多的


无损脑功能成像技术中,

近红外光谱术

functional Near-infrared Spectroscopy, fNIRS


作为一种新兴的功能成像技术已被广泛地应用于脑科学研究 中。该技术具有对环境


和被试限制小、使用方便、对人体无任何损伤等优势


[1-3]


,特别适合应用于儿童,老


年人及患者进 行脑功能活动检测。因此本文的

fNIRS-EEG

双模态检测技术研究工作


也是围绕

fNIRS

技术开展。

< p>
绪论部分首先对脑功能研究进行概述,包括脑活动相关的生理响应、脑功能检


测技术等。接着介绍本文的主要研究工具

近红外光谱术,以及在此基础上发展 的


fNIRS-EEG

双模态成像技术的研究概况和存在的问题。最后围 绕这些问题介绍本文


的主要研究内容。


1.2

脑功能研究概述


1.2.1

脑功能活动的生理机制


1





大脑是人体高级神经中枢所在地,位于半封闭式的颅腔内,从内到外分别被 软


脑膜、蛛网膜和硬脑膜包裹,并在膜内有脑脊液来缓冲外界震荡。人脑分为左右两


个半球,之间由粗大的神经纤维束

胼胝体将其连接。左右两 个半球的皮层又分为


若干不同的区域,

分别负责不同的功能。

< p>大脑在活动时需要血液输送氧来提供能量,


因此,大脑皮层表面和内部分布有复杂的 血管网络。人脑的功能使人类在进化中拥


有了独一无二的高级智慧。

< /p>


人的大脑有千亿数量的神经细胞和胶质细胞,在脑功能活动时,神经元会产生


动作电位。神经元的生物电活动是人脑调控各种生命活动的基础。神经细胞的胞内

和胞外存在大量的带电离子,这种内外的离子浓度梯度形成了膜电位。在没有外界


刺 激的情况下,膜电位维持在一个恒定的水平,一旦受到外部刺激,膜电位就会发


生改变从 而使神经元产生动作电位。不同神经元之间通过突触进行连接,并且通过


突触来进行信息 的传导。

神经胶质细胞包裹在神经细胞表面,

支持、

保护神经细胞 ,


并维持钾离子的浓度。另外,部分神经胶质细胞的一端贴附在毛细血管壁上,为神


经细胞供给氧和能量


[4]



当大脑某个区域的神经活动增加时,局部的血流、氧代谢、糖代谢会产生相应


的变化。因此,通过测量局部血流和代谢等神经元活动的间接反应也可以了解该部


位脑功能活动情况,相应为血液动力学参数的变化。大脑皮层的某个部位被激活需

要消耗能量来产生(或恢复)动作电位和细胞膜两侧的离子浓度梯度,因此该区域



Oxy-hemoglobin

HbO< /p>


2


)浓度的增加

,

而去氧血红蛋白 (

Deoxy-hemoglobin

Hb

)< /p>


浓度小幅度降低


[1]


大多数研究认为造成这种血氧响应的原因是神经元直接或通过神


经胶质细胞间接地向血管 释放血管活性因子来调节血流


[5]



1.2.2

脑功能成像技术


脑功能成像技术可以动态地检测活体脑的生理活动,对当代认知神经科学的发

展产生深刻而巨大的影响。根据所测量的内容,可以把无创的脑功能成像技术分为


两 大类。第一类是非侵入性电生理技术,第二类是通过测量神经元活动引起的次级


2





反应(如局部葡 萄糖代谢和血流、血氧变化等)研究与行为相关的脑局部神经元的


活动情况


[6]



1

)脑电图(

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>但是由于

MEG

的仪


器非常复杂和精密,并且需要建 造一个严密的电磁场屏蔽室,因此造价高昂,只有


为数不多的单位安装有此设备。


3

)功能磁共振成像

functional Magnetic Resonance Imaging

fMRI

< br>核磁共振现象是由原子核的磁矩与自旋角动量这两者在外磁场的联合作用下形


成的 。功能核磁共振成像一般是依据血氧含量对比度的原理:神经元的活动会引起


一系列血液 动力学参数的变化,

比如脑血流、

脑血容、

血红蛋白浓度等


[8]


其中

HbO

< br>2


是反磁性物质,

Hb

是顺磁性 物质。当

Hb

HbO


2


浓度比降低时,磁化率引起的散


相会减少,从而表现为

MR

信号的增强。这种现象就是血氧水平依赖效应

(Blood


Oxygenation Level Dependent, BOLD)

< p>。

fMRI

最大的优势是空间分辨率较高

2 ~4 mm


但是时间分辨率较低,为秒量级。功能磁 共振成像系统体积庞大,需要建立专门的


电磁屏蔽室来放置仪器。在实验中被试需要躺在 狭小的腔体内,扫描过程中不能发


生肢体运动


[9]

< p>
,因此不适合老人、小孩以及患有幽闭恐惧症的人群。


4

)正电子发射断层显像

Positron Emission Tomography

PET


PET

利用示踪剂测量局部糖代谢、

血流和血氧变化对各种脑活动进行 功能定位。


在成像之前首先注射带有放射性核素标记的代谢生物分子,比如葡萄糖。当带 有放


射性标记的代谢物聚集在相关区域时便开始进行成像


[8]


该区域的功能活动时由于代


谢需要就会聚集带 有标记的葡萄糖。具有很短半衰期的放射性核素开始衰变后,质


子转化为一个中子和一个 正电子,而正电子几乎在瞬间与组织内的电子结合发生湮


灭。根据质能守恒,湮灭后产生 一对飞行方向相反的

γ

光子。这一对光子被环绕在


生物体 外的探测器探测到,根据探测到光子的位置、数量、飞行时间,可以得到这


些代谢生物分 子在体内的分布情况并绘制分布图像。


PET

用于脑功 能检测时,由于人脑在活动时需要消耗能量,活动区域的葡萄糖


分子浓度增加,表现为该 区域图像信号增强。

PET

不仅可以测出大脑活动的区域,


并且还可以定量地测出激活的强度。

PET

利用示踪剂检测组织糖代谢和血液动 力学


参数变化的时间分辨率和空间分辨率都不如

fMRI

由于

PET

具有极高的化学选择性


和灵 敏度,

能够分辨出

10


-9


< /p>


mol

浓度的化合物,

目前主要用于神经递质等物质的测量 。


4





fMRI

类似,

PET

对于 环境和被试的限制也非常大。


5

)功能近红外光谱术

functional Near-infrared Spectroscopy

fNIRS


1977

年,

Jobsis

发现了生物组织对 于近红外波段(

780-1100 nm

)的光具有低吸


收、高散射的光学窗特性


[10]


。人体组织中水的含量占了很 大比例,而水对近红外波


段的光吸收非常少。因此,近红外光照射人脑时,可以穿透头皮 和颅骨进入大脑皮


层,并被脑组织散射出来。散射出来的光相比于入射光的衰减程度和脑 组织中相关


吸收色团的浓度相关。具体关系可以通过修正的比尔朗博定律来表示。入射到 大脑


的近红外光在具有不同光学属性的多层组织中传播,光子的分布近似

香蕉形

。在


近红外光学窗波段,

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]



能够检测多种血氧参数,

并且具有对环境和被试限制小,

制作和运行成本低等优势。


5





fNIRS

技术的局限在于近红外光的穿透深度有限,只能探测大脑灰质部分,而无


法检 测深层脑组织活动


[14]


。因此该技术主要应用于认知心理学 的研究中。


1.3

脑功能的近红外研究


1.3.1

近红外脑功能成像仪器


近红外脑功能成像技术是根据大 脑组织在光学窗波段低吸收高散射的特性来检


测血氧动力学参数,

从而反 映大脑的功能活动。

根据使用不同的近红外光检测方法,


用于脑功能成像 的仪器有三类:时域仪器、频域仪器和连续光仪器。三种仪器的工


作原理如图

< p>1.2

所示。



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]


。由于时域仪器和频域仪器的复杂程度较高,大部分仍处于实验样


机 阶段,商业化程度低,因此较少应用于脑功能的研究中。


连续光仪器使 用直流稳态或者低频调制的光源照射大脑,通过探测大脑出射光


信号的光强变化来反映血 液动力学参数的相对变化,

从而间接地评估大脑功能活动。


连续光仪器采 集到的光强数据通过修正的比尔朗伯定律转化为血氧参数,而且只能


6





对血氧参数的相 对变化进行测量


[19]


。一般情况下脑功能研究关注的是大脑 活动也就


是变化过程,因此连续光技术适用于脑功能活动研究。目前,包括哈佛大学、华 中


科技大学和中科院自动化所在内的多家研究单位都在研究连续光近红外技术和相关


应用


[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

个探测通道的近红外脑功能成像仪器


[26]


并利用


该仪器对视觉和运动皮层的脑功能活动进行检测


[4, 26-28]


该仪器在对指运动范式下测

< br>量到运动皮层的响应结果和

fMRI

测量结果一致,

Science

期刊以

The brain lights up


为标题进行了报道


[29]



在脑功能成像领域,

fNIRS

技术主要应用于认知 心理学的研究以及临床精神类疾


病的研究与诊断中


[12, < /p>

30-32]


。对于认知心理学的研究主要关注前额叶大脑皮层。比如


人们在进行记忆、情感、判断以及逻辑思维类的任务时,前额叶对应的大脑皮层会

< p>
出现血氧激活现象,即

HbO


2


浓 度的显著上升,任务结束后

HbO


2


浓度会逐渐 下降恢


复至基线状态


[33]


。由于前 额叶区域的颅骨较薄,探测难度较低,因此

fNIRS

在该领

域的应用非常广泛。而大脑在认知任务中的表现和许多神经和精神类疾病(比如精


神 分裂、阅读障碍、癫痫等)密切相关


[34,

35]

< br>。通过研究正常人和患者在认知心理学


任务中大脑活动的差异能够帮助研究人员深 入了解这类疾病的神经机制,从而指导


这类疾病的预防、诊断和治疗。随着光纤式

fNIRS

检测技术的发展,目前在临床上


8





fNIRS

技术还可以应用于躯体运动皮层神经活动的研究中,

指导运动功能受损的患者


进行康复训练


[36]


。近几年,

fNIRS

在躯体运动与感觉皮层的脑功能检测开始应用于


脑机接口技 术的研究中,作为辅助定位的手段来提高脑机接口的识别准确率


[37,

38]



另外,随着

fNIRS

系统通道数目的增加和脑连接分析技术的发展,全脑的

fNIRS

< p>
测技术开始应用于神经科学的脑网络研究中,帮助人们深入了解大脑的功能分化和

< br>功能整合


[39-41]


< p>
尽管

fNIRS

技术已经广泛应用于脑功能的研究中,目前仍存在制 约

fNIRS

技术


发展的三个主要问题:

1

由于近红外光需要穿透头皮和颅骨,

fNI RS

检测到的血氧


信号受表层血流影响严重,头皮、颅骨等浅表的信号会 对大脑皮层信号产生干扰,


降低脑功能信号的可靠程度


[42- 44]


2

)由于光纤式

fNIRS

成像技术发展较晚,系统


的长时间稳定性较差,

< p>上位机软件的可视化程度较低;

3

fNI RS

数据分析方法较为


单一。

目前的

fN IRS

研究一般只分析

HbO


2


在任务段的激活情况,

而激活情况只能反


映大脑的功能分化情况而不能反 映大脑的功能整合情况


[45]


。研究并解决以上问题是


提高

fNIRS

技术在脑功能研究中可靠性,进一步发掘

fNIRS

应用潜力的关键,也是


fNIRS- EEG

双模态联合检测技术的基础。


1.4

近红外

-

脑电双模态联合检测的研究

< p>
fNIRS

技术是通过检测神经活动引发的氧代谢来反映脑功能活动,

对神经活动的


响应速度较慢,因此不能确定处于毫秒量级的神经活动的时间特征。大脑 是一个极


其复杂的动态系统,大脑的活动涉及到从分子到大脑皮层整体的多个层次的变化 ,


因此需要更加多元化的脑功能检测技术。随着各种脑功能检测技术的不断发展,近


几年发展出了将多种技术联合在一起检测脑功能的多模态成像技术。将

fNI RS

与其


他脑功能成像技术进行联合检测,可以弥补

fN IRS

技术的缺陷,实现优势互补,帮


助研究人员更加系统、深入地了解 认知任务中的神经加工过程和大脑功能机制。


目前已有一些国外的课题 组将

fNIRS

EEG

联合,

f NIRS

fMRI

联合以及


fNIRS

PET

联合进行认知科学和临床方面的研究,

并 取得了一系列有价值的成果


[46,


47]

< br>。其中

EEG

作为一种发展了多年的脑功能检测技术,商业化仪器和信号处 理方法


9





都较为成熟,具有很高的时间分辨率,直接记录大脑皮层的神经电活动可以在毫秒

< p>
量级区分认知任务不同阶段的神经活动


[48]


fNIRS

EEG

都具有对被试限制小、探


头佩戴方便的优势,因此应用对象可以是儿童或者老年人。它们还都具有良好的环

< p>
境兼容性可以用于床边检测(

fMRI

PET

CT

等只能在特定房间内使用)

,能够提

< p>
供更加贴近生活场景的检测环境,因此适合应用于认知心理学和社会心理学的研究

< br>中


[49]


另外,

fN IRS

检测的内源光信号与

EEG

检测的电生理信号之间不会产生 串扰,


因此这两种技术非常适合联合在一起检测脑功能活动。将这两种技术联合在一起对


脑功能进行同步检测可以综合

fNIRS

EEG

的优势,

同时具备很高的时间分辨率以


及可以接 受的空间分辨率。另外,这两种技术采集到的脑功能信号可以互相补充、


互相印证,在神 经科学和临床研究中具有非常大的应用潜力


[50]



21

世纪初开始一些国外的研究小组使用

fNIRS-E EG

双模态检测技术来研究认知


心理学任务中的大脑活动,比如语言逻辑 、工作记忆、静息与觉醒状态等


[51-53]


。早

< p>
期研究目的主要是通过

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]


然而该技术目前并没有实现真正的光电联合脑机接口,


只是将

fNIRS

和眼电信号提取出来进行分类,而且缺乏两种信号的反馈环节,降低

< br>了被试的训练效果从而影响识别准确率。


总体来说,目前的

fNIRS-EEG

双模态检测技术仍 处于起步阶段,尚无成熟的系


统联合技术以及配套的采集软件,缺少可靠的双模态物理和 信号层面联合手段,存


在采集通道数量少、实验实施困难以及系统非闭环等问题。在双模 态检测技术中


fNIRS

检测技术由于发展时间较晚,

还 额外存在系统稳定性差、

信号受表层干扰严重


等问题,需要开展研究工作 来补齐该技术短板。另外双模态检测技术应用于认知心


理学研究时缺乏合适的数据分析方 法,需要开展更多的工作进行方法建立和验证。


这些问题影响了双模态联合检测的可靠性 ,限制了

fNIRS-EEG

双模态检测技术在认


知神经 科学和脑机接口领域中的应用。本文围绕以上影响双模态联合检测性能的因


素展开研究, 并将全脑双模态联合检测技术应用于大脑功能活动的研究中。


1.5

本文的主要内容


如前所述,

fNIRS-EEG

双模态成像技术对于脑功 能的研究有着非常重要的意义。


本文主要基于自主研发的

fNIRS

系统研究光电联合脑功能检测的方法,

在信号质量、


系统稳定 性、联合方式以及数据分析方法等方面对影响双模态检测性能的问题展开


研究。


11





第一章,绪论。首先说明了脑功能研究的意义,并介绍了常用的几种脑功能检


测技术。基于这些检测技术引出本文的主要研究对象:

fNIRS-EEG

双模态联合成像


技术。介绍了双模态联合成像技术的意义和目前存在的一些问题,然后 引出本文的


主要研究内容。


第二章,提高全脑

fNIRS

检测性能的研究。首先研究了影响

fNIRS

信号的表层


干扰问题,使用小间距配合自适应滤波的方式进行去除。围绕影响系统 稳定性的问


题研究并解决了光源发热问题、重新设计光纤探头。在光纤式

fNIRS

系统基础上开


发了一套功能完善、使用方便的采集软件以及配 套的可视化人机交互界面。最后通


过一个

N-back

实 验验证

fNIRS

技术的在体检测有效性。

第三章,

提高全脑双模态联合检测性能的研究。

首先研制了光电一体式探头帽 ,


通过特制的布帽基底和探头夹具将多通道的

fNIRS

探头以及

EEG

电极固定在一起。


使用光电打标的方式实 现高精度信号同步。围绕双模态脑机接口的应用研究了基于


fNIRS

EEG

信号的在线反馈功能并将这两种信号的采集与显示融合在一个软件平


台中。最后通过一个

Stroop

实验验证双模态联合技术的在 体检测有效性。


第四章,工作记忆中脑功能性别差异的双模态成像研究

。本章使用全脑


fNIRS-EEG

双模态成像技术研究 了语言类工作记忆任务下脑功能的性别差异。然后


使用包括功能连接、

效 应连接和小世界网络特性等方法来分析数据。

最后结合

fNIRS


EEG

的结果对于性别差异的神经机制进行了讨论。


第五章,总结与展望。对本文工作进行归纳总结,介绍本文的研究内容、研究

< p>
成果和主要创新点,并对接下来的研究工作进行展望。


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





集模式的

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


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

表层干扰的去除方法开展研究。这些方


法通常可以分为纯软件的方法和引入参考信 号的方法。纯软件的方法是通过离线的


数据处理方法,包括低通滤波、主成分分析和小波 滤波等去除包含表层干扰在内的


生理噪声,并将任务相关的血氧信号提取出来

< p>
[60,

63-65]


。其中常用的滤波方法是根据


干扰信号的频谱特性将噪声滤除。然而由于表层干扰的成分复杂,其频率范围和大

< p>
脑皮层中血氧信号(目标信号)的频率范围很大一部分是混叠在一起的


[6 6]


,因此只


通过软件滤波的方法很难在保留目标信号的前提下 去除表层干扰。


考虑到干扰信号与目标信号的来源不同,引入参考信号 的方法更加有效。该方


法首先需要加入额外的探测通道采集干扰信号,然后将干扰信号作 为参考信号通过


软件算法进行滤除。根据表层干扰信号来源于浅表的头皮和颅骨,可以仅 提取表层


组织中不包含大脑功能活动的血氧信号,然后将该信号作为参考来去除

fNIRS

信号


中的表层干扰成分


[67 ]


。根据课题组之前使用蒙特卡洛模拟光子在大脑皮层中的传输


分布,

fNIRS

探测深度和光源与探测器的距离有关,当间距只有

1 cm

时该通道的探


测深度为

~0.5 cm

,因此只能接收到来自头皮和颅骨的血氧信号


[68]


。这种比正常间距


15





短的通道称为小间距探测通道。小间距探测通道和正常间距探测通道的探测 原理如


2.2

所示。

其中黄色线条代表 头皮,

红色香蕉形曲线代表不同间距通道探测的光子


分布。通过放置小间 距和正常间距探测器可以分别记录到表层血氧信号和包含表层


血氧信号和脑功能活动的混 合信号。由于小间距通道和正常间距通道一样可以探测


人体组织的血氧浓度变化,只是无 法探测大脑活动的信号。因此小间距通道不仅可


以探测到表层的干扰,还可以接收到来自 系统或者人体的全局噪声,比如光源的温


度漂移、脉搏和呼吸波动等。这类噪声同样可以 通过将小间距通道作为参考的方式


进行滤除。



2.2

小间距通道去除表层干扰原理图。


Fig. 2.2 Schematic diagram of short separation channels removing superficial interface.


2.2.1

构建小间距探测通道


相比于单独使用滤波等信号处理方法,使用小间距通道引入参考信号的方法能


够在不影响目标信号的前提下对表层干扰进行去除或抑制。

本节通过设计新的光纤、


研究信号处理算法、改进采集软件等工作将该方法引入到光纤式

fNIRS< /p>

系统中,并


对其去除表层干扰的效果进行测试。小间距去除表层干扰的方法 首先需要一个或多


个短间距通道将该区域的表层干扰提取出来,然后通过构建数学模型, 将表层信号


作为参考来去除正常通道信号中的表层干扰。


为了小间距去除表层干扰方法应用于光纤式

fNIRS

系统中,首先需 要在系统中


构建小间距探测通道。当光源和探测器间距只有

1 cm

时光子在组织中的平均光程长


度将大大缩短,

因此组织的散射 和吸收对入射光造成衰减相比于正常通道大大减小。


由于光纤式

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>。根据探测面积比例,在理论上小间距光纤接收光强幅值


为正常探测光纤(芯径

2.5 mm

)接收光强的

156

分之一,可以满足探 测光强在

APD


动态范围内的要求。


SMA 905


90

°


18mm


Φ

1mm



2.3

小间距探头示意图。


Fig. 2.3 Sketch map of short separation detection probe.


小间距光纤末端设计为

SMA90 5

接头,与正常间距探测光纤相同,因此可以根


据实际需求灵活地替换机 箱上的正常探测光纤。制作好的小间距光纤首先要测试它


是否能够工作在正常的动态范围 ,测试方法是将光源探头和正常探测探头固定在组


织光学仿体上,间距设置为

< p>3 cm

。将光源调到正常工作状态,然后在距离光源

1 cm



17




处放置小间距光纤,并观察它和正常通道检测到的光信号。测试结果如图

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%

,滤波后为

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

包含


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

图中可以看出自适应


滤波处理后的信号(

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)


y


1


(n),w(n),e(n)



2.10

实时自适应滤波程序框图。


Fig. 2.10 Block diagram of real-time adaptive filter.


25





实时的表层干扰去除第一步是要在程序中嵌入自适应滤波器。出于简化处理过


程,提高程序运行速度的目的,使用一阶的自适应滤波器。根据公式

2.1

2.2

,只


需要当前时间点采集的数据和上一个时间 点采集的数据即可进行得出去除表层干扰


后的数据。由于自适应滤波在执行中需要设置初 始值,因此在程序中加入一个采集


循环运行次数的计数器,将计数器的值作为后续自适应 滤波的执行条件。实时自适


应滤波的程序框图如图

2.10

所示。在第一次采集循环对数据进行初始化。之后的每


一次循环调用上个循环的小间距 通道数据



,滤波系数

w

以及去除表层干扰后的数


e

值,自适应滤波 处理完后更新这个三个值并保存到一个数组中供下一个时间点


的自适应滤波程序来调用。 自适应滤波运算涉及较多矩阵乘法运算,如果使用


LabVIEW

自带的 乘法和

for

循环来实现会拖慢程序的运行速度,从而降低表层干扰


去除的实时性,因此在线的自适应滤波部分使用

Matlab

脚本实现 。每个时间点的滤


波处理之后,将数组

e

按照通道次序分 别送到每个通道的波形图中显示。


使用图

2.7

排布的探头帽进行测试。

测试时让被试保持静息状态。

在静息态中主


要的可见生理噪声为脉搏波


[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






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




而且

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






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

温度越高,发光效率


越低。

3.1

中可以看出,

在开机大约

1

小时 之后温度和输出功率仍在缓慢变化这种


变化一直持续到测试结束。

LD< /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

发 热


面的热量散到外部。同时使用温度传感器和上位机构成一个闭环的温度反馈模块来


控制

TEC

和风扇的工作情况,使实验中光源机箱的温度稳定 在提前设定好的温度范


围内。温度反馈模块可以接受上位机发送的设定温度的指令,并通 过电磁继电器控


TEC

和风扇供电电源的通断。整个主 动散热模块,包括

TEC

,风扇以及温度反馈


模块,只需 要系统自带的

12

V

电源进行供电,不需要加入额外的控制或电 源电路。


测试加入温控后的光强和温度随时间变化曲线。测试时先关闭温控模块让系统正 常


运行,

然后在大约

1

小时后打开温控,

并将温度设定为

32

℃,

记录光强和温度信号< /p>


的变化。二者的长时间变化曲线如图

2.15

所示,打开温 控后

LD

温度开始下降,输


出光强信号幅值随之上升,最 终温度温度稳定在

32

℃左右。打开温控到

LD

温度稳


32





定需要一段时间 ,因此在实验中也需要先打开

LD

和温控,等待

5-10

分钟再进行正


式地实验数据采集。


综上 ,本节将被动散热与主动降温结合来控制系统光源部分的温度,提高了光


源工作的稳定性 ,由于温度和发光功率线性相关,同时也提高了光源的工作效率。


加入温度显示还可以在 实验中实时观测仪器的运行状态,避免系统过热运行。


2.3.2

探测光纤设计

< p>
光纤式

fNIRS

系统的光源和探测器通过光纤与被试的头皮接触, 可以探测头发


覆盖的位置。作为信号获取的关键部件,光纤的传光特性和制作工艺直接影 响到


fNIRS

检测信号的质量


[76]


。制作光纤的材料一般可以选择塑料、玻璃或者石英。塑


料光纤尽管制 作工艺简单,但是传光效率较低,不适用于微弱生理信号探测。综合


考虑传光效率、

柔韧性、

制作成本等因素,

fNIRS

系统的光源光纤 材质选择石英光纤,


探测器选择玻璃光纤束。由于

fNIRS

< p>系统光纤需要传输两个波长的光信号,且传输


距离相对较短

1.5~5 m

因此光源光纤和探测光纤均选择

62.5/125

的多模光纤,


值孔径为< /p>

0.6

。目前使用的光纤探头在使用时固定不方便,容易晃动,而且长时间使


用时包层容易磨损,影响

fNIRS

探测信号的质量。光纤的结 构和材料直接影响到被


试在实验中的体验,本节对已有的光纤进行了重新设计,设计的基 本原则是舒适、


轻便和耐用。在实验中光源光纤的使用数量较少,因此本节的工作主要围 绕探测光


纤开展。


探测用的玻璃光纤束由数千根直径

50

μm

< p>的玻璃丝组成。光纤束的一端通过


FC/APC

接头与

APD

模块耦合在一起,另一端熔融为镜面接收光信号。细微的光纤


束 封装在护套中,使其能够弯曲而不至于断裂。护套一般选择黑色材料实现遮光的


目的。护 套材料的选择在很大程度上影响光纤的柔韧性。为了让被试更加舒适地进


行实验,希望光 纤在保证传光效率的基础上尽可能柔软。原光纤的包裹材料使用的


是热缩管,然而在使用 中发现这种材料导致光纤过于僵硬,对于固定在被试头皮的


光纤有很大的拉扯力。这种拉 扯不仅会降低信号质量,缩短光纤的使用寿命,还会


降低被试在实验中的舒适程度,从而 影响实验结果。因此,本节更换了更加柔软而


33





有韧性的橡胶管来包裹光纤。


探测光纤的末端 有一个

90

°的弯角,弯折后伸出一段接触到头皮。原光纤的伸

< br>出部分较长,不仅增加光纤的重量,使被试佩戴不舒服,而且使探头在实验中容易


发生倾斜。新版光纤通过缩减弯折的长度和倒角的体积以减轻重量。原光纤与被试


头皮接 触的部分为圆形平面,边缘比较尖锐,压迫被试头皮时会产生刺痛,从而影


响被试的配合 程度和脑功能信号。为了解决尖锐边缘产生的压迫感,新版光纤中在


末端用玻璃材料做了 一个弧状倒角处理。另外,探测光纤的末端需要长时间与头皮


或毛发接触,人体分泌的汗 液会对其产生腐蚀效果。为了提高探测光纤的耐用性,


探头末端部分使用重量更轻、硬度 更低、并且抗酸碱腐蚀的聚四氟乙烯作为包裹材


料。

< br>热缩管包层


橡胶管包层


10mm


7mm


普通塑料


22mm


13mm


聚四氟乙烯


原探测光纤


新探测光纤



2.16

两版光纤对比图。


Fig. 2.16 Compared graph of two probe fibers.


新旧两 版光纤对比如图

2.16

所示,左侧为原版光纤,右侧为新版光纤。底部的


黄色曲线代表与探头接触的头皮部分。新版光纤更加轻便、耐用、贴附性更好,而


且佩戴更加舒适。被试在舒适的实验环境下能够更好地配合实验,从而能够得到更

可靠的实验数据。


为了测试新版光纤的性能,

在光 学仿体上分别固定光源探头和两版探测器探头,


源探间距均为

3

cm

,然后将光源打开并设定为正常工作状态,观察两个通道检测到

的解调前的波形。结果表明在相同的光源和源探间距下,新版光纤探测到的信号与


原 光纤相比在数量级上保持一致,而且波形稳定。在长时间的实验中新版光纤能够



34




更加稳定地贴合 在头皮表面。说明本节设计的新版探测光纤可以替代原探测光纤用


于血氧参数的检测。< /p>


在实验中需要用到多根探测光纤和光源光纤,光纤数量一般是左右侧脑区 对称


分布。这些两米长的光纤需要搭在支撑架上,光纤之间会相互挤压,在下面的光纤< /p>


玻璃束有可能会折断。

而且无规则的放置也为实验中安放对应编号的光纤带 来麻烦。


为了方便光纤的安放以及保护柔软的探测光纤束不被折断,使用黑色套管包裹多 根


光纤。每一个套管中包含左侧或者右侧脑区用到的所有探测光纤和光源光纤,然后


将两个包含光纤的套管搭放在仪器左右两侧的支架上。

以上的研究工作能够良好地抑制系统发热以及探头与头皮接触不良带来的系统


不稳定 噪声,但是在实际实验中仍有可能因为市电电压的跳动或者外部触碰光纤而


导致探测信号 发生波动,这种波动由于偶然性和不可预测性而难以通过特定的方式


来避免。由于这些波 动一般都具有全局性,即对于每个探测通道的影响是一致的,


针对这类系统噪声本节给出 的解决方案是使用一个额外的探测通道来专门记录这类


异常波动。具体做法是将一路探测 光纤进行遮光处理,在暗环境下记录信号。正常


情况下此通道信号应该稳定在

< p>0

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

< p>脚本可以看做

LabVIEW

中的子函数,脚本使用

Matl ab

语言编写,通过脚本


的输入输出接口与

LabVIE W

主程序互联,脚本中调用

Matlab

软件对输入数据进行


36





处理。使用混合编程可以直接移植编写好的

Matlab

信号 处理程序,弥补

LabVIEW


难以进行复杂数据操作(比如矩阵乘法、 二维插值等)的缺陷,使软件使用更加灵


活,功能更加完善。


光纤式

fNIRS

系统的人机交互软件作为一种生理信号检测软件 ,除了应具备基


本的数据采集、显示与保存之外,还需要满足以下五项需求:

< p>(

1

)软件中可以实时


采集并记录每通道的数据以及 任务起止信号,保存路径可以使用默认或者人为设置;


2

软件中可以实时显示采集到的光强信号以及经过修正比尔朗博定律转换后的血


氧信号,二者可以根据需求实时切换;

3

)系统运行的 基本参数,包括采样率、光


源选择、调制频率等可以在界面中进行配置;

4

)软件中可以实现一些基本的数据


处理功能,

比如带通滤波和叠加平均,

能够在软件中看到任务中的血氧激活情况;

5


为了满足光电联合检测的应用需求,软件需要有同步 触发以及与其他上位机通信的


功能。另外,该软件还需要具有较强的可移植性,对于一些 关键性的功能提供二次


开发的接口。


2.4.2

可视化

fNIRS

系统软件设计与实现


信号显示


初始参数设置

信号采集


数据保存


2.17 fNIRS

系统软件框架。


Fig. 2.17 Software framework of fNIRS system.



根据上一节的软件设计,

fNIRS

系统软件的基本框架如图

2.17

所示,主要包含


参数设置、信号采集、信号显示 与数据保存四项功能。


1

初始参数设置


上位机软件在运行之前需要设置一些初始 的运行参数。包括采样率、使用的光


源、被试信息、同步信号等。由于光纤系统使用调制 光源,调制频率在

1 kHz

100


kHz

之间,在程序中需要对采集到的原始信号进行数字锁定运算。为了保证数字锁



37

山东大学贾磊-山东大学贾磊


我的大学在线阅读-我的大学在线阅读


南京航空航天大学分数-南京航空航天大学分数


华中农业大学学工处-华中农业大学学工处


扬州职业大学宿舍-扬州职业大学宿舍


上海财经大学会计专业-上海财经大学会计专业


南开大学管理科学与工程-南开大学管理科学与工程


天津外国语大学几本-天津外国语大学几本



本文更新与2020-11-27 23:08,由作者提供,不代表本网站立场,转载请注明出处:https://bjmy2z.cn/daxue/11055.html

华中科技大学博士论文的相关文章