一种大脑暂时性异常态的检测方法、装置和系统与流程
本发明涉及一种脑检测的系统和方法，具体涉及一种大脑暂时性异常态的检测方法、装置和系统。
背景技术：
目前，吸毒管控和毒驾检查等禁毒工作的主要问题在于发现的滞后性。年轻人通常在多次吸毒成瘾后才会被家人发觉，极容易错过成瘾前的最佳干预期。毒驾基本都是在造成恶性交通事故后，在事故起因调查过程中才查出是由于司机吸毒所致，无法事前避免悲剧的发生。因此，应对化解毒品问题带来的风险危害和保障人民群众健康福祉，最为有效的举措是实现吸毒诱发案件的前置性预防或者吸毒行为的及时发现介入。然而，现有的吸毒检测手段还无法满足前置性预防或者及时发现介入的需求。
现有技术中，发检和血检是最精准的技术手段，但是需要精密的检测设备，设备价格昂贵且需要一定的操作专业知识，使其在国内一线的公安基层难以得到广泛的采用，通常是设置在各地司法鉴定中心，由专业人员完成检测。“尿检板”和“唾液板”是通过检测体液(尿液和唾液)中的毒品或其代谢物来实现吸毒检测，操作上对比发检和血检来说相对简便。但由于尿检需要私密性，并且需要女性工作人员检查女性，但是基层单位不能保证女性工作人员总是在岗，使其在路边抽检或者公共场所抽检的可操作性不强。同时，采集体液(尿液和唾液)很容易造成疾病的传播，吸毒人群多数患有各种传染疾病(如肝病、艾滋病、梅毒等)，而一线民警或交警在现场执勤时很难有良好的安全防护措施，因此，严重威胁执勤警察的健康安全。此外，“尿检板”和“唾液板”的使用会产生大量的一次性消耗，并且检测时间也要长达五分钟左右。
功能近红外光谱成像技术(functionalnear-infraredspectroscopy，fnirs)作为一种无创光学脑成像技术，其具有无创、便携、低成本、可长时程连续监测、对适用人群无限制等特点，可以友好的在自然环境下进行脑功能活动的实时监测。fnirs技术主要利用生物组织在近红外波长为650nm～950nm的“光学窗”内不同的光学特性，通过测量入射光强及出射光强实现脑活动的检测。光源发出的光经过头皮、颅骨、脑脊液以及大脑皮层的吸收、散射后，可通过与光源一定距离的探测器检测到，近红外光在组织中传播的轨迹呈“香蕉形”。
传统的近红外光技术假设大脑皮层的血氧状态与局部的神经元活动之间存在联系,因此可以通过测量脑组织血红蛋白浓度变化，来观察大脑局部区域血流量的动态变化，间接获取到大脑的神经活动。近红外光在组织的散射传播过程中，光子被吸收导致光衰减，fnirs利用这一特性计算出血氧变化情况。
然而，传统的fnirs技术需要基于大脑皮层的血氧状态或者浓度的变化以得出大脑各部分之间的神经元活动联系(功能连接量)的检测方法需要相对较长的观测信号时间，同时由于需要长时间精细化采集大脑各局部特征，易受外界输入或噪声的影响，效率低且检测成本高。不利于在公路、站卡等场景的快速准确排查。
技术实现要素：
为解决上述至少一个问题，本发明提供的一种大脑暂时性异常态的检测方法、装置和系统。
本发明提供一种大脑暂时性异常态的检测方法，该方法包括以下步骤：
s100：基于出射光强信号获取受试者大脑动态系统的观测空间i；
s200：基于所述观测空间i重构所述大脑动态系统的相空间，得到本征分量构造；
s300：基于所述本征分量构造识别所述大脑动态系统的演化轨迹；
s400：判定所述大脑动态系统的活动态模式；
s500：基于所述判定结果确定所述受试者是否异常并输出结果。
进一步的，其中，所述步骤s100包括，
s110：采集短时高维信号；
s120：对所述短时高维信号预处理；
s130：利用预处理后的所述短时高维信号组成大脑动态系统的所述观测空间i；
其中，所述短时高维信号构成为多路检测通道同时采集的短时出射光强信号。
进一步的，所述重构实施为对所述大脑动态系统的观测空间i做奇异值分解，即：
i＝uσv*
其中，所述大脑动态系统的动力学相空间mv(t)可以通过右奇异矩阵v的前r个本征分量构造。
进一步的，所述演化轨迹基于线性解析分析得出，即
其中a和b分别是基本分量v(t)和残余力r(t)的线性动力学的回归系数。
进一步的，其中，所述步骤s400包括，
步骤s410：基于所述基本分量v(t)和残余力r(t)获得特征量；
步骤s420：根据所述特征量和分类模型确定大脑动态系统的活动态模式。
进一步的，其中，所述特征量包括，
第一个特征量，为基本分量v(t)；
第二个特征量，为残余力r(t)；
第三个特征量，为基本分量v(t)与残余力r(t)之和的希尔伯特变换。
进一步的，所述分类模型基于如下步骤预先构建：
步骤a)：获取大脑正常活动态演化轨迹的正常样本；
步骤b)：获取大脑暂时性异常态演化轨迹的异常样本；
步骤c)：建立大脑动态系统的活动态模式的分类模型，将所述正常样本与异常样本的数据点分成具有最大间距的不同类别。
进一步的，其中，
所述正常样本基于采集未摄入药品人员的出射光强信号，按照所述步骤s100至s300，对每个所述未摄入药品人员分别获得其大脑动态系统的演化轨迹，作为所述正常样本；
所述异常样本基于采集摄入药品人员的出射光强信号，按照所述步骤s100至s300，对每个所述摄入药品人员分别获得其大脑动态系统的演化轨迹，作为所述异常样本。
进一步的，所述预处理包括对所述出射光强信号归一化处理和滤波，其中，
所述归一化处理通过减均值除方差实现，所述滤波采用滤波器滤除低频和高频噪声。
进一步的，所述预处理实施为通过将对应于相同光源的参考信号来拟合每个探测信号，然后用每个探测信号减去相应的参考信号得到的拟合信号。
进一步的，所述观测空间i表示为：
其中，ij表示所述短时高维信号中的第j路出射光强信号，1≤j≤k，所述短时高维信号总共包括k路出射光强信号，m表示所述出射光强信号的总采样点数为m。
其中，所述步骤s400包括，
步骤s410：基于共模空间模式提取大脑相空间的特征获得特征量；
步骤s420：根据所述特征量和分类模型确定大脑动态系统的活动态模式。
本方案通过采用短时高维信号组成大脑动态系统的观测空间，可以避免要求观测信号时间足够长的需求，即缩短信号采集时间，提高检测效率，便于实战使用。此外，本发明的方法通过从大脑动态系统的演化轨迹识别关键全局特征。由于全局特征不同于局部特征，全局特征不易受外界输入或噪声的影响，从而可以实现从复杂的非平稳动态产生的噪声响应中推断出大脑暂时性异常态的关键特征信息。
本发明还提供一种大脑暂时性异常态的检测方法，该方法包括以下步骤：
步骤a100：基于出射光强信号构建受试者大脑第一类相空间；
步骤a200：提取所述受试者大脑第一类相空间特征；
步骤b100：基于所述受试者大脑高维hbo2信号构建所述大脑动态系统的第二相空间；
步骤b200：基于共模空间模式提取所述第二相空间的特征获得第二特征量集；
步骤c：将所述第一特征量和第二特征量进行整合，形成更高维的特征集；
步骤d：至少部分地基于所述特征集判定所述大脑动态系统的活动态模式；
步骤e：基于所述判定结果确定所述受试者是否异常。
其中，所述步骤a100包括：
步骤a110：基于出射光强信号获取受试者大脑动态系统的第一观测空间；
步骤a120：基于所述第一观测空间重构所述大脑动态系统的第一相空间，得到第一本征分量构造。
其中，所述步骤a200包括：
步骤a210：基于所述第一本征分量构造识别所述大脑动态系统的第一演化轨迹，获得第一基本分量v1(t)和第一残余力r1(t)；
步骤a220：基于所述第一基本分量v1(t)和第一残余力r1(t)获得第一特征量。
其中，所述第二相空间构建方法包括：
步骤b110，采集探头中的每个光源辐射的两种波长的光波，通过修正的lambert-beer定律，计算获得hbo2信号；
步骤b120，利用低通滤波器对获得的hbo2信号进行滤波去除噪声；
步骤b130，用滤波处理后的hbo2信号生成大脑相空间矩阵。
本方案同时分别用出射光强信号和hbo2信号构建两个状态空间，通过融合两类特征，可以形成高维的特征集，有助于提高判定的准确性。
本发明还提供一种大脑暂时性异常态检测装置，包括主机和探头，
所述主机包括处理器、存储器、以及通信接口；
所述存储器用于存储至少包括代码、数据和结果；
所述通信接口用于跟所述探头进行通信；
所述探头采集受试者大脑相应区域的出射光强信号并传输给主机；
其中，所述处理器用于执行如前述大脑暂时性异常态的检测方法以鉴定大脑是否处于药品影响下的大脑暂时性异常态，并反馈结果。
进一步的，所述主机还包括显示屏用于显示检测结果。
进一步的，所述探头包括多个光源和多个探测器，其中，所述多个光源可以辐射至少一个以上的光波长，所述探测器的探测光谱范围涵盖所述光源的辐射波长。
进一步的，所述多个探测器可以接收同一个光源发出的光线经大脑的出射光强信号，以形成多个探测信号通道。
进一步的，所述探头还包括陀螺仪和/或加速度器传感器，用于检测探头的姿态变化或者抖动，用于去除所述探测器输出的出射光强信号中的运动伪迹。
本发明还提供一种大脑暂时性异常态检测系统，包括远程数据管理平台和多个如前述大脑暂时性异常态检测装置，所述检测装置包括，
通信电路，用于与数据管理平台无线通信。
进一步的，所述检测装置通过所述通信电路将检测结果实时传输给数据管理平台，或者，
所述检测装置对所述检测结果进行本地存储，定时或在网络通畅时将所述检测结果传输给数据管理平台，
所述数据管理平台包含数据库，将收到的检测结果存储在所述数据库中。
进一步的，所述检测装置还包括身份识别单元，所述身份识别单元至少部分包括：
身份证读取单元，用于支持身份证信息读取并进行与数据管理平台联网查询比对功能，和/或
指纹采集单元格，用于支持指纹采集录入并进行与数据管理平台联网查询比对功能，和/或
摄像头单元，用于支持现场拍照实时人脸识别比对，或者是与数据管理平台联网比对功能。
附图说明
为了更清楚地说明本发明的实施方式或现有技术中的技术方案，下面将对实施方式或现有技术描述中所需要使用的附图作简单地介绍。显而易见地，下面描述中的附图仅仅是示例性的，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据提供的附图引伸获得其它的实施附图。
本说明书所绘示的结构、比例、大小等，均仅用以配合说明书所揭示的内容，以供熟悉此技术的人士了解与阅读，并非用以限定本发明可实施的限定条件，故不具技术上的实质意义，任何结构的修饰、比例关系的改变或大小的调整，在不影响本发明所能产生的功效及所能达成的目的下，均应仍落在本发明所揭示的技术内容得能涵盖的范围内。
图1为多巴胺传输路径示意图；
图2为一种大脑暂时性异常态检测方法的示意图；
图3为获取大脑动态系统的观测空间的示意图；
图4为大脑相空间矩阵示意图；
图5为大脑正常活动态演化轨迹示意图；
图6为大脑暂时性异常态演化轨迹示意图；
图7为判定大脑动态系统的活动态模式示意图；
图8为构建分类模型的方法示意图；
图9为构建的分类模型示意图；
图10为一种大脑暂时性异常态的检测方法示意图；
图11为应用如图10所示方法的结果示例图；
图12为大脑暂时性异常态检测装置的使用方式示意图；
图13为大脑暂时性异常态的检测装置示意图；
图14为探头中光源和光电探测器的一种排布示意图；
图15为探头中光源和光电探测器的另一种排布示意图；
图16为大脑暂时性异常态的检测系统示意图。
具体实施方式
以下由特定的具体实施例说明本发明的实施方式，熟悉此技术的人士可由本说明书所揭露的内容轻易地了解本发明的其他优点及功效，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。需要说明的是，在不冲突的情况下，本发明中的实施例及实施例中的特征可以相互组合。为使本发明的上述目的、特征和优点能够更为明显易懂，下面结合附图对本发明的具体实施例做详细的说明。
大脑动态系统具有混沌、分岔和奇异吸引子等动力学行为，可以用高度非线性的动力学系统来描述。即大脑系统在不接受外界输入情况下靠自身动力学就能维持的非静息的稳定状态(activestationarystate)，为此，需要两个基本条件：
1)神经元之间有兴奋性的互馈连接，没有外界输入的情况下，靠神经元之间的正反馈，网络就能维持住稳定活动。
2)系统同时还有抑制性的连接，避免系统活动因反复的正反馈而“爆炸”。大脑动态系统在正常活动态下，兴奋和抑制作用的比例总体上是平衡的，无外界输入情况下，大脑靠自身动力学维持获得非静息的稳定状态。
在某些对大脑产生影响的药品作用下，上述稳定状态被打破，大脑处于异常兴奋状态并且持续约为十几个小时后，大脑再恢复到正常活动态。大脑在药物影响下的异常兴奋状态被称为“大脑暂时性异常态”。此类兴奋类药物包括，但不限于，苯丙胺类、麦司卡林、可卡因、甲基苯丙胺、苯丙胺类兴奋剂等。如图1所示，人体在摄入兴奋类药物后，在中脑腹侧被盖区(ventraltagmentalarea，vta)产生超大量的多巴胺(dopamine)兴奋性神经递质。多巴胺两条主要传输路径分别是中脑边缘系统通道，即从vta通往伏隔核(nucleusaccumbens)，以及中脑皮层通道，从vta通往前额叶。多巴胺促使前额叶以及伏隔核的神经元兴奋性连接激增，破坏兴奋和抑制比例的平衡，使得大脑的动力学系统脱离非静息稳定态，形成异常的高兴奋活动状态，其效应的平均半衰期约为十几个小时。
基于此，如图2所示，本发明的一种大脑暂时性异常态的检测方法具体如下：
步骤s100：获取大脑动态系统的观测空间
大脑可以被认为是一个动态系统，每一段时刻或者每一个空间位置所探测到的出射光强信号都反映了脑动态系统中的一个特定状态，是脑动态系统的相空间的一个点。通过获取不同空间位置的出射光强信号，可以构建出脑动态系统的观测空间。
如图3所示，获取大脑动态系统的观测空间包括以下步骤：
步骤s110采集短时高维信号
通过包含多个近红外光源110和光电探测器120的探头100，形成多路检测通道。近红外光源110的发射光经过头皮、颅骨、脑脊液以及大脑皮层的吸收和散射后，出射光可通过与光源110一定距离的光电探测器120检测到。出射光的光强变化反映了对应空间部位的大脑活动状态。
在较短的采集时长内(例如10秒内)，多路检测通道同时采集出射光强信号，即获得多路的短时出射光强信号。多路的短时出射光强信号被称为短时高维信号。
步骤s120短时高维信号预处理
预处理可以包括对短时高维信号的归一化和滤波。将出射光强信号ik归一化可通过减均值mean(ik)除方差σ(ik)实现。滤波可采用滤波器滤除低频和高频噪声。例如，采用三阶butterworth滤波器，低通和高通的滤波截止频率分别设定为0.01hz和0.8hz。
近红外光在组织中传播的轨迹呈“香蕉形”，穿透深度相当于光电探测器120与光源110的距离的1/4。因此，通过光电探测器120与光源110的不同距离，可以获得不同深度反射的出射光强信号。当光电探测器120与光源110距离较近时，只能获得脑组织比较浅层的出射光强信号，只有光电探测器120与光源110距离较远时，才能获得脑组织比较深层的出射光强信号。但由于光子一定是先穿过浅层再抵达深层，因此，深层信号中会叠加浅层信号，并且呼吸、心脏活动和低频振荡等全身生理活动对浅层和深层信号都有影响。因此，将浅层信号作为参考信号，将深层信号作为探测信号。例如，对于相同的光源110，离此光源110距离最近的光电探测器120检测的出射光信号作为参考信号，其余接收此光源110出射光的光电探测器120检测的出射光信号作为探测信号。
因此，一些实施例中，通过拟合降噪进行预处理。拟合降噪是通过将对应于相同光源110的参考信号(浅层信号)来拟合每个探测信号(深层信号)，然后用每个探测信号减去对应于相同光源110的参考信号得到的拟合信号。由此，可以消除浅层干扰以及系统级干扰(例如，呼吸，心脏活动和低频振荡等全身生理活动)。
例如，用is表示探测信号，ir表示参考信号，if表示拟合信号。则拟合降噪步骤如下：
步骤s121：通过线性拟合获得if，if(t)＝a×ir(t)+b.其中t表示离散时间指数。可以通过在整个时间间隔内最小化每个时间点的拟合信号if和探测信号is之间的绝对差之和来确定参数a和b，即
步骤s122：拟合降噪后的出射光强信号为探测信号is减去对应于相同光源110的参考信号得到的拟合信号if，即
i(t)＝is(t)-if(t)
在一些实施例中，可以将短时高维信号中的每一路出射光强信号进行低通滤波，以去除生理噪声。例如：呼吸(0.2-0.3hz)，心跳信号(0.8-1.2hz)和血压震荡的mayer波(约0.1hz)所产生的生理噪声。
在一些实施例中，采用三阶切比雪夫(chebyshev)ii型滤波器，截止频率0.1hz，阻带频率0.5hz，通带损耗不超过6db且阻带衰减至少为50db。
在一些实施例中，采用小波滤波器，例如daubechies-12小波进行十级分解。
在另一些实施例中，采用三阶巴特沃斯(butterworth)滤波器，截止频率取为0.8hz和0.01hz。
步骤s130利用预处理后的短时高维信号组成大脑动态系统的观测空间
如图4所示，由预处理后的短时高维信号组成的观测空间i可表示为：
其中，ij代表短时高维信号中的第j路出射光强信号，1≤j≤k，短时高维信号总共包括k路出射光强信号。m代表出射光强信号的总采样点数为m。由此，可以分别生成大脑正常活动态下的，即未受毒品影响的状态，正常相空间矩阵in，和大脑暂时性异常态下，即受毒品影响的暂时性异常态，异常相空间矩阵ia。
步骤s200：重构大脑动态系统的相空间
假设大脑动力学相空间是由具有动力状态的隐藏非线性系统所产生的。存在从观测空间i构建的大脑动态系统的动力学相空间mv(t)能够逼近原始非线性动力状态。
对大脑动态系统的观测空间i做奇异值分解(singularvaluedecomposition,svd)，即下式：
i＝uσv*
大脑动态系统的动力学相空间
步骤s300：识别大脑动态系统的演化轨迹
为了识别大脑动态系统的整体演变的模式，将线性化应用于非线性动力学。采用线性解析分析(linearresolventanalysis，lra)，即
其中a和b分别是基本分量v(t)和残余力r(t)的线性动力学的回归系数。
图5和图6分别显示了没有毒品作用下的大脑正常活动态演化轨迹示例和毒品作用下的大脑异常态演化轨迹示例。可以看出两者存在明显差异，得到的演化轨迹可以分为两个区域，即裙环区域31和中央磁盘区域32。正常活动态的轨迹主要集中在裙环区域31，而暂时性异常态的轨迹主要集中在中央磁盘区域32。
步骤s400：判定大脑动态系统的活动态模式
活动态模式分为正常活动态和暂时性异常态。如图7所示，判定大脑动态系统的活动态模式包含如下步骤：
步骤s410：获得特征量
利用基本分量v(t)和残余力r(t)可以得到，例如三个特征量：第一个特征量为v(t)本身；第二个特征量为r(t)本身；第三个特征量为r(t)与v(t)之和的希尔伯特变换，即h[v(t)+r(t)]。
在一些实施例中，可以采用共模空间模式(commonspatialpatternfiltering,csp)提取大脑相空间的特征。cps的目标是最大化区分大脑正常活动态(对应正常相空间矩阵in)和暂时性异常态(对应异常相空间矩阵ia)。
假定coνn和coνa分别是正常相空间矩阵in和异常相空间矩阵ia的平均协方差矩阵。csp滤波器可以通过rayleigh商来获得，即
其中w表示空间滤波器，t代表转置。
其中nn代表正常类的数据总数，na代表异常类的数据总数。
一系列的空间滤波器w＝[w1,w2,…,wn]，可以通过解决本征值问题获得
coνnw＝λcoνaw
λ代表本征值并按照降序排列。可以选择前三个和后三个特征值，并且相应的特征向量表示两个类的空间滤波器。
估计的源信号可表示为：
y＝wti
计算在检测时间窗内信号的均值，斜率和方差，并将其作为特征。
在另一些实施例中，也可以采用主成分分析(principlecomponentanalysis,pca)，独立成分分析(independentcomponentanalysis,ica)，或者共平均参考(commonaveragereferencing,car)等方法来提取大脑相空间的特征。
步骤s420：根据特征量和分类模型确定大脑动态系统的活动态模式。
分类模型是预先通过大量数据训练的出的，也可以不断将新建测得数据增加到样本库中，更新训练集，优化分类模型。构建分类模型的方法如图8所示：
步骤a)：获取大脑正常活动态演化轨迹的样本。通过采集未摄入毒品人员的出射光强信号，按照上述步骤s100至s300，对每个人员分别获得其大脑动态系统的演化轨迹，作为大脑正常活动态演化轨迹的样本。
步骤b)：获取大脑暂时性异常态演化轨迹的样本。通过采集摄入毒品人员的出射光强信号，按照上述步骤s100至s300，对每个涉毒人员分别获得其大脑动态系统的演化轨迹，作为大脑暂时性异常态演化轨迹的样本。
步骤c)：建立大脑动态系统的活动态模式的分类模型。
通过线性判定分析(lineardiscriminantanalysis,lda)，确定决策超平面，将数据点分成具有最大间距的不同类别。其中，线性判定分析(lda)是对费舍尔的线性鉴别方法的归纳，使用统计学，模式识别和机器学习方法，试图找到两类物体或事件的特征的一个线性组合，以能够特征化或区分它们。所得的组合可用来作为一个线性分类器，或者，为后续的分类做降维处理。在对自变量每一次观察测量值都是连续量的时候，lda能有效的起作用。下面对lda具体计算方法进行描述：
考虑在已知类y中每一个对象或事件的一组观察量
lda的方法是，通过假设条件概率密度函数
lda做出额外简化的方差齐性假设(即，不同类之间协方差是相同的，所以∑0＝∑1＝∑，并且协方差是满秩的。在这种情况下，可以消掉一些项：
如此，上面的判断准则就变成判断一个点积的阈值
对于某些阈值常数c，当
这意味着，一个输入
从几何学的角度来看：判断一个输入
在另一些实施例中，也可以用支持向量机(supportvectormachine,svm)，深度学习(ann)等。
图9为从18人获得异常类实测数据和14人获得正常类实测数据，按照上述方法获得的分类模型。其中粗虚线代表分界面。
其中正常类数据全部落入正常态分类中。15个异常类实测数据落入异常态分类中，有3个异常类实测数据落在两类之外，成为不能正确归类的不确定点。
再次参见图7，根据特征量和分类模型确定大脑动态系统的活动态模式，即为将获得的特征量输入分类模型中，即可根据所在投影空间位置确定当前测试是正常活动态还是暂时性异常态。
s500：基于所述判定结果确定所述受试者是否异常并输出结果。
当判定受试者处于正常活动态时，输出其目前正常的指示；当判定受试者处于暂时性异常态时，输出其目前异常的指示；当判定受试者处于不确定状态时，输出风险提示，从而可对其进行进一步检测，以排除其滥用药物导致精神状态异常的可能。
本发明的方法通过采用短时高维信号组成大脑动态系统的观测空间，可以避免要求观测信号时间足够长的需求，即缩短信号采集时间，提高检测效率，便于实战使用。此外，本发明的方法通过从大脑动态系统的演化轨迹识别关键全局特征。由于全局特征不同于局部特征，全局特征不易受外界输入或噪声的影响，从而可以实现从复杂的非平稳动态产生的噪声响应中推断出大脑暂时性异常态的关键特征信息。
图10为大脑暂时性异常态的检测方法的另一个示例。
如图10所示，其中检测方法如下：
步骤a100：基于出射光强信号构建受试者大脑第一类相空间。具体包括：
步骤a110：基于出射光强信号获取受试者大脑动态系统的第一观测空间；
步骤a120：基于所述第一观测空间重构所述大脑动态系统的第一相空间，得到第一本征分量构造；
步骤a200：提取受试者大脑第一类相空间特征。具体包括：
步骤a210：基于所述第一本征分量构造识别所述大脑动态系统的第一演化轨迹，获得第一基本分量v1(t)和第一残余力r1(t)；
步骤a220：基于所述第一基本分量v1(t)和第一残余力r1(t)获得第一特征量；
步骤b100：基于所述受试者大脑高维hbo2信号构建所述大脑动态系统的第二相空间。其中，hb是指血红蛋白，hbo2是指携带了氧的血红蛋白即氧合血红蛋白。
第二相空间构建方法包括：
步骤b110，采集探头100中的每个光源110辐射的两种波长的光波，通过修正的lambert-beer定律，计算获得hbo2信号。
步骤b120，利用低通滤波器对获得的hbo2信号进行滤波去除噪声；例如：可以采用20阶fir汉明窗滤波器，截止频率为0.1hz。
步骤b130，用滤波处理后的hbo2信号生成大脑相空间矩阵。
步骤b200：基于共模空间模式提取所述第二相空间的特征获得第二特征量集。此步骤类似于步骤s410，此处不再累述。
步骤c：将所述第一特征量和第二特征量进行整合，形成更高维的特征集。
步骤d：至少部分地基于所述特征集判定所述大脑动态系统的活动态模式；分类大脑暂时性异常态和正常态。与步骤s420类似，此处不再累述。
步骤e：基于所述判定结果确定所述受试者是否异常并输出结果。
当判定受试者处于正常活动态时，输出其目前正常的指示；当判定受试者处于暂时性异常态时，输出其目前异常的指示；当判定受试者处于不确定状态时，输出风险提示，从而可对其进行进一步检测，以排除其滥用药物导致精神状态异常的可能。
图11为从40人获得异常类实测数据和50人获得正常类实测数据，按照上述方法，从第一类相空间获取包含三个特征的第一特征量集，从第二类相空间获取包含三个特征的第二特征量级，整合为六维特征集，获得的分类模型。为了显示分类结果，图11选取其中三位特征进行显示，灰色阴影代表代表分界面。其中正常类数据全部落入正常态分类中。39个异常类实测数据落入异常态分类中，有1个异常类实测数据落正常类数据中。
由此可知，图10所示方法与图2所示方法不同之处在于，同时分别用出射光强信号和hbo2信号构建两个状态空间，通过融合两类特征，可以形成高维的特征集，有助于提高判定的准确性。
本发明还提供一种可实施上述暂时性异常态检测方法的大脑暂时性异常态检测装置，其使用方式如图12所示。检测装置可以分为主机200和探头100两部分，测试过程中，可将探头100贴合在被测人的额头或者太阳穴等位置。由图1中的多巴胺传输路径可知，额头对应的前额叶脑区或者太阳穴对应伏隔核脑区均在多巴胺传输路径内。
图13为大脑暂时性异常态的检测装置。主机200包括处理器210、存储器220、以及通信接口230。处理器210用于执行大脑暂时性异常态的检测方法。处理器210可以是单片机、cpu等。存储器220用于存储代码、数据、结果等。通信接口230用于跟探头100的通信接口130进行通信。存储器220可以是闪存、云盘、硬盘等。主机200和探头100之间可以有线连接(如usb、串口等)或者无线连接(如蓝牙、wifi等)。主机200和探头100的通信包括主机200向探头100传输的控制指令(如检测指令)以及探头100向主机200传输数据(如采集到的出射光强信号)等。探头100采集出射光强信号并传输给主机200，主机200执行大脑暂时性异常态的检测方法，鉴定大脑是否处于毒品影响下的大脑暂时性异常态。进一步主机200还可以包括显示装置用于显示检测结果，显示装置可以是触摸感应显示屏、手持智能终端、道路站卡警示屏、可交互头戴设备等。
探头100包括多个光源110和多个光电探测器120。多个光源110可以辐射至少一个以上的光波长，辐射的光波长在700-900nm之间。在一些实施例中，每个光源110可以辐射单个波长的光波，例如840nm。在一些实施例中，每个光源110可以辐射两种波长的光波，可以是一个波长在氧合血红蛋白和脱氧血红蛋白的等吸收点之上(即大于等吸收点的波长)，另一个波长在等吸收点之下(即小于等吸收点的波长)，例如770nm和840nm。光电探测器120的探测光谱范围涵盖光源110的辐射波长。
可以采用光源110复用的方式(即多个光电探测器120可以接收同一个光源110发出的光线经大脑的出射光强信号，以形成多个探测信号通道)，由此，在获取高维信号的同时减小探头100尺寸，有助于提高设备的便携性。
图14为探头100中多个光源110和多个光电探测器120的一种示例性排布。探头100中包括四个光源110(图中由字母s表示)和四个光电探测器120(图中由字母d表示)。光源110与光电探测器120的间距范围可为5mm-30mm之间，光源110间距可以在5-10mm之间。光源110和光电探测器120间的每一条连线代表一个探测信号通道，每一条虚线连线代表一个参考信号通道。一个光源110与周围多个光电探测器120形成多个探测信号通道，与一个距离最近的光电探测器120形成一个对应的参考信号通道。例如，光源s2分别与d1、d2、d4、d5和d6形成s2-d1、s2-d2、s2-d4、s2-d5和s2-d6五个探测信号通道，以及与d3形成一个对应的参考信号通道s2-d3。一个光源与周围多个光电探测器120形成多个探测信号通道，由此，利用四个光源和四个光电探测器120可以形成16个空间维度检测通道。
需要注意的是，图14仅为示例性的，探头100中光源110和光电探测器120的数量以及排布方式不限于此，只要能形成高维探测信号通道即可。例如，可以采用三个光源110和六个光电探测器120，即去掉图14中的任意一列光源110和光电探测器120。再如，也可以采用四个光源110和四个光电探测器120，即去掉图14中第一排光电探测器120。
在一些实施例中，探头100中多个光源110和多个探测器120的排布可以如图15所示，探头100中包括四个光源110(图中由字母s表示)和四个探测器120(图中由字母d表示)。光源110与探测器120的间距范围可为5mm-30mm之间，光源110间距可以在5-10mm之间。光源110和探测器120间的每一条连线代表一个探测信号通道。一个光源110与周围多个探测器120形成多个探测信号通道，由此，利用四个光源110和四个探测器120可以形成10个空间维度检测通道。
在一些实施例中，探头100还可以包括陀螺仪或者加速度器等传感器，用于检测探头100的姿态变化或者抖动，用于去除近红外光电探测器120输出的信号中的运动伪迹。
本发明还公开了一种大脑暂时性异常态检测系统。如图16所示，检测系统包括远程数据管理平台和多个所述大脑暂时性异常态的检测装置。检测装置进一步包含通信电路用于与数据管理平台无线通信。检测装置通过通信电路可以将检测结果实时传输给数据管理平台。也可以是本地存储数据，定时或者检测有网络时将检测结果传输给数据管理平台。数据管理平台包含数据库，将收到的检测结果存储在数据库中。此外，检测装置还可以进一步包括身份证读取单元，用于支持身份证信息读取并进行与数据管理平台联网查询比对功能。检测装置还可以进一步包括指纹采集单元，用于支持指纹采集录入并进行与数据管理平台联网查询比对功能。检测装置还可以进一步包括摄像头单元，用于支持现场拍照实时人脸识别比对，或者是与数据管理平台联网比对功能。
虽然，上文中已经用一般性说明及具体实施例对本发明作了详尽的描述，但在本发明基础上，可以对之作一些修改或改进，这对本领域技术人员而言是显而易见的。因此，在不偏离本发明精神的基础上所做的这些修改或改进，均属于本发明要求保护的范围。
- 还没有人留言评论。精彩留言会获得点赞！