CN113226174B - 一种心脏舒张功能评估方法、设备和系统 - Google Patents
一种心脏舒张功能评估方法、设备和系统 Download PDFInfo
- Publication number
- CN113226174B CN113226174B CN201980074575.9A CN201980074575A CN113226174B CN 113226174 B CN113226174 B CN 113226174B CN 201980074575 A CN201980074575 A CN 201980074575A CN 113226174 B CN113226174 B CN 113226174B
- Authority
- CN
- China
- Prior art keywords
- information
- subject
- diastolic function
- hemodynamic
- characteristic value
- 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
- 230000003205 diastolic effect Effects 0.000 title claims abstract description 79
- 238000011156 evaluation Methods 0.000 title claims abstract description 36
- 230000000004 hemodynamic effect Effects 0.000 claims abstract description 52
- 230000000747 cardiac effect Effects 0.000 claims abstract description 20
- 238000007781 pre-processing Methods 0.000 claims abstract description 9
- 230000006870 function Effects 0.000 claims description 73
- 238000000034 method Methods 0.000 claims description 36
- 230000002861 ventricular Effects 0.000 claims description 34
- 239000000835 fiber Substances 0.000 claims description 21
- 230000017531 blood circulation Effects 0.000 claims description 17
- 210000004369 blood Anatomy 0.000 claims description 16
- 239000008280 blood Substances 0.000 claims description 16
- 210000005240 left ventricle Anatomy 0.000 claims description 14
- 230000033001 locomotion Effects 0.000 claims description 14
- 210000001765 aortic valve Anatomy 0.000 claims description 13
- 210000000038 chest Anatomy 0.000 claims description 11
- 238000004590 computer program Methods 0.000 claims description 10
- 238000012545 processing Methods 0.000 claims description 10
- 238000001914 filtration Methods 0.000 claims description 8
- 210000003205 muscle Anatomy 0.000 claims description 8
- 230000003116 impacting effect Effects 0.000 claims description 4
- 210000000115 thoracic cavity Anatomy 0.000 claims description 4
- 238000000605 extraction Methods 0.000 claims description 3
- 238000009499 grossing Methods 0.000 claims description 3
- 230000000630 rising effect Effects 0.000 claims description 3
- 210000001991 scapula Anatomy 0.000 claims description 3
- 238000000718 qrs complex Methods 0.000 claims description 2
- 238000013210 evaluation model Methods 0.000 claims 2
- 238000010801 machine learning Methods 0.000 claims 2
- 238000012549 training Methods 0.000 claims 2
- 238000012935 Averaging Methods 0.000 claims 1
- 239000013307 optical fiber Substances 0.000 abstract description 18
- 238000012544 monitoring process Methods 0.000 abstract description 8
- 206010019280 Heart failures Diseases 0.000 description 19
- 238000010586 diagram Methods 0.000 description 14
- 230000008859 change Effects 0.000 description 8
- 238000005259 measurement Methods 0.000 description 7
- 230000006866 deterioration Effects 0.000 description 6
- 238000006073 displacement reaction Methods 0.000 description 6
- 230000001746 atrial effect Effects 0.000 description 5
- 238000004891 communication Methods 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 230000004872 arterial blood pressure Effects 0.000 description 3
- 210000004204 blood vessel Anatomy 0.000 description 3
- 230000000737 periodic effect Effects 0.000 description 3
- 230000029058 respiratory gaseous exchange Effects 0.000 description 3
- 230000001360 synchronised effect Effects 0.000 description 3
- 238000012360 testing method Methods 0.000 description 3
- 230000009466 transformation Effects 0.000 description 3
- 206010052337 Diastolic dysfunction Diseases 0.000 description 2
- 230000001133 acceleration Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 2
- 210000002376 aorta thoracic Anatomy 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 210000002837 heart atrium Anatomy 0.000 description 2
- 230000010365 information processing Effects 0.000 description 2
- 230000010354 integration Effects 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 210000001147 pulmonary artery Anatomy 0.000 description 2
- 230000010349 pulsation Effects 0.000 description 2
- 230000000241 respiratory effect Effects 0.000 description 2
- 230000035945 sensitivity Effects 0.000 description 2
- 206010007558 Cardiac failure chronic Diseases 0.000 description 1
- 201000001943 Tricuspid Valve Insufficiency Diseases 0.000 description 1
- 206010044640 Tricuspid valve incompetence Diseases 0.000 description 1
- 230000001154 acute effect Effects 0.000 description 1
- 206010000891 acute myocardial infarction Diseases 0.000 description 1
- 230000032683 aging Effects 0.000 description 1
- 238000003491 array Methods 0.000 description 1
- 238000010009 beating Methods 0.000 description 1
- 230000002457 bidirectional effect Effects 0.000 description 1
- 230000002146 bilateral effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000036772 blood pressure Effects 0.000 description 1
- 210000000748 cardiovascular system Anatomy 0.000 description 1
- 230000001684 chronic effect Effects 0.000 description 1
- 208000035850 clinical syndrome Diseases 0.000 description 1
- 230000006378 damage Effects 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000035487 diastolic blood pressure Effects 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 230000002526 effect on cardiovascular system Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 239000011521 glass Substances 0.000 description 1
- 230000004217 heart function Effects 0.000 description 1
- 210000003709 heart valve Anatomy 0.000 description 1
- 230000002427 irreversible effect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
- 238000005192 partition Methods 0.000 description 1
- 230000008506 pathogenesis Effects 0.000 description 1
- 230000002685 pulmonary effect Effects 0.000 description 1
- 210000003102 pulmonary valve Anatomy 0.000 description 1
- 230000004083 survival effect Effects 0.000 description 1
- 208000024891 symptom Diseases 0.000 description 1
- 230000035488 systolic blood pressure Effects 0.000 description 1
- 238000000844 transformation Methods 0.000 description 1
- 230000002792 vascular Effects 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
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/021—Measuring pressure in heart or blood vessels
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/103—Detecting, measuring or recording devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
- A61B5/11—Measuring movement of the entire body or parts thereof, e.g. head or hand tremor, mobility of a limb
-
- 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/0205—Simultaneously evaluating both cardiovascular conditions and different types of body conditions, e.g. heart and respiratory condition
-
- 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/0261—Measuring blood flow using optical means, e.g. infrared light
-
- 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/318—Heart-related electrical modalities, e.g. electrocardiography [ECG]
- A61B5/346—Analysis of electrocardiograms
- A61B5/349—Detecting specific parameters of the electrocardiograph cycle
- A61B5/366—Detecting abnormal QRS complex, e.g. widening
-
- 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
- A61B5/7239—Details of waveform analysis using differentiation including higher order derivatives
-
- 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
- A61B5/7264—Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems
- A61B5/7267—Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems involving training the classification device
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Cardiology (AREA)
- Physics & Mathematics (AREA)
- Public Health (AREA)
- Surgery (AREA)
- Veterinary Medicine (AREA)
- Physiology (AREA)
- General Health & Medical Sciences (AREA)
- Animal Behavior & Ethology (AREA)
- Molecular Biology (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Artificial Intelligence (AREA)
- Psychiatry (AREA)
- Signal Processing (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Fuzzy Systems (AREA)
- Evolutionary Computation (AREA)
- Mathematical Physics (AREA)
- Dentistry (AREA)
- Oral & Maxillofacial Surgery (AREA)
- Hematology (AREA)
- Vascular Medicine (AREA)
- Pulmonology (AREA)
- Measuring And Recording Apparatus For Diagnosis (AREA)
- Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
Abstract
一种心脏舒张功能评估方法，适用于心脏监测领域，包括：通过一个或多个光纤传感器获取对象胸腔体表的振动信息(101)；对振动信息进行预处理生成血流动力学相关信息(102)；对血流动力学相关信息进行高频分量提取生成高频分量信息(103)；基于血流动力学相关信息和高频分量信息确定第一特征值和第二特征值(104)；基于第一特征值和第二特征值生成指示参数，并基于指示参数评估对象的心脏舒张功能(105)。
Description
技术领域
本发明属于心脏监测领域，尤其涉及一种非侵入式心脏舒张功能评估方法、设备和系统。
背景技术
心力衰竭(简称心衰)是一个包含多种病因和发病机制的临床综合征。随着人口老龄化以及急性心肌梗死患者存活率的升高，慢性心衰患者的数量快速增长。心衰患者从慢性状态过度到急性的恶化状态时，伴随着心脏充盈压的升高。其中高充盈会使心脏功能进入快速的恶性循环，但是患者本身要在充盈压持续升高约20天左右，才感觉到症状从而需要紧急入院，这时心脏的损伤已经发生，且不可逆转。当识别到患者处于高充盈压状态时，需要及时干预，从而可以避免患者进一步恶化，这一点已经成为临床医生的共识。
目前有植入性产品用于评估心脏舒张功能，但价格较高，而且若只为监测所需，患者接受度较差。因此需要一种对人体更友好、使用方便的产品来评估心脏舒张功能。
发明内容
本发明的目的在于提供一种能评估测量对象的心脏舒张功能的评估方法、设备、系统和计算机可读存储介质，旨在实现非侵入式评估心脏舒张功能。
第一方面，本发明提供了一种心脏舒张功能评估方法，所述方法包括：
通过一个或多个光纤传感器获取对象胸腔体表的振动信息；其中，所述光纤传感器被配置为置于所述对象的身体下方；
对所述振动信息进行预处理生成血流动力学相关信息；
对所述血流动力学相关信息进行高频分量提取生成高频分量信息；
基于所述血流动力学相关信息和所述高频分量信息确定第一特征值和第二特征值；其中，第一特征值被认为表征了包括心室射血期主动脉瓣打开事件，第二特征值被认为表征了心室舒张早期血液流入左心室后冲击左室内壁的事件；
基于所述第一特征值和第二特征值生成指示参数，并基于所述指示参数评估所述对象的心脏舒张功能。
第二方面，本发明提供了一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，所述计算机程序被处理器执行时实现如上述的心脏舒张功能评估方法的步骤。
第三方面，本发明提供了一种心脏舒张功能评估设备，包括：一个或多个处理器；存储器；以及一个或多个计算机程序，其中所述一个或多个计算机程序被存储在所述存储器中，并且被配置成由所述一个或多个处理器执行，所述处理器执行所述计算机程序时实现如上述的心脏舒张功能评估方法的步骤。
第四方面，本发明提供了一种心脏舒张功能评估系统，所述系统包括：
一个或多个光纤传感器，被配置为放置于对象身体下方以获取对象的振动信息；和
与光纤传感器连接的，如上述的心脏舒张功能评估设备。
本发明是通过采集使用者的振动信息来监测心脏的舒张功能，无需侵入人体，被动测量，而且可以实现连续监测，使用者只需要躺在测量设备上即可进行测量，无需专业人员辅助，并且具有测量精度高、操作简单的优点，能提高测试者的舒适性，可以适用于医院和家庭等场景。本发明提供的心脏舒张功能评估系统，可以评估使用者的心脏舒张功能，进而在出现恶化迹象时提前预警，帮助使用者避免恶化后果。
附图说明
图1是依据本发明实施例一提供的心脏舒张功能评估方法的流程图；
图2是光纤传感器采集得到的对象A的振动信息的波形示意图；
图3是血流动力学相关信息的时域波形示意图；
图4是血流动力学相关信息和高频分量信息位于同一时间轴的时域波形示意图；
图5是一个心动周期内血流动力学相关信息和高频分量信息的第一波群和第二波群的示意图；
图6是使用心电信息、心音信息、振动能量信息曲线划分第一波群、第二波群的示意图；
图7是第一特征值的示意图；
图8是第二特征值的示意图；
图9A是基于对象B的振动信息的第一特征值和第二特征值示意图；
图9B是指示参数的ROC曲线图；
图10是依据本发明实施例三提供的心脏舒张功能评估设备的结构框图；
图11是依据本发明实施例四提供的心脏舒张功能评估系统的结构框图。
具体实施方式
为了使本发明的目的、技术方案及有益效果更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。
为了说明本发明所述的技术方案，下面通过具体实施例来进行说明。
实施例一：
请参阅图1，本发明实施例一提供的一种心脏舒张功能评估方法100包括以下步骤：需注意的是，若有实质上相同的结果，本发明的心脏舒张功能评估方法并不以图1所示的流程顺序为限。
S101、通过一个或多个光纤传感器获取对象胸腔体表的振动信息。
通过光纤传感器采集对象胸腔体表的振动信息时，光纤传感器可被放置于对象身体下方。例如对象可以呈仰卧、俯卧、侧卧等姿势，光纤传感器可放置于床上，对象仰卧(俯卧或侧卧)于其上。以对象呈仰卧姿势为例，较佳的测量位置是光纤传感器被配置为置于对象的背部下方，例如第一胸椎骨到第十二胸椎骨对应的背部体表区域之下，具体的可以是右肩肩胛骨下方。为描述方便，将对象的肩部分为左肩、中肩和右肩，其中左肩包括对象左肩胛骨对应的体表区域，右肩包括对象右肩胛骨对应的体表区域，左肩和右肩之间的区域是中肩。本领域普通人员可以理解的是，当对象呈俯卧姿势时，与仰卧姿势时的背部对应的测量位置对应的是对象的胸部。此外，还可以将光纤传感器放置于一定倾斜角仰卧人体背后的接触面、轮椅或其它可倚靠物体的倚卧人体背后的接触面等进行振动信息的采集。
光纤传感器至少一个，其感应区域面积至少是20平方厘米，例如其尺寸可以是5厘米*4厘米，此处的感应区域是指指光纤传感器中光纤分布的区域。当光纤传感器置于对象右肩肩胛骨下方时，其感应区域面积覆盖对象的右肩肩胛骨对应的体表区域。
图2所示为光纤传感器采集得到的对象A的振动信息的波形示意图。其中，曲线21的横轴表示时间，纵轴表示归一化处理后的振动信息，无量纲。光纤传感器采集的振动信息包含被测对象呼吸信号成分、血流动力学信号成分、以及环境微震动、被测对象体动引起的干扰和电路自身的噪声信号。此时的信号大轮廓即为人体呼吸产生的信号包络，而血流动力学信号与其它干扰噪声则叠加在呼吸包络曲线上。
S102、对振动信息进行预处理生成血流动力学相关信息。
不同的传感器获得的振动信息包含的信息量不同，有的包含的信息量比较丰富，因此需要对其进行预处理来捕获相关信号。例如，振动传感器采用光纤传感器时获得的振动信息中还包含被测对象的呼吸信号、体动信号、血流动力学信号、传感器固有的一些噪声等信号。
在本发明实施例一中，S102具体可以包括：
对振动信息进行滤波、去噪、信号缩放中的至少一种，得到血流动力学相关信息；具体可以为：根据对滤波后信号特征的需求采用IIR滤波器、FIR滤波器、小波滤波器、零相位双向滤波器、多项式拟合平滑滤波器、积分变换、微分变换中的一种或多种组合，对振动信息进行滤波去噪。例如对所述振动信息滤除2Hz以下的信息可以滤除呼吸信号和体动信号。预处理还可以包括：判断振动信息是否携带工频干扰信号，如果有，则通过工频陷波器滤除工频噪声。还可以对一些高频噪声(例如45Hz以上)进行去噪处理，处理后的信息可以根据情况进行信号缩放后可得到血流动力学相关信息。也可以直接设置滤波区间，例如滤波区间可以是1Hz-50Hz之间的区间。
如图3所示为对图2所示的光纤传感器获取的振动信息进行预处理后的血流动力学相关信息的时域波形示意图，曲线31滤波区间选择的是2Hz-45Hz。曲线31每个波形特征明显且一致性良好、周期规律、轮廓清晰、基线平稳，也即信号质量较佳。
S103、对血流动力学相关信息进行高频分量提取生成高频分量信息。
心脏的周期性搏动会导致各种变化的周期现象，例如心内压和心血管压、心房与心室的容积、心内瓣膜(包括二尖瓣、三尖瓣、主动脉瓣、肺动脉瓣)的启闭、血流速度等的周期变化。这些变化驱使血液在血管中沿着一定的方向流动。血流动力学(hemodynamics)研究的是血液在心血管系统中流动的力学，是以血液与血管的变形和流动为研究对象。本发明描述的“血流动力学相关信息”指任何与血流动力学相关的信息，可以包括但不限于，与血流产生相关的信息(例如心脏的收缩舒张导致射血)、与血流流动相关的信息(例如心排量CO(cardiac output)、左心室射血冲击主动脉弓)、与血流压力相关的信息(例如动脉收缩压、舒张压、平均动脉压)、与血管相关的信息(例如血管弹性)中的一种或几种。心脏的周期性搏动可以维持血液循环，因此与心脏搏动相关的各种参数，例如心内瓣膜的启闭、心房与心室容积的变化、心房与心室压力的变化、心房与心室中血流的流速和方向等，都是与血流动力学相关的信息。
光纤传感器采集的振动信息本质上对应的是位移变化信息，位移变化信息较为平滑，一些加速度或速度的变化细节在位移变化信息上较难识别。例如速度从零逐渐增加到某一峰值，再从该峰值逐渐递减到零，速度变化曲线形成先上升后下降的波形，但是位移变化曲线为单调波形。因此相较于与位移对应的信号分量而言，与速度、加速度对应的信号分量对应的峰谷时间宽度更窄，称之为高频分量信息。高频分量提取方法可以是多项式拟合平滑滤波的方式，还可以对血流动力学相关信息进行微分处理生成高频分量信息，例如，S103可以是对血流动力学相关信息进行二阶微分处理生成高频分量信息。
如图4所示，曲线42是高频分量信息的时域波形曲线。横轴表示时间，纵轴无量纲。选取图3所示的血流动力学相关信息即曲线31中的一个心动周期内的波形曲线，进行二阶微分处理后生成波形曲线42。为方便对照，将曲线31和曲线42置于同一时间轴进行同步显示。
S104、基于血流动力学相关信息和高频分量信息确定第一特征值和第二特征值。第一特征值用于表征包括心室射血期主动脉瓣打开事件。第二特征值用于表征心室舒张早期血液流入左心室后冲击左室内壁的事件。
在实施例一中，S104具体可以包括：
S1041、将血流动力学相关信息和高频分量信息置于同一时间轴上进行同步，划分心拍。
在一些例子中，振动信息是连续获取时，基于振动信息处理生成的血流动力学相关信息、高频分量信息、振动能量信息同样也是连续数据，此时需要对上述信息进行心拍划分。其中，心拍划分可以依据血流动力学相关信息或高频分量信息波形信号中重复出现的特征来进行划分。由于心脏活动具有明显的周期性，有一些明显的特征重复性高，例如，正常人的心动周期大约是0.6s至1秒之间，可以据此设定搜索区间，搜索最高峰，将最高峰作为心拍划分特征。类似的，最低谷也同样可以作为心拍划分特征。
在获取对象振动信息的同时，可以通过心电传感器获取对象的心电信息，由于ECG信号噪声小，信号干净，用来划分心拍准确度高，因此可以基于与振动信息同步获取的ECG信号来对血流动力学相关信息和高频分量信息进行心拍划分。
在另一些例子中，振动信息是以心动周期为单位离散获取时，并不需要进行心拍划分，S1041可以被省略。在本发明实施例一中，后续流程可以是对逐个心拍内的血流动力学相关信息的高频分量信息进行处理，还可以是将预设的一段时间内(例如30分钟)的血流动力学相关信息和高频分量信息按照心拍进行数据叠加平均处理后生成对应的平均信息，再对各个平均信息进行后续处理。因此下文描述的步骤中，血流动力学相关信息和高频分量信息可以是指一个心拍的数据，也可以是预设一段时间内的按心拍进行叠加平均处理后的数据。
S1042、在高频分量信息上划分第一波群和第二波群。
首先，基于血流动力学相关信息生成振动能量信息。具体的，对血流动力学相关信息逐点计算指定时间窗口的能量积分，生成振动能量信息。积分窗口的时间宽度可以是10ms，50ms，100ms等合适宽度，能量积分可以是取均值后的绝对值、平方、平方根等计算方式。图5中曲线51即是一个心动周期的振动能量信息的波形曲线。振动能量曲线有两个能量包络带，被认为表征了心脏收缩过程和心脏舒张过程产生的能量积聚.
其次，确定血流动力学相关信息曲线在同一个心动周期内的最高峰，如图5中点511所示。最高峰被认为表征了主动脉射血之后血液冲击主动脉弓所产生的震荡，将振动能量信息中包含点511的能量带确定为第一能量带，另一个确定为第二能量带。第一能量带对应的时间窗口确定为第一时间窗，第二能量带对应的时间窗口确定为第二时间窗。
最后，将高频分量信息位于第一时间窗内的波簇确定为第一波群，位于第二时间窗内的波簇确定为第二波群。如图5所示，是以一个心动周期为例的波群划分示意图，其中，曲线31、曲线42同图4所示，曲线51是振动能量信息的波形曲线，曲线42在方框区域502中的波簇即是第一波群，在方框区域504之内的波簇即是第二波群。
在一些例子中，波群划分还可以是基于心电信息进行辅助划分。在获取对象振动信息的同时，通过心电传感器同步获取对象的心电信息，并基于心电信息辅助划分第一波群和第二波群。例如，同一心动周期内的第一波群与心电信息的QRS波群时间几乎对应，第一波群通常处于心电信息的Q波之后80毫秒-200毫秒之间，可以用心电信息来辅助验证第一波群。
在一些例子中，波群划分还可以是基于心音信息进行划分。在获取对象振动信息的同时，通过心音传感器获取对象的心音信息，并基于心音信息划分第一波群和第二波群。例如，将与第一心音持续时间同步的波群作为第一波群，与第二心音持续时间同步的波群作为第二波群。
如图6所示，将同步获取的心电信息、心音信息、振动能量信息与图5的三条曲线置于同一时间轴同步，曲线61是心电信息示意图，曲线62是心音信息示意图，曲线42在方框区域602的波簇即是第一波群，在方框区域604之内的波簇即是第二波群。
S1043、在第一波群中搜索“W形”，并将其中点后下降沿幅度作为第一特征值。
“W形”的中点是W左右对称的交点。如图7所示，点711即为“W形”的中点。取点711后的下降沿幅度作为第一特征值，即点711与其后的第一个波谷点712间的幅度作为第一特征值L1。
S1044、在第二波群中搜索“W形”，并将其第二个波谷在其后预设时间段内的上升沿幅度作为第二特征值。
首先，在第二波群中搜索“W形”并确定其第二个波谷，第二个波谷被认为可以表征肺动脉瓣关闭这一事件。如图8所示，点902即为第二个波谷。其次，确定第二个波谷之后预设时间段内的最高峰，预设时间段可以是20毫秒至110毫秒之间任意值，还可以根据所指定的人群调整。如图8所示，点902之后预设时间段内的最高波峰标记为点904。最后计算第二个波谷和最高峰之间的幅度，即点902在预设时间段内的上升沿，作为第二特征值L2。
心室射血期主动脉瓣打开事件和心室舒张早期血液流入左心室后冲击左室内壁的事件可以通过不同的传感器来获取不同维度的信息，例如电生理传感器可以获取该事件的电信号，振动传感器可以获取该事件的振动信号。具体地，可以通过振动传感器获取对象的胸腔体表运动，进而从中提取出该对象的心室射血期主动脉瓣打开事件和心室舒张早期血液流入左心室后冲击左室内壁的事件。心室射血期主动脉瓣打开事件包括心室射血期主动脉瓣打开造成的肌肉和血流运动在体表形成的振动；心室舒张早期血液流入左心室后冲击左室内壁的事件包括心室舒张早期血液流入左心室后冲击左室内壁造成的肌肉和血流运动在体表形成的振动。本发明实施例一中，我们选取第一特征值用来表征心室射血期主动脉瓣打开造成的肌肉运动和血流运动在对象体表形成的振动幅度，选取第二特征值用来表征心室舒张早期血液流入左心室后冲击左室内壁造成的肌肉和血流运动在体表形成的振动幅度。可以理解的是，除了振动幅度我们还可以选取表征振动能量、振动频率、或振动时间等参数，用于表征心室射血期主动脉瓣打开事件和心室舒张早期血液流入左心室后冲击左室内壁的事件。
S105、基于第一特征值和第二特征值生成指示参数，并基于指示参数评估所述对象的心脏舒张功能。
在本发明实施例一中，可以将第二特征值与第一特征值的比值作为指示参数。当指示参数大于阈值时，判定对象心脏充盈压状态为高充盈压状态。心脏高充盈压状态是指满足三尖瓣反流流速>2.8m/s、E/e'>14、和E/A>1条件时的心脏状态。
图9A所示是依据对象B的胸腔体表振动信息的第一特征值和第二特征值的示意图。对象B是一位处于高充盈压状态的心衰患者。图9A中，曲线90是振动能量信息曲线，曲线91是血流动力学相关信息曲线，曲线92是高频分量信息曲线。L11是第一特征值，L22是第二特征值。同图7和图8相比，可见该心衰患者的第二特征值变化较大，这里选取第二特征值和第一特征值的比值来表征变化。
本领域普通技术人员可以据此得到将第一特征值与第二特征值的比值作为指示参数时心脏充盈压状态评估的方法，也包含在本发明的保护范围内。另外，本领域普通技术人员容易想到可以将第一特征值和第二特征值进行其他运算后生成指示参数，包括但不限于：加、减、乘、除、指数等操作，也均在本发明的保护范围之内。
选取25名心衰患者作为测试对象入组进行临床试验，25名心衰患者中包括12名高充盈压状态的患者(标记为阳性)，13名非高充盈压状态患者(标记为阴性)。按照上述的心脏舒张功能评估方法100计算25位测试对象的指示参数，并对25位测试对象的指示参数进行敏感性、特异性分析，构建ROC曲线如图9B所示，指示参数的准确率为86.5％，敏感性为76.9％，特异性91.7％。其中，阈值是0.4015，所述阈值是基于心衰人群确定的阈值。在本发明一些实施例中，阈值还可以是绝对阈值，用来区分正常人和患有心脏舒张功能障碍的人。阈值还可以是基于对象自身的阈值，例如，基于对监测对象的个人历史数据分析可以得到其心脏舒张功能恶化时的相对阈值。
在本发明实施例一中，心脏舒张功能是用心室充盈压状态来表征的，例如高充盈压状态表征了重度舒张功能障碍。此外，心脏舒张功能还可以用心房压来表征。心脏结构使得左心室充盈压与左房压、肺动脉压均具有相关性，因此在一些实施例中，指示参数除了可以用来评价充盈压状态，还可以经过一系列变换后用作间接评估左心房压力状态、肺动脉压力状态以及心衰程度等，同样在本发明的保护范围之内。
实施例二：
本发明实施例二提供了一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，所述计算机程序被处理器执行时实现如本发明实施例一提供的心脏舒张功能评估方法的步骤。
实施例三：
本发明实施例三提供了一种心脏舒张功能评估设备，图10所示是一种心脏舒张功能评估设备200的具体结构框图。该心脏舒张功能评估设备200可以是专门设计用于处理光纤传感器的振动信息的专用计算机设备。
例如，心脏舒张功能评估设备200可以包括连接到与其连接的网络的通信端口201，以便于数据通信。心脏舒张功能评估设备200还可以包括处理器203，处理器203以一个或多个处理器的形式，用于执行计算机指令。计算机指令可以包括例如执行本文描述的心脏舒张功能评估方法的例程、程序、对象、组件、数据结构、过程、模块和功能。例如，所述处理器203可以获得光纤传感器的振动信息，并对振动信息进行预处理生成血流动力学相关信息等。
在一些例子中，处理器203可以包括一个或多个硬件处理器，例如微控制器、微处理器、精简指令集计算机(RISC)、专用集成电路(ASIC)、图形处理单元(GPU)、中央处理单元(CPU)、数字信号处理器(DSP)、现场可编程门阵列(FPGA)、高级RISC机器(ARM)、可编程逻辑器件(PLD)等能够执行一个或多个功能的任何电路或处理器等，或其任何组合。
心脏舒张功能评估设备200可以包括内部通信总线205、存储器207用于由计算机处理和/或发送的各种数据文件，以及存储在存储器207中由处理器203执行的其他类型的非暂时性存储介质中的程序指令。本申请的方法和/或过程可以作为程序指令实现。心脏舒张功能评估设备200还包括输入/输出组件209，支持计算机和其他组件(例如，用户界面元件)之间的输入/输出。
应当理解的是，为了描述方便在本申请中心脏舒张功能评估设备200中仅描述了一个处理器。然而，应当注意，本申请中的心脏舒张功能评估设备200还可以包括多个处理器，因此，本申请中披露的操作和/或方法步骤可以如本申请所述的由一个处理器执行，也可以由多个处理器联合执行。例如，如果在本申请中心脏舒张功能评估设备200的处理器203执行步骤A和步骤B，则应该理解，步骤A和步骤B也可以由信息处理中的两个不同处理器联合或分开执行(例如，第一处理器执行步骤A，第二处理器执行步骤B，或者第一和第二处理器共同执行步骤A和B)。
实施例四：
本发明实施例四提供了一种心脏舒张功能评估系统，包括：
一个或多个光纤传感器；和
本发明实施例三提供的心脏舒张功能评估设备。
如图11所示是一种心脏舒张功能评估系统300的结构框图。一种心脏舒张功能评估系统300可以包括一个或多个光纤传感器301，一个或多个心脏舒张功能评估设备303，一个或多个存储装置305。
其中，光纤传感器301可被放置于对象身体下方。例如对象可以呈仰卧、俯卧、侧卧等姿势，光纤传感器可放置于床上，对象仰卧(俯卧或侧卧)于其上。以对象呈仰卧姿势为例，较佳的测量位置是光纤传感器被配置为置于对象的背部下方，较佳的测量状态是光纤传感器被配置为置于对象的左肩胛骨和右肩胛骨之间的区域。一般地，为描述方便，将对象左肩胛骨和右肩胛骨之间所对应的体表区域定义为中肩。本领域普通人员可以理解的是，当对象呈俯卧姿势时，与仰卧姿势时的背部对应的测量位置对应的是对象的胸部。此外，还可以振动传感器放置于一定倾斜角仰卧人体背后的接触面、轮椅或其它可倚靠物体的倚卧人体背后的接触面等进行振动信息的采集。
心脏舒张功能评估设备303如本发明实施例三所述。心脏舒张功能评估设备303可以通过网络320与光纤传感器301相连接。网络320可以是单一网络，例如有线网络或无线网络，还可以是多种网络的组合。网络320可以包括但不限于局域网、广域网、共用网络、专用网络等。网络320可以包括多种网络接入点，例如无线或有线接入点、基站或网络接入点，通过以上接入点使心脏舒张功能评估系统300的其他组成部分可以连接网络103并通过网络传送信息。
存储装置305可以被配置为存储数据和指令。存储装置305可以包括但不限于随机存储器、只读存储器、可编程只读存储器等。存储装置305可以是利用电能方式、磁能方式、光学方式等存储信息的设备，例如硬盘、软盘、磁芯存储器、CD、DVD等。以上提及的存储设备只是列举了一些例子，存储装置305使用的存储设备并不局限于此。
在一些例子中，该心脏舒张功能评估系统300还可以包括输出装置307，输出装置307被配置为将心脏舒张功能评估结果输出，输出方式包括但不限于图形、文字、数据、语音等，例如图形显示、数字显示、语音播报、盲文显示等中的一种或多种。输出装置307可以是显示器、手机、平板电脑、投影仪、可穿戴设备(手表、耳机、眼镜等)、盲文显示器等中的一种或多种。在一些例子中，输出装置307可以实时显示对象102的心脏舒张功能评估结果，在另一些例子中，输出装置307可以非实时显示一份报告，该报告是对象在预设时间段内的测量结果，例如用户在入睡时间段内的心脏充盈压监测结果。当监测对象是心衰患者时，如果心脏舒张功能评估设备将其心脏舒张功能评估为高充盈压状态，那心衰患者此时会面临心衰状况恶化，需要接受住院治疗，该监测系统的输出装置可以给该心衰患者发送提示信息，例如发送短信、电子邮件、电话、微信等即时聊天信息，还可以给该心衰患者的家庭医生发送信息，提示该患者有可能面临心衰恶化以帮助医生做决策。该系统还可以进一步包括医患交流平台，当医生接收到系统推送的患者可能面临心衰恶化时，及时与患者沟通。
再如，输出装置307还可以实现预警功能，例如通过语音预警，当心脏舒张功能评估设备将心衰患者心脏舒张功能评估为高充盈压状态，那心衰患者此时会面临心衰状况恶化，通过声音来提醒患者及时就诊。
本发明是通过采集使用者的振动信息来监测心脏的舒张功能，无需侵入人体，被动测量，而且可以实现连续监测，使用者只需要躺在测量设备上即可进行测量，无需专业人员辅助，并且具有测量精度高、操作简单的优点，能提高测试者的舒适性，可以适用于医院和家庭等场景。本发明提供的心脏舒张功能评估系统，可以评估使用者的心脏舒张功能状态，进而在出现恶化迹象时提前预警，帮助使用者避免恶化后果。
本领域普通技术人员可以理解上述实施例的各种方法中的全部或部分步骤是可以通过程序来指令相关的硬件来完成，该程序可以存储于一计算机可读存储介质中，存储介质可以包括：只读存储器(ROM，Read Only Memory)、随机存取记忆体(RAM，RandomAccess Memory)、磁盘或光盘等。
以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。
Claims (20)
1.一种心脏舒张功能评估方法，其特征在于，所述方法包括：
通过一个或多个光纤传感器获取对象胸腔体表的振动信息,其中，所述光纤传感器被配置为置于所述对象的身体下方；
对所述振动信息进行预处理生成血流动力学相关信息；
对所述血流动力学相关信息进行高频分量提取生成高频分量信息；
基于所述血流动力学相关信息和所述高频分量信息确定第一特征值和第二特征值，其中，第一特征值用于表征心室射血期主动脉瓣打开事件，第二特征值用于表征心室舒张早期血液流入左心室后冲击左室内壁的事件；
基于所述第一特征值和第二特征值生成指示参数，并基于所述指示参数评估所述对象的心脏舒张功能；
所述基于所述指示参数评估所述对象的心脏舒张功能包括评估所述对象的心脏充盈压状态；所述基于所述指示参数评估所述对象的心脏充盈压状态，具体是：
确定第二特征值与第一特征值的比值作为指示参数；
当所述指示参数大于阈值时，判定所述对象心脏为高充盈压状态。
2.如权利要求1所述的方法，其特征在于，所述阈值是0.4015。
3.如权利要求1所述的方法，其特征在于，所述阈值包括基于人群的阈值。
4. 如权利要求1所述的方法，其特征在于，所述血流动力学相关信息是：
一个心动周期内的数据；或
预设时间段内的以心动周期为单位进行叠加和平均后的数据。
5.如权利要求1所述的方法，其特征在于，所述光纤传感器被配置为置于所述对象的第一胸椎体到第十二胸椎体对应的体表区域下方。
6.如权利要求5所述的方法，其特征在于，所述光纤传感器被配置为置于所述对象的右肩肩胛骨下方。
7.如权利要求1所述的方法，其特征在于，所述光纤传感器的感应区域的面积至少是二十平方厘米。
8.如权利要求6所述的方法，其特征在于，所述光纤传感器的感应区域的面积覆盖所述对象的右肩肩胛骨对应的体表区域。
9.如权利要求1所述的方法，其特征在于，所述对象的身体姿势之一是仰卧。
10.如权利要求1所述的方法，其特征在于，所述预处理包括滤波、去噪、信号缩放中的至少一种。
11.如权利要求1所述的方法，其特征在于，所述心室射血期主动脉瓣打开事件是心室射血期主动脉瓣打开造成的肌肉和血流运动在体表形成的振动；所述心室舒张早期血液流入左心室后冲击左室内壁的事件是心室舒张早期血液流入左心室后冲击左室内壁造成的肌肉和血流运动在体表形成的振动。
12.如权利要求11所述的方法，其特征在于，所述第一特征值包括心室射血期主动脉瓣打开造成的肌肉和血流运动在体表形成的振动幅度，所述第二特征值包括心室舒张早期血液流入左心室后冲击左室内壁造成的肌肉和血流运动在体表形成的振动幅度。
13. 如权利要求12所述的方法，其特征在于，所述对所述血流动力学相关信息进行高频分量提取生成高频分量信息，具体是：
对所述血流动力学相关信息进行二阶微分处理生成高频分量信息；或者是
对所述血流动力学相关信息进行多项式拟合平滑滤波生成高频分量信息。
14. 如权利要求1所述的方法，其特征在于， 所述基于所述血流动力学相关信息和所述高频分量信息确定第一特征值和第二特征值，包括：
基于所述血流动力学相关信息生成振动能量信息，所述振动能量信息在一个心动周期内包含两个能量包络带；
将所述血流动力学相关信息和所述振动能量信息置于同一时间轴同步，确定所述血流动力学相关信息在同一个心动周期内的最高峰；
将所述振动能量信息两个能量包络带中包括所述血流动力学相关信息的最高峰的能量包络带的持续时间确定为第一时间窗，另一个能量包络带的持续时间确定为第二时间窗；
将所述高频分量信息位于第一时间窗内的波簇确定为第一波群，位于第二时间窗内的波簇确定为第二波群；
基于所述第一波群和第二波群确定第一特征值和第二特征值。
15.如权利要求14所述的方法，其特征在于，所述方法还包括：
通过心电传感器与所述振动信息同步获取所述对象的心电信息；
基于心电信息的QRS波群辅助确定所述高频分量信息的第一波群。
16.如权利要求14所述的方法，其特征在于，所述基于所述第一波群和第二波群确定第一特征值和第二特征值，包括：
在所述第一波群中搜索“W形”，并将其中点后下降沿幅度作为第一特征值；
在所述第二波群中搜索“W形”，并将其第二个波谷在其后预设时间段的上升沿幅度作为第二特征值。
17.一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，其特征在于，所述计算机程序被处理器执行时实现如权利要求1至16任一项所述的心脏舒张功能评估方法的步骤。
18.一种用于心脏舒张功能评估的设备，包括：
一个或多个处理器；
存储器；以及
一个或多个计算机程序，其中所述一个或多个计算机程序被存储在所述存储器中，并且被配置成由所述一个或多个处理器执行，其特征在于，所述处理器执行所述计算机程序时实现如权利要求1至16任一项所述的心脏舒张功能评估方法的步骤。
19. 一种用于心脏舒张功能评估的系统，其特征在于，所述系统包括：
一个或多个光纤传感器，被配置为放置于对象的身体下方以获取所述对象胸腔体表的振动信息；和
与光纤传感器连接的，如权利要求18所述的用于心脏舒张功能评估的设备。
20.一种基于机器学习的心脏舒张功能评估系统，其特征在于，所述系统包括：
一个或多个处理器，所述处理器被集中或各自编程为实现：
接收对象的胸腔体表振动信息作为训练输入信息；
通过机器学习对训练输入信息进行分析建立评估模型；
接收待评价对象的胸腔体表振动信息，所述评估模型对所述待评价对象的心脏舒张功能做出评估；
其中，评估模型执行以下操作：
对所述振动信息进行预处理生成血流动力学相关信息；
对所述血流动力学相关信息进行高频分量提取生成高频分量信息；
基于所述血流动力学相关信息和所述高频分量信息确定第一特征值和第二特征值，其中，第一特征值用于表征心室射血期主动脉瓣打开事件，第二特征值用于表征心室舒张早期血液流入左心室后冲击左室内壁的事件；
基于所述第一特征值和第二特征值生成指示参数，并基于所述指示参数评估所述对象的心脏舒张功能；
所述基于所述指示参数评估所述对象的心脏舒张功能包括评估所述对象的心脏充盈压状态；
所述基于所述指示参数评估所述对象的心脏充盈压状态，具体是：
确定第二特征值与第一特征值的比值作为指示参数；
当所述指示参数大于阈值时，判定所述对象心脏为高充盈压状态。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410568666.0A CN118383756A (zh) | 2019-05-20 | 2019-05-20 | 一种心脏舒张功能评估方法、设备和系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| PCT/CN2019/087640 WO2020232606A1 (zh) | 2019-05-20 | 2019-05-20 | 一种心脏舒张功能评估方法、设备和系统 |
Related Child Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202410568666.0A Division CN118383756A (zh) | 2019-05-20 | 2019-05-20 | 一种心脏舒张功能评估方法、设备和系统 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113226174A CN113226174A (zh) | 2021-08-06 |
| CN113226174B true CN113226174B (zh) | 2024-06-07 |
Family
ID=73459530
Family Applications (2)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201980074575.9A Active CN113226174B (zh) | 2019-05-20 | 2019-05-20 | 一种心脏舒张功能评估方法、设备和系统 |
| CN202410568666.0A Pending CN118383756A (zh) | 2019-05-20 | 2019-05-20 | 一种心脏舒张功能评估方法、设备和系统 |
Family Applications After (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202410568666.0A Pending CN118383756A (zh) | 2019-05-20 | 2019-05-20 | 一种心脏舒张功能评估方法、设备和系统 |
Country Status (3)
| Country | Link |
|---|---|
| US (1) | US20220218208A1 (zh) |
| CN (2) | CN113226174B (zh) |
| WO (1) | WO2020232606A1 (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2020232607A1 (zh) * | 2019-05-20 | 2020-11-26 | 深圳市大耳马科技有限公司 | 一种心脏舒张功能评估方法、设备和系统 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5743267A (en) * | 1995-10-19 | 1998-04-28 | Telecom Medical, Inc. | System and method to monitor the heart of a patient |
| CN108056769A (zh) * | 2017-11-14 | 2018-05-22 | 深圳市大耳马科技有限公司 | 一种生命体征信号分析处理方法、装置和生命体征监测设备 |
| CN109310371A (zh) * | 2016-06-16 | 2019-02-05 | 阿克瑞克公司 | 定量心震描记法 |
Family Cites Families (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7846104B2 (en) * | 2007-02-08 | 2010-12-07 | Heart Force Medical Inc. | Monitoring physiological condition and detecting abnormalities |
| US8428698B2 (en) * | 2009-03-04 | 2013-04-23 | Pacesetter, Inc. | Systems and methods for monitoring DP, IVRT, DiFT, diastolic function and/or HF |
| US20140288442A1 (en) * | 2010-02-25 | 2014-09-25 | Tonino Bombardini | Method and apparatus for quantification and monitoring of cardiovascular function during induced stress or physical activity and at rest |
| CA2834102A1 (en) * | 2011-05-03 | 2012-11-08 | Heart Force Medical Inc. | Method and apparatus for estimating myocardial contractility using precordial vibration signals |
| WO2014015300A1 (en) * | 2012-07-19 | 2014-01-23 | Regents Of The University Of Minnesota | Cardiac assist device with pulse wave analysis |
| US20180153420A1 (en) * | 2013-09-30 | 2018-06-07 | Elfi-Tech Ltd. | Apparatus and method for optical measurement of cardiovascular fitness, stress and physiological parameters |
| JP2017535316A (ja) * | 2014-09-30 | 2017-11-30 | 深▲せん▼市大耳馬科技有限公司Ｓｈｅｎｚｈｅｎ Ｄａｒｍａ Ｔｅｃｈｎｏｌｏｇｙ Ｃｏ．，Ｌｔｄ． | 姿勢と生命徴候用モニターリングシステム及び方法 |
| US20170188978A1 (en) * | 2016-01-04 | 2017-07-06 | AventuSoft, LLC | System and method of measuring hemodynamic parameters from the heart valve signal |
| US11241170B2 (en) * | 2017-04-07 | 2022-02-08 | University Of Maryland, College Park | Monitor for blood pressure and other arterial properties |
-
2019
- 2019-05-20 CN CN201980074575.9A patent/CN113226174B/zh active Active
- 2019-05-20 US US17/613,475 patent/US20220218208A1/en active Pending
- 2019-05-20 WO PCT/CN2019/087640 patent/WO2020232606A1/zh active Application Filing
- 2019-05-20 CN CN202410568666.0A patent/CN118383756A/zh active Pending
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5743267A (en) * | 1995-10-19 | 1998-04-28 | Telecom Medical, Inc. | System and method to monitor the heart of a patient |
| CN109310371A (zh) * | 2016-06-16 | 2019-02-05 | 阿克瑞克公司 | 定量心震描记法 |
| CN108056769A (zh) * | 2017-11-14 | 2018-05-22 | 深圳市大耳马科技有限公司 | 一种生命体征信号分析处理方法、装置和生命体征监测设备 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN118383756A (zh) | 2024-07-26 |
| US20220218208A1 (en) | 2022-07-14 |
| WO2020232606A1 (zh) | 2020-11-26 |
| CN113226174A (zh) | 2021-08-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US9451921B2 (en) | System and method for indicating coronary artery disease risk based on low and high frequency bands | |
| JP6131404B2 (ja) | 心臓の機能不全と異常を表す情報を測定するための方法及び器械 | |
| EP2840955B1 (en) | Method and apparatus for determining information indicative of cardiac malfunctions and abnormalities | |
| EP3010407B1 (en) | Method and apparatus for determining information indicative of cardiac malfunctions | |
| US20170071564A1 (en) | Heart rate detection method and device using heart sound acquired from ausculation positions | |
| TWI595860B (zh) | 心震圖譜特徵點量測方法 | |
| AU2019350718A1 (en) | Model-based sensor technology for detection of cardiovascular status | |
| CA3122115C (en) | Method of predicting fluid responsiveness in patients | |
| CN113226170B (zh) | 一种心脏舒张功能评估方法、设备和系统 | |
| RU2346653C2 (ru) | Способ и система для анализа сердечно-сосудистых звуков | |
| CN113226174B (zh) | 一种心脏舒张功能评估方法、设备和系统 | |
| TW201739415A (zh) | 心臟搏動生理之量測系統 | |
| Granados et al. | Acoustic heart. Interpretation of phonocardiograms by computer | |
| CN113164124B (zh) | 一种心脏舒张功能评估方法、设备和系统 | |
| US20220304631A1 (en) | Multisensor pulmonary artery and capillary pressure monitoring system | |
| Forouzanfar et al. | Model-based oscillometric blood pressure estimation | |
| Munck et al. | Body surface mapping of the mechanical cardiac activity | |
| Cosenza | Design of an Inertial Multi-Sensor Network for the Monitoring of the Heart Rate During Sleep using Ballistocardiographic Signals |
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 | ||
| TA01 | Transfer of patent application right | ||
| TA01 | Transfer of patent application right |
Effective date of registration: 20240325 Address after: 518000, 5th Floor, Building B, Baoshi Science and Technology Park, Baoshi Road, Tangtou Community, Shiyan Street, Bao'an District, Shenzhen City, Guangdong Province Applicant after: Oufeikang Technology (Shenzhen) Co.,Ltd. Country or region after: China Address before: 518052 501, Shangmei technology building, 2009 Shahe West Road, Yuehai street, Nanshan District, Shenzhen City, Guangdong Province Applicant before: SHENZHEN DARMA TECHNOLOGY Co.,Ltd. Country or region before: China |
|
| GR01 | Patent grant |