CN115381400A - 一种基于特征选择和svm的脑部疾病患者分类系统和方法 - Google Patents
一种基于特征选择和svm的脑部疾病患者分类系统和方法 Download PDFInfo
- Publication number
- CN115381400A CN115381400A CN202210999624.3A CN202210999624A CN115381400A CN 115381400 A CN115381400 A CN 115381400A CN 202210999624 A CN202210999624 A CN 202210999624A CN 115381400 A CN115381400 A CN 115381400A
- Authority
- CN
- China
- Prior art keywords
- svm
- feature selection
- value
- classifier
- calculating
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/0033—Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room
- A61B5/004—Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room adapted for image acquisition of a particular organ or body part
- A61B5/0042—Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room adapted for image acquisition of a particular organ or body part for the brain
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/05—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves
- A61B5/055—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves involving electronic [EMR] or nuclear [NMR] magnetic resonance, e.g. magnetic resonance imaging
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/16—Devices for psychotechnics; Testing reaction times ; Devices for evaluating the psychological state
- A61B5/165—Evaluating the state of mind, e.g. depression, anxiety
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/40—Detecting, measuring or recording for evaluating the nervous system
- A61B5/4076—Diagnosing or monitoring particular conditions of the nervous system
- A61B5/4082—Diagnosing or monitoring movement diseases, e.g. Parkinson, Huntington or Tourette
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/40—Detecting, measuring or recording for evaluating the nervous system
- A61B5/4076—Diagnosing or monitoring particular conditions of the nervous system
- A61B5/4088—Diagnosing of monitoring cognitive diseases, e.g. Alzheimer, prion diseases or dementia
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7235—Details of waveform analysis
- A61B5/7264—Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7235—Details of waveform analysis
- A61B5/7264—Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems
- A61B5/7267—Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems involving training the classification device
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/70—Arrangements for image or video recognition or understanding using pattern recognition or machine learning
- G06V10/72—Data preparation, e.g. statistical preprocessing of image or video features
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/70—Arrangements for image or video recognition or understanding using pattern recognition or machine learning
- G06V10/764—Arrangements for image or video recognition or understanding using pattern recognition or machine learning using classification, e.g. of video objects
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/70—Arrangements for image or video recognition or understanding using pattern recognition or machine learning
- G06V10/77—Processing image or video features in feature spaces; using data integration or data reduction, e.g. principal component analysis [PCA] or independent component analysis [ICA] or self-organising maps [SOM]; Blind source separation
- G06V10/771—Feature selection, e.g. selecting representative features from a multi-dimensional feature space
Landscapes
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Medical Informatics (AREA)
- Artificial Intelligence (AREA)
- Veterinary Medicine (AREA)
- Animal Behavior & Ethology (AREA)
- Surgery (AREA)
- Public Health (AREA)
- Neurology (AREA)
- Molecular Biology (AREA)
- Heart & Thoracic Surgery (AREA)
- Biomedical Technology (AREA)
- Pathology (AREA)
- Biophysics (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Psychiatry (AREA)
- Evolutionary Computation (AREA)
- Physiology (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Multimedia (AREA)
- Developmental Disabilities (AREA)
- Databases & Information Systems (AREA)
- Computing Systems (AREA)
- Neurosurgery (AREA)
- General Physics & Mathematics (AREA)
- Radiology & Medical Imaging (AREA)
- Software Systems (AREA)
- Mathematical Physics (AREA)
- Signal Processing (AREA)
- Child & Adolescent Psychology (AREA)
- Hospice & Palliative Care (AREA)
- Psychology (AREA)
- Fuzzy Systems (AREA)
- High Energy & Nuclear Physics (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
- Educational Technology (AREA)
Abstract
本发明公开了一种基于特征选择和SVM的脑部疾病患者分类系统和方法，其中方法包括基于患者的fMRI脑影像数据，计算ALFF，ReHo和FC值、根据解剖自动贴标模板提取ALFF，ReHo和FC值作为原始特征数据、基于原始特征数据，进一步进行特征选择，将选择后的特征作为最终的分类特征、选择相应的分类器和结果计算与输出的步骤；系统包括特征选择模块，分类器选择模块和结果计算与输出模块。本发明通过对患者组和健康对照组的rs‑fMRI数据进行分析处理，分别计算出ALFF、ReHo和ALFF值，选取当中具有显著差异的有效特征作为分类特征，进而利用SVM搭建分类器进行模式识别。最后将分类器组装成一种基于特征选择和SVM的脑部疾病患者分类系统，以此帮助医生对疾病的诊断治疗。
Description
技术领域
本发明涉及脑成像数据处理与分析技术领域，尤其涉及一种基于特征选择和 SVM的脑部疾病患者分类系统和方法。
背景技术
帕金森、阿尔茨海默症、精神分裂症、抑郁症和脑白质病变等是一种以精神障碍为主要特征的一类疾病，此类疾病的治愈几率比较小，一般需要早期治疗。静息状态功能磁共振成像(rs-fMRI)通过检测血氧水平依赖(BOLD)信号的自发波动，揭示神经元在静息状态下的自发活动现象。对于rs-fMRI最常用的研究方法有低频振荡振幅(ALFF)、局部一致性(ReHo)和功能连接(FC)等。ALFF 能够可靠的测量自发BOLD信号强度，反映的是静息状态下局部脑区的自发活动强度；ReHo描述的是给定体素与周围体素之间的脑活动同步性，可直接观察到局部脑区神经活动的时间同步性；FC描述的是不同脑区的功能性信号在时间维度上的相关程度，反映不同脑区之间的连接强度。目前，已有研究表明帕金森、阿尔茨海默症、精神分裂症、抑郁症和脑白质病变等存在ALFF、ReHo和FC的异常，并且对于帕金森和抑郁症有着潜在的诊断潜力。因此，发明出一种快速精准的自动分类系统能够有效的帮助医生诊断治疗疾病。
随着机器学习技术的快速发展，将机器学习应用于神经成像数据的研究引起了研究者的兴趣。近年来，机器学习方法尤其是多变量模式分析(MVPA)被用于提取神经精神疾病的结构或功能特征，为探索一种基于神经影像数据的神经精神疾病高效诊断方法提供了一个有前景的方向。MVPA是一种数据驱动技术，能够获得大脑活动的空间分辨模式，在fMRI数据分析中起着至关重要的作用。此外，MVPA使我们能够识别出在个体分类中贡献最大的特征，这可能为研究群体差异提供新的见解。支持向量机(SVM)是机器学习方法之一，与其他机器学习方法相比，SVM在识别复杂数据集中的细微模式方面非常强大。
发明内容
基于现有技术中存在的不足之处，本发明提供了一种基于特征选择和SVM 的脑部疾病患者分类系统和方法，具体技术方案如下：
一种基于特征选择和SVM的脑部疾病患者分类方法，包括以下步骤：
步骤1：基于患者的fMRI脑影像数据，计算ALFF，ReHo和FC值；
步骤2：根据解剖自动贴标模板提取ALFF，ReHo和FC值，作为原始特征数据；
步骤3：基于原始特征数据，进一步进行特征选择；
步骤4：将选择后的特征作为最终的分类特征，然后进行SVM分类器的选择；
步骤5：结果计算与输出。
具体的，所述计算ALFF，ReHo和FC值还包括预处理的步骤，包括去剔除前 10个时间点、头动校正、回归协变量、归一化、平滑处理、去线性趋势和去除协变量。
具体的，所述ALFF值的计算包括：
将fMRI测量所得的时间序列信号进行傅里叶变换，得到0.01-0.08Hz之间的功率谱；
计算功率谱中每个频率的平方根，并在0.01–0.08Hz内获得的均方根即为 ALFF值；
计算公式为：
其中，ak和bk均为不同频率下所对应的系数，N为体素总数，f为时间序列信号；
所述ReHo值的计算包括：
采用肯德尔一致性系数来度量给定体素的时间序列与周围体素的时间序列的相似性；
计算公式为：
其中，其中W表示肯德尔一致性系数，取值范围(0，1)，n表示时间点个数，K表示特定体素与周围相邻体素的个数，其取值有7、19和27，Ri是第i 个时间点的各个体素点的体素值的等级总数，
是R的平均值；
所述FC值的计算包括：
计算两个感兴趣区域之间时间序列的皮尔逊相关性。
具体的，所述根据解剖自动贴标模板提取ALFF，ReHo和FC值包括以下步骤：
步骤21：将解剖自动贴标(AAL)作为掩码，将大脑区域分为90个大脑区域和26个小脑区域
步骤22：将每个受试者的ALFF值、ReHo值分别与AAL相乘，求取每个大脑区域内的均值即为每个脑区的ALFF值、ReHo值的M×N的矩阵
其中，M为受试者数目，特征数目N为116；
步骤23：提取N×N的FC矩阵的上三角元素作为该受试者的FC特征，将所有受试者的FC特征组合为一个M×H的矩阵，即为FC值的原始特征数据；
其中N为116，M为受试者数目，特征数目H为6670。
具体的，所述特征选择包括一种基于F-分数特征选择方法，利用F-分数数值对特征进行排序，F-分数数值越大，表示此特征的辨别能力越强；
计算公式为：
其中
分别是整个数据集、病人组数据集和正常人组数据集的第i个特征的平均值，xk,i (+)表示病人组第i个特征的第K个数据；xk,i (-)表示正常人组第i个特征的第K个数据。
具体的，所述特征选择还包括一种基于T-分数特征选择方法：
T-分数使用双样本t检验返回原假设的检验决策，根据检验统计量信息对特征从高到低进行排序。
具体的，所述特征选择还包括一种基于相关性特征选择方法：
使用corr函数计算两条特征之间的皮尔逊相关系数，并根据皮尔逊相关系数对特征由高到低进行排序。
具体的，所述SVM分类器选择基于libsvm工具包搭建，使用svmtrain和svmpredict函数进行SVM训练和测试；所述SVM分类器根据核函数的不同，可分为核函数为linear的线性SVM和核函数为高斯径向基核函数的RBF SVM；
所述RBF SVM还包括：
使用SVMcgForClass函数在(-8，8)之间以0.2的步长对惩罚因子C和不敏感损失函数的参数g进行优化；
所述分类器的选择还包括使用留一法交叉验证(LOOCV)的步骤，对分类器的性能进行优化。
具体的，所述结果计算与输出包括ROC曲线绘制及AUC计算、计算平均正确率和计算分类器预测结果的敏感性和特异性的子步骤；
所述ROC曲线绘制及AUC计算子步骤包括：
计算并绘制真阳性率(Y轴)和假阳性率(X轴)的ROC曲线；
根据ROC曲线以下的面积，得到AUC值，用于评估整个分类过程的性能；
所述计算平均正确率子步骤包括：
求取每次LOOCV返回正确率的平均值，作为SVM分类器的平均正确率；
所述计算分类器预测结果的敏感性和特异性子步骤，计算公式为：
其中，Accuracy表示每次LOOCV的正确率，Sensitivity表示将患者正确识别为患者的概率，Specificity表示将正常人正确识别为正常人的概率；TP为预测的正确正类数量，FN表示预测的错误正类数量；TN为预测的正确负类数量， FP为预测的错误负类数量。
一种基于特征选择和SVM的脑部疾病患者分类系统，用以实现权利要求1～9 任意一项所述的一种基于特征选择和SVM的脑部疾病患者分类方法，包括特征选择模块，分类器选择模块和结果计算与输出模块，其中，
所述特征选择模块包括F-分数特征选择子模块、T-分数特征选择子模块和相关性特征选择子模块；
所述F-分数特征选择子模块利用F-分数数值对特征从高到低进行排序；
所述T-分数特征选择子模块根据检验统计量信息对特征从高到低进行排序；
所述相关性特征选择子模块根据皮尔逊相关系数将特征由高到低进行排序；
所述分类器选择模块包括线性SVM子模块、RBF SVM子模块和留一法交叉验证子模块；
所述线性SVM子模块选定核函数为linear，对模型进行训练；
所述RBF SVM子模块选定核函数为高斯径向基核函数，对模型进行训练；
所述留一法交叉验证子模块对分类器的性能进行优化；
所述结果计算与输出模块包括ROC曲线绘制及AUC计算子模块、计算平均正确率子模块和计算分类器预测结果的敏感性和特异性子模块；
所述ROC曲线绘制及AUC计算子模块用于评估分类器的分类性能，输出ROC 曲线图和AUC值；
所述计算平均正确率子模块使用LOOCV返回正确率的平均值；
所述计算分类器预测结果的敏感性和特异性子模块输出Accuracy、 Sensitivity和Specificity三个数值，体现对患者的诊断结果。
本发明的有益效果：
本发明通过对患者组和健康对照组的rs-fMRI数据进行分析处理，分别计算出ALFF、ReHo和ALFF值，选取当中具有显著差异的有效特征作为分类特征，进而利用SVM搭建分类器进行模式识别。最后将分类器组装成一种基于特征选择和SVM的脑部疾病患者分类系统，以此帮助医生对疾病的诊断治疗。
附图说明
图1所示是本发明的系统模型示意图。
图2所示为本发明的线性SVM分类流程图。
图3所示为本发明的RBF SVM分类流程图。
具体实施方式
为了对本发明的技术特征、目的和效果有更加清楚的理解，现对照附图说明本发明的具体实施方式。
如图1-3所示，本发明提出一种基于特征选择和SVM的脑部疾病患者分类方法，包括以下步骤：
步骤1：基于患者的fMRI脑影像数据，计算ALFF，ReHo和FC值；
步骤2：根据解剖自动贴标模板提取ALFF，ReHo和FC值，作为原始特征数据；
步骤3：基于原始特征数据，进一步进行特征选择；
步骤4：将选择后的特征作为最终的分类特征，然后进行SVM分类器的选择；
步骤5：结果计算与输出。
具体的，所述计算ALFF，ReHo和FC值还包括预处理的步骤，包括去剔除前 10个时间点、头动校正、回归协变量、归一化、平滑处理、去线性趋势和去除协变量。
具体的，所述ALFF值的计算包括：
将fMRI测量所得的时间序列信号进行傅里叶变换，得到0.01-0.08Hz之间的功率谱；
计算功率谱中每个频率的平方根，并在0.01–0.08Hz内获得的均方根即为 ALFF值；
计算公式为：
其中，ak和bk均为不同频率下所对应的系数，N为体素总数，f为时间序列信号；
所述ReHo值的计算包括：
采用肯德尔一致性系数来度量给定体素的时间序列与周围体素的时间序列的相似性；
计算公式为：
其中，其中W表示肯德尔一致性系数，取值范围(0，1)，n表示时间点个数，K表示特定体素与周围相邻体素的个数，其取值有7、19和27，Ri是第i 个时间点的各个体素点的体素值的等级总数，
是R的平均值；
所述FC值的计算包括：
计算两个感兴趣区域之间时间序列的皮尔逊相关性。
具体的，所述根据解剖自动贴标模板提取ALFF，ReHo和FC值包括以下步骤：
步骤21：将解剖自动贴标(AAL)作为掩码，将大脑区域分为90个大脑区域和26个小脑区域
步骤22：将每个受试者的ALFF值、ReHo值分别与AAL相乘，求取每个大脑区域内的均值即为每个脑区的ALFF值、ReHo值的M×N的矩阵
其中，M为受试者数目，特征数目N为116；
步骤23：提取N×N的FC矩阵的上三角元素作为该受试者的FC特征，将所有受试者的FC特征组合为一个M×H的矩阵，即为FC值的原始特征数据；
其中N为116，M为受试者数目，特征数目H为6670。
具体的，所述特征选择包括一种基于F-分数特征选择方法，利用F-分数数值对特征进行排序，F-分数数值越大，表示此特征的辨别能力越强；
计算公式为：
其中
分别是整个数据集、病人组数据集和正常人组数据集的第i个特征的平均值，xk,i (+)表示病人组第i个特征的第K个数据；xk,i (-)表示正常人组第i个特征的第K个数据。
具体的，所述特征选择还包括一种基于T-分数特征选择方法：
T-分数使用双样本t检验返回原假设的检验决策，根据检验统计量信息对特征从高到低进行排序。
具体的，所述特征选择还包括一种基于相关性特征选择方法：
使用corr函数计算两条特征之间的皮尔逊相关系数，并根据皮尔逊相关系数对特征由高到低进行排序。
具体的，所述SVM分类器选择基于libsvm工具包搭建，使用svmtrain和svmpredict函数进行SVM训练和测试；所述SVM分类器根据核函数的不同，请着重参阅图2，可分为核函数为linear的线性SVM和核函数为高斯径向基核函数的RBF SVM；
请着重参阅图3，所述RBF SVM还包括：
使用SVMcgForClass函数在(-8，8)之间以0.2的步长对惩罚因子C和不敏感损失函数的参数g进行优化；
所述分类器的选择还包括使用留一法交叉验证(LOOCV)的步骤，对分类器的性能进行优化。
具体的，所述结果计算与输出包括ROC曲线绘制及AUC计算、计算平均正确率和计算分类器预测结果的敏感性和特异性的子步骤；
所述ROC曲线绘制及AUC计算子步骤包括：
计算并绘制真阳性率(Y轴)和假阳性率(X轴)的ROC曲线；
根据ROC曲线以下的面积，得到AUC值，用于评估整个分类过程的性能；
所述计算平均正确率子步骤包括：
求取每次LOOCV返回正确率的平均值，作为SVM分类器的平均正确率；
所述计算分类器预测结果的敏感性和特异性子步骤，计算公式为：
其中，Accuracy表示每次LOOCV的正确率，Sensitivity表示将患者正确识别为患者的概率，Specificity表示将正常人正确识别为正常人的概率；TP为预测的正确正类数量，FN表示预测的错误正类数量；TN为预测的正确负类数量， FP为预测的错误负类数量。
请着重参阅图1，一种基于特征选择和SVM的脑部疾病患者分类系统，用以实现权利要求1～9任意一项所述的一种基于特征选择和SVM的脑部疾病患者分类方法，包括特征选择模块，分类器选择模块和结果计算与输出模块，其中，
所述特征选择模块包括F-分数特征选择子模块、T-分数特征选择子模块和相关性特征选择子模块；
所述F-分数特征选择子模块利用F-分数数值对特征从高到低进行排序；
所述T-分数特征选择子模块根据检验统计量信息对特征从高到低进行排序；
所述相关性特征选择子模块根据皮尔逊相关系数将特征由高到低进行排序；
所述分类器选择模块包括线性SVM子模块、RBF SVM子模块和留一法交叉验证子模块；
所述线性SVM子模块选定核函数为linear，对模型进行训练；
所述RBF SVM子模块选定核函数为高斯径向基核函数，对模型进行训练；
所述留一法交叉验证子模块对分类器的性能进行优化；
所述结果计算与输出模块包括ROC曲线绘制及AUC计算子模块、计算平均正确率子模块和计算分类器预测结果的敏感性和特异性子模块；
所述ROC曲线绘制及AUC计算子模块用于评估分类器的分类性能，输出ROC 曲线图和AUC值；
所述计算平均正确率子模块使用LOOCV返回正确率的平均值；
所述计算分类器预测结果的敏感性和特异性子模块输出Accuracy、 Sensitivity和Specificity三个数值，体现对患者的诊断结果。
以上显示和描述了本发明的基本原理和主要特征和本发明的优点。本行业的技术人员应该了解，本发明不受上述实施例的限制，上述实施例和说明书中描述的只是说明本发明的原理，在不脱离本发明精神和范围的前提下，本发明还会有各种变化和改进，这些变化和改进都落入要求保护的本发明范围内。
Claims (10)
1.一种基于特征选择和SVM的脑部疾病患者分类方法，其特征在于，包括以下步骤：
步骤1：基于患者的fMRI脑影像数据，计算ALFF，ReHo和FC值；
步骤2：根据解剖自动贴标模板提取ALFF，ReHo和FC值，作为原始特征数据；
步骤3：基于原始特征数据，进一步进行特征选择；
步骤4：将选择后的特征作为最终的分类特征，然后进行SVM分类器的选择；
步骤5：结果计算与输出。
2.根据权利要求1所述的一种基于特征选择和SVM的脑部疾病患者分类方法，其特征在于，所述计算ALFF，ReHo和FC值还包括预处理的步骤，包括去剔除前10个时间点、头动校正、回归协变量、归一化、平滑处理、去线性趋势和去除协变量。
3.根据权利要求1所述的一种基于特征选择和SVM的脑部疾病患者分类方法，其特征在于，所述ALFF值的计算包括：
将fMRI测量所得的时间序列信号进行傅里叶变换，得到0.01-0.08Hz之间的功率谱；
计算功率谱中每个频率的平方根，并在0.01–0.08Hz内获得的均方根即为ALFF值；
计算公式为：
其中，ak和bk均为不同频率下所对应的系数，N为体素总数，f为时间序列信号；
所述ReHo值的计算包括：
采用肯德尔一致性系数来度量给定体素的时间序列与周围体素的时间序列的相似性；
计算公式为：
其中，其中W表示肯德尔一致性系数，取值范围(0，1)，n表示时间点个数，K表示特定体素与周围相邻体素的个数，其取值有7、19和27，Ri是第i个时间点的各个体素点的体素值的等级总数，
是R的平均值；
所述FC值的计算包括：
计算两个感兴趣区域之间时间序列的皮尔逊相关性。
4.根据权利要求3所述的一种基于特征选择和SVM的脑部疾病患者分类方法，其特征在于，所述根据解剖自动贴标模板提取ALFF，ReHo和FC值包括以下步骤：
步骤21：将解剖自动贴标(AAL)作为掩码，将大脑区域分为90个大脑区域和26个小脑区域
步骤22：将每个受试者的ALFF值、ReHo值分别与AAL相乘，求取每个大脑区域内的均值即为每个脑区的ALFF值、ReHo值的M×N的矩阵；
其中，M为受试者数目，特征数目N为116；
步骤23：提取N×N的FC矩阵的上三角元素作为该受试者的FC特征，将所有受试者的FC特征组合为一个M×H的矩阵，即为FC值的原始特征数据；
其中N为116，M为受试者数目，特征数目H为6670。
5.根据权利要求1所述的一种基于特征选择和SVM的脑部疾病患者分类方法，其特征在于，所述特征选择包括一种基于F-分数特征选择方法，利用F-分数数值对特征进行排序，F-分数数值越大，表示此特征的辨别能力越强；
计算公式为：
其中
分别是整个数据集、病人组数据集和正常人组数据集的第i个特征的平均值，xk,i (+)表示病人组第i个特征的第K个数据；xk,i (-)表示正常人组第i个特征的第K个数据。
6.根据权利要求1所述的一种基于特征选择和SVM的脑部疾病患者分类方法，其特征在于，所述特征选择还包括一种基于T-分数特征选择方法：
T-分数使用双样本t检验返回原假设的检验决策，根据检验统计量信息对特征从高到低进行排序。
7.根据权利要求1所述的一种基于特征选择和SVM的脑部疾病患者分类方法，其特征在于，所述特征选择还包括一种基于相关性特征选择方法：
使用corr函数计算两条特征之间的皮尔逊相关系数，并根据皮尔逊相关系数对特征由高到低进行排序。
8.根据权利要求1所述的一种基于特征选择和SVM的脑部疾病患者分类方法，其特征在于，所述SVM分类器选择基于libsvm工具包搭建，使用svmtrain和svmpredict函数进行SVM训练和测试；所述SVM分类器根据核函数的不同，可分为核函数为linear的线性SVM和核函数为高斯径向基核函数的RBF SVM；
所述RBF SVM还包括：
使用SVMcgForClass函数在(-8，8)之间以0.2的步长对惩罚因子C和不敏感损失函数的参数g进行优化；
所述分类器的选择还包括使用留一法交叉验证(LOOCV)的步骤，对分类器的性能进行优化。
9.根据权利要求1所述的一种基于特征选择和SVM的脑部疾病患者分类方法，其特征在于，所述结果计算与输出包括ROC曲线绘制及AUC计算、计算平均正确率和计算分类器预测结果的敏感性和特异性的子步骤；
所述ROC曲线绘制及AUC计算子步骤包括：
计算并绘制真阳性率(Y轴)和假阳性率(X轴)的ROC曲线；
根据ROC曲线以下的面积，得到AUC值，用于评估整个分类过程的性能；
所述计算平均正确率子步骤包括：
求取每次LOOCV返回正确率的平均值，作为SVM分类器的平均正确率；
所述计算分类器预测结果的敏感性和特异性子步骤，计算公式为：
其中，Accuracy表示每次LOOCV的正确率，Sensitivity表示将患者正确识别为患者的概率，Specificity表示将正常人正确识别为正常人的概率；TP为预测的正确正类数量，FN表示预测的错误正类数量；TN为预测的正确负类数量，FP为预测的错误负类数量。
10.一种基于特征选择和SVM的脑部疾病患者分类系统，用以实现权利要求1～9任意一项所述的一种基于特征选择和SVM的脑部疾病患者分类方法，其特征在于，包括特征选择模块，分类器选择模块和结果计算与输出模块，其中，
所述特征选择模块包括F-分数特征选择子模块、T-分数特征选择子模块和相关性特征选择子模块；
所述F-分数特征选择子模块利用F-分数数值对特征从高到低进行排序；
所述T-分数特征选择子模块根据检验统计量信息对特征从高到低进行排序；
所述相关性特征选择子模块根据皮尔逊相关系数将特征由高到低进行排序；
所述分类器选择模块包括线性SVM子模块、RBF SVM子模块和留一法交叉验证子模块；
所述线性SVM子模块选定核函数为linear，对模型进行训练；
所述RBF SVM子模块选定核函数为高斯径向基核函数，对模型进行训练；
所述留一法交叉验证子模块对分类器的性能进行优化；
所述结果计算与输出模块包括ROC曲线绘制及AUC计算子模块、计算平均正确率子模块和计算分类器预测结果的敏感性和特异性子模块；
所述ROC曲线绘制及AUC计算子模块用于评估分类器的分类性能，输出ROC曲线图和AUC值；
所述计算平均正确率子模块使用LOOCV返回正确率的平均值；
所述计算分类器预测结果的敏感性和特异性子模块输出Accuracy、Sensitivity和Specificity三个数值，体现对患者的诊断结果。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210999624.3A CN115381400A (zh) | 2022-08-19 | 2022-08-19 | 一种基于特征选择和svm的脑部疾病患者分类系统和方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210999624.3A CN115381400A (zh) | 2022-08-19 | 2022-08-19 | 一种基于特征选择和svm的脑部疾病患者分类系统和方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN115381400A true CN115381400A (zh) | 2022-11-25 |
Family
ID=84120151
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210999624.3A Pending CN115381400A (zh) | 2022-08-19 | 2022-08-19 | 一种基于特征选择和svm的脑部疾病患者分类系统和方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115381400A (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116936035A (zh) * | 2023-07-14 | 2023-10-24 | 南京脑科医院 | 一种ai引导下的抑郁症认知行为治疗效果评估系统 |
-
2022
- 2022-08-19 CN CN202210999624.3A patent/CN115381400A/zh active Pending
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116936035A (zh) * | 2023-07-14 | 2023-10-24 | 南京脑科医院 | 一种ai引导下的抑郁症认知行为治疗效果评估系统 |
| CN116936035B (zh) * | 2023-07-14 | 2024-04-12 | 南京脑科医院 | 一种ai引导下的抑郁症认知行为治疗效果评估系统 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Ghazal et al. | Alzheimer disease detection empowered with transfer learning | |
| Gunduz | An efficient dimensionality reduction method using filter-based feature selection and variational autoencoders on Parkinson's disease classification | |
| Wang et al. | Detection analysis of epileptic EEG using a novel random forest model combined with grid search optimization | |
| Das et al. | Schizophrenia detection technique using multivariate iterative filtering and multichannel EEG signals | |
| Rubin et al. | Recognizing abnormal heart sounds using deep learning | |
| Miranda et al. | A survey of medical image classification techniques | |
| Kaur et al. | Diagnosis of Parkinson’s disease using deep CNN with transfer learning and data augmentation | |
| CN110188836B (zh) | 一种基于变分自编码器的脑功能网络分类方法 | |
| Pourbabaee et al. | Feature leaning with deep convolutional neural networks for screening patients with paroxysmal atrial fibrillation | |
| Rehman et al. | Notice of Retraction: Infrared Sensing Based Non-Invasive Initial Diagnosis of Chronic Liver Disease Using Ensemble Learning | |
| CN111513717A (zh) | 一种大脑功能状态提取方法 | |
| Peixoto et al. | Automatic classification of pulmonary diseases using a structural co-occurrence matrix | |
| Rajinikanth et al. | Hand-sketchs based Parkinson's disease screening using lightweight deep-learning with two-Fold training and fused optimal features | |
| Carvajal-Castaño et al. | Effective detection of abnormal gait patterns in Parkinson's disease patients using kinematics, nonlinear, and stability gait features | |
| Houssein et al. | EEG signals classification for epileptic detection: a review | |
| CN115381400A (zh) | 一种基于特征选择和svm的脑部疾病患者分类系统和方法 | |
| Kanse et al. | HG-SVNN: harmonic genetic-based support vector neural network classifier for the glaucoma detection | |
| Pradhan | A Novel Threshold based Method for Vessel Intensity Detection and Extraction from Retinal Images | |
| Shukla et al. | An improved decision support system for identification of abnormal EEG signals using a 1D convolutional neural network and savitzky-golay filtering | |
| Cimr et al. | Enhancing EEG signal analysis with geometry invariants for multichannel fusion | |
| Shekerbek et al. | APPLICATION OF MATHEMATICAL METHODS AND MACHINE LEARNING ALGORITHMS FOR CLASSIFICATION OF X-RAY IMAGES. | |
| Li et al. | Diagnosis of Alzheimer's disease by feature weighted-LSTM: a preliminary study of temporal features in brain resting-state fMRI | |
| Omar et al. | Enhancing EEG signals classification using LSTM‐CNN architecture | |
| Samiei et al. | A complex network approach to time series analysis with application in diagnosis of neuromuscular disorders | |
| Yang et al. | Functional Connectivity Based Classification for Autism Spectrum Disorder Using Spearman’s Rank Correlation |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination |