一种基于瞬时频率稳定性参数的帕金森病实时监测方法技术

本发明专利技术提出了一种基于瞬时频率稳定性参数的帕金森病实时监测方法，涉及帕金森病诊断、监测。目前帕金森病还无法实现实时的判断、监测，而且判别的花费较高。本发明专利技术通过分析三轴加速度计采集到的手指运动信号，计算模信号，利用计算震颤的瞬时频率稳定性，提出了评价帕金森病严重程度的新指标，既能诊断震颤的有无，也能诊断震颤的严重程度，从而对于PD的诊断及治疗效果评估都具有重要意义。

【技术实现步骤摘要】
一种基于瞬时频率稳定性参数的帕金森病实时监测方法
本专利技术涉及帕金森病诊断、监测。
技术介绍
帕金森病(PD)是一种常见的神经系统变性疾病，其临床特征主要是包括静止性震颤、运动迟缓、肌肉僵硬、姿势障碍在内的运动症状以及包括睡眠障碍、自主神经功能障碍、嗅觉减退在内的非运动症状[1]。对于帕金森病的诊断，运动障碍专家的诊断准确率估计为80％，如果该疾病不是运动障碍专家的神经科医生诊断的，则为74％[2]。即使使用UKBrainBank帕金森病的标准，大约有20％的帕金森病患者仍然会被误诊[3]。随着影像学的发展以及成像技术和设备的进步，以MRI、PET和SPECT为代表的影像学技术对帕金森病与帕金森叠加综合症的鉴别诊断效果良好，但影像学手段成本过高，难以推广[4]，同时，也无法胜任对于PD治疗后的症状实时评价。而临床神经生理学信号例如手指震颤可被广泛获得，目前已被尝试作为PD的诊断辅助。例如，利用帕金森病患者的手指震颤频率在4-8Hz这一特征，国外有研究提出了震颤稳定性指数(TSI)和平均谐波功率(MHP)，对于区分PD震颤和原生性震颤(ET)，或判断PD有无发作，有一定的效果[4]。但利用手指震颤进行PD诊断的现有方法仍存在不足。第一，由于手的转动自由度，手指的运动会在不同方向转换。目前众多研究或者仅分析三个方向中幅值变化最大的方向，或者只分析主成分分解后的主方向。单一方向信号分析会造成信息缺失，而三维矢量分析或主成分分解还会带来大的计算复杂度，不利于算法的实时性。第二，手指运动信号具有非平稳性，忽略非平稳性会造成结果的偏差[5]，而常规用时频分析解决非平稳性，又使得计算量加大，不利于实时实现。第三，现有方法只能判断震颤的有或无，而缺乏对于震颤严重程度的诊断。鉴于上述问题，本专利技术提出了一种易于实时实现的PD诊断方法，方法通过分析三轴加速度计采集到的运动信号，既能诊断震颤的有无，也能诊断震颤的严重程度，从而对于PD的诊断及治疗效果评估都具有重要意义。
技术实现思路
专利技术目的本专利技术的目的在于利用加速度计采集到的患者手指震颤数据，实时计算相关指标，基于指标判断出是否发病及和发病程度，为医生的诊断和治疗提供可评估依据。技术方案本专利技术的目的是这样实现的：读取加速度计得到患者手指运动数据。为获得不受转动影响的瞬时频率，计算三轴加速度的模，每次取60s数据，先计算判断信号为震颤的阈值，再对模信号序列每秒钟幅值最大值进行判断，未达阈值则认为未发作，不做处理。超过阈值部分，每10s数据采取减去这10s数据均值的方法达到去直流的目的，然后计算这10s震颤的瞬时频率。对上述去除直流量的信号，定位极小值位置，计算相邻极小值之间的时间间隔，获得瞬时周期T，瞬时频率f＝1/T。对瞬时频率序列求差分，记为Δf。Δf与相应阈值比较，小于阈值可诊断为帕金森震颤。进一步，对于已判断为PD发作的震颤信号，计算60s数据内帕金森引起的震颤持续时间D、模信号强度A，利用程度指标L＝D*A作为震颤的严重程度指标，越大提示越严重。附图说明附图1本专利技术的操作流程图。附图2判断信号为震颤的流程图。附图3对信号求取、判断瞬时频率稳定性的流程图。附图4三轴加速度、主成分分解后的主方向信号、模信号的瞬时频率求解示意图。附图5去除非震颤信号的效果图。具体实施方式为了更了解本专利技术的
技术实现思路
，特配合所附图示说明如下。图1是本专利技术一种利用加速度计诊断、监测帕金森病的方法的总体操作流程图。图2是判断信号为震颤的流程图。图3对信号求取瞬时频率稳定性的流程图。图4是三分立轴、PCA信号、模信号的瞬时频率求解示意图，图4(a)-(e)分别画出了三个分立方向、主成分分解后的主方向和模信号的同步时间序列图(g为重力加速度)，其中三个分立方向和主成分分解后的主方向分别画出正梯度的过零点位置，而模信号画出了极小点的位置。由图中可见，由于手的转动自由度，会导致不同方向间的运动相互转化，其结果就是从单一方向无法。，而模信号体现的运动稳定性要优于另外四个单方向信号。可以证明，模信号的极小点即对应主方向的过零点，因此，我们基于模信号极小点进行瞬时频率的估计。图5是去除非震颤信号的效果图。一种利用加速度计诊断、监测帕金森病的方法，其步骤包括：S1：用加速度计采集患者手指运动信号；S2：计算三轴加速度的模，S3：取60秒数据，计算每秒钟内信号幅值最大值，如第i秒钟内幅值最大值，记为max_i，并计算这些每秒钟幅值最大值的平均值后续每秒钟信号幅值最大值max_i与threshold＝C*mean_max_1s比较，C为设置的某常数。若小于threshold，认为这段数据不是震颤信号，则跳过这段数据。图5(a)是原始模序列，(b)是跳过不达阈值后的模序列；S4：对10s震颤信号，计算其平均值。然后这10s信号减去这个平均值，达到去直流量的效果；S5：对求取、判断去除直流量后的信号瞬时频率稳定性，如图3中的流程图。首先计算模序列极小点位置，可以得到图4(e)示意图；利用相邻极小点的时间间隔求得瞬时周期T，进一步获得瞬时频率f＝1/T。利用后一个频率减去前一个频率，记为Δf。对Δf进行判断，小于阈值，认为是PD震颤；S6：若震颤信号未计算完，重复步骤S4-S5，直至60秒内震颤数据取完；S7：计算60秒数据内PD震颤持续时间D和模信号强度A。模信号A计算流程：先求出PD震颤信号每秒钟幅值最大值，然后求得这些最大值的平均值，即为A；S8：计算两者的积，即本专利技术提出的指标L＝D*A，评价患者患病严重程度，该指标越大，则发作越严重；S9：重复步骤S3-S8，直至所有的采集数据取完。下面以本方法的实际应用，结合附图对本专利技术作进一步说明。参考图1，是本专利技术的操作流程图。参考图2，是判断信号为震颤的流程图。参考图3，是对信号求取、判断瞬时频率稳定性的流程图。参考图4，是三轴加速度、主成分分解后的主方向信号、模信号的频率图。参考图5，是去除非震颤信号的效果图。上述图示结果说明，本专利技术提出的一种利用加速度计诊断、监测帕金森病的方法，可以有效判断是否为帕金森病引起的震颤，及观察患者患病程度情况。虽然本专利技术已以较佳实施例揭露如上，然其并非用以限定本专利技术。本专利技术所属
中具有通常知识者，在不脱离本专利技术的精神和范围内，当可作各种的更动与润饰。因此，本专利技术的保护范围当视权利要求书所界定者为准。参考文献[1]侯莹，刘丽华，江钟立.帕金森病运动症状的评估与康复治疗进展[J].中国康复医学杂志，2018，33(11)：1356-1360.[2]RizzoG，CopettiM，ArcutiS，etal.AccuracyofclinicaldiagnosisofParkinsondisease：Asystematicreviewandmeta-analysis[J].Neurology，2016：WNL.0000000000002350.[3]SelikhovaM，KempsterPA，ReveszT，HoltonJL，LeesAJ.Neuropathologicalfindingsinbenigntremulousparkinsonism.MovDisord2013；28：145-52.[4]DiBiaseL，Britta本文档来自技高网...

【技术保护点】
1.一种基于瞬时频率稳定性参数的帕金森病实时监测方法，其特征在于：1，以三轴加速度的模信号为分析对象；2，利用模信号的幅值作为判断是否为PD震颤初步筛选依据；3，初筛后的模序列，以其极小点为依据估计震颤瞬时频率，并计算瞬时频率稳定性参数；4，对判断为PD震颤的信号段，提出了新的帕金森病PD程度的评价指标。

【技术特征摘要】
1.一种基于瞬时频率稳定性参数的帕金森病实时监测方法，其特征在于：1，以三轴加速度的模信号为分析对象；2，利用模信号的幅值作为判断是否为PD震颤初步筛选依据；3，初筛后的模序列，以其极小点为依据估计震颤瞬时频率，并计算瞬时频率稳定性参数；4，对判断为PD震颤的信号段，提出了新的帕金森病PD程度的评价指标。2.如权利要求1所述的以三轴加速度的模信号为分析对象，其特征在于：以计算公式得到的模序列为分析对象，能有效避免手的转动引起运动在三个方向间转换带来的影响，为获得准确的运动瞬时频率信息提供保障。3.如权利要求1所述的利用模信号的幅值作为判断是否为PD震颤初步筛选依据，其特征在...

【专利技术属性】
技术研发人员：黄晓林，徐东旭，王程，余强，葛云，陈颖，
申请(专利权)人：南京大学，
类型：发明
国别省市：江苏,32