张英，王骏，鲍国强，张春香，王士同

（江南大学 数字媒体学院，江苏 无锡 214122）

自闭症又叫孤独症，是一种广泛性神经失调的发展障碍疾病，其临床表现主要为社会交往障碍、言语交流障碍和动作刻板重复等[1-3]。自1943年利奥·坎纳医生首次对自闭症做出经典性描述以来，该疾病因其病因不明确，又集多种障碍于一身，已经得到了广泛的关注。目前，在美国平均每68名儿童中就有1名自闭症儿童；在世界范围内，自闭症患者超过1%。据统计，约2/3孤独症患儿成年后无法独立生活，需要终生照顾和看护。为此，人们迫切希望寻找科学而高效的手段对自闭症进行准确诊断以达到早发现早治疗的目的。

已有研究表明，自闭症患者行为和认知缺陷与潜在的脑功能异常有关。静息态fMRI技术[4-9]通过使用血氧水平依赖(blood oxygenation level dependent, BOLD)信号可以有效地反映患者在静息状态下脑部血流和代谢活动等功能性变化，因此在自闭症的辅助诊断中起着重要的作用。基于静息态的医学影像诊断大致可分为特征提取和分类器设计这两个阶段。在特征提取阶段，通过使用相关性特征选择等方法从高维医学影像特征中得到有鉴别意义的特征；在分类器设计阶段，通过使用线性或非线性分类器进行分类。但是目前多采用线性特征学习方法，虽然具有简单高效等优点，但是尚不能有效挖掘出医学影像数据中的复杂信息。

本文基于静息态fMRI提出了一种新型的无监督模糊特征学习方法，并在此基础上构造多视角分类器从而实现自闭症的辅助诊断。其技术路线可概括为：首先，将静息态fMRI经过预处理后，基于Pearson相关系数构造低阶功能连接矩阵，并使用文献[7]的方法构造高阶功能连接矩阵；然后，提出新型的无监督模糊特征学习方法，将2种功能连接矩阵分别映射到不同的低维嵌入空间中；最后，以低阶和高阶功能连接在脑疾病诊断中的互补性为出发点，设计多视角SVM[10]分类器进行分类。

1 实验数据集及处理

本文使用的数据集来自ABIDE(autism brain imaging datae xchange)数据库(http://preprocessedconnectomes-project.org/abide/download.html)的NYU子集和USM子集，采集对象的FIQ(full-scale intelligence quotient)范围为65～150。实验所处理对象的相关信息如表1和表2所示。

表1 NYU实验对象信息Table 1 The information of NYU experimental object

表2 USM实验对象信息Table 2 The information of USM experimental object

每个对象的fMRI数据使用DPARSF(data processing assistant for resting-state fMRI)[11]工具进行处理。处理步骤主要包括：1）去除功能磁共振图像序列前10个时间点的数据；2）时间层校正和头部运动校正；3）数据通过统一使用T1加权图像分割并将其归一化到MNI152(montreal neurological institute 152)标准空间中；4）使用AAL (anatomical automatic labeling)模板将大脑划分成116个脑区，每个区域重新采样3×3×3 mm3的体素；5）使用半高全宽高斯核进行空间平滑处理；6）应用带通滤波(0.01～0.1 Hz)去除噪声；7）去线性漂移并进行全局信号校正去干扰变量；8）计算每个脑区的平均时间序列。

2 基于fMRI数据的无监督模糊特征学习方法

2.1 方法流程

基于静息态fMRI进行自闭症辅助诊断包括训练和测试两个阶段。本文所采用的训练过程如图1所示，各步骤包括：

1）根据预处理后的静息态功能磁共振成像数据，提取出各脑区的平均时间序列信号，计算脑区之间的Pearson系数，得到低阶功能连接矩阵。

2）将低阶功能连接矩阵的每一行作为各脑区的特征描述，再次计算脑区之间的Pearson系数，得到高阶功能连接矩阵。

4）使用无监督模糊特征映射方法对低阶和高阶功能连接的特征矩阵进行特征学习，得到相应的变换矩阵，并将和变换到低维空间中得到相应的嵌入矩阵和。

6）将2种不同的核矩阵进行线性组合得到复合核矩阵。

7）构造相应的SVM分类器。

图1 无监督模糊特征学习方法的框架流程图Fig. 1 The framework of an unsupervised fuzzy feature learning method

2.2 低阶和高阶功能连接矩阵的生成

虽然大脑各脑区在空间上相互隔离，但它们之间的神经活动相互影响，通常使用功能连接矩阵来表示脑区之间的时间相关性[12]。本文使用文献[8]中的方法，基于每个脑区fMRI时间序列来计算每个对象各脑区之间的Pearson相关系数，得到低阶功能连接矩阵。将低阶功能连接矩阵的每一行作为一个脑区的特征描述，使用同样的方法再次计算脑区之间的Pearson相关系数，从而得到高阶功能连接矩阵。

2.3 基于流形正则化约束的模糊特征学习

为了在无监督环境中将高维空间中的样本点映射到S维的低维空间，本文引入多输出TSK模糊建模方法。TSK模糊系统是由形式为“IFTHEN”的模糊规则组成的，对于多输出TSK模糊系统，给定任意一个D维输入向量

本文采用高斯函数表示模糊规则的模糊集，即

并使用式(2)、式(3)来计算每个样本相对于每条模糊规则的点火强度：

令，1为全1列向量，多输出TSK模糊系统的建模过程可用式(7)～ (13)的计算方法表述[15]：

式(7)～(13)表述的多输出TSK模糊系统建模过程可以看作将维空间中的数据通过变换矩阵降维到S维数据空间的过程，通常。为了学习得到适合的，本文构造的目标学习准则为

无监督模型模糊特征学习的过程如图2所示。

图2 无监督模糊特征学习流程图Fig. 2 Flowchart of unsupervised fuzzy feature learning

2.4 基于功能连接矩阵的多视角SVM

使用式(21)对核矩阵进行线性组合，得到最终的复合核矩阵：

3 实验结果与分析

3.1 性能指标及参数设置

本节通过实验验证本文方法的有效性并与已有方法进行比较。本文选取精确度(ACC)、AUC值、敏感度(SEN)、特异性(SPE)作为实验的评价指标，实验的参数设置使用嵌套五折交叉验证来进行寻优。相关方法的表述及各参数的详细设置等信息如表3所示。

表3 各算法说明及参数详细设置Table 3 Details and parameter settings of algorithms

本文采用式(22)～(24)计算性能指标，以此进行算法性能评估：

式中：TP和TN分别代表正确分类为自闭症患者的人数和正确分类为正常对照组的人数；FP和FN分别代表正常个体被分类为自闭症患者的人数以及自闭症患者被分类为正常对照组的人数[17]；灵敏度SEN作为衡量该方法准确预测出自闭症患者的能力的指标；特异性SPE是衡量该方法准确预测出健康个体的能力指标。本文实验使用LIBSVM[18]工具包中的高斯核SVM。

3.2 分类结果及分析

从表4和表5中可以看到，UFFM-M2SVC使用无监督模糊特征映射算法UFFM进行特征学习，在此基础上构造分类器M2SVC，可以取得比LOFCC、HOFCC及M2SVC等传统分类器更好的性能。

M2SVC算法联合使用两个视角的数据进行分类，其分类结果要明显优于基于单视角的分类器LOFCC和HOFCC；而本文方法采用新型的模糊特征学习技术进行特征学习，在此基础上使用M2SVC进行分类，从而取得了比仅使用M2SVC更好的效果，这表明本文所使用的模糊特征学习技术能够得到更有鉴别能力的特征，这对自闭症的诊断具有更优的判别能力。

表4 各算法在自闭症分类中的性能比较(NYU)Table 4 The performance comparison of each algorithm in the classification of autism(NYU)

表5 各算法在自闭症分类中的性能比较(USM)Table 5 The performance comparison of each algorithm in the classification of autism(USM)

3.3 特征数目对实验结果的影响

本节研究UFFM所提取的特征数目S对M2SVC性能的影响，图3和图4中给出了算法在数据集NYU和USM上S取不同值时M2SVC的分类结果。从图3和图4可以看出，UFFMM2SVC算法的分类准确率随着S取值的不同而发生变化，图3中最高达到78.38%，图4中最高达到78.95%，均比M2SVC算法提高了约5%。可知，在特征提取数S合理取值的情况下，本文使用UFFM进行无监督模糊特征学习，不仅可以有效保留原始有效的特征信息，而且去除了其中包含的冗余信息，使得自闭症患者和正常对照组之间的特征差异更有区分性。

图3 特征数目对分类结果的影响(NYU)Fig. 3 The influence of the number of feature learning on the classification results (NYU)

图4 特征数目对分类结果的影响(USM)Fig. 4 The influence of the number of feature learning on the classification results (USM)

4 结束语

自闭症患者的脑区神经活动异于常人。以此为出发点，本文基于静息态fMRI提出了一种面向自闭症辅助诊断的无监督模糊特征学习方法UFFM，并与M2SVC相结合，从而实现了自闭症的辅助诊断。相对于目前传统的分类算法，UFFM-M2SVC算法训练出的分类模型具有优越且稳定的分类性能。实验结果表明，本文所使用的无监督模糊特征学习方法能够有效地提取出fMRI数据中的重要特征，从而有助于自闭症的临床辅助诊断。但是，本文仅基于静息态fMRI数据进行实验，针对任务态fMRI数据的分类效果还未可知，有待进一步的研究。