视网膜电图振荡电位的特征值提取方法及其应用技术

本发明专利技术公开了一种视网膜电图振荡电位的特征值提取方法及其应用，包括采集视网膜电图；从视网膜电图中利用滤波器提取振荡电位；将提取的振荡电位在时域上经Hilbert变换后进行包络函数拟合，对包络函数拟合后的波形进行高斯拟合，获得作为时域特征值的时域高斯参数；同时；取振荡电位幅值与视网膜电图中b波幅值之比，获得作为时域特征值的比值参数，根据时域特征值和频域特征值的有效提取和分析，扩大对振荡电位的研究范围，同时在临床应用可通过对比时域特征值对GlyR、GABAc和GABAa受体抑制环路功能进行检测，提高检测效率。

Method for extracting characteristic value of electroretinogram oscillatory potential and its application

The invention discloses a electroretinogram extraction method and its application, including the acquisition of electroretinogram oscillatory potentials of electroretinogram; extracted from the filter will be extracted; oscillatory potentials in time domain by Hilbert transform for envelope function fitting, Gauss fitting of the envelope function after fitting the waveform obtained in time domain Gauss parameters as the time-domain eigenvalue; at the same time; the b wave amplitude of oscillatory potentials than amplitude and ERG, obtained the ratio of parameters as the time-domain eigenvalue, according to the effective extraction and analysis of time-domain eigenvalue and frequency domain characteristic value, expand the scope of the study of oscillatory potentials at the same time, the GlyR, GABAc and GABAa receptor suppression loop function were tested by comparing the time-domain characteristic value in the clinical application, improve the detection efficiency.

【技术实现步骤摘要】
视网膜电图振荡电位的特征值提取方法及其应用
本专利技术涉及一种，具体涉及一种视网膜电图振荡电位的特征值提取方法及其应用。
技术介绍
视网膜电图(electroretinoram，ERG)是光诱导视网膜产生的综合电活动，因其能够客观、定量、无创地反映视网膜内各层不同神经元的功能，使得ERG的临床研究和应用的日益广泛应用。ERG主要成份包括起源于光感受器的a波、起源于ON型双极细胞的b波以及一组叠加在b波上升支上具有较高频率且较低幅值的节律性小波的振荡电位(oscillatorypotentials,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为高频段高斯函数，低频高斯参数至少包括参数a1，高频高斯参数至少包括参数a2。进一步，在步骤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受体抑制环路的药理学实验数据如下：表1：不同药物注射后参数a值(刺激光强为0dB，幅值单位uV)药物实验组a值对照组a值P值GABAa受体阻断剂138.9711±36.25027545.3304±152.81111.96*10-7GlyR受体阻断剂155.8389±23.12676253.628±87.437740.02GABAc受体阻断剂250.4941±110.823409.8047±28.759770.01(P<0.05说明有统计学差异)表2：不同药物注射后OPs/b-wave比值(刺激光强为0dB，幅值单位uV)药物实验组对照组P值GABAa受体阻断剂0.153288±0.0785640.469871±0.1536850.017GlyR受体阻断剂0.253777±0.0515130.454788±0.0958390.001GABAc受体阻断剂0.604921±0.0101680.457399±0.0215170.01(P<0.05说明有统计学差异)表3：不同药物注射后参数m(刺激光强为0dB)药物实验组/对照组GABAc受体阻断剂1.038143±0.016542**GABAa受体阻断剂1.002234±0.002143GlyR受体阻断剂0.997563±0.003682(**P<0.01；P<0.05)P值为统计学根据显著性检验方法所得，OPs/b-wave比值为取振荡电位幅值与视网膜电图中b波幅值之比，其中GABAa受体阻断剂为Gabazine(CAS104104-50-9)，GABAc受体阻断剂为TPMPA((1,2,5,6-tetrahydropyridin-4-yl)methylphosphinicacid)，GlyR受体阻断剂为士的宁(strychnine)。根据表1和表2可知，GABAa、GABAc和GlyR受体抑制环路受阻会引起a值的降低，GABAc受体抑制环路受阻会引起OPs/b比值升高；根据a值和OPs/b值的提取和分析，扩大对振荡电位的研究范围，同时在临床应用可通过对比时域特征值的提取值与正常值，对GABAc受体抑制环路功能进行检测，提高检测效率。本实施例中，在步骤b中，所述滤波器的阶数大于4阶，滤波带宽为60-300Hz，利于提取特征值提取的准确性。本实施例中，还包括步骤d，步骤d为对步骤b提取的振荡电位进行频域特征值提取，具体包括对步骤b提取的振荡电位进行傅立叶变换并对傅立叶变换后的波形在高频段和低频段进行双高斯函数拟合，获得作为频域特征值的高频高斯参数和低频高斯参数。本实施例中，在步骤d中，所述双高斯函数为i＝1为低频段高斯函数，i＝2为高频段高斯函数，低频高斯参数至少包括参数a1，高频高斯参数至少包括参数a2。关于对GlyR、GABAa和GABAc受体抑制环路的药理学实验数据如下：表4：实验组值与对应对照组值的比值频域特征值GABAa受体阻断剂GABAc受体阻断剂GlyR受体阻断剂a10.48262±0.150314*0.810048±0.1903890.852429±0.139084m10.891001±0.088本文档来自技高网...

【技术保护点】
一种视网膜电图振荡电位的特征值提取方法，其特征在于：包括下列步骤：a.采集视网膜电图；b.从步骤a获得的视网膜电图中利用滤波器提取振荡电位；c.将步骤b提取的振荡电位在时域上经Hilbert变换后进行包络函数拟合，对包络函数拟合后的波形进行高斯拟合，获得作为时域特征值的时域高斯参数；同时；取振荡电位幅值与视网膜电图中b波幅值之比，获得作为时域特征值的比值参数。

【技术特征摘要】
1.一种视网膜电图振荡电位的特征值提取方法，其特征在于：包括下列步骤：a.采集视网膜电图；b.从步骤a获得的视网膜电图中利用滤波器提取振荡电位；c.将步骤b提取的振荡电位在时域上经Hilbert变换后进行包络函数拟合，对包络函数拟合后的波形进行高斯拟合，获得作为时域特征值的时域高斯参数；同时；取振荡电位幅值与视网膜电图中b波幅值之比，获得作为时域特征值的比值参数。2.根据权利要求1所述的视网膜电图振荡电位的特征值提取方法，其特征在于：在步骤c中，所述高斯拟合函数为时域高斯参数至少包括参数a。3.根据权利要求1所述的视网膜电图振荡电位的特征值提取方法，其特征在于：在步骤b中，所述滤波器的阶数大于4阶，滤波带宽为60-300Hz。4.根据权利要求1所述的视网膜电图振荡电位的特征值提取方法，其特征在于：还包括步骤d，所述步骤d为对步骤b提取的...

【专利技术属性】
技术研发人员：戴加满，李世迎，阴正勤，
申请(专利权)人：中国人民解放军第三军医大学第一附属医院，
类型：发明
国别省市：重庆,50