CN108830865A - 一种用于动态脑电图像的稳定上下边界的确定方法 - Google Patents
一种用于动态脑电图像的稳定上下边界的确定方法 Download PDFInfo
- Publication number
- CN108830865A CN108830865A CN201810433468.8A CN201810433468A CN108830865A CN 108830865 A CN108830865 A CN 108830865A CN 201810433468 A CN201810433468 A CN 201810433468A CN 108830865 A CN108830865 A CN 108830865A
- Authority
- CN
- China
- Prior art keywords
- boundary
- less
- coboundary
- microvolts
- window
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 25
- 230000006641 stabilisation Effects 0.000 title claims abstract description 14
- 238000011105 stabilization Methods 0.000 title claims abstract description 14
- 230000011218 segmentation Effects 0.000 claims abstract description 27
- 238000001914 filtration Methods 0.000 claims abstract description 7
- 238000005070 sampling Methods 0.000 claims description 6
- 238000012360 testing method Methods 0.000 claims description 6
- 238000012545 processing Methods 0.000 claims description 5
- 241000208340 Araliaceae Species 0.000 claims description 3
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 claims description 3
- 235000003140 Panax quinquefolius Nutrition 0.000 claims description 3
- 235000008434 ginseng Nutrition 0.000 claims description 3
- 210000004556 brain Anatomy 0.000 abstract description 3
- 238000012544 monitoring process Methods 0.000 abstract description 3
- 230000009286 beneficial effect Effects 0.000 abstract description 2
- 238000012216 screening Methods 0.000 abstract description 2
- 230000003925 brain function Effects 0.000 description 2
- 230000007812 deficiency Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000007831 electrophysiology Effects 0.000 description 2
- 238000002001 electrophysiology Methods 0.000 description 2
- 208000027418 Wounds and injury Diseases 0.000 description 1
- 230000002490 cerebral effect Effects 0.000 description 1
- 230000006378 damage Effects 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 208000014674 injury Diseases 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/12—Edge-based segmentation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/70—Denoising; Smoothing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20024—Filtering details
- G06T2207/20032—Median filtering
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Abstract
本发明公开了一种用于动态脑电图像的稳定上下边界的确定方法,其特征是,通过对动态脑电图像的上边界和下边界的变化程度进行分段后,对每一段进行滤波确定稳定的上下边界,最后进行背景模式分类。本发明所达到的有益效果:本方法通过对动态脑电图进行处理,采用逐步分段筛选后,并加入滤波手段的方法,使得处理后的图像,能够充分反馈监护对象的状态,其所进行的处理效率高,且对原始数据的影响低,减弱了动态脑电图产生过程中所附带的随机噪声,对于脑电波监控领域有着很强的意义。
Description
技术领域
本发明涉及一种用于动态脑电图像的稳定上下边界的确定方法,属于脑电图图像处理技术领域。
背景技术
振幅整合EEG(aEEG)是临床上脑功能监护和评估的重要方法,是一种广泛使用的脑电趋势图。通过aEEG上下边界的确定来对aEEG进行背景分类是临床上确定脑功能状态和损伤严重度的最重要的方法之一。目前临床上主要通过有经验的电生理医生人工读图来进行分类。本专利提出的自动上下边界确定方法在临床上解决专业电生理医生不足和规模使用脑功能监护上具有重大意义。图一是人工读图的示意。
发明内容
为解决现有技术的不足,本发明的目的在于提供一种用于动态脑电图像的稳定上下边界的确定方法,通过对得到的脑电图经过分段筛选后进行滤波处理,使得脑电图像的上下边界的随机噪声降低,形成能够良好反馈监护对象状态的脑电图。
为了实现上述目标,本发明采用如下的技术方案:
一种用于动态脑电图像的稳定上下边界的确定方法,其特征是,包括如下步骤:
1)采集动态脑电图图像,基于该图像的下边界值进行分段;
2)对于分段结果进行筛选,时长小于t1的段落将不确定上边界和下边界,且不做背景模式分类,时长不小于t1的段落进行下一步骤;
3)根据步骤2)的分段筛选结果,在每一段内,对段落对应的原始图像的上边界和下边界进行中值滤波,得到滤波后的上边界和下边界;
4)在每一段内将滤波后的上边界和下边界每t2时长分成一个窗口,取每个窗口内上边界和下边界的平均值作为该窗口的上边界和下边界,并根据每个窗口的上边界和下边界进行背景模式分类;
进行窗口分割时:若最后的窗口小于则将其与之前的t2时长的窗口合并处理;若最后的窗口大于分钟但小于t2分钟,则作为独立窗口处理。
进一步地,所述步骤1)中分段的具体步骤如下:
101)取开始n点的动态脑电图图像的下边界值作为参考边界集合,取集合的中值作为参考边界;
102)从第(n+1)点开始,比较当前动态脑电图图像的下边界值与参考边界的差别,并进行分段;
103)不断重复101)和102)的处理直至分段结束。
进一步地,所述步骤102)中分段准则是:
I)若两者相差小于或等于a微伏,则将当前下边界值加入参考边界集合;
II)若新的参考边界集合元素超过b个,则去除最前面的那个元素;然后取新的参考边界集合的中值作为新的参考边界,继续检测后续采样直到新的分段出现或信号记录结束;
III)若两者相差大于a微伏,如果有连续n0个下边界值与参考边界相差都大于a微伏,则确定为新的平稳段的开始并进行分段;如果少于n0点,则作为异常值不做分段,也不更新参考边界,继续检测后续采样直到新的分段出现或信号记录结束。
优选地,所述n=10,a=3,b=20,n0=4。
进一步地,所述步骤4)中的背景模式分类准则为:
1)连续正常电压:下边界大于等于5微伏但小于10微伏,上边界大于等于10微伏;
2)不连续正常电压:下边界小于5微伏,上边界大于等于10微伏;
3)爆发抑制:下边界小于等于2微伏,上边界大于等于25微伏;
4)低电压:下边界小于5微伏,上边界小于10微伏但大于5微伏;
5)Inactive Flat:上下边界都小于5微伏。
优选地,所述步骤2)中t1=15min。
优选地,所述步骤3)中中值滤波采用20点的窗口范围。
优选地,所述步骤4)中t2=30min。
本发明所达到的有益效果:本方法通过对动态脑电图进行处理,采用逐步分段筛选后,并加入滤波手段的方法,使得处理后的图像,能够充分反馈监护对象的状态,其所进行的处理效率高,且对原始数据的影响低,减弱了动态脑电图产生过程中所附带的随机噪声,对于脑电波监控领域有着很强的意义。
附图说明
图1是本发明自动确定上下边界进行aEEG背景分类示例图;
图2是临床上人工分类aEEG背景示例图。
具体实施方式
下面对本发明作进一步描述。以下实施例仅用于更加清楚地说明本发明的技术方案,而不能以此来限制本发明的保护范围。
首先对于本发明所涉及的动态脑电图像上下边界确定方法,其主要是针对目前动态脑电图像中往往伴随很强的随机噪声的现象,为了要找到图像中稳定的上边界和下边界,才使用本方法中的内容来进行处理,结合具体的例子,进行说明。
本方法包括如下步骤:
步骤1)采集动态脑电图图像,基于该图像的下边界值进行分段:分段的具体步骤如下:
101)取开始10点的动态脑电图图像的下边界值作为参考边界集合,取集合的中值作为参考边界;
102)从第11点开始,比较当前动态脑电图图像的下边界值与参考边界的差别,并进行分段,分段准则是:
I)若两者相差小于或等于3微伏,则将当前下边界值加入参考边界集合;
II)若新的参考边界集合元素超过20个,则去除最前面的那个元素;然后取新的参考边界集合的中值作为新的参考边界,继续检测后续采样直到新的分段出现或信号记录结束;
III)若两者相差大于3微伏,如果有连续4个下边界值与参考边界相差都大于3微伏,则确定为新的平稳段的开始并进行分段;如果少于4点,则作为异常值不做分段,也不更新参考边界,继续检测后续采样直到新的分段出现或信号记录结束。
103)不断重复101)和102)的处理直至分段结束。
步骤2)对于分段结果进行筛选,时长小于15分钟的段落将不确定上边界和下边界,且不做背景模式分类,时长不小于15分钟的段落进行下一步骤。
步骤3)根据步骤2)的分段筛选结果,在每一段内,对段落对应的原始图像的上边界和下边界进行20点的中值滤波,得到滤波后的上边界和下边界。
步骤4)在每一段内将滤波后的上边界和下边界每30分钟分成一个窗口,取每个窗口内上边界和下边界的平均值作为该窗口的上边界和下边界,并根据每个窗口的上边界和下边界进行背景模式分类。
进行窗口分割时:若最后的窗口小于15分钟,则将其与之前的30分钟的窗口合并处理;若最后的窗口大于15分钟但小于15分钟,则作为独立窗口处理。
背景模式分类准则为:
1)连续正常电压:下边界大于等于5微伏但小于10微伏,上边界大于等于10微伏;
2)不连续正常电压:下边界小于5微伏,上边界大于等于10微伏;
3)爆发抑制:下边界小于等于2微伏,上边界大于等于25微伏;
4)低电压:下边界小于5微伏,上边界小于10微伏但大于5微伏;
5)Inactive Flat:上下边界都小于5微伏。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明技术原理的前提下,还可以做出若干改进和变形,这些改进和变形也应视为本发明的保护范围。
Claims (8)
1.一种用于动态脑电图像的稳定上下边界的确定方法,其特征是,包括如下步骤:
1)采集动态脑电图图像,基于该图像的下边界值进行分段;
2)对于分段结果进行筛选,时长小于t1的段落将不确定上边界和下边界,且不做背景模式分类,时长不小于t1的段落进行下一步骤;
3)根据步骤2)的分段筛选结果,在每一段内,对段落对应的原始图像的上边界和下边界进行中值滤波,得到滤波后的上边界和下边界;
4)在每一段内将滤波后的上边界和下边界每t2时长分成一个窗口,取每个窗口内上边界和下边界的平均值作为该窗口的上边界和下边界,并根据每个窗口的上边界和下边界进行背景模式分类;
进行窗口分割时:若最后的窗口小于则将其与之前的t2时长的窗口合并处理;若最后的窗口大于分钟但小于t2分钟,则作为独立窗口处理。
2.根据权利要求1所述的一种用于动态脑电图像的稳定上下边界的确定方法,其特征是,所述步骤1)中分段的具体步骤如下:
101)取开始n点的动态脑电图图像的下边界值作为参考边界集合,取集合的中值作为参考边界;
102)从第(n+1)点开始,比较当前动态脑电图图像的下边界值与参考边界的差别,并进行分段;
103)不断重复101)和102)的处理直至分段结束。
3.根据权利要求2所述的一种用于动态脑电图像的稳定上下边界的确定方法,其特征是,所述步骤102)中分段准则是:
I)若两者相差小于或等于a微伏,则将当前下边界值加入参考边界集合;
II)若新的参考边界集合元素超过b个,则去除最前面的那个元素;然后取新的参考边界集合的中值作为新的参考边界,继续检测后续采样直到新的分段出现或信号记录结束;
III)若两者相差大于a微伏,如果有连续n0个下边界值与参考边界相差都大于a微伏,则确定为新的平稳段的开始并进行分段;如果少于n0点,则作为异常值不做分段,也不更新参考边界,继续检测后续采样直到新的分段出现或信号记录结束。
4.根据权利要求3所述的一种用于动态脑电图像的稳定上下边界的确定方法,其特征是,所述n=10,a=3,b=20,n0=4。
5.根据权利要求1所述的一种用于动态脑电图像的稳定上下边界的确定方法,其特征是,所述步骤4)中的背景模式分类准则为:
1)连续正常电压:下边界大于等于5微伏但小于10微伏,上边界大于等于10微伏;
2)不连续正常电压:下边界小于5微伏,上边界大于等于10微伏;
3)爆发抑制:下边界小于等于2微伏,上边界大于等于25微伏;
4)低电压:下边界小于5微伏,上边界小于10微伏但大于5微伏;
5)Inactive Flat:上下边界都小于5微伏。
6.根据权利要求1所述的一种用于动态脑电图像的稳定上下边界的确定方法,其特征是,所述步骤2)中t1=15min。
7.根据权利要求1所述的一种用于动态脑电图像的稳定上下边界的确定方法,其特征是,所述步骤3)中中值滤波采用20点的窗口范围。
8.根据权利要求1所述的一种用于动态脑电图像的稳定上下边界的确定方法,其特征是,所述步骤4)中t2=30min。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810433468.8A CN108830865B (zh) | 2018-05-08 | 2018-05-08 | 一种用于动态脑电图像的稳定上下边界的确定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810433468.8A CN108830865B (zh) | 2018-05-08 | 2018-05-08 | 一种用于动态脑电图像的稳定上下边界的确定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108830865A true CN108830865A (zh) | 2018-11-16 |
CN108830865B CN108830865B (zh) | 2021-06-15 |
Family
ID=64148550
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810433468.8A Active CN108830865B (zh) | 2018-05-08 | 2018-05-08 | 一种用于动态脑电图像的稳定上下边界的确定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108830865B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109300544A (zh) * | 2018-12-10 | 2019-02-01 | 南京伟思医疗科技股份有限公司 | 一种新生儿脑功能标准化七步分析方法 |
CN113476059A (zh) * | 2021-06-02 | 2021-10-08 | 南京伟思医疗科技股份有限公司 | 一种基于振幅整合脑电图判断左右脑对称性的方法 |
Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7299088B1 (en) * | 2002-06-02 | 2007-11-20 | Nitish V Thakor | Apparatus and methods for brain rhythm analysis |
CN103070683A (zh) * | 2012-12-20 | 2013-05-01 | 中山大学 | 基于生物电阻抗的睡眠呼吸模式识别方法及装置 |
CN103228209A (zh) * | 2010-09-21 | 2013-07-31 | 皮层动力学有限公司 | 复合脑功能监视和显示系统 |
CN103429150A (zh) * | 2011-03-11 | 2013-12-04 | 皇家飞利浦有限公司 | 用于监测生理信号的监测装置 |
CN103976740A (zh) * | 2014-05-23 | 2014-08-13 | 清华大学 | 一种面向网络环境的脑电信号身份识别系统及识别方法 |
CN104173046A (zh) * | 2014-09-01 | 2014-12-03 | 清华大学 | 一种色彩标记振幅整合脑电图的提取方法 |
CN104644160A (zh) * | 2015-03-18 | 2015-05-27 | 电子科技大学 | 心电图伪差信号识别方法及装置 |
US20160081616A1 (en) * | 2014-09-23 | 2016-03-24 | Boe Technology Group Co., Ltd. | Apparatus and method for processing electroencephalogram, and sleep monitoring wearable device |
CN105708451A (zh) * | 2016-01-29 | 2016-06-29 | 中山衡思健康科技有限公司 | 一种脑电信号质量实时判定方法 |
CN105962929A (zh) * | 2015-03-13 | 2016-09-28 | 胡渐佳 | 脑波数据分段记录方法以及记录和显示装置 |
CN106618559A (zh) * | 2016-12-06 | 2017-05-10 | 华东师范大学 | 一种aEEG信号图像特征重构表示和提取方法 |
CN106889986A (zh) * | 2017-03-08 | 2017-06-27 | 浙江大学 | 一种基于卡诺模型的脑电信号标定方法 |
CN107788976A (zh) * | 2017-09-22 | 2018-03-13 | 复旦大学 | 基于振幅整合脑电图的睡眠监测系统 |
-
2018
- 2018-05-08 CN CN201810433468.8A patent/CN108830865B/zh active Active
Patent Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7299088B1 (en) * | 2002-06-02 | 2007-11-20 | Nitish V Thakor | Apparatus and methods for brain rhythm analysis |
CN103228209A (zh) * | 2010-09-21 | 2013-07-31 | 皮层动力学有限公司 | 复合脑功能监视和显示系统 |
CN103429150A (zh) * | 2011-03-11 | 2013-12-04 | 皇家飞利浦有限公司 | 用于监测生理信号的监测装置 |
CN103070683A (zh) * | 2012-12-20 | 2013-05-01 | 中山大学 | 基于生物电阻抗的睡眠呼吸模式识别方法及装置 |
CN103976740A (zh) * | 2014-05-23 | 2014-08-13 | 清华大学 | 一种面向网络环境的脑电信号身份识别系统及识别方法 |
CN104173046A (zh) * | 2014-09-01 | 2014-12-03 | 清华大学 | 一种色彩标记振幅整合脑电图的提取方法 |
US20160081616A1 (en) * | 2014-09-23 | 2016-03-24 | Boe Technology Group Co., Ltd. | Apparatus and method for processing electroencephalogram, and sleep monitoring wearable device |
CN105962929A (zh) * | 2015-03-13 | 2016-09-28 | 胡渐佳 | 脑波数据分段记录方法以及记录和显示装置 |
CN104644160A (zh) * | 2015-03-18 | 2015-05-27 | 电子科技大学 | 心电图伪差信号识别方法及装置 |
CN105708451A (zh) * | 2016-01-29 | 2016-06-29 | 中山衡思健康科技有限公司 | 一种脑电信号质量实时判定方法 |
CN106618559A (zh) * | 2016-12-06 | 2017-05-10 | 华东师范大学 | 一种aEEG信号图像特征重构表示和提取方法 |
CN106889986A (zh) * | 2017-03-08 | 2017-06-27 | 浙江大学 | 一种基于卡诺模型的脑电信号标定方法 |
CN107788976A (zh) * | 2017-09-22 | 2018-03-13 | 复旦大学 | 基于振幅整合脑电图的睡眠监测系统 |
Non-Patent Citations (3)
Title |
---|
L G M VAN ROOIJ 等: "Recovery of amplitude integrated electroencephalographic background patterns within 24hours of perinatal asphyxia", 《WWWARCHDISCHILD.COM》 * |
张爱巧: "足月小于胎龄儿危险因素及振幅整合脑电图对其监测意义的研究", 《中国优秀硕士学位论文全文数据库 医药卫生科技辑》 * |
王愈: "新生儿振幅整合脑电自动识别研究", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109300544A (zh) * | 2018-12-10 | 2019-02-01 | 南京伟思医疗科技股份有限公司 | 一种新生儿脑功能标准化七步分析方法 |
CN109300544B (zh) * | 2018-12-10 | 2024-04-26 | 南京伟思医疗科技股份有限公司 | 一种新生儿脑功能标准化七步分析方法 |
CN113476059A (zh) * | 2021-06-02 | 2021-10-08 | 南京伟思医疗科技股份有限公司 | 一种基于振幅整合脑电图判断左右脑对称性的方法 |
CN113476059B (zh) * | 2021-06-02 | 2023-06-27 | 南京伟思医疗科技股份有限公司 | 一种基于振幅整合脑电图判断左右脑对称性的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN108830865B (zh) | 2021-06-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Sharma et al. | QRS complex detection in ECG signals using locally adaptive weighted total variation denoising | |
CN108830865A (zh) | 一种用于动态脑电图像的稳定上下边界的确定方法 | |
CN108511055B (zh) | 基于分类器融合及诊断规则的室性早搏识别系统及方法 | |
US11147502B2 (en) | Signal trimming and false positive reduction of post-segmentation swallowing accelerometry data | |
Agurto et al. | Detection of neovascularization in the optic disc using an AM-FM representation, granulometry, and vessel segmentation | |
Davies | On the noise suppression and image enhancement characteristics of the median, truncated median and mode filters | |
Zhang et al. | The detection of crackles based on mathematical morphology in spectrogram analysis | |
CN109300544A (zh) | 一种新生儿脑功能标准化七步分析方法 | |
CN107230208A (zh) | 一种高斯噪声的图像噪声强度估计方法 | |
CN108109159A (zh) | 一种基于hessian矩阵和区域增长相结合的视网膜血管分割系统 | |
EP4108167A1 (en) | Method and device for generating photoplethysmography signal | |
CN108805826A (zh) | 改善去雾效果的方法 | |
CN108703753B (zh) | 一种基于动态脑电图的新生儿睡眠觉醒周期检测方法 | |
US11179109B1 (en) | Apparatus and method for cardiac signal noise detection and disposition based on physiologic relevance | |
Datta et al. | A new contrast enhancement method of retinal images in diabetic screening system | |
DE102018127776A1 (de) | Detektieren von Atmungsfrequenzen in Audio unter Verwendung eines adaptiven Tiefpassfilters | |
CN108765316A (zh) | 雾气浓度自适应判断方法 | |
Huang et al. | An effective tooth isolation method for bitewing dental X-ray images | |
Jadhav et al. | Detection of exudates for diabetic retinopathy using wavelet transform | |
CN110233682A (zh) | 一种实时门限信号检测方法 | |
Teijeiro et al. | A noise robust QRS delineation method based on path simplification | |
Kopman et al. | An analysis of the TOF-watch algorithm for modifying the displayed train-of-four ratio. | |
Raju Maher et al. | A Decision Support System for Automatic Screening of Non-proliferative Diabetic Retinopathy | |
Riebold et al. | Self-adapting classification system for swallow intention detection in dysphagia therapy | |
TWI696191B (zh) | 降低雜訊影響之qrs波偵測演算法 |
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 |