CN106618559B - 一种aEEG信号图像特征重构表示和提取方法 - Google Patents
一种aEEG信号图像特征重构表示和提取方法 Download PDFInfo
- Publication number
- CN106618559B CN106618559B CN201611106785.6A CN201611106785A CN106618559B CN 106618559 B CN106618559 B CN 106618559B CN 201611106785 A CN201611106785 A CN 201611106785A CN 106618559 B CN106618559 B CN 106618559B
- Authority
- CN
- China
- Prior art keywords
- signal
- aeeg
- amplitude
- contour map
- contour
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
- A61B5/316—Modalities, i.e. specific diagnostic methods
- A61B5/369—Electroencephalography [EEG]
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
- A61B5/316—Modalities, i.e. specific diagnostic methods
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7235—Details of waveform analysis
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Surgery (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Veterinary Medicine (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Public Health (AREA)
- Molecular Biology (AREA)
- Physics & Mathematics (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Psychiatry (AREA)
- Psychology (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physiology (AREA)
- Signal Processing (AREA)
- Image Analysis (AREA)
Abstract
本发明公开了一种aEEG信号图像特征重构表示和提取方法,包括如下步骤:信号分段步骤:获取aEEG信号并分成多段信号段;等高线图绘制步骤:统计各所述信号段中不同幅值出现的频率,连接频率相同的点,绘制等高线图;图像特征提取步骤:提取所述等高线图中的图像特征。通过本发明方法绘制等高线图后,aEEG信号的幅值频率等高图中,色泽区分明显,可以直观地显示不同幅值的分布情况,从而帮助医生更好地界定幅值上下边界,在重构的幅值频率等高图基础上,可以进一步提取其图像特征。
Description
技术领域
本发明专利涉及振幅整合脑电图的特征表示和提取方法,特别是一种基于图像信息的aEEG信号特征表示和提取方法。
背景技术
振幅整合脑电图(amplitude-integrated electroencephalography,aEEG)技术是一种新生儿监测脑功能状态的有效手段(Hellstrom–Westas et al.,2008),已广泛作为新生儿脑损伤的辅助检查手段,为医生对新生儿缺氧缺血脑病、窒息、癫痫、颅内出血和脑积水,以及脑病治疗对象筛选提供了客观的评价指标。通常aEEG信号呈条带状轨迹,一般幅值在10~40μV之间。它的宽度与新生儿脑功能状态有关,通常在新生儿深度睡眠期间轨迹较宽,而在轻度睡眠状态或清醒状态下aEEG信号轨迹较窄。
目前国际上对aEEG信号分析主要依靠人工实现,即医生对aEEG的背景图进行目测判读,主要看其背景活动形式(主要关注其上下边界)、睡眠-清醒周期是否存在、痫样波是否存在等。目前对aEEG信号有两种分类方法(De Vries,L.S.et al.,2005)。一类是根据新生儿aEEG的背景活动分为三种:(1)振幅正常:振幅波谱带上边界>10μV,下边界>5μV;(2)振幅轻度异常:波谱带上边界>10μV,下边界<5μV;(3)振幅重度异常:波谱带上边界<10μV,下边界<5μV。三种形式aEEG均可伴有痫样活动,表现为幅值突然增高伴波谱带变窄,随后短暂抑制。根据aEEG背景活动及有无痫样活动,将aEEG结果分为三种:振幅正常为正常aEEG,振幅轻度异常及振幅正常伴痫样活动为轻度异常aEEG,其余均为重度异常aEEG。另一类是将新生儿aEEG的背景活动描述为六种:(1)连续正常电压:连续电活动,幅值为5-10/10-50uV;(2)不连续正常电压:不连续电活动,幅值>5uV;(3)连续低电压:连续电活动,幅值<5uV;(4)爆发抑制:不连续电活动,间歇期幅值极低,间有高幅爆发;(5)平台:<5uV的极低电压。(6)痫样活动。亦根据aEEG背景活动及有无痫样活动,将aEEG结果分为三种:连续正常电压为正常aEEG,连续不正常电压为轻度异常aEEG,其余为重度异常aEEG。
近年来,开始了对新生儿aEEG信号的自动分析研究工作。其基本方法是,首先提取aEEG信号的特征,然后送入分类模型对aEEG信号自动分类判别,其中关键是特征提取这一环节。目前对于aEEG信号提取特征方法主要是针对信号的幅值,提取信号的线性特征包括幅值均值、方差、最大最小值以及上下边界等。如Bowen通过计算aEEG信号的最大最小值、均值和标准差来量化分析aEEG(Bowen,Paradisis et al.2010),其中幅值均值、最大最小幅值给出了aEEG数据的大致分布,方差则是描述aEEG信号幅值的离散程度。此外,也出现了一些非线性特征的使用,例如李雷、陈伟婷等将近似熵理论应用到新生儿aEEG信号分析上,该研究将非线性动力学理论引入aEEG特征提取中(Li,Chen et al.2010);王愈等将线性特征以及幅值分布直方图结合用于aEEG的特征分析(Chen,Wang et al.2014)。
现有的aEEG人工判读中,医生对幅值上下边界根据监护仪上波带颜色浓淡而决定,把颜色较浓处所对应的高低幅值作为上下边界。但背景图的浓淡边界通常不明显,不同人对颜色浓淡的判断不同,导致不同医生对同一aEEG信号的判读结果有异。而在aEEG自动分析中,目前工作基本为基于信号幅值的特征提取,未出现图像特征的使用。
发明内容
本发明提出了一种aEEG信号图像特征重构表示和提取方法,包括如下步骤:
信号分段步骤:获取aEEG信号,对aEEG信号幅值半对数化处理及取整运算,并分成多段信号段;
等高线图绘制步骤:统计各所述信号段中不同幅值出现的频率,连接频率相同的点,绘制各信号段的等高线图,合并得到所述aEEG信号的幅值频率等高图;
图像特征提取步骤:提取所述等高线图中的图像特征。
本发明提出的所述aEEG信号图像特征重构表示和提取方法中,所述等高线图中,统计信号段中不同幅值出现的频率,以直线连接频率相同的点形成等值线,对等值线所围成的等值区域以颜色填充。
本发明提出的所述aEEG信号图像特征重构表示和提取方法中,从所述等高线图中提取的图像特征包括纹理特征和颜色特征。
本发明提出的所述aEEG信号图像特征重构表示和提取方法中,所述信号分段步骤中,所述aEEG信号按时长等分。
本发明提出的所述aEEG信号图像特征重构表示和提取方法中,在所述信号分段步骤之前进一步包括:
信号预处理步骤:剔除aEEG信号中的无用信号;剔除方式包括手动剔除或根据脑功能状态监测仪阻值大小自动剔除过大阻值所对应的aEEG信号。
本发明的有益效果在于:通过本发明方法绘制等高线图后,aEEG信号的幅值频率等高图中,色泽区分明显,可以直观地显示不同幅值的分布情况,从而帮助医生更好地界定幅值上下边界。且根据幅值频率等高图,可以进一步提取其图像特征,包括纹理特征和颜色特征等。图像特征可以更好地提高自动分类效果,其分类结果在正确率、敏感度以及特异性上都有不同程度的提升。
附图说明
图1是本发明aEEG信号图像特征重构表示和提取方法的流程示意图。
图2为正常的新生儿原始aEEG信号图。
图3为异常的新生儿原始aEEG信号图。
图4为正常的新生儿aEEG信号的幅值频率等高图,段长取100秒。
图5为异常的新生儿aEEG信号的幅值频率等高图,段长取100秒。
具体实施方式
结合以下具体实施例和附图,对本发明作进一步的详细说明。实施本发明的过程、条件、实验方法等,除以下专门提及的内容之外,均为本领域的普遍知识和公知常识,本发明没有特别限制内容。
如图1所示,本发明aEEG信号图像特征重构表示和提取方法,包括如下步骤:
信号分段步骤:获取aEEG信号,对aEEG信号幅值半对数化处理及取整运算,并分成多段信号段;
等高线图绘制步骤:统计各所述信号段中不同幅值出现的频率,连接频率相同的点,绘制各信号段的等高线图,合并得到所述aEEG信号的幅值频率等高图;
图像特征提取步骤:提取所述等高线图中的图像特征。
本发明提出的所述aEEG信号图像特征重构表示和提取方法中,所述等高线图中,统计信号段中不同幅值出现的频率,以直线连接频率相同的点形成等值线,对等值线所围成的等值区域以颜色填充。
本发明提出的所述aEEG信号图像特征重构表示和提取方法中,从所述等高线图中提取的图像特征包括纹理特征和颜色特征。
本发明提出的所述aEEG信号图像特征重构表示和提取方法中,所述信号分段步骤中,所述aEEG信号按时长等分。
本发明提出的所述aEEG信号图像特征重构表示和提取方法中,在所述信号分段步骤之前进一步包括:
信号预处理步骤:剔除aEEG信号中的无用信号;剔除方式包括手动剔除或根据脑功能状态监测仪阻值大小自动剔除过大阻值所对应的aEEG信号。
以下结合新生儿aEEG信号为例,具体阐述本发明特征表示和提取方法的实施步骤。
对由脑功能状态监测仪记录的有效新生儿aEEG数据(如图2和图3),经半对数化和取整处理,以固定时长进行信号分段,统计每段信号中各幅值出现的频率,据此画出各段幅值频率等高线图,加以合并后得到整段信号的幅值频率等高图(如图4和图5)。进一步可提取等高图图像特征,如:等高线图灰度化后提取其图像的LBP直方图特征等。以matlab2012b环境下实施为例,步骤如下:
1、数据清洗和筛选,剔除无效信号。
脑功能状态监测仪在监测新生儿脑电时,可能因为电极接触状况以及新生儿自身头部的转动,产生一些电阻过大的无效信号。剔除这些电阻过大的信号,从而保证实验数据的有效性。剔除方式可采用手动剔除,或根据脑功能状态监测仪阻值将过大阻值(如,阻值大于10千欧)所对应的aEEG信号自动剔除。
2、将aEEG信号进行幅值半对数化处理并取整。
对aEEG幅值进行半对数化处理,即幅值小于10μV保持不变,而对幅值大于10μV的进行对数化处理。对半对数化处理后的数据进行取整运算。
3、将信号进行无交叠等长分段。
可按不同时长对信号进行等长分段。以3小时时长的aEEG信号为例。由于该实施例中所用的脑功能状态监测仪的aEEG信号采样率为100Hz,则3小时时长的aEEG包含1080000个数据点。若按100秒每段进行等长无重叠分段,则3小时时长的aEEG可分为108段,每段包含10000个点。通过matlab中tabulate函数统计每段中各个幅值出现的频率。
4、统计各幅值的频率并绘制幅值频率等高图。
采用matlab中contourf函数对各分段aEEG信号画出幅值频率等高图,然后将各分段幅值频率等高图进行合并保存。如图4、图5的幅值频率等高图所示。取不同分段长度时,幅值频率等高图显示形式不同。
5、等高图灰度化。
采用matlab中rgb2gray函数灰度化图像。
6、提取等高图图像特征。
该实施例中通过LBP(local binary pattern)算法对图像提取LBP直方图特征。采用3*3的邻域,共产生58维特征向量,然后送入支持向量机SVM进行分类。
表1所示为支持向量机SVM对幅值特征(最大幅值、最小幅值、均值、方差和下边界)和根据等高图所得LBP图像特征的分类结果。
表1幅值特征和等高图LBP特征的SVM分类结果对比
特征组 | 正确率 | 敏感度 | 特异性 |
幅值特征 | 89.21% | 92.99% | 73.31% |
LBP特征 | 91.12% | 93.09% | 82.87% |
幅值&LBP特征 | 93.64% | 97.63% | 82.07% |
表1中,aEEG信号的LBP图像特征整体的分类能力比通常的幅值特征好,在正确率、敏感度以及特异性上都有不同程度的提升。特别地,以幅值特征和LBP图像特征组合在一起的特征,可使SVM的分类能力有更大的提升,相较于单纯的幅值特征,分类正确率提高超过4个百分点。
本发明的保护内容不局限于以上实施例。在不背离发明构思的精神和范围下,本领域技术人员能够想到的变化和优点都被包括在本发明中,并且以所附的权利要求书为保护范围。
Claims (1)
1.一种aEEG信号图像特征重构表示和提取方法,其特征在于,包括如下步骤:
信号预处理步骤:剔除aEEG信号中的无用信号;剔除方式包括手动剔除或根据脑功能状态监测仪阻值大小自动剔除过大阻值所对应的aEEG信号;
信号分段步骤:获取aEEG信号,对aEEG信号幅值半对数化处理及取整运算,并按时长等分成多段信号段;
等高线图绘制步骤:统计各所述信号段中不同幅值出现的频率,以直线连接频率相同的点形成等值线,对等值线所围成的等值区域以颜色填充,绘制各信号段的等高线图,合并得到所述aEEG信号的幅值频率等高图;
图像特征提取步骤:提取所述等高线图中的图像特征,从所述等高线图中提取的图像特征包括纹理特征和颜色特征。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611106785.6A CN106618559B (zh) | 2016-12-06 | 2016-12-06 | 一种aEEG信号图像特征重构表示和提取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611106785.6A CN106618559B (zh) | 2016-12-06 | 2016-12-06 | 一种aEEG信号图像特征重构表示和提取方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106618559A CN106618559A (zh) | 2017-05-10 |
CN106618559B true CN106618559B (zh) | 2020-08-25 |
Family
ID=58819762
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611106785.6A Active CN106618559B (zh) | 2016-12-06 | 2016-12-06 | 一种aEEG信号图像特征重构表示和提取方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106618559B (zh) |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108830865B (zh) * | 2018-05-08 | 2021-06-15 | 南京伟思医疗科技股份有限公司 | 一种用于动态脑电图像的稳定上下边界的确定方法 |
CN108703753B (zh) * | 2018-05-08 | 2020-10-27 | 南京伟思医疗科技股份有限公司 | 一种基于动态脑电图的新生儿睡眠觉醒周期检测方法 |
CN109009089B (zh) * | 2018-05-08 | 2021-06-15 | 南京伟思医疗科技股份有限公司 | 一种适用于新生儿的脑电信号爆发抑制检测方法 |
CN109009090B (zh) * | 2018-05-08 | 2021-02-26 | 南京伟思医疗科技股份有限公司 | 一种基于直方图的脑信号图显示方法 |
CN109106365A (zh) * | 2018-09-04 | 2019-01-01 | 杭州航弈生物科技有限责任公司 | 基于脑电信号处理的癫痫发作预警设备 |
CN111513675B (zh) * | 2020-04-07 | 2023-09-22 | 南京伟思医疗科技股份有限公司 | 一种基于振幅整合脑电图的新生儿睡眠觉醒周期检测方法 |
CN112244871A (zh) * | 2020-09-25 | 2021-01-22 | 吉林大学 | 一种基于机器学习的振幅整合脑电图分类识读系统 |
CN112545535B (zh) * | 2020-12-07 | 2023-01-06 | 杭州沃维医疗科技有限公司 | 一种基于振幅整合脑电图的睡眠觉醒周期的分析方法 |
CN113476059B (zh) * | 2021-06-02 | 2023-06-27 | 南京伟思医疗科技股份有限公司 | 一种基于振幅整合脑电图判断左右脑对称性的方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101627909A (zh) * | 2009-05-05 | 2010-01-20 | 复旦大学附属儿科医院 | 一种数字化振幅整合脑功能监护仪 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3581361B1 (ja) * | 2003-02-17 | 2004-10-27 | 株式会社脳機能研究所 | 脳活動測定装置 |
EP2492690A1 (en) * | 2011-02-22 | 2012-08-29 | BIOCRATES Life Sciences AG | Method and use of metabolites for the diagnosis of inflammatory brain injury in preterm born infants |
US11096619B2 (en) * | 2013-07-12 | 2021-08-24 | Innara Health, Inc. | Neural analysis and treatment system |
-
2016
- 2016-12-06 CN CN201611106785.6A patent/CN106618559B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101627909A (zh) * | 2009-05-05 | 2010-01-20 | 复旦大学附属儿科医院 | 一种数字化振幅整合脑功能监护仪 |
Non-Patent Citations (1)
Title |
---|
基于小波神经网络的脑电信号数据压缩与棘波识别研究;郁阿丽等;《中国医疗器械杂志》;19981231;第22卷(第5期);正文第252页,图4-6 * |
Also Published As
Publication number | Publication date |
---|---|
CN106618559A (zh) | 2017-05-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106618559B (zh) | 一种aEEG信号图像特征重构表示和提取方法 | |
Islam et al. | A wavelet-based artifact reduction from scalp EEG for epileptic seizure detection | |
Bairagi | EEG signal analysis for early diagnosis of Alzheimer disease using spectral and wavelet based features | |
KR101400362B1 (ko) | 마취 심도 측정 방법 및 장치 | |
Mirzaei et al. | EEG analysis based on wavelet-spectral entropy for epileptic seizures detection | |
Petersen et al. | Generic single-channel detection of absence seizures | |
Lopes et al. | Automatic electroencephalogram artifact removal using deep convolutional neural networks | |
CN110418604B (zh) | 用于检测电生理诱发电位变化的医疗系统和方法 | |
CN107693004A (zh) | 基于hilbert变换的胎儿心电提取与胎儿心率识别方法 | |
Guarascio et al. | Automatic minimization of ocular artifacts from electroencephalogram: A novel approach by combining Complete EEMD with Adaptive Noise and Renyi's Entropy | |
Stevenson et al. | Artefact detection in neonatal EEG | |
Turnip | Automatic artifacts removal of EEG signals using robust principal component analysis | |
CN113974653A (zh) | 基于约登指数优化棘波检测方法和装置、存储介质和终端 | |
US11944445B2 (en) | Method for detecting elements of interest in electrophysiological signals and detector | |
Zhang et al. | High-speed image analysis reveals chaotic vibratory behaviors of pathological vocal folds | |
De Moraes et al. | Automated Detection of Interictal Spikes in EEG: A literature review | |
Fathima et al. | Wavelet based features for classification of normal, ictal and interictal EEG signals | |
Das et al. | Statistical parameters in the dual tree complex wavelet transform domain for the detection of epilepsy and seizure | |
CN113208629A (zh) | 一种基于eeg信号的阿尔兹海默症筛查方法及系统 | |
CN112022151B (zh) | 一种脑电棘慢波的处理及识别方法 | |
Zammouri et al. | Intelligentocular artifacts removal in a noninvasive singlechannel EEG recording | |
Shete et al. | Comparison of sub-band decomposition and reconstruction of EEG signal by daubechies9 and symlet9 wavelet | |
Deepa et al. | EEG feature extraction and classification of Alzheimer’s disease using support vector machine classifier | |
Mayeli et al. | An automatic ICA-based method for removing artifacts from EEG data acquired during fMRI in real time | |
CN113855047A (zh) | 一种基于异常片段识别的fvep信号增强方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CP02 | Change in the address of a patent holder | ||
CP02 | Change in the address of a patent holder |
Address after: 200241 No. 500, Dongchuan Road, Shanghai, Minhang District Patentee after: EAST CHINA NORMAL University Address before: 200062 No. 3663, Putuo District, Shanghai, Zhongshan North Road Patentee before: EAST CHINA NORMAL University |