CN105433936B - 心输出量获得方法及装置 - Google Patents

心输出量获得方法及装置 Download PDF

Info

Publication number
CN105433936B
CN105433936B CN201610061495.8A CN201610061495A CN105433936B CN 105433936 B CN105433936 B CN 105433936B CN 201610061495 A CN201610061495 A CN 201610061495A CN 105433936 B CN105433936 B CN 105433936B
Authority
CN
China
Prior art keywords
data segment
crest
waveform
wave
wave crest
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
Application number
CN201610061495.8A
Other languages
English (en)
Other versions
CN105433936A (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.)
Zhejiang Maomao Technology Co ltd
Original Assignee
Beijing Xinliang Technology Co Ltd
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 Beijing Xinliang Technology Co Ltd filed Critical Beijing Xinliang Technology Co Ltd
Priority to CN201610061495.8A priority Critical patent/CN105433936B/zh
Publication of CN105433936A publication Critical patent/CN105433936A/zh
Application granted granted Critical
Publication of CN105433936B publication Critical patent/CN105433936B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/02—Detecting, 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
    • A61B5/026—Measuring blood flow
    • A61B5/0295—Measuring blood flow using plethysmography, i.e. measuring the variations in the volume of a body part as modified by the circulation of blood therethrough, e.g. impedance plethysmography
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/02—Detecting, 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
    • A61B5/02028—Determining haemodynamic parameters not otherwise provided for, e.g. cardiac contractility or left ventricular ejection fraction
    • 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/7225—Details of analog processing, e.g. isolation amplifier, gain or sensitivity adjustment, filtering, baseline or drift compensation

Abstract

本发明提供了一种心输出量获得方法及装置，该心输出量获得方法包括获得时长为第二时间片段的多个数据段；对每个数据段修正基线漂移后，获得每个数据段所对应的心跳间隔的多个波峰，随后获得多个波峰相邻两个波峰中低波峰的波峰高度与高波峰的波峰高度的比值作为与数据段分别对应的多个波峰高度比值；获得波峰高度比值属于预定阈值范围的数据段，并在上述的数据段中获得多个波峰高度比值的平均值最大的数据段的每个完整波形；根据从每个完整波形获得的平均值、波峰以及波谷分别获得对应的多个初选心输出量，取多个初选心输出量的中值，便可以获得心输出量。本发明实施例提供的方法能改善现有技术中测量方法过于复杂，测量门槛较高的不足。

Description

心输出量获得方法及装置
技术领域
本发明涉及心脏测量技术领域，具体而言，涉及一种心输出量获得方法及装置。
背景技术
现在广泛采用的测量心输出量的方法为热稀释法。热稀释法是将漂浮导管经右心房输入肺动脉，然后经该导管向右心房注入冷生理盐水，该导管的前端放置有温度传感器，当该生理盐水与血液混合后发生温度变化，混合后的血液进入肺动脉时，将被温度传感器感知，根据冷生理盐水注入的时刻和混合后温度的变化情况，分析心输出量。
采用热稀释法时，由于血液和冷生理盐水在右心房混合，而在肺动脉处测量血液和冷生理盐水混合后的温度变化，进而测量心脏射血的能力。在混合的同时，心脏也在射血，整个测量过程是一个动态的过程，故注射速度和均匀性会直接影响冷生理盐水和血液混合的充分性，因此对心输出量的测量结果有直接的影响，并且该测量方法过于复杂，测量门槛较高。
发明内容
有鉴于此，本发明实施例提供了一种心输出量获得方法及装置，以改善现有技术中测量方法过于复杂，测量门槛较高的不足。
为实现上述目的，本发明实施例提供如下技术方案：
一种心输出量获得方法，所述方法包括：获得初选脉搏信号所对应的波形，每间隔第一时间片段，对所述波形进行时长为第二时间片段的截取以得到多个数据段，对所述多个数据段中的每个数据段修正基线漂移；获得修正基线漂移后的每个数据段对应的心跳间隔的多个波峰，获得所述多个波峰的相邻两个波峰中的低波峰的波峰高度与高波峰的波峰高度的比值作为对应数据段的多个波峰高度比值；当修正基线漂移后的所述多个数据段中的数据段的所述多个波峰高度比值均属于预定阈值范围时，获得该数据段；获得所述多个波峰高度比值均属于预定阈值范围的多个数据段中的所述多个波峰高度比值的平均值最大的数据段的每个完整波形；根据所述每个完整波形的平均值、波峰以及波谷获得所述每个完整波形分别对应的多个初选心输出量，对所述多个初选心输出量取中值，获得心输出量。
相应的，本发明实施例还提供了一种心输出量获得装置，所述装置包括：第一获取模块，用于获得初选脉搏信号所对应的波形，每间隔第一时间片段，对所述波形进行时长为第二时间片段的截取以得到多个数据段，对所述多个数据段中的每个数据段修正基线漂移；第二获取模块，用于获得修正基线漂移后的每个数据段对应的心跳间隔的多个波峰，获得所述多个波峰的相邻两个波峰中的低波峰的波峰高度与高波峰的波峰高度的比值作为对应数据段的多个波峰高度比值；第三获取模块，用于在修正基线漂移后的所述多个数据段中的数据段的所述多个波峰高度比值均属于预定阈值范围时，获得该数据段；第四获取模块，用于获得所述多个波峰高度比值均属于预定阈值范围的多个数据段中的所述多个波峰高度比值的平均值最大的数据段的每个完整波形；第五获取模块，用于根据所述每个完整波形的平均值、波峰以及波谷获得所述每个完整波形分别对应的多个初选心输出量，对所述多个初选心输出量取中值，获得心输出量。
本发明实施例提供的心输出量获得方法及装置，具有如下有益效果：
本发明实施例提供的心输出量获得方法及装置获得时长为第二时间片段的多个数据段；对每个数据段修正基线漂移后，获得每个数据段所对应的心跳间隔的多个波峰，随后获得多个波峰相邻两个波峰中低波峰的波峰高度与高波峰的波峰高度的比值作为与数据段分别对应的多个波峰高度比值；获得波峰高度比值属于预定阈值范围的数据段，并在上述的数据段中获得多个波峰高度比值的平均值最大的数据段的每个完整波形；根据从每个完整波形获得的平均值、波峰以及波谷分别获得对应的多个初选心输出量，取多个初选心输出量的中值，便可以获得心输出量。本发明实施例提供的方法能改善现有技术中测量方法过于复杂，测量门槛较高的不足。
附图说明
为了更清楚的说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单的介绍，显而易见地，下面描述中的附图仅仅是本发明实施例的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1示出了本发明实施例提供的一种心输出量获得方法及装置的应用环境示意图；
图2示出了本发明实施例提供的数据分析终端的方框示意图；
图3是本发明第一实施例提供的心输出量获得方法的流程图；
图4是本发明第二实施例提供的心输出量获得方法的流程图；
图5是本发明第三实施例提供的心输出量获得装置的结构示意图；
图6是本发明第四实施例提供的心输出量获得装置的结构示意图。
具体实施方式
下面将结合本发明实施例中附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。以下对在附图中提供的本发明的实施例的详细描述并非旨在限制要求保护的本发明的范围，而是仅仅表示本发明的选定实施例。基于本发明的实施例，本领域技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例，都属于本发明保护的范围。
应注意到：相似的标号和字母在下面的附图中表示类似项，因此，一旦某一项在一个附图中被定义，则在随后的附图中不需要对其进行进一步定义和解释。同时，在本发明实施例的描述中，术语“第一”、“第二”等仅用于区分描述，而不能理解为指示或暗示相对重要性。
如图1所示，为本发明实施例提供的一种心输出量获得方法及装置的应用环境示意图，包括心跳间隔采集终端101和数据分析终端102，心跳间隔采集终端101为采集人体等生物的心跳间隔，具体可以是心电监护仪或者采用采集人体的脉搏图像，然后再通过光电容积描记法(Photoplethysmography，PPG)对脉搏图像处理后获得心跳间隔。数据分析终端102为具有数据分析能力的终端，可以是计算机、智能手机或者其他由具有数据分析能力的芯片构成的集成电路装置。
如图2所示，为所述数据分析终端102的方框示意图。所述数据分析终端102包括：心输出量获得装置、存储器202、存储控制器203、处理器204、外设接口201、输入输出单元206和显示单元205。
所述存储器202、存储控制器203、处理器204、外设接口201、输入输出单元206和显示单元205各元件相互之间直接或间接地电性连接，以实现数据的传输或交互。例如，这些元件相互之间可通过一条或多条通讯总线或信号线实现电性连接。所述心输出量获得装置包括至少一个可以软件或固件(firmware)的形式存储于所述存储器202中或固化在所述数据分析终端102的操作系统(operating system，OS)中的软件功能模块。所述处理器204用于执行存储器202中存储的可执行模块，例如所述目标检测装置包括的软件功能模块或计算机程序。
其中，存储器202可以是，但不限于，随机存取存储器(Random Access Memory，RAM)，只读存储器(Read Only Memory，ROM)，可编程只读存储器(Programmable Read-OnlyMemory，PROM)，可擦除只读存储器(Erasable Programmable Read-Only Memory，EPROM)，电可擦除只读存储器(Electric Erasable Programmable Read-Only Memory，EEPROM)等。其中，存储器202用于存储程序，所述处理器204在接收到执行指令后，执行所述程序，前述本发明实施例任一实施例揭示的流过程定义的服务器所执行的方法可以应用于处理器204中，或者由处理器204实现。
处理器204可能是一种集成电路芯片，具有信号的处理能力。上述的处理器204可以是通用处理器，包括中央处理器(Central Processing Unit，简称CPU)、网络处理器(Network Processor，简称NP)等；还可以是数字信号处理器(DSP)、专用集成电路(ASIC)、现成可编程门阵列(FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件。可以实现或者执行本发明实施例中的公开的各方法、步骤及逻辑框图。通用处理器204可以是微处理器或者该处理器204也可以是任何常规的处理器等。
所述外设接口201将各种输入/输出装置耦合至处理器204以及存储器202。在一些实施例中，外设接口201，处理器204以及存储控制器203可以在单个芯片中实现。在其他一些实例中，他们可以分别由独立的芯片实现。
显示单元205在所述数据分析终端102与用户之间提供一个交互界面(例如用户操作界面)或用于显示图像数据给用户参考。在本实施例中，所述显示单元205可以是液晶显示器或触控显示器。若为触控显示器，其可为支持单点和多点触控操作的电容式触控屏或电阻式触控屏等。支持单点和多点触控操作是指触控显示器能感应到来自该触控显示器上一个或多个位置处同时产生的触控操作，并将该感应到的触控操作交由处理器204进行计算和处理。
参阅图3，是本发明第一实施例提供的应用于图2所示的数据分析终端的一种心输出量获得方法的流程图，所述方法包括：
步骤S101，获得初选脉搏信号所对应的波形，每间隔第一时间片段，对所述波形进行时长为第二时间片段的截取以得到多个数据段，对所述多个数据段中的每个数据段修正基线漂移。
初选脉搏信号指的是未经过滤波处理的含有噪声的脉搏信号。初选脉搏信号具体可以由以下步骤获得：
获得多个连续的图像，所述图像包括反映人体脉搏跳动的区域的图像；
获得所述多个连续的图像中的每一个图像的像素的红色通道数值，对所述多个连续的图像中的每一个图像的像素的红色通道数值分别进行数值累加，得到初选脉搏信号。
获得初选脉搏信号对应的波形后，可以以时长为第二时间片段的时间对上述的波形截取以获得多个数据段，截取上一个数据段与截取下一个数据段之间的时间间隔为第一时间间隔。分析时可以针对数据段进行分析。
以第一时间片段为2秒，所述第二时间片段为5秒为例。获得初选脉搏信号所对应的波形后，每隔2秒，截取时长为5秒的数据段。
可以理解，所述第一时间片段和第二时间片段的时长也可以根据实际情况设定为其他数值，例如，第一时间片段为2.3秒，第二时间片段为5.4秒，所述第一时间片段和第二时间片段的具体时长不应该理解为是对本发明的限制。
截取多个数据段后，对截取后的多个数据段中的每个数据段进行基线漂移修正，基线漂移修正具体可以通过进行小波变换来修正基线漂移，也可以采用RC高通滤波器来补偿基线漂移。
步骤S102，获得修正基线漂移后的每个数据段对应的心跳间隔的多个波峰，获得所述多个波峰的相邻两个波峰中的低波峰的波峰高度与高波峰的波峰高度的比值作为对应数据段的多个波峰高度比值。
对初选脉搏信号修正基线漂移后，可以获得去除基线漂移的脉搏信号；然后还需要进行低通滤波以去除高频噪声，获得去除高频噪声后的脉搏信号；对滤波后的信号去除时间异常波峰以及高度异常波峰后就可以根据筛选后的波峰获得心跳间隔的波形。心跳间隔的波形由于经过了滤波，所以比所述初选脉搏信号的波形更加精准。当然，除了修正基线漂移，为了获得更精确的心跳间隔还可以对初选脉搏信号进行其他处理，本发明的具体实现方式并不以此为限。
每个数据段所对应的心跳间隔的波形均具有多个波峰，获得所述多个波峰的相邻两个波峰中的低波峰的波峰高度与高波峰的波峰高度的比值以作为对应数据段的多个波峰高度比值。
例如，所述每个数据段的心跳间隔的多个波峰可以包括五个波峰。当所述多个波峰为五个波峰时，相邻两个波峰的所述波峰高度比值的数量为四个。
可以理解，所述心跳间隔的波峰以及波峰高度比值的具体数量不应该理解为是对本发明的限制。
步骤S103，当修正基线漂移后的所述多个数据段中的数据段的所述多个波峰高度比值均属于预定阈值范围时，获得该数据段。
优选的，预定阈值范围为大于0.9并且小于或等于1，也就是说当一个数据段的多个波峰高度比值均大于0.9且小于或等于1时，获得该数据段。当属于上述预定阈值范围时，多个波峰的高度较为均匀，起伏波动较小。
步骤S104，获得所述多个波峰高度比值均属于预定阈值范围的多个数据段中的所述多个波峰高度比值的平均值最大的数据段的每个完整波形。
从属于预定阈值范围的多个数据段中，获得数据段对应的多个波峰高度比值的平均值最大的数据段，进而获得该数据段的每个完整波形。
步骤S105，根据所述每个完整波形的平均值、波峰以及波谷获得所述每个完整波形分别对应的多个初选心输出量，对所述多个初选心输出量取中值，获得心输出量。
优选地，可以根据来获得参数K，其中平均值为M，波峰为H以及波谷为L。获得参数K后，可以根据 来获得初选心输出量CO，然后取多个初选心输出量CO的中值，即可以获得心输出量。
中值指的是将一组数据按照从小到大排列或者从大到小排列，当该组数据的个数为奇数时，取该组数据中间的数字；当该组数据的个数为偶数时，取该组数据中间两个数的平均值。
故先将所述多个初选心输出量按照从大到小或者从小到大的顺序排列，当所述多个初选心输出量的个数为奇数时，取多个初选心输出量中间的数值；当所述多个初选心输出量的个数为偶数时，取多个初选心输出量中间两个数值的平均值。
本发明实施例提供的方法能改善现有技术中测量方法过于复杂，测量门槛较高的不足。
参阅图4，是本发明第二实施例提供的应用于图2所示的数据分析终端的一种心输出量获得方法的流程图，所述方法包括：
步骤S201，获得多个连续的图像，所述图像包括反映人体脉搏跳动的区域的图像。
具体可以通过手机的摄像头来获得多个连续的图像，所述多个连续的图像具体可以为视频。反映人体脉搏跳动的区域具体可以为人的手指。
也就是说，具体可以通过手机的摄像头来获得人的手指部位的视频信息。
具体地，反映人体脉搏跳动的区域可以为手指，也可以为其他的区域如脖颈处，反映人体脉搏跳动的具体的区域不应该理解为是对本发明的限制。
步骤S202，获得所述多个连续的图像中的每一个图像的像素的红色通道数值，对所述多个连续的图像中的每一个图像的像素的红色通道数值分别进行数值累加，得到初选脉搏信号。
该步骤是对光电容积描记法(Photoplethysmography,PPG)测量原理的应用。PPG方法需要特定波长的光源，常用红光或近红光，所以进行红色通道数值的累加。
对所述多个连续的图像中的每一个图像的像素的红色通道数值分别进行数值累加，由于每一个图像的像素的红色通道数值的累加值不同，故可以获得多个互不相同的红色通道数据的累加值。
上述的多个互不相同的红色通道数据的累加值具体可以反映皮肤表面光学特性的变化，即通过探测皮肤表面光学特性的变化可以获取血液容积脉搏信号。
步骤S203，获得初选脉搏信号所对应的波形，每间隔第一时间片段，对所述波形进行时长为第二时间片段的截取以得到多个数据段。
获得初选脉搏信号对应的波形后，可以以时长为第二时间片段的时间对上述的波形截取以获得多个数据段，分析时可以针对数据段进行分析。
所述第一时间片段具体可以为2秒，所述第二时间片段具体可以为5秒。
可以理解，所述第一时间片段的时长可以为2秒，也可以为其他数值如2.3秒，所述第一时间片段的具体时长不应该理解为是对本发明的限制。
可以理解，所述第二时间片段的时长可以为5秒，也可以为其他数值如5.4秒，所述第二时间片段的具体时长不应该理解为是对本发明的限制。
步骤S204，对所述初选脉搏信号进行滤波，获得滤波后的脉搏信号。
具体地，可以通过非递归型滤波器对所述初选脉搏信号进行低通滤波以去除高频噪声，获得所述去除高频噪声后的脉搏信号。
然后再通过离散小波滤波器对所述去除高频噪声后的脉搏信号进行滤波以去除基线漂移，获得所述滤波后的脉搏信号。
在使用离散小波滤波器对所述去除高频噪声后的脉搏信号进行滤波时，具体可以使用离散小波滤波器对上述的脉搏信号进行六层分解。
步骤S205，根据所述滤波后的脉搏信号，获得所述滤波后的脉搏信号的波峰以及波谷，并根据所述波峰以及波谷获得归一化波峰。
具体地，可以对经过滤波后的脉搏信号利用二阶导数找出滤波后的脉搏信号的每个波峰以及每个波谷的位置。
对上述每个波峰中的一个波峰做垂线。
连接与该波峰相邻的两个波谷获得所述两个波谷的连线。获得波峰的垂线与所述两个波谷的连线的交点，随后取所述波峰到所述交点的距离作为该波峰的高度。
通过上述的方法可以获得归一化波峰。
步骤S206，去除相邻两个波峰的时间差值超过第一预定阈值范围的时间异常波峰以及所述归一化波峰的波峰高度超过第二预定阈值范围的高度异常波峰，得到筛选后的波峰。
具体地，所述第一预定阈值范围可以为150毫秒，所述第二预定阈值范围可以为所述波峰高度的均值的1.5倍。
可以将波峰的时间差值超过150毫秒的时间异常波峰去除，并可以将波峰高度超过波峰高度的均值的1.5倍的高度异常波峰去除，获得经过筛选后的波峰。
步骤S207，获得所述筛选后的波峰中相邻的两个波峰的时间间隔，所述相邻的两个波峰的时间间隔为所述心跳间隔。
经过时间长度的筛选以及波峰高度的筛选后获得筛选后的波峰。可以根据上述的筛选后的波峰来获得相邻两个波峰的时间间隔，该相邻两个波峰的时间间隔即为心跳间隔。
步骤S208，获得每个数据段对应的心跳间隔的多个波峰，获得所述多个波峰的相邻两个波峰中的低波峰的波峰高度与高波峰的波峰高度的比值作为对应数据段的多个波峰高度比值。
对初选脉搏信号修正基线漂移以及其他的滤波处理后可以获得心跳间隔的波形。心跳间隔的波形由于经过了滤波，所以比所述初选脉搏信号的波形更加精准。
每个数据段所对应的心跳间隔的波形均具有多个波峰，获得所述多个波峰的相邻两个波峰中的低波峰的波峰高度与高波峰的波峰高度的比值以作为对应数据段的多个波峰高度比值。
具体地，所述每个数据段的心跳间隔的多个波峰可以包括五个波峰。当所述多个波峰为五个波峰时，相邻两个波峰的所述波峰高度比值可以包括四个波峰高度比值。
步骤S209，当修正基线漂移后的所述多个数据段中的数据段的所述多个波峰高度比值均属于预定阈值范围时，获得该数据段。
所述多个波峰高度比值均属于预定阈值范围具体可以包括所述多个波峰高度比值均大于0.9的同时，小于或等于1。当多个波峰高度均属于上述预定阈值范围时，多个波峰高度较为均匀，起伏波动较小。
数据段的所述多个波峰高度比值均属于预定阈值范围具体指的是该数据段的所有波峰高度比值均满足大于0.9且小于或等于1的条件时，才获得该数据段。
步骤S210，对所述多个波峰高度比值均满足预定阈值范围的数据段的所述多个波峰高度比值取平均值作为该数据段的波形质量。
步骤S211，按照所述波形质量的数值大小进行排序，获得数值最大的波形质量所对应的数据段。
步骤S212，获得所述数值最大的波形质量所对应的数据段的每个完整波形。
步骤S210至步骤S212与第一实施例中的步骤S104实现的效果相同，在此便不做赘述。
步骤S213，对所述完整波形的数值做积分，获得所述完整波形的平均值M，波峰H以及波谷L。
步骤S214，根据第一公式获得参数K。
所述第一公式具体可以为根据公式可以算出参数K。
步骤S215，根据第二公式、参数K，获得初选心输出量CO。
所述第二公式具体可以为根据公式可以获得初选心输出量。
步骤S216，对所述多个初选心输出量取中值，获得心输出量。
因为一个数据段中具有多个完整波形，故可以获得与所述多个完整波形分别对应的多个初选心输出量。所述多个初选心输出量的数量与所述多个完整波形的数量相同。
然后对该数据段的多个初选心输出量取中值，便可以获得心输出量。
图5示出了本发明第三实施例提供的一种心输出量获得装置500，包括：
第一获取模块510，用于获得初选脉搏信号所对应的波形，每间隔第一时间片段，对所述波形进行时长为第二时间片段的截取以得到多个数据段，对所述多个数据段中的每个数据段修正基线漂移；
第二获取模块520，用于获得修正基线漂移后的每个数据段对应的心跳间隔的多个波峰，获得所述多个波峰的相邻两个波峰中的低波峰的波峰高度与高波峰的波峰高度的比值作为对应数据段的多个波峰高度比值；
第三获取模块530，用于在修正基线漂移后的所述多个数据段中的数据段的所述多个波峰高度比值均属于预定阈值范围时，获得该数据段；
第四获取模块540，用于获得所述多个波峰高度比值均属于预定阈值范围的多个数据段中的所述多个波峰高度比值的平均值最大的数据段的每个完整波形；
第五获取模块550，用于根据所述每个完整波形的平均值、波峰以及波谷获得所述每个完整波形分别对应的多个初选心输出量，对所述多个初选心输出量取中值，获得心输出量。
图6示出了本发明第四实施例提供的一种心输出量获得装置600，包括：
图像获取模块610，用于获得多个连续的图像，所述图像包括反映人体脉搏跳动的区域的图像。
初选脉搏信号获取模块620，用于获得所述多个连续的图像中的每一个图像的像素的红色通道数值，对所述多个连续的图像中的每一个图像的像素的红色通道数值分别进行数值累加，得到初选脉搏信号。
第一获取模块510，用于获得初选脉搏信号所对应的波形，每间隔第一时间片段，对所述波形进行时长为第二时间片段的截取以得到多个数据段。
滤波模块630，用于对所述初选脉搏信号进行滤波，获得滤波后的脉搏信号。
归一化模块640，用于根据所述滤波后的脉搏信号，获得所述滤波后的脉搏信号的波峰以及波谷，并根据所述波峰以及波谷获得归一化波峰。
波峰筛选模块650，用于去除相邻两个波峰的时间差值超过第一预定阈值范围的时间异常波峰以及所述归一化波峰的波峰高度超过第二预定阈值范围的高度异常波峰，得到筛选后的波峰。
心跳间隔获取模块660，用于获得所述筛选后的波峰中相邻的两个波峰的时间间隔，所述相邻的两个波峰的时间间隔为所述心跳间隔。
第二获取模块520，还用于获得每个数据段对应的心跳间隔的多个波峰，获得所述多个波峰的相邻两个波峰中的低波峰的波峰高度与高波峰的波峰高度的比值作为对应数据段的多个波峰高度比值。
第三获取模块530，用于在修正基线漂移后的所述多个数据段中的数据段的所述多个波峰高度比值均属于预定阈值范围时，获得该数据段。
所述第四获取模块540具体包括：
波形质量获取模块541，用于对所述多个波峰高度比值均属于预定阈值范围的数据段的所述多个波峰高度比值取平均值作为该数据段的波形质量。
排序模块542，用于按照所述波形质量的数值大小进行排序，获得数值最大的波形质量所对应的数据段。
完整波形获取模块543，用于获得所述数值最大的波形质量所对应的数据段的每个完整波形。
所述第五获取模块550具体包括：
第一积分模块551，用于对所述完整波形的数值做积分，获得所述完整波形的平均值M，波峰H以及波谷L。
第一运算模块552，用于根据第一公式获得参数K。
第二运算模块553，用于根据第二公式、参数K，获得初选心输出量CO。
中值获取模块554，用于对所述多个初选心输出量取中值，获得心输出量。
本发明实施例提供的心输出量获得方法及装置获得时长为第二时间片段的多个数据段；对每个数据段修正基线漂移后，获得每个数据段所对应的心跳间隔的多个波峰，随后获得多个波峰相邻两个波峰中低波峰的波峰高度与高波峰的波峰高度的比值作为与数据段分别对应的多个波峰高度比值；获得波峰高度比值属于预定阈值范围的数据段，并在上述的数据段中获得多个波峰高度比值的平均值最大的数据段的每个完整波形；根据从每个完整波形获得的平均值、波峰以及波谷分别获得对应的多个初选心输出量，取多个初选心输出量的中值，便可以获得心输出量。本发明实施例提供的方法能改善现有技术中测量方法过于复杂，测量门槛较高的不足。
在本申请所提供的几个实施例中，应该理解到，所揭露的装置和方法，也可以通过其它的方式实现。以上所描述的装置实施例仅仅是示意性的，例如，附图中的流程图和框图显示了根据本发明的多个实施例的装置、方法和计算机程序产品的可能实现的体系架构、功能和操作。在这点上，流程图或框图中的每个方框可以代表一个模块、程序段或代码的一部分，所述模块、程序段或代码的一部分包含一个或多个用于实现规定的逻辑功能的可执行指令。也应当注意，在有些作为替换的实现方式中，方框中所标注的功能也可以以不同于附图中所标注的顺序发生。例如，两个连续的方框实际上可以基本并行地执行，它们有时也可以按相反的顺序执行，这依所涉及的功能而定。也要注意的是，框图和/或流程图中的每个方框、以及框图和/或流程图中的方框的组合，可以用执行规定的功能或动作的专用的基于硬件的系统来实现，或者可以用专用硬件与计算机指令的组合来实现。
另外，在本发明各个实施例中的各功能模块可以集成在一起形成一个独立的部分，也可以是各个模块单独存在，也可以两个或两个以上模块集成形成一个独立的部分。
所述功能如果以软件功能模块的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括：U盘、移动硬盘、只读存储器(ROM，Read-Only Memory)、随机存取存储器(RAM，Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。需要说明的是，在本文中，诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来，而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。
以上所述仅为本发明的优选实施例而已，并不用于限制本发明，对于本领域的技术人员来说，本发明可以有各种更改和变化。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。应注意到：相似的标号和字母在下面的附图中表示类似项，因此，一旦某一项在一个附图中被定义，则在随后的附图中不需要对其进行进一步定义和解释。
以上所述，仅为本发明的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应所述以权利要求的保护范围为准。

Claims (10)

1.一种心输出量获得方法，其特征在于，所述方法包括：
获得初选脉搏信号所对应的波形，每间隔第一时间片段，对所述波形进行时长为第二时间片段的截取以得到多个数据段，对所述多个数据段中的每个数据段修正基线漂移；
获得修正基线漂移后的每个数据段对应的心跳间隔的多个波峰，获得所述多个波峰的相邻两个波峰中的低波峰的波峰高度与高波峰的波峰高度的比值作为对应数据段的多个波峰高度比值；
当修正基线漂移后的所述多个数据段中的数据段的所述多个波峰高度比值均属于预定阈值范围时，获得该数据段；
获得所述多个波峰高度比值均属于预定阈值范围的多个数据段中的所述多个波峰高度比值的平均值最大的数据段的每个完整波形；
根据所述每个完整波形的平均值、波峰以及波谷获得所述每个完整波形分别对应的多个初选心输出量，对所述多个初选心输出量取中值，获得心输出量。
2.根据权利要求1所述的方法，其特征在于，所述获得所述多个波峰高度比值均属于预定阈值范围的多个数据段中的所述多个波峰高度比值的平均值最大的数据段的每个完整波形，包括：
对所述多个波峰高度比值均属于预定阈值范围的数据段的所述多个波峰高度比值取平均值作为该数据段的波形质量；
按照所述波形质量的数值大小进行排序，获得数值最大的波形质量所对应的数据段；
获得所述数值最大的波形质量所对应的数据段的每个完整波形。
3.根据权利要求2所述的方法，其特征在于，所述根据所述每个完整波形的平均值、波峰以及波谷获得所述每个完整波形分别对应的多个初选心输出量，包括：
对所述完整波形的数值做积分，获得所述完整波形的平均值M，波峰H以及波谷L；
根据公式获得参数K；
根据公式获得初选心输出量CO。
4.根据权利要求1所述的方法，其特征在于，所述获得初选脉搏信号，包括：
获得多个连续的图像，所述图像包括反映人体脉搏跳动的区域的图像；
获得所述多个连续的图像中的每一个图像的像素的红色通道数值，对所述多个连续的图像中的每一个图像的像素的红色通道数值分别进行数值累加，得到初选脉搏信号。
5.根据权利要求4所述的方法，其特征在于，所述对所述多个数据段中的每个数据段修正基线漂移以及获得修正基线漂移后的每个数据段对应的心跳间隔的多个波峰，包括：
对所述初选脉搏信号进行滤波，获得滤波后的脉搏信号；
根据所述滤波后的脉搏信号，获得所述滤波后的脉搏信号的波峰以及波谷，并根据所述波峰以及波谷获得归一化波峰；
去除相邻两个所述滤波后的脉搏信号的波峰的时间差值超过第一预定阈值范围的时间异常波峰以及所述归一化波峰的波峰高度超过第二预定阈值范围的高度异常波峰，得到筛选后的波峰；
获得所述筛选后的波峰中相邻的两个波峰的时间间隔，所述相邻的两个波峰的时间间隔为所述心跳间隔。
6.根据权利要求5所述的方法，其特征在于，所述对所述多个数据段中的每个数据段修正基线漂移，包括：
通过对所述多个数据段中的每个数据段做小波变换的方式以修正基线漂移。
7.根据权利要求1所述的方法，其特征在于，所述多个比值均属于预定阈值范围，包括：
所述多个比值均大于0.9且所述多个比值均小于或等于1。
8.一种心输出量获得装置，其特征在于，所述装置包括：
第一获取模块，用于获得初选脉搏信号所对应的波形，每间隔第一时间片段，对所述波形进行时长为第二时间片段的截取以得到多个数据段，对所述多个数据段中的每个数据段修正基线漂移；
第二获取模块，用于获得修正基线漂移后的每个数据段对应的心跳间隔的多个波峰，获得所述多个波峰的相邻两个波峰中的低波峰的波峰高度与高波峰的波峰高度的比值作为对应数据段的多个波峰高度比值；
第三获取模块，用于在修正基线漂移后的所述多个数据段中的数据段的所述多个波峰高度比值均属于预定阈值范围时，获得该数据段；
第四获取模块，用于获得所述多个波峰高度比值均属于预定阈值范围的多个数据段中的所述多个波峰高度比值的平均值最大的数据段的每个完整波形；
第五获取模块，用于根据所述每个完整波形的平均值、波峰以及波谷获得所述每个完整波形分别对应的多个初选心输出量，对所述多个初选心输出量取中值，获得心输出量。
9.根据权利要求8所述的装置，其特征在于，所述第四获取模块还包括：
波形质量获取模块，用于对所述多个波峰高度比值均属于预定阈值范围的数据段的所述多个波峰高度比值取平均值作为该数据段的波形质量；
排序模块，用于按照所述波形质量的数值大小进行排序，获得数值最大的波形质量所对应的数据段；
完整波形获取模块，用于获得所述数值最大的波形质量所对应的数据段的每个完整波形。
10.根据权利要求9所述的装置，其特征在于，所述第五获取模块还包括：
第一积分模块，用于对所述完整波形的数值做积分，获得所述完整波形的平均值M，波峰H以及波谷L；
第一运算模块，用于根据公式获得参数K；
第二运算模块，用于根据公式获得初选心输出量CO。
CN201610061495.8A 2016-01-29 2016-01-29 心输出量获得方法及装置 Active CN105433936B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610061495.8A CN105433936B (zh) 2016-01-29 2016-01-29 心输出量获得方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610061495.8A CN105433936B (zh) 2016-01-29 2016-01-29 心输出量获得方法及装置

Publications (2)

Publication Number Publication Date
CN105433936A CN105433936A (zh) 2016-03-30
CN105433936B true CN105433936B (zh) 2018-06-05

Family

ID=55544902

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610061495.8A Active CN105433936B (zh) 2016-01-29 2016-01-29 心输出量获得方法及装置

Country Status (1)

Country Link
CN (1) CN105433936B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106175742A (zh) * 2016-07-19 2016-12-07 北京心量科技有限公司 一种心脏体征获取方法以及装置
CN110840430B (zh) * 2018-08-21 2022-09-13 北京万生人和科技有限公司 腹内压数据筛选方法、计算机可读存储介质、腹内压数据筛选装置
CN113093194B (zh) * 2019-12-23 2023-12-29 佛山市顺德区美的饮水机制造有限公司 超声波波峰的标注方法与装置、存储介质、检测方法
CN111581599B (zh) * 2020-04-29 2023-10-03 四川虹美智能科技有限公司 重量值输出方法和数字称重变送器

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN100344257C (zh) * 2004-06-17 2007-10-24 肖行贯 心血管动力学参数的检测方法
CN101006919A (zh) * 2007-01-26 2007-08-01 北京工业大学 一种用于大压差下的心输出量检测方法及其装置
CN101176663B (zh) * 2007-12-06 2010-10-27 山东大学 基于脉搏波的无创心输出量检测方法和装置
GB201119982D0 (en) * 2011-11-18 2012-01-04 Maarek Albert Medical device and measurement of cardiac output
CN102499669B (zh) * 2011-10-26 2014-12-24 中国科学院深圳先进技术研究院 心脏参数的测量装置
CN104244814B (zh) * 2012-05-15 2017-08-11 皇家飞利浦有限公司 心输出量的监测
CN104323768B (zh) * 2012-11-20 2017-03-29 深圳市理邦精密仪器股份有限公司 一种心输出量连续监测的参数校准方法

Also Published As

Publication number Publication date
CN105433936A (zh) 2016-03-30

Similar Documents

Publication Publication Date Title
CN105433936B (zh) 心输出量获得方法及装置
CN105266772B (zh) 一种生理参数的测量方法
CN105550530B (zh) 精神散乱程度获得装置
CN103561661B (zh) 超声图像处理装置和程序
Chittaro et al. Data mining on temporal data: a visual approach and its clinical application to hemodialysis
CN109480808A (zh) 一种基于ppg的心率检测方法、系统、设备和存储介质
CN105433917A (zh) 心跳间隔获得方法及装置
CN106236049A (zh) 基于视频图像的血压测量方法
Karperien et al. Automated detection of proliferative retinopathy in clinical practice
CN110352035A (zh) 经由信号变化放大的静脉穿刺和动脉线引导
CN110522463B (zh) 一种基于脑功能连接分析的抑郁症辅助诊断系统
CN107550488A (zh) 心电伪迹消除方法及装置
CN110464333A (zh) 一种心电图数据的存储方法及装置
JP2016127923A (ja) 被験者の映像から生理学的パラメータを抽出するための対象領域を選択する方法
CN108937957A (zh) 检测方法、装置及检测设备
CN105025786A (zh) 对从电阻抗断层成像所获得信息进行简化的方法和设备
CN106491114A (zh) 一种心率检测方法及装置
CN106667439B (zh) 一种心电信号的降噪方法及装置
CN106974637A (zh) 心率实时计算方法及装置
CN109350031A (zh) 一种光电容积脉搏波信号特征点检测方法及装置
CN105595991B (zh) 一种心脏相干指数的测量方法及装置
Dumoulin et al. Contour extracting networks in early extrastriate cortex
CN113712526B (zh) 一种脉搏波提取方法、装置、电子设备及存储介质
CN106485287A (zh) 扫描区域的获取方法、数字切片的获取方法及装置
CN106096544B (zh) 基于二阶盲辨识的非接触式眨眼与心率联合检测系统及方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
TR01 Transfer of patent right

Effective date of registration: 20190102

Address after: 610000 1, 3, 1, 366 north section of lakeside road, Tianfu New District, Chengdu, Sichuan

Patentee after: Rent Seeking (Chengdu) Technology Co.,Ltd.

Address before: 100000 Beijing Haidian District Qinghe Anningzhuang East Road, No. 18, 20, 4-storey 4051

Patentee before: BEIJING XINLIANG TECHNOLOGY Co.,Ltd.

TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20231225

Address after: 318000 Medical Equipment Town Incubation Park, Xianju County Economic Development Zone, Taizhou City, Zhejiang Province (No. 8 Yong'an Avenue, Baita Block, Xianju Economic Development Zone, Baita Town, Xianju County), Zone A, West Side of Building 1, 2nd Floor

Patentee after: Zhejiang Maomao Technology Co.,Ltd.

Address before: 610000 1, 3, 1, 366 north section of lakeside road, Tianfu New District, Chengdu, Sichuan

Patentee before: Rent Seeking (Chengdu) Technology Co.,Ltd.

TR01 Transfer of patent right