CN113598794A - 一种冰毒成瘾者检测模型的训练方法和系统 - Google Patents
一种冰毒成瘾者检测模型的训练方法和系统 Download PDFInfo
- Publication number
- CN113598794A CN113598794A CN202110925526.0A CN202110925526A CN113598794A CN 113598794 A CN113598794 A CN 113598794A CN 202110925526 A CN202110925526 A CN 202110925526A CN 113598794 A CN113598794 A CN 113598794A
- Authority
- CN
- China
- Prior art keywords
- picture
- ice
- layer
- network model
- electroencephalogram
- 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
- 238000001514 detection method Methods 0.000 title claims abstract description 84
- 238000012549 training Methods 0.000 title claims abstract description 64
- 229940079593 drug Drugs 0.000 title claims abstract description 63
- 239000003814 drug Substances 0.000 title claims abstract description 63
- 238000000034 method Methods 0.000 title claims abstract description 35
- 206010012335 Dependence Diseases 0.000 claims abstract description 26
- 230000000638 stimulation Effects 0.000 claims abstract description 23
- 238000012360 testing method Methods 0.000 claims abstract description 16
- 230000000763 evoking effect Effects 0.000 claims abstract description 13
- 230000000007 visual effect Effects 0.000 claims abstract description 12
- 238000007781 pre-processing Methods 0.000 claims abstract description 6
- 238000013527 convolutional neural network Methods 0.000 claims description 9
- 238000011176 pooling Methods 0.000 claims description 9
- 208000006379 syphilis Diseases 0.000 claims description 8
- 238000001914 filtration Methods 0.000 claims description 6
- 238000012880 independent component analysis Methods 0.000 claims description 6
- 238000002474 experimental method Methods 0.000 claims description 5
- 230000008569 process Effects 0.000 claims description 5
- 238000012795 verification Methods 0.000 claims description 5
- 230000001988 toxicity Effects 0.000 claims description 4
- 231100000419 toxicity Toxicity 0.000 claims description 4
- 238000002790 cross-validation Methods 0.000 description 8
- 206010013663 drug dependence Diseases 0.000 description 7
- 208000011117 substance-related disease Diseases 0.000 description 7
- 230000006870 function Effects 0.000 description 6
- 238000002599 functional magnetic resonance imaging Methods 0.000 description 6
- 229960001252 methamphetamine Drugs 0.000 description 5
- MYWUZJCMWCOHBA-VIFPVBQESA-N methamphetamine Chemical compound CN[C@@H](C)CC1=CC=CC=C1 MYWUZJCMWCOHBA-VIFPVBQESA-N 0.000 description 5
- 230000001965 increasing effect Effects 0.000 description 4
- 210000002700 urine Anatomy 0.000 description 4
- 230000003321 amplification Effects 0.000 description 3
- 210000004556 brain Anatomy 0.000 description 3
- 230000001149 cognitive effect Effects 0.000 description 3
- 230000007547 defect Effects 0.000 description 3
- 238000003199 nucleic acid amplification method Methods 0.000 description 3
- 238000010521 absorption reaction Methods 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 238000004497 NIR spectroscopy Methods 0.000 description 1
- 208000005374 Poisoning Diseases 0.000 description 1
- 230000004913 activation Effects 0.000 description 1
- 230000003044 adaptive effect Effects 0.000 description 1
- 238000013473 artificial intelligence Methods 0.000 description 1
- 230000003190 augmentative effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000009535 clinical urine test Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- PCHJSUWPFVWCPO-UHFFFAOYSA-N gold Chemical compound [Au] PCHJSUWPFVWCPO-UHFFFAOYSA-N 0.000 description 1
- 238000010801 machine learning Methods 0.000 description 1
- 229960002532 methamphetamine hydrochloride Drugs 0.000 description 1
- TWXDDNPPQUTEOV-FVGYRXGTSA-N methamphetamine hydrochloride Chemical compound Cl.CN[C@@H](C)CC1=CC=CC=C1 TWXDDNPPQUTEOV-FVGYRXGTSA-N 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 210000000653 nervous system Anatomy 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
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/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
- A61B5/316—Modalities, i.e. specific diagnostic methods
- A61B5/369—Electroencephalography [EEG]
- A61B5/377—Electroencephalography [EEG] using evoked responses
- A61B5/378—Visual stimuli
-
- 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/7203—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal
-
- 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
-
- 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/725—Details of waveform analysis using specific filters therefor, e.g. Kalman or adaptive filters
-
- 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
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Artificial Intelligence (AREA)
- Heart & Thoracic Surgery (AREA)
- Molecular Biology (AREA)
- Psychiatry (AREA)
- Veterinary Medicine (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Biomedical Technology (AREA)
- Public Health (AREA)
- Medical Informatics (AREA)
- Signal Processing (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physiology (AREA)
- Psychology (AREA)
- Evolutionary Computation (AREA)
- Fuzzy Systems (AREA)
- Mathematical Physics (AREA)
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Abstract
本发明公开了一种冰毒成瘾者检测模型的训练方法，包括：获取多个待测试者在看到视觉诱发刺激图片后的脑电图信号，对每个待检测者对应的脑电图信号进行预处理，以得到对应于每个待检测者的预处理后的脑电图信号，针对对应于每个待检测者的预处理后的脑电图信号而言，将其中每个对应于标准刺激图片的子脑电图信号的电压值分别与每个对应于成瘾偏差图片的子脑电图信号的电压值相减，从而得到对应于该待检测者的数据集，将所有待检测者的数据集划分为训练集和测试集，将训练集输入冰毒成瘾者检测网络模型中，并使用反向传播算法对冰毒成瘾者检测网络模型中每层的权重参数和偏置参数进行更新和优化，以得到更新后的冰毒成瘾者检测网络模型。
Description
技术领域
本发明属于机器学习和特征检测技术领域，更具体地，涉及一种冰毒成瘾者检测模型的训练方法和系统。
背景技术
冰毒成瘾者检测的金标准是进行尿液检测，主要特点是快速、方便、便于携带和准确率高。但该方法受到毒品半衰期影响，冰毒用药后尿检为阳性的时间段为吸毒后1小时—3天，即只能检测出吸毒后不久的成瘾者。鉴于冰毒成瘾者神经系统受损，有较大的潜在社会危害，因此如何从认知神经科学的角度，成功识别冰毒成瘾者，弥补尿液检测的短时性缺点，已经成为业界的研究热点。
现有从认知神经科学角度的冰毒成瘾者检测方法主要包括使用脑电图(Electroencephalogram，简称EEG)、使用功能核磁共振(Functional magnetic resonanceimaging，简称fMRI)、以及使用功能性近红外光谱((Functional near-infraredspectroscopy，简称fNIRS)。
然而，上述冰毒成瘾者检测方法均存在一些不可忽略的技术问题，对于 EEG方法而言，其虽然简单便携、时间分辨率高，但其空间分辨率较低；对于fMRI方法而言，其能直观地反映脑区激活程度，但设备昂贵；对于 fNIRS方法而言，其时空分辨率和设备价格处于EEG方法和fMRI方法之间，但待检测者的头发会影响检测的准确性。
发明内容
针对现有技术的以上缺陷或改进需求，本发明提供了一种冰毒成瘾者检测模型的训练方法和系统，其目的在于，弥补尿液检测的短时性缺点，解决现有基于EEG实现的冰毒成瘾者检测方法存在的空间分辨率较低的技术问题，以及基于fMRI的冰毒成瘾者检测方法存在的设备昂贵的技术问题，以及基于fNIRS的冰毒成瘾者检测方法存在的待检测者的头发会影响检测准确性的技术问题。
为实现上述目的，按照本发明的一个方面，提供了一种冰毒成瘾者检测模型的训练方法，包括以下步骤：
(1)获取多个待测试者在看到视觉诱发刺激图片后的脑电图信号；
(2)对步骤(1)得到的每个待检测者对应的脑电图信号进行预处理，以得到对应于每个待检测者的预处理后的脑电图信号；
(3)针对步骤(2)得到的对应于每个待检测者的预处理后的脑电图信号而言，将其中每个对应于标准刺激图片的子脑电图信号的电压值a1、 a2、…、am分别与每个对应于成瘾偏差图片的子脑电图信号的电压值b1、 b2、…、bn相减，从而得到对应于该待检测者、并由多个电压差值构成的数据集，即{a1-b1，a1-b2，…，a1-bn，a2-b1，a2-b2，...，a2-bn，..，am-b1， am-b2，…am-bn}，将所有待检测者的数据集按照4:1的比例划分为训练集和测试集，其中n表示对应于成瘾偏差图片的子脑电图信号的总数，m表示对应于标准刺激图片的子脑电图信号的总数。
(4)将步骤(3)得到的训练集输入冰毒成瘾者检测网络模型中，并使用反向传播算法对冰毒成瘾者检测网络模型中每层的权重参数和偏置参数进行更新和优化，以得到更新后的冰毒成瘾者检测网络模型。
(5)对步骤(4)更新后的冰毒成瘾者检测网络模型进行迭代训练，直到该冰毒成瘾者检测网络模型的损失函数达到最小为止；
(6)使用步骤(3)得到的测试集对迭代训练后的冰毒成瘾者检测网络模型进行迭代验证，直到得到的分类精度达到最优为止，从而得到训练好的冰毒成瘾者检测网络模型。
优选地，步骤(1)具体为，首先在Oddball实验模式下，采集待测试者先后看到不同视觉诱发刺激图片后的64通道脑电图信号，然后，在预定的时间间隔后继续重复该过程，从而最终得到最终的脑电图信号。
优选地，视觉诱发刺激图片包括标准刺激图片和成瘾偏差图片，其中标准刺激图片是与冰毒完全无关的图片，成瘾偏差图片是与冰毒高度相关的图片。
优选地，步骤(2)具体为，首先，针对每个待检测者对应的脑电图信号而言，选择TP9、TP10通道的平均电压作为参考电压，然后将其他所有 62路通道的电压分别减去该参考电压，从而得到去参考后的脑电图信号，然后，对该脑电图信号进行1到30Hz的带通滤波，随后，去除带通滤波后的脑电图信号中的伪迹、信号漂移和非生物信号，随后，采用独立成分分析 (ICA)方法去除脑电图信号中的眼电干扰，以得到去除干扰后的脑电图信号；然后，从去除干扰后的脑电图信号中，截取每张标准刺激图片被呈现给待测试者的时刻前200ms、后1000ms的脑电图信号，从而生成对应于每张标准刺激图片的子脑电图信号，并截取每张成瘾偏差图片被呈现给待测试者的时刻前200ms、后1000ms的脑电图信号，从而生成对应于每张成瘾偏差图片的子脑电图信号，所有对应于所有标准刺激图片的子脑电图信号和所有对应于所有成瘾偏差图片的子脑电图信号共同构成对应于该待检测者的预处理后的脑电图信号。
优选地，步骤(4)中权重参数的初始值是使用标准差为0.1的截断式正态分布输出的随机值，偏置参数的初始值设为0；
反向传播算法是采用Adam-Optimizer作为优化器。
优选地，冰毒成瘾者检测网络模型的损失函数L为：
其中N表示训练集的样本总数，K表示训练集中类别的数量，ti，k表示第 k类第i个训练样本输入冰毒成瘾者检测网络模型后的预测结果，yi,k表示第k类第i个训练样本对应的真实结果，k∈[1，K]，λ表示正则化程度，其取值为0.01，Wi,k表示第k类第i个训练样本输入冰毒成瘾者检测网络模型时的权重参数，其随着CNN冰毒成瘾者检测网络模型的训练自行变化。
优选地，冰毒成瘾者检测网络模型采用卷积神经网络，其具体结构如下：
第一层是输入层，其输入大小为300×62的脑电图信号；
第二层是一维卷积层，其卷积核大小为1×100，卷积核个数为64，其输入为输入层的300×62的脑电信号，其输出维度为201×64。
第三层是最大池化层，其步距为2，其输入维度为201×64，其输出维度为100×64。
第四层是卷积层，其卷积核大小为1×32，卷积核个数为128，其输入维度为100×64，其输出维度为69×128。
第五层是最大池化层，其步距为2，其输入维度为69×128，其输出维度为34×128。
第六层是卷积层，其卷积核大小为1×32，卷积核个数为256，其输入维度为34×128，其输出维度为3×256。
第七层是最大池化层，其步距为2，其输入维度为3×256，其输出维度为1×256。
第八层是Dropout层，其比率为0.3，其输出维度为1×256。
第九层是全连接层，其输入维度为1×256，其输出维度为1×2；
按照本发明的另一方面，提供了一种冰毒成瘾者检测模型的训练系统，包括：
第一模块，用于获取多个待测试者在看到视觉诱发刺激图片后的脑电图信号。
第二模块，用于对第一模块得到的每个待检测者对应的脑电图信号进行预处理，以得到对应于每个待检测者的预处理后的脑电图信号；
第三模块，用于针对第二模块得到的对应于每个待检测者的预处理后的脑电图信号而言，将其中每个对应于标准刺激图片的子脑电图信号的电压值a1、a2、…、am分别与每个对应于成瘾偏差图片的子脑电图信号的电压值b1、b2、…、bn相减，从而得到对应于该待检测者、并由多个电压差值构成的数据集，即{a1-b1，a1-b2，…，a1-bn，a2-b1，a2-b2，...，a2-bn，..， am-b1，am-b2，…am-bn}，将所有待检测者的数据集按照4:1的比例划分为训练集和测试集，其中n表示对应于成瘾偏差图片的子脑电图信号的总数， m表示对应于标准刺激图片的子脑电图信号的总数；
第四模块，用于将第三模块得到的训练集输入冰毒成瘾者检测网络模型中，并使用反向传播算法对冰毒成瘾者检测网络模型中每层的权重参数和偏置参数进行更新和优化，以得到更新后的冰毒成瘾者检测网络模型；
第五模块，用于对第四模块更新后的冰毒成瘾者检测网络模型进行迭代训练，直到该冰毒成瘾者检测网络模型的损失函数达到最小为止；
第六模块，用于使用第三模块得到的测试集对迭代训练后的冰毒成瘾者检测网络模型进行迭代验证，直到得到的分类精度达到最优为止，从而得到训练好的冰毒成瘾者检测网络模型。
总体而言，通过本发明所构思的以上技术方案与现有技术相比，能够取得下列有益效果：
1、由于本发明采用了步骤(1)到步骤(6)，其采用脑电信号，从认知科学角度对冰毒成瘾者检测，因此能够解决尿液检测的时间有效性问题；脑电采集设备的成本较低，因此能够解决现有基于fMRI的冰毒成瘾者检测方法存在的设备昂贵的技术问题；脑电采集过程中，脑电信号不受头发等因素的影响，因此能够解决现有基于fNIRS的冰毒成瘾者检测方法存在的待检测者的头发会影响检测准确性的技术问题；
2、由于本发明采用了步骤(1)到步骤(6)，其采用10-20标准的64 通道脑电信号，因此能够解决现有基于EEG实现的冰毒成瘾者检测方法存在的空间分辨率较低的技术问题；
3、由于本发明采用了步骤(1)到步骤(6)，其采用人工智能技术，跨被试的预测平均正确率能够达到86％，因此能够解决现有基于EEG实现的冰毒成瘾者检测方法的识别的可行性问题；
4、由于本发明采用了步骤(3)，其采用样本扩增技术，将样本量扩增了30倍，同时减少了个体脑电的特异性影响，因此能够解决现有深度卷积网络模型样本量不够、特征不够明显，导致的识别效率较低的问题。
附图说明
图1是本发明冰毒成瘾者检测模型的训练方法的流程图；
图2是本发明冰毒成瘾者检测网络模型的网络结构图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。此外，下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。
如图1所示，本发明提供了一种冰毒成瘾者检测模型的训练方法，包括以下步骤：
(1)获取多个待测试者在看到视觉诱发刺激图片后的脑电图信号；
具体而言，本步骤具体为，首先是在Oddball实验模式下，采集待测试者先后看到不同视觉诱发刺激图片(在本示例中是总计85张)后的64通道(10-20国际电极标准系统)脑电图信号(前后两张图片的间隔是2.8到 3.8秒)，然后，间隔某个预定的时间间隔(其可以为5分钟到10分钟之间) 后继续重复该过程，从而最终得到最终的脑电图信号。
其中视觉诱发刺激图片被分成两种类型，即标准刺激图片(其是与冰毒完全无关的图片，在本示例中其是70张篮球图片)、以及成瘾偏差图片(其是与冰毒高度相关的图片，例如吸毒装置、冰毒本身等的图片，在本示例中是从45张预设的成瘾偏差图片集中随机获取15张)。
需要注意的是，本发明呈现给待测试者的85张图片，并不是以固定的顺序进行排列，可以以任何顺序采集以获得最终的多个脑电图信号。
(2)对步骤(1)得到的每个待检测者对应的脑电图信号进行预处理，以得到对应于每个待检测者的预处理后的脑电图信号；
本步骤具体为，首先，针对每个待检测者对应的脑电图信号而言，选择 TP9、TP10通道的平均电压作为参考电压，然后将其他所有62路通道的电压分别减去该参考电压，从而得到去参考后的脑电图信号，然后，对该脑电图信号进行1到30Hz的带通滤波(将采样频率下采样到250Hz)，随后，去除带通滤波后的脑电图信号中的伪迹、信号漂移和非生物信号，随后，采用独立成分分析(independent component analysis，简称ICA)方法去除脑电图信号中的眼电干扰，以得到去除干扰后的脑电图信号；然后，从去除干扰后的脑电图信号中，截取每张标准刺激图片被呈现给待测试者的时刻前 200ms、后1000ms的脑电图信号，从而生成对应于每张标准刺激图片的子脑电图信号，并截取每张成瘾偏差图片被呈现给待测试者的时刻前200ms、后1000ms的脑电图信号，从而生成对应于每张成瘾偏差图片的子脑电图信号，所有对应于所有标准刺激图片的子脑电图信号和所有对应于所有成瘾偏差图片的子脑电图信号共同构成对应于该待检测者的预处理后的脑电图信号。
本步骤的优点在于，流程化去除脑电伪迹，批量提取特定刺激诱发的脑电信号。
(3)针对步骤(2)得到的对应于每个待检测者的预处理后的脑电图信号而言，将其中每个对应于标准刺激图片的子脑电图信号的电压值a1、 a2、…、am分别与每个对应于成瘾偏差图片的子脑电图信号的电压值b1、 b2、…、bn相减，从而得到对应于该待检测者、并由多个电压差值构成的数据集，即{a1-b1，a1-b2，…，a1-bn，a2-b1，a2-b2，...，a2-bn，..，am-b1，am-b2，…am-bn}，其中n表示对应于成瘾偏差图片的子脑电图信号的总数， m表示对应于标准刺激图片的子脑电图信号的总数，共得到353319个数据。将所有待检测者的数据集按照4:1的比例划分为训练集和测试集；
本步骤的优点在于，扩增了用于深度卷积神经网络学习的训练样本，其数量提升了约30倍，因此提高了深度卷积神经网络的泛化能力。同时，标准刺激诱发电位与成瘾偏差刺激诱发电位的差值，一定程度上弱化了不同被试的个体差异，更有效地反映出两类被试间的固有差异，使数据集具有更好的分类特征。
(4)将步骤(3)得到的训练集输入冰毒成瘾者检测网络模型中，并使用反向传播算法对冰毒成瘾者检测网络模型中每层的权重参数和偏置参数进行更新和优化，以得到更新后的冰毒成瘾者检测网络模型；
具体而言，权重参数的初始值是使用标准差为0.1的截断式正态分布输出的随机值，偏置参数的初始值设为0；
具体而言，为了加快损失值的下降速度，本步骤中的反向传播算法是采用Adam-Optimizer作为优化器，这是由于该优化器具有动量以及自适应学习率的优势。
(5)对步骤(4)更新后的冰毒成瘾者检测网络模型进行迭代训练，直到该冰毒成瘾者检测网络模型的损失函数达到最小为止；
冰毒成瘾者检测网络模型的损失函数L为：
其中N表示训练集的样本总数(在本实施方式中，其为282655)，K表示训练集中类别的数量(在本实施方式中，其为2)，ti,k表示第k类第i个训练样本输入冰毒成瘾者检测网络模型后的预测结果，yi,k表示第k类第i 个训练样本对应的真实结果，k∈[1，K]，λ表示正则化程度，其取值为0.01，Wi,k表示第k类第i个训练样本输入冰毒成瘾者检测网络模型时的权重参数，其随着CNN冰毒成瘾者检测网络模型的训练自行变化。
(6)使用步骤(3)得到的测试集对迭代训练后的冰毒成瘾者检测网络模型进行迭代验证，直到得到的分类精度达到最优为止，从而得到训练好的冰毒成瘾者检测网络模型。
如图2所示，本发明中的冰毒成瘾者检测网络模型采用卷积神经网络(Convolutional neural network，简称CNN)，其具体结构如下：
第一层是输入层，其输入大小为300×62的脑电图信号；
第二层是一维卷积层，其卷积核大小为1×100，卷积核个数为64，其输入为输入层的300×62的脑电信号，其输出维度为201×64。
第三层是最大池化层，其步距为2，其输入维度为201×64，其输出维度为100×64。
第四层是卷积层，其卷积核大小为1×32，卷积核个数为128，其输入维度为100×64，其输出维度为69×128。
第五层是最大池化层，其步距为2，其输入维度为69×128，其输出维度为34×128。
第六层是卷积层，其卷积核大小为1×32，卷积核个数为256，其输入维度为34×128，其输出维度为3×256。
第七层是最大池化层，其步距为2，其输入维度为3×256，其输出维度为1×256。
第八层是Dropout层，其比率为0.3，其输出维度为1×256。
第九层是全连接层，其输入维度为1×256，其输出维度为1×2。
实验结果
实验进行了70次标准试次，15次成瘾偏差试次。实验中，采用的刺激为视觉刺激，被试人员在电脑屏幕前观察图片，并迅速做出反应。其中标准刺激图片为70张篮球图片，记为S1，45张冰毒相关图片，记为S2。
本发明旨在识别数据是来源于成瘾被试还是正常被试，将训练集分为 S1，S2两种。预测集随机选择。由于一次运算具有偶然性，又考虑到运算的速度，本发明选择五重交叉验证，对五种可能训练集和两种可能预测集分别测试，结果如下。
表1预测集随机选择交叉验证结果
表1的结果表明预测结果较好，然而根据实际情况，本发明需要以被试人数划分数据集，将五分之一的被试所有数据作为预测集，剩余被试的数据作为训练集，进行五重交叉验证，循环五次。结果如下：
表2按被试人数划分交叉验证结果
如表2中所示，按被试人数划分预测集的正确率普遍低于随机选择的正确率。这是因为随机选择数据作为预测集，训练集中会包含该被试的其他数据，同一被试的脑电波形是相似的，所以预测的准确率更高。按照被试人数划分数据，只能依靠其他被试的数据建立模型，预测正确率自然偏低。
为了提高冰毒成瘾者检测网络模型的性能，本发明将训练数据进行了扩增。如201号被试中，S1有206个数据，记为a1，a2，a3…a206，S2有 44个数据，记为b1，b2，b3…b44。取a1-b1，a1-b2…a1-b44，a2-b1，a2- b2…a2-b44等插值作为新的输入特征。S1-S2差值过后，数据总数为353319 个，数据格式为mat格式。比三种刺激的总数11864扩大了近三十倍。在扩充后的数据上进行五重交叉验证，结果如表3所示。
表3扩增数据后交叉验证结果
从表3中可以观察到，数据扩增的方法对提高争取率是有效的，模型的准确率进一步得到了提高。今后，可以继续完善数据集，增加被试的样本数，采集更多数据作为训练集。
本发明利用成瘾被试与正常被试在面对刺激时脑电波形不同的特征，完成成瘾与正常被试的脑电数据的二分类问题。在脑电数据处理时，首先需要考虑电脑相关配置的型号与功能。
表4硬件配置
本发明选择的电脑内存大小和显卡型号基本可以实现项目数据运算的要求，运算快速，运行流畅。
表5交叉验证运行时间
表5可知，按被试人数划分训练集与预测集时，五重交叉验证时间最短的为S3刺激下的数据训练，运行一次需要的时间为370.565s。耗时最长的为差值运算S1-S2的数据运算。这是由于S1-S2差值运算下扩充的数据量最大，所以耗时长。但总体看来，模型的运算时间都在可以接受的范围内，该模型是可行的。
本领域的技术人员容易理解，以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。
Claims (8)
1.一种冰毒成瘾者检测模型的训练方法，其特征在于，包括以下步骤：
(1)获取多个待测试者在看到视觉诱发刺激图片后的脑电图信号；
(2)对步骤(1)得到的每个待检测者对应的脑电图信号进行预处理，以得到对应于每个待检测者的预处理后的脑电图信号；
(3)针对步骤(2)得到的对应于每个待检测者的预处理后的脑电图信号而言，将其中每个对应于标准刺激图片的子脑电图信号的电压值a1、a2、…、am分别与每个对应于成瘾偏差图片的子脑电图信号的电压值b1、b2、…、bn相减，从而得到对应于该待检测者、并由多个电压差值构成的数据集，即{a1-b1，a1-b2，…，a1-bn，a2-b1，a2-b2，...，a2-bn，..，am-b1，am-b2，…am-bn}，将所有待检测者的数据集按照4:1的比例划分为训练集和测试集，其中n表示对应于成瘾偏差图片的子脑电图信号的总数，m表示对应于标准刺激图片的子脑电图信号的总数。
(4)将步骤(3)得到的训练集输入冰毒成瘾者检测网络模型中，并使用反向传播算法对冰毒成瘾者检测网络模型中每层的权重参数和偏置参数进行更新和优化，以得到更新后的冰毒成瘾者检测网络模型。
(5)对步骤(4)更新后的冰毒成瘾者检测网络模型进行迭代训练，直到该冰毒成瘾者检测网络模型的损失函数达到最小为止；
(6)使用步骤(3)得到的测试集对迭代训练后的冰毒成瘾者检测网络模型进行迭代验证，直到得到的分类精度达到最优为止，从而得到训练好的冰毒成瘾者检测网络模型。
2.根据权利要求1所述的冰毒成瘾者检测模型的训练方法，其特征在于，步骤(1)具体为，首先在Oddball实验模式下，采集待测试者先后看到不同视觉诱发刺激图片后的64通道脑电图信号，然后，在预定的时间间隔后继续重复该过程，从而最终得到最终的脑电图信号。
3.根据权利要求1或2所述的冰毒成瘾者检测模型的训练方法，其特征在于，视觉诱发刺激图片包括标准刺激图片和成瘾偏差图片，其中标准刺激图片是与冰毒完全无关的图片，成瘾偏差图片是与冰毒高度相关的图片。
4.根据权利要求1至3中任意一项所述的冰毒成瘾者检测模型的训练方法，其特征在于，步骤(2)具体为，首先，针对每个待检测者对应的脑电图信号而言，选择TP9、TP10通道的平均电压作为参考电压，然后将其他所有62路通道的电压分别减去该参考电压，从而得到去参考后的脑电图信号，然后，对该脑电图信号进行1到30Hz的带通滤波，随后，去除带通滤波后的脑电图信号中的伪迹、信号漂移和非生物信号，随后，采用独立成分分析(ICA)方法去除脑电图信号中的眼电干扰，以得到去除干扰后的脑电图信号；然后，从去除干扰后的脑电图信号中，截取每张标准刺激图片被呈现给待测试者的时刻前200ms、后1000ms的脑电图信号，从而生成对应于每张标准刺激图片的子脑电图信号，并截取每张成瘾偏差图片被呈现给待测试者的时刻前200ms、后1000ms的脑电图信号，从而生成对应于每张成瘾偏差图片的子脑电图信号，所有对应于所有标准刺激图片的子脑电图信号和所有对应于所有成瘾偏差图片的子脑电图信号共同构成对应于该待检测者的预处理后的脑电图信号。
5.根据权利要求1至4中任意一项所述的冰毒成瘾者检测模型的训练方法，其特征在于，
步骤(4)中权重参数的初始值是使用标准差为0.1的截断式正态分布输出的随机值，偏置参数的初始值设为0；
反向传播算法是采用Adam-Optimizer作为优化器。
6.根据权利要求1至5中任意一项所述的冰毒成瘾者检测模型的训练方法，其特征在于，冰毒成瘾者检测网络模型的损失函数L为：
其中N表示训练集的样本总数，K表示训练集中类别的数量，t_(i,k)表示第k类第i个训练样本输入冰毒成瘾者检测网络模型后的预测结果，y_(i,k)表示第k类第i个训练样本对应的真实结果，k∈[1，K]，λ表示正则化程度，其取值为0.01，W_(i,k)表示第k类第i个训练样本输入冰毒成瘾者检测网络模型时的权重参数，其随着CNN冰毒成瘾者检测网络模型的训练自行变化。
7.根据权利要求1所述的冰毒成瘾者检测模型的训练方法，其特征在于，冰毒成瘾者检测网络模型采用卷积神经网络，其具体结构如下：
第一层是输入层，其输入大小为300×62的脑电图信号；
第二层是一维卷积层，其卷积核大小为1×100，卷积核个数为64，其输入为输入层的300×62的脑电信号，其输出维度为201×64。
第三层是最大池化层，其步距为2，其输入维度为201×64，其输出维度为100×64。
第四层是卷积层，其卷积核大小为1×32，卷积核个数为128，其输入维度为100×64，其输出维度为69×128。
第五层是最大池化层，其步距为2，其输入维度为69×128，其输出维度为34×128。
第六层是卷积层，其卷积核大小为1×32，卷积核个数为256，其输入维度为34×128，其输出维度为3×256。
第七层是最大池化层，其步距为2，其输入维度为3×256，其输出维度为1×256。
第八层是Dropout层，其比率为0.3，其输出维度为1×256。
第九层是全连接层，其输入维度为1×256，其输出维度为1×2。
8.一种冰毒成瘾者检测模型的训练系统，其特征在于，包括：
第一模块，用于获取多个待测试者在看到视觉诱发刺激图片后的脑电图信号。
第二模块，用于对第一模块得到的每个待检测者对应的脑电图信号进行预处理，以得到对应于每个待检测者的预处理后的脑电图信号；
第三模块，用于针对第二模块得到的对应于每个待检测者的预处理后的脑电图信号而言，将其中每个对应于标准刺激图片的子脑电图信号的电压值a1、a2、…、am分别与每个对应于成瘾偏差图片的子脑电图信号的电压值b1、b2、…、bn相减，从而得到对应于该待检测者、并由多个电压差值构成的数据集，即{a1-b1，a1-b2，…，a1-bn，a2-b1，a2-b2，...，a2-bn，..，am-b1，am-b2，…am-bn}，将所有待检测者的数据集按照4:1的比例划分为训练集和测试集，其中n表示对应于成瘾偏差图片的子脑电图信号的总数，m表示对应于标准刺激图片的子脑电图信号的总数；
第四模块，用于将第三模块得到的训练集输入冰毒成瘾者检测网络模型中，并使用反向传播算法对冰毒成瘾者检测网络模型中每层的权重参数和偏置参数进行更新和优化，以得到更新后的冰毒成瘾者检测网络模型；
第五模块，用于对第四模块更新后的冰毒成瘾者检测网络模型进行迭代训练，直到该冰毒成瘾者检测网络模型的损失函数达到最小为止；
第六模块，用于使用第三模块得到的测试集对迭代训练后的冰毒成瘾者检测网络模型进行迭代验证，直到得到的分类精度达到最优为止，从而得到训练好的冰毒成瘾者检测网络模型。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110925526.0A CN113598794A (zh) | 2021-08-12 | 2021-08-12 | 一种冰毒成瘾者检测模型的训练方法和系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110925526.0A CN113598794A (zh) | 2021-08-12 | 2021-08-12 | 一种冰毒成瘾者检测模型的训练方法和系统 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN113598794A true CN113598794A (zh) | 2021-11-05 |
Family
ID=78308318
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110925526.0A Pending CN113598794A (zh) | 2021-08-12 | 2021-08-12 | 一种冰毒成瘾者检测模型的训练方法和系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113598794A (zh) |
Citations (24)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1745701A (zh) * | 2005-07-21 | 2006-03-15 | 高春平 | 高真实度职业测试方法及其装置 |
| CN100998503A (zh) * | 2006-12-28 | 2007-07-18 | 山东大学 | 自动识别与消除脑电信号中眼动干扰的方法 |
| CN101474070A (zh) * | 2009-01-21 | 2009-07-08 | 电子科技大学 | 一种脑电信号中眼电伪迹的去除方法 |
| CN102697493A (zh) * | 2012-05-03 | 2012-10-03 | 北京工业大学 | 一种快速的脑电信号中眼电伪迹自动识别和去除的方法 |
| CN109002798A (zh) * | 2018-07-19 | 2018-12-14 | 大连理工大学 | 一种基于卷积神经网络的单导视觉诱发电位提取方法 |
| CN109300529A (zh) * | 2018-09-12 | 2019-02-01 | 阿呆科技（北京）有限公司 | 基于人工智能和虚拟现实/增强现实的毒瘾渴求靶向干预矫治系统 |
| CN109325410A (zh) * | 2018-08-16 | 2019-02-12 | 广东工业大学 | 一种基于卷积神经网络的脑电信号eeg特征提取方法 |
| CN110059565A (zh) * | 2019-03-20 | 2019-07-26 | 杭州电子科技大学 | 一种基于改进卷积神经网络的p300脑电信号识别方法 |
| CN110123342A (zh) * | 2019-04-17 | 2019-08-16 | 西北大学 | 一种基于脑电波的网瘾检测方法及系统 |
| CN110222643A (zh) * | 2019-06-06 | 2019-09-10 | 西安交通大学 | 一种基于卷积神经网络的稳态视觉诱发电位信号分类方法 |
| CN110580450A (zh) * | 2019-08-12 | 2019-12-17 | 西安理工大学 | 一种基于卷积神经网络的交通标志识别方法 |
| CN111166327A (zh) * | 2020-01-06 | 2020-05-19 | 天津大学 | 基于单通道脑电信号和卷积神经网络的癫痫诊断装置 |
| CN111358475A (zh) * | 2020-03-16 | 2020-07-03 | 国网安徽省电力有限公司 | 电力系统高危作业人员精神状态评估方法及系统 |
| CN111460892A (zh) * | 2020-03-02 | 2020-07-28 | 五邑大学 | 一种脑电图模式分类模型的训练方法、分类方法及系统 |
| CN111466909A (zh) * | 2020-04-14 | 2020-07-31 | 清华大学 | 基于脑电特征的目标检测方法和系统 |
| CN111714089A (zh) * | 2020-06-11 | 2020-09-29 | 阿呆科技（北京）有限公司 | 一种基于多刺激短视频事件相关电位的毒瘾评估系统 |
| CN111724896A (zh) * | 2020-06-11 | 2020-09-29 | 阿呆科技（北京）有限公司 | 一种基于多刺激图片或视频erp的毒瘾评估系统 |
| CN112089398A (zh) * | 2020-08-17 | 2020-12-18 | 上海大学 | 一种毒瘾程度检测方法 |
| CN112381008A (zh) * | 2020-11-17 | 2021-02-19 | 天津大学 | 一种基于并行序列通道映射网络的脑电情感识别方法 |
| CN112528819A (zh) * | 2020-12-05 | 2021-03-19 | 西安电子科技大学 | 基于卷积神经网络的p300脑电信号分类方法 |
| CN112667080A (zh) * | 2020-12-28 | 2021-04-16 | 西安电子科技大学 | 基于深度卷积对抗网络的脑电信号无人平台智能控制方法 |
| CN112932486A (zh) * | 2021-01-20 | 2021-06-11 | 安徽建筑大学 | 基于多模态的大学生情绪压力检测系统与方法 |
| CN112957014A (zh) * | 2021-02-07 | 2021-06-15 | 广州大学 | 一种基于脑电波与神经网络的疼痛检测定位方法及系统 |
| CN113040772A (zh) * | 2021-03-10 | 2021-06-29 | 效隆神思(厦门)科技发展有限公司 | 一种基于脑电信号的冰毒成瘾者和正常人的分类方法 |
-
2021
- 2021-08-12 CN CN202110925526.0A patent/CN113598794A/zh active Pending
Patent Citations (24)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1745701A (zh) * | 2005-07-21 | 2006-03-15 | 高春平 | 高真实度职业测试方法及其装置 |
| CN100998503A (zh) * | 2006-12-28 | 2007-07-18 | 山东大学 | 自动识别与消除脑电信号中眼动干扰的方法 |
| CN101474070A (zh) * | 2009-01-21 | 2009-07-08 | 电子科技大学 | 一种脑电信号中眼电伪迹的去除方法 |
| CN102697493A (zh) * | 2012-05-03 | 2012-10-03 | 北京工业大学 | 一种快速的脑电信号中眼电伪迹自动识别和去除的方法 |
| CN109002798A (zh) * | 2018-07-19 | 2018-12-14 | 大连理工大学 | 一种基于卷积神经网络的单导视觉诱发电位提取方法 |
| CN109325410A (zh) * | 2018-08-16 | 2019-02-12 | 广东工业大学 | 一种基于卷积神经网络的脑电信号eeg特征提取方法 |
| CN109300529A (zh) * | 2018-09-12 | 2019-02-01 | 阿呆科技（北京）有限公司 | 基于人工智能和虚拟现实/增强现实的毒瘾渴求靶向干预矫治系统 |
| CN110059565A (zh) * | 2019-03-20 | 2019-07-26 | 杭州电子科技大学 | 一种基于改进卷积神经网络的p300脑电信号识别方法 |
| CN110123342A (zh) * | 2019-04-17 | 2019-08-16 | 西北大学 | 一种基于脑电波的网瘾检测方法及系统 |
| CN110222643A (zh) * | 2019-06-06 | 2019-09-10 | 西安交通大学 | 一种基于卷积神经网络的稳态视觉诱发电位信号分类方法 |
| CN110580450A (zh) * | 2019-08-12 | 2019-12-17 | 西安理工大学 | 一种基于卷积神经网络的交通标志识别方法 |
| CN111166327A (zh) * | 2020-01-06 | 2020-05-19 | 天津大学 | 基于单通道脑电信号和卷积神经网络的癫痫诊断装置 |
| CN111460892A (zh) * | 2020-03-02 | 2020-07-28 | 五邑大学 | 一种脑电图模式分类模型的训练方法、分类方法及系统 |
| CN111358475A (zh) * | 2020-03-16 | 2020-07-03 | 国网安徽省电力有限公司 | 电力系统高危作业人员精神状态评估方法及系统 |
| CN111466909A (zh) * | 2020-04-14 | 2020-07-31 | 清华大学 | 基于脑电特征的目标检测方法和系统 |
| CN111714089A (zh) * | 2020-06-11 | 2020-09-29 | 阿呆科技（北京）有限公司 | 一种基于多刺激短视频事件相关电位的毒瘾评估系统 |
| CN111724896A (zh) * | 2020-06-11 | 2020-09-29 | 阿呆科技（北京）有限公司 | 一种基于多刺激图片或视频erp的毒瘾评估系统 |
| CN112089398A (zh) * | 2020-08-17 | 2020-12-18 | 上海大学 | 一种毒瘾程度检测方法 |
| CN112381008A (zh) * | 2020-11-17 | 2021-02-19 | 天津大学 | 一种基于并行序列通道映射网络的脑电情感识别方法 |
| CN112528819A (zh) * | 2020-12-05 | 2021-03-19 | 西安电子科技大学 | 基于卷积神经网络的p300脑电信号分类方法 |
| CN112667080A (zh) * | 2020-12-28 | 2021-04-16 | 西安电子科技大学 | 基于深度卷积对抗网络的脑电信号无人平台智能控制方法 |
| CN112932486A (zh) * | 2021-01-20 | 2021-06-11 | 安徽建筑大学 | 基于多模态的大学生情绪压力检测系统与方法 |
| CN112957014A (zh) * | 2021-02-07 | 2021-06-15 | 广州大学 | 一种基于脑电波与神经网络的疼痛检测定位方法及系统 |
| CN113040772A (zh) * | 2021-03-10 | 2021-06-29 | 效隆神思(厦门)科技发展有限公司 | 一种基于脑电信号的冰毒成瘾者和正常人的分类方法 |
Non-Patent Citations (3)
| Title |
|---|
| 毛薇;韩春玉;孙伟;王红星;甄文凤;王玉平;: "人脑对视觉信息变化的处理过程的事件相关电位研究", 中国心理卫生杂志, no. 11, 10 November 2009 (2009-11-10) * |
| 王湘;陈斌;刘鼎;王晓燕;姚树桥;: "不同记忆负荷水平下执行控制的ERP效应", 中国心理卫生杂志, no. 08, 10 August 2008 (2008-08-10) * |
| 邹军;屠嘉衡;苑建齐;吕爽;张丽;董苗淼;沈学军;陈养胜;: "散打运动员赛前与冬训期间视觉事件相关电位P300的比较", 中国运动医学杂志, no. 06, 1 November 2009 (2009-11-01) * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Hu et al. | Real-time sensing of trust in human-machine interactions | |
| Abootalebi et al. | A comparison of methods for ERP assessment in a P300-based GKT | |
| KR101247748B1 (ko) | 사용자 집중도 분석장치 및 방법 | |
| JP7221693B2 (ja) | ヒトの脳の皮質機能の目録を作るための方法および磁気画像化デバイス | |
| CN108309329B (zh) | 一种基于脑网络分析的cnv脑电测谎方法 | |
| CN109255309B (zh) | 面向遥感图像目标检测的脑电与眼动融合方法及装置 | |
| Simbolon et al. | An experiment of lie detection based EEG-P300 classified by SVM algorithm | |
| US10085684B2 (en) | State identification in data with a temporal dimension | |
| KR20150029969A (ko) | 뇌파 측정을 이용한 감성상태 분류 방법 | |
| CN111616702A (zh) | 一种基于认知负荷增强的测谎分析系统 | |
| CN109684296A (zh) | 脑纹识别数据库构建方法、装置、计算机设备和存储介质 | |
| CN111714089B (zh) | 一种基于多刺激短视频事件相关电位的毒瘾评估系统 | |
| CN115439920B (zh) | 基于情绪视听刺激和面部表情的意识状态检测系统、设备 | |
| CN113598794A (zh) | 一种冰毒成瘾者检测模型的训练方法和系统 | |
| CN111081374A (zh) | 基于共同注意力范式的孤独症辅助诊断装置 | |
| CN113768520B (zh) | 一种脑电检测模型的训练方法及装置 | |
| Azhari et al. | Classification of concentration levels in adult-early phase using brainwave signals by applying K-nearest neighbor | |
| US11950919B2 (en) | System for real-time measurement of the activity of a cognitive function and method for calibrating such a system | |
| CN113855022A (zh) | 基于眼动生理信号的情绪评估方法及装置 | |
| CN112603337A (zh) | 一种脑电信号识别方法 | |
| Xu et al. | Unsupervised profiling of microglial arbor morphologies and distribution using a nonparametric Bayesian approach | |
| EP3598940A1 (en) | Method and system for post-traumatic stress disorder (ptsd) and mild traumatic brain injury (mtbi) diagnosis using magnetic resonance spectroscopy | |
| KR20210157229A (ko) | 인공지능을 기반으로 한 심리 검사 장치 및 그 방법 | |
| CN106333675B (zh) | 清醒状态下的脑电信号数据类型的标注方法和系统 | |
| CN116421202B (zh) | 基于脑电快速周期性视觉刺激奇异范式的大脑视觉功能快速检测方法、装置及存储介质 |
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 |