CN114711733A - 一种脉搏信号提取方法、装置、电子设备及存储介质 - Google Patents

一种脉搏信号提取方法、装置、电子设备及存储介质 Download PDF

Info

Publication number
CN114711733A
CN114711733A CN202210632256.9A CN202210632256A CN114711733A CN 114711733 A CN114711733 A CN 114711733A CN 202210632256 A CN202210632256 A CN 202210632256A CN 114711733 A CN114711733 A CN 114711733A
Authority
CN
China
Prior art keywords
pulse
period
signal
pulse signal
filtered
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
Application number
CN202210632256.9A
Other languages
English (en)
Other versions
CN114711733B (zh
Inventor
尚怡辰
王新安
李秋平
陈珙
杨婕
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Peking University Shenzhen Graduate School
Original Assignee
Peking University Shenzhen Graduate School
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Peking University Shenzhen Graduate School filed Critical Peking University Shenzhen Graduate School
Priority to CN202210632256.9A priority Critical patent/CN114711733B/zh
Publication of CN114711733A publication Critical patent/CN114711733A/zh
Application granted granted Critical
Publication of CN114711733B publication Critical patent/CN114711733B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/02Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/48Other medical applications
    • A61B5/4854Diagnosis based on concepts of traditional oriental medicine
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7203Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7225Details of analog processing, e.g. isolation amplifier, gain or sensitivity adjustment, filtering, baseline or drift compensation

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Medical Informatics (AREA)
  • Surgery (AREA)
  • Veterinary Medicine (AREA)
  • Pathology (AREA)
  • Physics & Mathematics (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Biophysics (AREA)
  • Molecular Biology (AREA)
  • Signal Processing (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Physiology (AREA)
  • Artificial Intelligence (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Psychiatry (AREA)
  • Alternative & Traditional Medicine (AREA)
  • Cardiology (AREA)
  • Power Engineering (AREA)
  • Manipulation Of Pulses (AREA)
  • Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)

Abstract

本发明涉及一种脉搏信号提取方法、装置、电子设备及存储介质,涉及脉搏信号技术领域。通过对在变化压力下采集得到的初始脉搏信号进行去噪,通过滤除多种噪声而得到的滤波脉搏信号,进而对滤波脉搏信号进行处理得到的中脉、浮脉以及沉脉,可通过本发明提供的脉搏提取方法为后续的脉搏信号的分析人员减少了对脉搏信号数据的分析处理量,通过将得到的浮脉、中脉、沉脉对脉搏信号更具体精细化处理,进一步提高脉搏信号提取的精确性以及效率性。

Description

一种脉搏信号提取方法、装置、电子设备及存储介质
技术领域
本发明涉及脉搏信号技术领域,尤其涉及一种脉搏信号提取方法、装置、电子设备及存储介质。
背景技术
中医脉诊是通过按触人体不同部位的脉搏,体察脉象变化的切诊方法。中医切脉时,施加不同的手指压力,按浮、中、沉的不同取法感知患者的脉律、深浅、强弱等综合印象,获得诊病开方的所需的重要信息。然而,中医号脉是以手指的感觉来获取脉象,根据医师对于脉搏的主观判断并结合自身经验,将手指感知到的脉形象化,称为“脉象”。通过人为对脉象进行判别,也就是人为对脉搏信号进行提取会造成误差性高的问题;通过现有的技术提取出来的脉搏信号还需要进行大量的数据处理量,不够便捷且人为提取容易造成误差。
发明内容
本发明提供了一种脉搏信号提取方法、装置、电子设备及存储介质,以解决现有技术中依靠人为提取脉搏信号,误差性高的问题。
为了解决上述问题,本发明采用以下的技术方案:
第一方面,本发明提供了一种脉搏信号提取方法,所述脉搏信号提取方法包括:
将初始脉搏信号中的噪声去除得到滤波脉搏信号,所述初始脉搏信号是在变化压力条件下采集到的;
提取所述滤波脉搏信号的包络信号;
根据所述包络信号对所述滤波脉搏信号进行周期分割,得到多个脉搏周期;
从所述包络信号中获取信号最大值,将所述信号最大值对应的脉搏周期作为极值脉搏周期;
将所述极值脉搏周期以及与所述极值脉搏周期相邻的两个脉搏周期作为中脉周期,并将所述中脉周期所对应的滤波脉搏信号作为中脉;
从多个所述脉搏周期中获取三个在所述中脉周期前的脉搏周期作为浮脉周期,并将所述浮脉周期所对应的滤波脉搏信号作为浮脉;
从多个所述脉搏周期中获取三个在所述中脉周期后的脉搏周期作为沉脉周期,并将所述沉脉周期所对应的滤波脉搏信号作为沉脉。
其进一步的技术方案为,将初始脉搏信号中的噪声去除得到滤波脉搏信号,包括:
通过带通滤波器将所述脉搏信号进行滤波,滤除频率为50赫兹的工频干扰、低频基线漂移及肌电干扰的高频噪声,得到所述滤波脉搏信号。
其进一步的技术方案为,所述提取所述滤波脉搏信号的包络信号,包括:
计算所述滤波脉搏信号的斜率;
检测所述滤波脉搏信号的波峰,根据所述斜率以及所述波峰,查找属于所述滤波脉搏信号范围内的波峰,得到波峰序列;
记录所述波峰序列内的波峰的横坐标以及纵坐标;
根据所述横坐标以及所述纵坐标,将所述波峰序列内的波峰进行拟合,得到所述包络信号。
其进一步的技术方案为,所述从所述包络信号中获取信号最大值,包括:
根据所述纵坐标,从所述波峰序列中获取最大值波峰;
将所述最大值波峰作为所述信号最大值。
其进一步的技术方案为,所述根据所述包络信号对所述滤波脉搏信号进行周期分割,得到多个脉搏周期,包括:
将所述波峰序列中的波峰作为周期分割点,所述波峰序列中的波峰分别对应单个脉搏周期;
根据所述周期分割点对所述滤波脉搏信号进行周期分割得到多个所述脉搏周期。
其进一步的技术方案为,所述从多个所述脉搏周期中获取三个在所述中脉周期前的脉搏周期作为浮脉周期,并将所述浮脉周期所对应的滤波脉搏信号作为浮脉,包括:
根据所述最大值波峰的横坐标乘以第一预设系数得到的第一预设结果;
根据所述第一预设结果获取所述浮脉周期对应的波峰;
将所述浮脉周期对应的波峰所对应的滤波脉搏信号作为浮脉。
其进一步的技术方案为,所述从多个所述脉搏周期中获取三个在所述中脉周期后的脉搏周期作为沉脉周期,并将所述沉脉周期所对应的滤波脉搏信号作为沉脉,包括:
根据所述最大值波峰的横坐标乘以第二预设系数得到的第二预设结果;
根据所述第二预设结果获取所述沉脉周期对应的波峰;
将所述沉脉周期对应的波峰所对应的滤波脉搏信号作为沉脉。
第二方面,本发明提供了一种脉搏信号提取装置,包括用于执行如第一方面所述的方法的单元。
第三方面,本发明提供了一种电子设备,所述电子设备包括控制器,所述控制器包括处理器、通信接口、存储器和通信总线,其中,处理器,通信接口,存储器通过通信总线完成相互间的通信;
存储器,用于存放计算机程序;
处理器,用于执行存储器上所存放的程序时,实现如第一方面所述的方法的步骤。
第四方面,本发明提供了一种计算机可读存储介质,其上存储有计算机程序,所述计算机程序被处理器执行时实现如第一方面所述的方法的步骤。
本发明的有益效果:
通过对在变化压力下采集得到的初始脉搏信号进行去噪,通过滤除多种噪声而得到的滤波脉搏信号,进而对滤波脉搏信号进行处理得到的中脉、浮脉以及沉脉,通过对中医脉诊中脉象的具体精细化,进而得到更标准化的信号处理结果,可通过本发明提供的脉搏提取方法为后续的脉搏信号的分析人员减少了对脉搏信号数据的分析处理量,通过将得到的浮脉、中脉、沉脉对脉搏信号更具体精细化处理,进一步提高脉搏信号提取的精确性以及效率性。
附图说明
此处的附图被并入说明书中并构成本说明书的一部分,示出了符合本发明的实施例,并与说明书一起用于解释本发明的原理。
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,对于本领域普通技术人员而言,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例1提供的一种脉搏信号提取方法的流程图;
图2为本发明实施例1提供的一种脉搏信号提取方法的子流程图;
图3为本发明实施例2提供的一种脉搏信号提取装置的结构图;
图4为本发明提供的一种控制器的结构图。
具体实施方式
为了更充分理解本发明的技术内容,下面结合具体实施例对本发明的技术方案进一步介绍和说明,但不局限于此。
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
应当理解,当在本说明书和所附权利要求书中使用时,术语“包括”和“包含”指示所描述特征、整体、步骤、操作、元素和/或组件的存在,但并不排除一个或多个其它特征、整体、步骤、操作、元素、组件和/或其集合的存在或添加。
还应当理解,在本发明说明书和所附权利要求书中使用的术语“和/或”是指相关联列出的项中的一个或多个的任何组合以及所有可能组合,并且包括这些组合。
实施例1
参见图1,本发明实施例1提供了一种脉搏信号提取方法,实施例1的脉搏信号提取方法包括步骤:S101-S107。
S101,将初始脉搏信号中的噪声去除得到滤波脉搏信号,所述初始脉搏信号是在变化压力条件下采集到的。
在本发明实施例中,初始脉搏信号可具体为变化压力下的脉搏信号,所使用的变化压力下的脉搏信号需要通过脉搏传感器获得,同时,通过脉搏传感器连续加压至200mmhg(毫米汞柱)后停止加压随后放压,变化压力下的脉搏信号为该过程中测量到的脉搏信号。通过脉搏传感器获得的变化压力下的脉搏信号会受到噪声的影响,包括:频率在50HZ附近的工频干扰,当使用交流信号供电或者受交流电的辐射,都会有较大幅度的工频干扰信号;频率小于1HZ的基线漂移,该噪声是由于呼吸时胸腔变化引起的低频噪声;有运动肌肉颤动产生的肌电信号,频率分布在高频;这些噪声的存在会影响后续对于脉搏信号的分析,本发明使用0.1-30HZ巴特沃斯带通滤波器,滤除工频干扰、基线漂移和肌电干扰等噪声,得到的滤波脉搏信号进行后续分析。
在一实施例中,将初始脉搏信号中的噪声去除得到滤波脉搏信号,包括:
通过带通滤波器将所述脉搏信号进行滤波,滤除频率为50赫兹的工频干扰、低频基线漂移及肌电干扰的高频噪声,得到所述滤波脉搏信号。
在本发明实施例中,通过使用0.1-30HZ巴特沃斯带通滤波器,滤除工频干扰、基线漂移和肌电干扰等噪声,得到的滤波脉搏信号进行后续分析。
S102,提取所述滤波脉搏信号的包络信号。
具体而言,通过对滤波脉搏信号进行拟合得到的包络信号,进而通过包络信号中的波峰序列可以进行周期分割,进而得到的多个脉搏周期,可以将整一个的滤波脉搏信号具体化,细分成多个脉搏周期形成的整个脉搏信号,进而提高脉搏信号提取的准确性。
在一实施例中,以上步骤S102包括步骤:S1021-S1024。
S1021,计算所述滤波脉搏信号的斜率。
具体而言,通过计算滤波脉搏信号中的斜率,判断斜率的正负,来找出滤波脉搏信号中的各个波峰,进而得到各个波峰的位置也就是波峰的横坐标以及纵坐标。
S1022,检测所述滤波脉搏信号的波峰,根据所述斜率以及所述波峰,查找属于滤波脉搏信号范围内的波峰,得到波峰序列。
具体而言,检测所述滤波脉搏信号的波峰,依次判断斜率的正负,若某点斜率为正,同时其相邻的下一点斜率为负,认为此处为极大值点即滤波脉搏信号的波峰所在处。同时记录所有波峰的横坐标X={X1,X2,...,Xm}与纵坐标Y={Y1,Y2,...,Ym}。
S1023,记录所述波峰序列内的波峰的横坐标以及纵坐标。
具体而言,通过得到的波峰序列的位置后,记录所有波峰的横坐标X={X1,X2,...,Xm}与纵坐标Y={Y1,Y2,...,Ym}。
S1024,根据所述横坐标以及所述纵坐标,将所述波峰序列内的波峰进行拟合,得到所述包络信号。
具体而言,得到的波峰序列为离散的点,要得到滤波脉搏信号的包络信号需要对其进行拟合,求连续两点即连续两个波峰之间的直线方程,将波峰之间利用直线进行拟合,得到包络信号与输入信号长度相同,便于后续处理,其中两点之间的直线方程公式为:
k=(Yp+1-Yp)/(Xp+1-Xp) 0<p<m
b=(Xp+1*Yp-Xp*Yp+1)/(Xp-Xp+1)0<p<m
y=k*x+b
其中,Xp和Xp+1为邻近两个波峰的横坐标,Yp和Yp+1分别为Xp和Xp+1对应的纵坐标,即波峰的幅值。
S103,根据所述包络信号对所述滤波脉搏信号进行周期分割,得到多个脉搏周期。
具体而言,通过得到的波峰序列,可以从波峰序列中查找得到包络信号中的信号最大值,所述信号最大值为所有波峰中纵坐标为最大值的那个波峰。首先将波峰序列中的第一个值记为临时最大值,然后依次将波峰序列中的值与这个临时的最大值进行比较,如果波峰序列中的值大于此临时最大值,则更新该波峰序列中的值为新的临时的最大值;如果波峰序列中的值小于或等于此临时最大值,则保持该临时最大值不变,继续与下一个波峰进行比较;直到波峰序列中所有值都比较完成后,此时的临时最大值就是最终的整个脉搏信号中的最大值,需要说明的是,上述所提到的波峰序列中的值都是参考波峰的纵坐标。记录信号最大值的横坐标Xm和纵坐标Ym。
在一实施例中,参见图2,以上步骤S103包括步骤S1031-S1032。
S1031,将所述波峰序列中的波峰作为周期分割点,所述波峰序列中的波峰分别对应单个脉搏周期。
具体而言,波峰序列中的波峰都是每个脉搏周期中的极大值点,所以通过将所述波峰序列中的波峰作为周期分割点,可以将整个滤波脉搏信号进行标准化的细分,更容易对细分下的滤波脉搏信号进行分析处理,以便中医诊脉更好地提高其精准性。
S1032,根据所述周期分割点对所述滤波脉搏信号进行周期分割得到多个所述脉搏周期。
具体而言,如上述步骤S1033所说,通过根据所述周期分割点对所述滤波脉搏信号进行周期分割得到多个所述脉搏周期,可以更容易对细分下的滤波脉搏信号进行分析处理,以便中医脉诊更进一步的发展。
S104,从所述包络信号中获取信号最大值,将所述信号最大值对应的脉搏周期作为极值脉搏周期。
具体而言,通过将信号最大值对应的脉搏周期作为极值脉搏周期,可以通过极值脉搏周期作为标准点,对中脉所对应的滤波脉搏信号进行划分,以此来得出本发明所说的中脉。
在一实施例中,以上步骤S104包括步骤S1041-S1042。
S1041,根据所述纵坐标,从所述波峰序列中获取最大值波峰。
具体而言,通过所述纵坐标来找出波峰序列中的最大值波峰,也就是找出纵坐标为最大值的那个波峰作为所述最大值波峰。
S1042,将所述最大值波峰作为所述信号最大值。
具体而言,将所述最大值波峰作为信号最大值,可以通过所述信号最大值来得出中脉所在的位置对应的滤波脉搏信号的位置,以此来得出中脉。
S105,将所述极值脉搏周期以及与所述极值脉搏周期相邻的两个脉搏周期作为中脉周期,并将所述中脉周期所对应的滤波脉搏信号作为中脉。
具体而言,如上述所说,将所述极值脉搏周期以及与所述极值脉搏周期相邻的两个脉搏周期作为中脉周期,也就是将所述信号最大值所在的脉搏周期以及其相邻的两个脉搏周期作为中脉周期,那么这个中脉周期包括了三个脉搏周期,这三个脉搏周期所对应的滤波脉搏信号就作为中脉。
取信号最大值前的第二个波峰作为起点,最大值后的第二个波峰作为终点,中间的脉搏信号切分为中医理论中脉对应的脉搏波,公式为:
Pulse_zhong = Pulse_bp[Xm-2:Xm+2]
其中,Pulse_zhong为中医中中脉对应的脉搏波,Pulse_bp为滤波后的脉搏信号,Xm-2为信号最大值的横坐标Xm前的第二个波峰的横坐标,Xm+2为信号最大值的横坐标Xm后的第二个波峰的横坐标;
S106,从多个所述脉搏周期中获取三个在所述中脉周期前的脉搏周期作为浮脉周期,并将所述浮脉周期所对应的滤波脉搏信号作为浮脉。
具体而言,浮脉是在所述中脉周期前的位置,通过在所述中脉周期前的位置获取三个脉搏周期作为浮脉周期,那么所述浮脉周期所对应的滤波脉搏信号作为浮脉。
在一实施例中,以上步骤S106包括以下步骤S1061-S1063。
S1061,根据所述最大值波峰的横坐标乘以第一预设系数得到的第一预设结果。
在本发明实施例中,第一预设系数可具体为0.5,第一预设结果可具体为Xf=β*Xm,其中Xf为浮脉对应的横坐标,Xm为信号最大值的横坐标,β为设定的系数为0.5,需要说明的是,β可以根据实际需要去设置,在此本发明并不作任何限定。
S1062,根据所述第一预设结果获取所述浮脉周期对应的波峰。
具体而言,基于浮脉所在位置将其前后3个周期的脉搏信号作为浮脉对应的脉搏波,取Xf前的第二个波峰作为起点,Xf后的第二个波峰作为终点,中间的脉搏信号切分为中医理论中的浮脉。
S1063,将所述浮脉周期对应的波峰所对应的滤波脉搏信号作为浮脉。
具体而言,通过将所述浮脉周期对应的波峰所对应的滤波脉搏信号作为浮脉可以将所述滤波脉搏信号进行整体标准化细分,可以使所提取得到的脉搏信号更为准确。
S107,从多个所述脉搏周期中获取三个在所述中脉周期后的脉搏周期作为沉脉周期,并将所述沉脉周期所对应的滤波脉搏信号作为沉脉。
具体而言,沉脉是在所述中脉周期后的位置,通过在所述中脉周期后的位置获取三个脉搏周期作为沉脉周期,那么所述沉脉周期所对应的滤波脉搏信号作为沉脉。
在一实施例中,以上步骤S107包括步骤S1071-S1073。
S1071,根据所述最大值波峰的横坐标乘以第二预设系数得到的第二预设结果。
在本发明实施例中,第二预设系数可具体为1.5,第二预设结果可具体为Xc=α*Xm,其中Xc为沉脉对应的横坐标,Xm为信号最大值的横坐标,α为设定的系数为1.5,需要说明的是,α可以根据实际需要去设置,在此本发明并不作任何限定。
S1072,根据所述第二预设结果获取所述沉脉周期对应的波峰。
具体而言,基于沉脉所在位置将其前后3个周期的脉搏信号作为沉脉对应的脉搏波,取Xc前的第二个波峰作为起点,Xc后的第二个波峰作为终点,中间的脉搏信号切分为中医理论中的沉脉。
S1073,将所述沉脉周期对应的波峰所对应的滤波脉搏信号作为沉脉。
具体而言,通过将所述沉脉周期对应的波峰所对应的滤波脉搏信号作为沉脉可以将所述滤波脉搏信号进行整体标准化细分,可以使所提取得到的脉搏信号更为准确。
实施例2
参见图3,本发明实施例2提供了一种脉搏信号提取装置400,该基于变化压力下脉搏波的浮中沉提取装置400包括:预处理单元401、提取单元402、获取单元403、极值单元404、中脉单元405、浮脉单元406、沉脉单元407。
预处理单元401,用于将初始脉搏信号中的噪声去除得到滤波脉搏信号,所述初始脉搏信号是在变化压力条件下采集到的。
提取单元402,用于提取所述滤波脉搏信号的包络信号。
获取单元403,用于根据所述包络信号对所述滤波脉搏信号进行周期分割,得到多个脉搏周期。
极值单元404,用于从所述包络信号中获取信号最大值,将所述信号最大值对应的脉搏周期作为极值脉搏周期。
中脉单元405,用于将所述极值脉搏周期以及与所述极值脉搏周期相邻的两个脉搏周期作为中脉周期,并将所述中脉周期所对应的滤波脉搏信号作为中脉。
浮脉单元406,用于从多个所述脉搏周期中获取三个在所述中脉周期前的脉搏周期作为浮脉周期,并将所述浮脉周期所对应的滤波脉搏信号作为浮脉。
沉脉单元407,用于从多个所述脉搏周期中获取三个在所述中脉周期后的脉搏周期作为沉脉周期,并将所述沉脉周期所对应的滤波脉搏信号作为沉脉。
在一实施例中,将初始脉搏信号中的噪声去除得到滤波脉搏信号,包括:
通过带通滤波器将所述脉搏信号进行滤波,滤除频率为50赫兹的工频干扰、低频基线漂移及肌电干扰的高频噪声,得到所述滤波脉搏信号。
在一实施例中,所述提取所述滤波脉搏信号的包络信号,包括:
计算所述滤波脉搏信号的斜率;
检测所述滤波脉搏信号的波峰,根据所述斜率以及所述波峰,查找属于滤波脉搏信号范围内的波峰,得到波峰序列;
记录所述波峰序列内的波峰的横坐标以及纵坐标;
根据所述横坐标以及所述纵坐标,将所述波峰序列内的波峰进行拟合,得到所述包络信号。
在一实施例中,所述从所述包络信号中获取信号最大值,包括:
根据所述纵坐标,从所述波峰序列中获取最大值波峰;
将所述最大值波峰作为所述信号最大值。
在一实施例中,所述根据所述包络信号对所述滤波脉搏信号进行周期分割,得到多个脉搏周期,包括:
将所述波峰序列中的波峰作为周期分割点,所述波峰序列中的波峰分别对应单个脉搏周期;
根据所述周期分割点对所述滤波脉搏信号进行周期分割得到多个所述脉搏周期。
在一实施例中,所述从多个所述脉搏周期中获取三个在所述中脉周期前的脉搏周期作为浮脉周期,并将所述浮脉周期所对应的滤波脉搏信号作为浮脉,包括:
根据所述最大值波峰的横坐标乘以第一预设系数得到的第一预设结果;
根据所述第一预设结果获取所述浮脉周期对应的波峰;
将所述浮脉周期对应的波峰所对应的滤波脉搏信号作为浮脉。
在一实施例中,所述从多个所述脉搏周期中获取三个在所述中脉周期后的脉搏周期作为沉脉周期,并将所述沉脉周期所对应的滤波脉搏信号作为沉脉,包括:
根据所述最大值波峰的横坐标乘以第二预设系数得到的第二预设结果;
根据所述第二预设结果获取所述沉脉周期对应的波峰;
将所述沉脉周期对应的波峰所对应的滤波脉搏信号作为沉脉。
实施例3
参见图4,本发明实施例提供的一种电子设备,所述电子设备包括控制器,控制器包括处理器111、通信接口112、存储器113和通信总线114,其中,处理器111,通信接口112,存储器113通过通信总线114完成相互间的通信。
存储器113,用于存放计算机程序;
处理器111,用于执行存储器113上所存放的程序,实现实施例1提供的基于变化压力下脉搏波的浮中沉提取方法。
本发明实施例3还提供了一种计算机可读存储介质,其上存储有计算机程序,所述计算机程序被处理器111执行时实现如实施例1提供的脉搏信号提取方法的步骤。
需要说明的是,在本文中,诸如“第一”和“第二”等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来,而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。
以上所述仅是本发明的具体实施方式,使本领域技术人员能够理解或实现本发明。对这些实施例的多种修改对本领域的技术人员来说将是显而易见的,本文中所定义的一般原理可以在不脱离本发明的精神或范围的情况下,在其它实施例中实现。因此,本发明将不会被限制于本文所示的这些实施例,而是要符合与本文所发明的原理和新颖特点相一致的最宽的范围。

Claims (10)

1.一种脉搏信号提取方法,其特征在于,包括:
将初始脉搏信号中的噪声去除得到滤波脉搏信号,所述初始脉搏信号是在变化压力条件下采集到的;
提取所述滤波脉搏信号的包络信号;
根据所述包络信号对所述滤波脉搏信号进行周期分割,得到多个脉搏周期;
从所述包络信号中获取信号最大值,将所述信号最大值对应的脉搏周期作为极值脉搏周期;
将所述极值脉搏周期以及与所述极值脉搏周期相邻的两个脉搏周期作为中脉周期,并将所述中脉周期所对应的滤波脉搏信号作为中脉;
从多个所述脉搏周期中获取三个在所述中脉周期前的脉搏周期作为浮脉周期,并将所述浮脉周期所对应的滤波脉搏信号作为浮脉;
从多个所述脉搏周期中获取三个在所述中脉周期后的脉搏周期作为沉脉周期,并将所述沉脉周期所对应的滤波脉搏信号作为沉脉。
2.根据权利要求1所述的脉搏信号提取方法,其特征在于,将初始脉搏信号中的噪声去除得到滤波脉搏信号,包括:
通过带通滤波器将所述脉搏信号进行滤波,滤除频率为50赫兹的工频干扰、低频基线漂移及肌电干扰的高频噪声,得到所述滤波脉搏信号。
3.根据权利要求1所述的脉搏信号提取方法,其特征在于,所述提取所述滤波脉搏信号的包络信号,包括:
计算所述滤波脉搏信号的斜率;
检测所述滤波脉搏信号的波峰,根据所述斜率以及所述波峰,查找属于所述滤波脉搏信号范围内的波峰,得到波峰序列;
记录所述波峰序列内的波峰的横坐标以及纵坐标;
根据所述横坐标以及所述纵坐标,将所述波峰序列内的波峰进行拟合,得到所述包络信号。
4.根据权利要求3所述的脉搏信号提取方法,其特征在于,所述从所述包络信号中获取信号最大值,包括:
根据所述纵坐标,从所述波峰序列中获取最大值波峰;
将所述最大值波峰作为所述信号最大值。
5.根据权利要求3所述的脉搏信号提取方法,其特征在于,所述根据所述包络信号对所述滤波脉搏信号进行周期分割,得到多个脉搏周期,包括:
将所述波峰序列中的波峰作为周期分割点,所述波峰序列中的波峰分别对应单个脉搏周期;
根据所述周期分割点对所述滤波脉搏信号进行周期分割得到多个所述脉搏周期。
6.根据权利要求1所述的脉搏信号提取方法,其特征在于,所述从多个所述脉搏周期中获取三个在所述中脉周期前的脉搏周期作为浮脉周期,并将所述浮脉周期所对应的滤波脉搏信号作为浮脉,包括:
根据所述最大值波峰的横坐标乘以第一预设系数得到的第一预设结果;
根据所述第一预设结果获取所述浮脉周期对应的波峰;
将所述浮脉周期对应的波峰所对应的滤波脉搏信号作为浮脉。
7.根据权利要求1所述的脉搏信号提取方法,其特征在于,所述从多个所述脉搏周期中获取三个在所述中脉周期后的脉搏周期作为沉脉周期,并将所述沉脉周期所对应的滤波脉搏信号作为沉脉,包括:
根据所述最大值波峰的横坐标乘以第二预设系数得到的第二预设结果;
根据所述第二预设结果获取所述沉脉周期对应的波峰;
将所述沉脉周期对应的波峰所对应的滤波脉搏信号作为沉脉。
8.一种脉搏信号提取装置,其特征在于,包括用于执行如权利要求1-7任一项所述的方法的单元。
9.一种电子设备,其特征在于,所述电子设备包括控制器,所述控制器包括处理器、通信接口、存储器和通信总线,其中,处理器,通信接口,存储器通过通信总线完成相互间的通信;
存储器,用于存放计算机程序;
处理器,用于执行存储器上所存放的程序时,实现如权利要求1-7任一项所述的方法的步骤。
10.一种计算机可读存储介质,其上存储有计算机程序,其特征在于,所述计算机程序被处理器执行时实现如权利要求1-7任一项所述的方法的步骤。
CN202210632256.9A 2022-06-07 2022-06-07 一种脉搏信号提取方法、装置、电子设备及存储介质 Active CN114711733B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210632256.9A CN114711733B (zh) 2022-06-07 2022-06-07 一种脉搏信号提取方法、装置、电子设备及存储介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210632256.9A CN114711733B (zh) 2022-06-07 2022-06-07 一种脉搏信号提取方法、装置、电子设备及存储介质

Publications (2)

Publication Number Publication Date
CN114711733A true CN114711733A (zh) 2022-07-08
CN114711733B CN114711733B (zh) 2023-07-28

Family

ID=82232534

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210632256.9A Active CN114711733B (zh) 2022-06-07 2022-06-07 一种脉搏信号提取方法、装置、电子设备及存储介质

Country Status (1)

Country Link
CN (1) CN114711733B (zh)

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH07299044A (ja) * 1994-05-10 1995-11-14 Sanyo Electric Co Ltd 脈拍検出装置
JP2006122377A (ja) * 2004-10-29 2006-05-18 Sanyo Electric Co Ltd 周期算出装置及びプログラム
CN103027669A (zh) * 2011-09-30 2013-04-10 Ge医疗系统环球技术有限公司 用于判断脉象浮沉度的方法和设备
CN105748056A (zh) * 2016-04-01 2016-07-13 广州视源电子科技股份有限公司 血压检测方法及其系统
US20170156609A1 (en) * 2015-06-19 2017-06-08 Boe Technology Group Co., Ltd. Pulse cycle detection device and method, and wearable electronic device
CN108403094A (zh) * 2018-03-23 2018-08-17 中国科学院微电子研究所 一种识别脉搏波波峰的方法
CN109464133A (zh) * 2018-11-15 2019-03-15 天津大学 一种便携自动脉象检测系统
CN113288075A (zh) * 2021-06-01 2021-08-24 山西恒朝科技有限责任公司 浮中沉脉位识别方法及系统
CN114431840A (zh) * 2022-04-08 2022-05-06 北京大学深圳研究生院 一种脉搏采集装置、脉搏采集方法和系统

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH07299044A (ja) * 1994-05-10 1995-11-14 Sanyo Electric Co Ltd 脈拍検出装置
JP2006122377A (ja) * 2004-10-29 2006-05-18 Sanyo Electric Co Ltd 周期算出装置及びプログラム
CN103027669A (zh) * 2011-09-30 2013-04-10 Ge医疗系统环球技术有限公司 用于判断脉象浮沉度的方法和设备
US20170156609A1 (en) * 2015-06-19 2017-06-08 Boe Technology Group Co., Ltd. Pulse cycle detection device and method, and wearable electronic device
CN105748056A (zh) * 2016-04-01 2016-07-13 广州视源电子科技股份有限公司 血压检测方法及其系统
CN108403094A (zh) * 2018-03-23 2018-08-17 中国科学院微电子研究所 一种识别脉搏波波峰的方法
CN109464133A (zh) * 2018-11-15 2019-03-15 天津大学 一种便携自动脉象检测系统
CN113288075A (zh) * 2021-06-01 2021-08-24 山西恒朝科技有限责任公司 浮中沉脉位识别方法及系统
CN114431840A (zh) * 2022-04-08 2022-05-06 北京大学深圳研究生院 一种脉搏采集装置、脉搏采集方法和系统

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
李慎安 主编: "《测量不确定度的简化评定》", 31 May 2004 *

Also Published As

Publication number Publication date
CN114711733B (zh) 2023-07-28

Similar Documents

Publication Publication Date Title
Wu et al. Research of fetal ECG extraction using wavelet analysis and adaptive filtering
CN109171712B (zh) 心房颤动识别方法、装置、设备及计算机可读存储介质
US9572504B2 (en) Continuous non-invasive monitoring of a pregnant human subject
CA2979135C (en) Systems, apparatus and methods for sensing fetal activity
Wang et al. A robust signal preprocessing framework for wrist pulse analysis
CN108388912A (zh) 基于多传感器特征优化算法的睡眠分期方法
KR101933338B1 (ko) 임신부의 연속적인 비침습적 모니터링
CN110680308B (zh) 基于改进emd与阈值法融合的心电信号去噪方法
Aston et al. Beyond HRV: attractor reconstruction using the entire cardiovascular waveform data for novel feature extraction
CN107233093B (zh) R波的检测方法、装置及电子设备
CN102835954A (zh) 一种心拍波形模板生成、室性早搏检测处理方法及装置
Zheng et al. Single-lead fetal electrocardiogram estimation by means of combining R-peak detection, resampling and comb filter
CN107997754B (zh) 一种智能床垫系统及人体生理特征数据提取方法
CN108509823A (zh) Qrs波群的检测方法及装置
CN110192853A (zh) 一种基于心电熵值图的心电异常评估方法
Jaswal et al. QRS detection using wavelet transform
CN107273827A (zh) 一种心电信号r波检测方法及装置
CN114732424A (zh) 基于表面肌电信号提取肌肉疲劳状态复杂网络属性的方法
CN103829944A (zh) 基于模式识别的胸阻抗信号处理方法
Martin et al. Detection of ECG significant waves for biventricular pacing treatment
CN107320097A (zh) 利用肌电信号边际谱熵提取肌肉疲劳特征的方法和装置
Elbuni et al. ECG parameter extraction algorithm using (DWTAE) algorithm
KR101048763B1 (ko) 신호 검출 장치 및 방법
CN114711733A (zh) 一种脉搏信号提取方法、装置、电子设备及存储介质
CN106725487A (zh) 一种呼吸率提取方法及装置

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