[发明专利]视网膜电图振荡电位的特征值提取方法及其应用

说明书

技术领域

本发明涉及一种，具体涉及一种视网膜电图振荡电位的特征值提取方法及其应用。

背景技术

视网膜电图(electroretinoram，ERG)是光诱导视网膜产生的综合电活动，因其能够客观、定量、无创地反映视网膜内各层不同神经元的功能，使得ERG的临床研究和应用的日益广泛应用。ERG主要成份包括起源于光感受器的a波、起源于ON型双极细胞的b波以及一组叠加在b波上升支上具有较高频率且较低幅值的节律性小波的振荡电位(oscillatory potentials,OPs)(如图1所示，A为视网膜电图ERG，B为振荡电位OPs)；振荡电位一般认为起源于内层视网膜，可能来自无长突细胞或抑制性反馈回路。而对其产生机制尚未完全明了，尽管振荡电位在多种眼部疾病的临床诊断中作用很明显，如青光眼、糖尿病视网膜病变及其他疾病；因OPs的起源及测量方法的不确定性，使得其应用受到一定限制。

因此，为解决以上问题，需要一种视网膜电图振荡电位的特征值提取方法及其应用，通过对振荡电位的不同成分进行特征值有效提取，进而扩大对振荡电位的研究范围，利用特征值提高振荡电位应用范围。

发明内容

有鉴于此，本发明的目的是克服现有技术中的缺陷，提供视网膜电图振荡电位的特征值提取方法，通过对振荡电位的不同成分进行特征值有效提取，进而扩大对振荡电位的研究范围，利用特征值提高振荡电位应用范围。

本发明的视网膜电图振荡电位的特征值提取方法，包括下列步骤：

a.采集视网膜电图；

b.从步骤a获得的视网膜电图中利用滤波器提取振荡电位；

c.将步骤b提取的振荡电位在时域上经Hilbert变换后进行包络函数拟合，对包络函数拟合后的波形进行高斯拟合，获得作为时域特征值的时域高斯参数；同时；取振荡电位幅值与视网膜电图中b波幅值之比，获得作为时域特征值的比值参数。

进一步，在步骤c中，所述高斯拟合函数为时域高斯参数至少包括参数a。

进一步，在步骤b中，所述滤波器的阶数大于4阶，滤波带宽为60-300Hz。

进一步，还包括步骤d，所述步骤d为对步骤b提取的振荡电位进行频域特征值提取，具体包括对步骤b提取的振荡电位进行傅立叶变换并对傅立叶变换后的波形在高频段和低频段进行双高斯函数拟合，获得作为频域特征值的高频高斯参数和低频高斯参数。

进一步，在步骤d中，所述双高斯函数为i＝1为低频段高斯函数，i＝2为高频段高斯函数，低频高斯参数至少包括参数a₁，高频高斯参数至少包括参数a₂。

进一步，在步骤d中，振荡电位在傅立叶变换后且高斯拟合前进行波形平滑处理，平滑指标为15-30。

本发明还公开了一种视网膜电图振荡电位的应用，利用上述的视网膜电图振荡电位的特征值提取方法提取特征值，根据提取的特征值分析GlyR、GABAa或/和GABAc受体抑制环路。

本发明的有益效果是：本发明公开的一种视网膜电图振荡电位的特征值提取方法及其应用，根据时域特征值和频域特征值的有效提取和分析，扩大对振荡电位的研究范围，同时在临床应用可通过对比时频域特征值对GlyR、GABAc和GABAa受体抑制环路功能进行检测，提高检测效率。

附图说明

下面结合附图和实施例对本发明作进一步描述：

图1为视网膜电图a、b波及OPs波的幅值和峰时标定；

图2为OPs波的时域分析流程图；

图3为OPs波的频域分析流程图。

具体实施方式

图1为视网膜电图a、b波及OPs波的幅值和峰时标定，图2为OPs波的时域分析流程图，图3为OPs波的频域分析流程图，如图所示，本实施例中的视网膜电图振荡电位的特征值提取方法；包括下列步骤：

a.采集视网膜电图；

b.从步骤a获得的视网膜电图中利用滤波器提取振荡电位；

c.将步骤b提取的振荡电位在时域上经Hilbert变换后进行包络函数拟合，对包络函数拟合后的波形进行高斯拟合，获得作为时域特征值的时域高斯参数；同时；取振荡电位幅值与视网膜电图中b波幅值之比，获得作为时域特征值的比值参数，其中a为幅值，m为峰时，t为敏感度。

在步骤c中，所述高斯拟合函数为时域高斯参数至少包括参数a。

关于对GlyR、GABAa和GABAc受体抑制环路的药理学实验数据如下：