一种辅助诊断脑瘫的婴儿运动姿态分析系统

本发明属于医学测量，具体涉及一种辅助诊断脑瘫的婴儿运动姿态分析系统。

背景技术：

1、脑瘫(cerebralpalsy)，全称脑性瘫痪，是指婴儿出生前到出生后一个月内脑发育早期，由于多种原因导致的非进行性脑损伤综合征。主要表现为中枢性运动障碍以及姿势异常，还可伴有智力低下、癫痫、感知觉障碍、语言障碍及精神行为异常等，是引起小儿机体运动残疾的主要疾病之一。

2、目前对于婴幼儿脑瘫的早期诊断，主要依靠gmfm粗大运动功能评估表、头颅影像学(如核磁共振、ct或b超的异常)、表面肌电信号semg、步态运动特征参数分析等，结合医生经验，得出脑瘫的诊断结果。采用以上方案，一般在幼儿三岁后确诊脑瘫，延误了脑瘫的早期干预治疗。此外，肌电信号检测信噪比低，干扰排除难；gmfm量表判断精确度不够，且依赖人的主观判断；考虑到辐射等因素，头颅影像学不适宜用于婴幼儿检测。

3、全身运动评估作为鉴别早期婴儿脑性瘫痪等神经发育障碍疾病的重要手段，可以在不接触婴儿对其造成影响的情况下对其健康发展状况提出有一定说服力的诊断结果。但是仅靠观察员来进行评估工作不仅需要长时间的训练模式，还可能因为个别性偏差导致误诊，因此借助视频、电磁跟踪等辅助设备来进全身运动评估具有一定的重要现实意义。

技术实现思路

1、有鉴于此，本发明的目的在于提供一种辅助诊断脑瘫的婴儿运动姿态分析系统。在婴幼儿早期的爬行阶段，基于运动姿态，判断待测婴儿是否患有脑瘫风险，从而更方便地实现对高危婴儿的早期治疗与康复。

2、为达到上述目的，本发明提供了一种辅助诊断脑瘫的婴儿运动姿态分析系统，包括视频处理模块、图像处理模块、数据处理模块、分析模块、判别模块和用户图形界面模块；所述视频处理模块向图像处理模块传送图像，所述图像处理模块向数据处理模块传送数据；所述数据处理模块向分析模块传输数据，所述分析模块向判别模块传送数据，所述数据处理模块、分析模块和判别模块均与用户图形界面模块交互数据；

3、视频处理模块：将婴儿视频每隔10帧拆分为图像，时间间隔为δt；

4、图像处理模块：对视频处理模块拆分的图像进行特征提取，每隔δt时间提取婴儿各特征关节点的平面坐标，并对拆分的图像进行高斯滤波；

5、所述婴儿特征关节点包括：鼻子、脖子、右肩膀、右手肘、右手腕、左肩膀、左手肘、左手腕、右髋、右膝、右脚踝、左髋、左膝、左脚踝、右眼、左眼、右耳、左耳18个关节点；

6、数据处理模块：从时域和频域两个方面对婴儿运动特征信息进行数据处理；

7、在时域上，计算婴儿上一帧的关节点像素坐标与下一帧的关节点像素坐标的距离变化，从而得到婴儿各特征关节点的速度序列与加速度序列；

8、

9、

10、在频域上，计算婴儿运动时的能量谱密度；

11、分析模块：

12、i.用于特征数据的降维处理：对于婴儿18个特征关节点的特征数据，采用pca主成分析法进行降维处理，并对降维后的信息保留量进行评估；

13、ii.相关性分析：从运动特征数据的大小和曲线趋势两个方面进行相关性分析，比较标准视频与样本视频的相关性，用相关系数表示样本视频与标准视频的相关程度；

14、判别模块：设有基于catboost建立的婴儿神经发育判别模型和时间窗阈值决策损失函数，用于辅助婴儿脑瘫诊断。

15、进一步，所述数据处理模块在频域上，计算婴儿运动时的能量谱密度的流程如下：

16、1).设一个运动信号s(t)的能量为e，运动信号的能量由式(3)决定：

17、

18、2).运动信号s(t)傅里叶变换，即频谱密度为s(f)，则由巴塞伐尔定理可知：

19、

20、式中，|s(f)|2为能量谱密度，表示在频率f处宽度为f的频带内的信号能量，或单位频带内的信号能量。

21、进一步，所述分析模块用于降维处理的流程如下：

22、s1.假设婴儿的特征关节点共涉及到p个关节，用x1,x2,…,xp分别表示，每一列表示对应关节每隔12帧运动的速度或加速度；

23、

24、s2.计算标准化样本的协方差矩阵；

25、首先对标准化样本进行协方差计算：

26、

27、得到样本相关系数矩阵为：

28、

29、s3.计算r的特征值和特征向量，得到婴儿肢体的特征值及特征向量；

30、特征值：λ1≥λ2≥…≥λp≥0；

31、特征向量：

32、将婴儿全身运动特征用p个主成分y表示，每个主成分都是p个关节的线性组合；

33、

34、s4.计算各关节部位的主成分y的贡献率及累计贡献率；

35、贡献率是特征值的归一化体现，累计贡献率则是对贡献率的叠加，计算公式如下：

36、

37、

38、进一步，所述相关性分析包括婴儿运动速度和加速度的大小与神经发育的关联分析、灰色关联分析。

39、进一步，所述婴儿运动速度和加速度与神经发育的关联分析的流程如下：

40、采用式(12)所示的距离度量来衡量样本与正常婴儿在运动特征数据上存在的差距，d1i距离越小说明样本视频与标准视频的运动数据值相差较小，表明个体间的差异小；

41、

42、进一步，所述灰色关联分析的流程如下：

43、首先标定母序列和子序列，定义子序列中各个指标与母序列的相关系数：

44、

45、式中，a代表|x0(k)-xi(k)|的两极最小差；b表示|x0(k)-xi(k)|的两极最大差；ρ1为分辨系数，取ρ1＝0.5；

46、计算灰色关联度：

47、

48、灰色相关度越大，曲线越接近，相应序列之间的关联度就越大，反之就越小。

49、进一步，所述判别模块辅助诊断脑瘫的流程如下：

50、①.提取初始层的运动姿态特征，并基于初始层前十个重要关节点进行预测分类，后续骨架特征层则对18个关节点信息进行处理，据此构建5s时间窗序列；

51、②.根据5s时间窗序列中各帧图像所获取特征与神经发育异常的相关性程度，按阈值进行判断划分：当检测出异常特征占比超过50％时，判别为具有潜在神经异常疾病的对象；当该占比低于50％时，判别为低风险对象；

52、③.进一步将异常特征占比划分为75％和25％两条标准，结合gms评估中，对于扭动运动阶段中对观察婴儿中出现的轻度异常所存在的单调性gms模式和混乱性gms模型进行补充评估预测；

53、单调性gms：存在一定异于正常婴儿的特征，但对临床诊断无参考价值，异常特征占比为25％～50％；

54、混乱性gms：在检测时特征不明显，在数周后发展明显异常，异常特征占比为50％～75％。

55、本发明的有益效果：

56、1.本发明一种辅助诊断脑瘫的婴儿运动姿态分析系统设有视频处理模块、图像处理模块、数据处理模块、分析模块、判别模块和用户图形界面模块，可通过婴儿全身运动视频，进行辅助诊断，判断待测婴儿是否患有脑瘫风险，从而更方便地实现对高危婴儿的早期治疗与康复。

57、2.本发明一种辅助诊断脑瘫的婴儿运动姿态分析系统，分析模块基于统计学原理，通过主成分分析、方差分析以及灰色相关度计算，获得了具有一定可靠性的预测分类特征向量，可定量研究婴儿运动特征与神经发育异常的相关性；判别模块基于关键点时间窗序列决策阈值策略和对称决策树算法，解决了典型预测分类算法在脑瘫预测中存在的梯度偏差和预测偏移的问题，提高了预测分类的准确性和泛化性。

58、本发明的其他优点、目标和特征在某种程度上将在随后的说明书中进行阐述，并且在某种程度上，基于对下文的考察研究，对本领域技术人员而言将是显而易见的，或者可以从本发明的实践中得到教导。本发明的目标和其他优点可以通过下面的说明书来实现和获得。