RU2013115110A - Способ, устройство и программа для автоматической обработки сигналов кровяного давления - Google Patents
Способ, устройство и программа для автоматической обработки сигналов кровяного давления Download PDFInfo
- Publication number
- RU2013115110A RU2013115110A RU2013115110/14A RU2013115110A RU2013115110A RU 2013115110 A RU2013115110 A RU 2013115110A RU 2013115110/14 A RU2013115110/14 A RU 2013115110/14A RU 2013115110 A RU2013115110 A RU 2013115110A RU 2013115110 A RU2013115110 A RU 2013115110A
- Authority
- RU
- Russia
- Prior art keywords
- pressure
- pressure signal
- point
- derivative
- points
- Prior art date
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/021—Measuring pressure in heart or blood vessels
- A61B5/02108—Measuring pressure in heart or blood vessels from analysis of pulse wave characteristics
- A61B5/02125—Measuring pressure in heart or blood vessels from analysis of pulse wave characteristics of pulse wave propagation time
-
- 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
- A61B5/02108—Measuring pressure in heart or blood vessels from analysis of pulse wave characteristics
-
- 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
- A61B5/0215—Measuring pressure in heart or blood vessels by means inserted into the body
-
- 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/725—Details of waveform analysis using specific filters therefor, e.g. Kalman or adaptive filters
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/74—Details of notification to user or communication with user or patient ; user input means
- A61B5/742—Details of notification to user or communication with user or patient ; user input means using visual displays
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Cardiology (AREA)
- Molecular Biology (AREA)
- General Health & Medical Sciences (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Veterinary Medicine (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Public Health (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Physics & Mathematics (AREA)
- Physiology (AREA)
- Vascular Medicine (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Psychiatry (AREA)
- Signal Processing (AREA)
- Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
- External Artificial Organs (AREA)
- Investigating Or Analysing Biological Materials (AREA)
Abstract
1. Способ автоматической обработки сигналов кровяного давления, содержащий по меньшей мере следующие этапы, на которых:А. дискретизируют обнаруженный сигнал давления P(t) для одного или более сердечных сокращений, причем каждое сердечное сокращение начинается в начальный момент, совпадающий с моментом диастолического давления, и оканчивается в последний момент, совпадающий с моментом следующего диастолического давления, и содержит дикротическую точку, причем каждое сердечное сокращение имеет систолическую фазу, заключающуюся между начальной диастолической точкой и дикротической точкой; иВ. автоматически анализируют и выделяют морфологию дискретизированного сигнала давления P(t) для каждого сердечного сокращения, определяют момент и значение давления в одной или более характеристических точках сигнала P(t), выбранных из группы, содержащей: точку начального диастолического кровяного давления; точку систолического давления; дикротическую точку; и одну или более резонансных точек, каждая из которых возникает, когда вторая производная dP/dtот сигнала давления P(t) имеет локальный максимум, в котором по меньшей мере характеристическая точка сигнала давления P(t) принадлежит к систолической фазе сердечного сокращения и отличается от точки начального диастолического давления;причем способ отличается тем, что он дополнительно содержит следующий этап, на котором:С. для каждого сердечного сокращения определяют значение энергетической эффективности посредством следующих подэтапов, на которых:C.1 определяют импеданс Z(t) прямой динамической волны давления для каждой из упомянутых одной или более характеристически�
Claims (14)
1. Способ автоматической обработки сигналов кровяного давления, содержащий по меньшей мере следующие этапы, на которых:
А. дискретизируют обнаруженный сигнал давления P(t) для одного или более сердечных сокращений, причем каждое сердечное сокращение начинается в начальный момент, совпадающий с моментом диастолического давления, и оканчивается в последний момент, совпадающий с моментом следующего диастолического давления, и содержит дикротическую точку, причем каждое сердечное сокращение имеет систолическую фазу, заключающуюся между начальной диастолической точкой и дикротической точкой; и
В. автоматически анализируют и выделяют морфологию дискретизированного сигнала давления P(t) для каждого сердечного сокращения, определяют момент и значение давления в одной или более характеристических точках сигнала P(t), выбранных из группы, содержащей: точку начального диастолического кровяного давления; точку систолического давления; дикротическую точку; и одну или более резонансных точек, каждая из которых возникает, когда вторая производная d2P/dt2 от сигнала давления P(t) имеет локальный максимум, в котором по меньшей мере характеристическая точка сигнала давления P(t) принадлежит к систолической фазе сердечного сокращения и отличается от точки начального диастолического давления;
причем способ отличается тем, что он дополнительно содержит следующий этап, на котором:
С. для каждого сердечного сокращения определяют значение энергетической эффективности посредством следующих подэтапов, на которых:
C.1 определяют импеданс Zd-D(t) прямой динамической волны давления для каждой из упомянутых одной или более характеристических точек, принадлежащих систолической фазе рассматриваемого сердечного сокращения, за исключением точки начального диастолического давления, причем упомянутый импеданс Zd-D(t) прямой волны давления задается соотношением между значением сигнала давления P(t) в характеристической точке и расстоянием момента от начального момента рассматриваемого сердечного сокращения, и определяют импеданс zd прямой волны давления путем сложения с чередующимися знаками значений динамических прямых импедансов Zd-D(t), упорядоченных согласно прямому временному порядку, начиная с начального момента рассматриваемого сердечного сокращения, до дикротического момента, применяют положительный знак к первому динамическому прямому импедансу Zd-D(t) согласно прямому временному порядку;
С.2 определяют для каждой из упомянутых одной или более характеристических точек динамический отраженный импеданс Zd_R(t), который задается отношением между давлением в характеристических точках и расстоянием соответствующего момента от момента последнего сокращения, и определяют значение импеданса ZR отраженных волн давления, получая его путем сложения с чередующимися знаками динамических импедансов второго набора точек, определенных таким образом, отсортированных в соответствии с их обратным временным порядком, начиная с момента последнего сокращения до момента начального диастолического давления, применяя положительный знак к первому динамическому импедансу Zd_R(t)в соответствии с обратным временным порядком;
С.3 определяют упомянутую энергетическую эффективность как соотношение между импедансом zd давления прямой волны и импедансом ZR отраженных волн: RES=ZD/ZR.
2. Способ по п.1, отличающийся тем, что дополнительно содержит следующие этапы, на которых:
D. проверяют в отношении упомянутой энергетической эффективности, определенной на этапе С, является ли первая производная dP/dt сигнала давления P(t) более низкой, чем первое максимальное пороговое значение Td на всем протяжении рассматриваемого сердечного сокращения, и является ли вторая производная d2P/dt2 более низкой, чем второе максимальное пороговое значение Td2 на всем протяжении рассматриваемого сердечного сокращения, и в отрицательном случае выполняют этап Е, в противном случае выполняют этап F;
Е. выбирают частоту отсечки фильтра нижних частот на основе энергетической эффективности, определенной на этапе С, из первой производной dP/dt и второй производной d2P/dt2, и применяют упомянутый фильтр нижних частот к сигналу давления P(t) для того, чтобы получить новый дискретизированный сигнал давления и выполнить предыдущие этапы, начиная с этапа В;
F. обеспечивают в качестве выходного сигнала сигнал P(t), при котором этап В был выполнен в последний раз.
3. Способ по п.1, отличающийся тем, что одну или более упомянутых резонансных точек определяют на этапе В посредством следующих подэтапов, на которых:
В.2 определяют общее число Ndp-max точек локального максимума первой производной dP/dt сигнала давления (дискретизированного) в диапазоне одного сердечного сокращения;
В.3 определяют точки относительного максимума второй производной d2P/dt2 сигнала давления в диапазоне одного сердечного сокращения, и
В.4 определяют Ndp-max точек относительного максимума второй производной d2P/dt2, имеющих самые высокие значения, и определяют Ndp-max моментов времени td2p_max (i) в которых они возникают, рассматривая точки сигнала давления в упомянутые Ndp-max моментов времени td2p-max (i) как резонансные точки.
4. Способ по п.1 или 3, отличающийся тем, что на этапе В определяют следующие характеристические точки сигнала давления P(t):
- точку начального диастолического давления,
- точку систолического давления,
- дикротическую точку, и
- одну или более резонансных точек.
5. Способ по п.2, отличающийся тем, что упомянутые первое максимальное пороговое значение Td и второе максимальное пороговое значение Td2 зависят от энергетической эффективности, определенной на этапе С.
6. Способ по п.2, отличающийся тем, что на упомянутом этапе D проверяют, принадлежит ли упомянутая энергетическая эффективность, определенная на этапе С, к одному из трех или более, предпочтительно четырех, смежных интервалов изменчивости, при этом первое максимальное пороговое значение Td и второе максимальное пороговое значение Td2 предпочтительно зависят от интервала, к которому принадлежит энергетическая эффективность, определенная на этапе С.
7. Способ по п.6, отличающийся тем, что на упомянутом этапе Е частоту отсечки выбирают следующим образом:
- значения RES выделяют в три или более, предпочтительно в четыре, смежных интервала изменчивости;
- для каждого из упомянутых трех или более смежных интервалов изменчивости энергетической эффективности, определенной на этапе С, значения первой производной dP/dt сигнала давления на протяжении всего сердечного сокращения выделяют в три или более, предпочтительно шесть, смежных интервалов изменчивости;
- и для каждого из упомянутых интервалов изменчивости первой производной dP/dt сигнала давления P(t) на протяжении всего сердечного сокращения значения второй производной d2P/dt2 сигнала давления выделяют в три или более, предпочтительно четыре, неперекрывающихся интервала изменчивости, которым соответствует соответствующее значение упомянутой частоты отсечки.
8. Способ по п.2, отличающийся тем, что упомянутая частота отсечки возрастает по мере того, как возрастает первая производная dP/dt сигнала давления P(t), причем энергетическая эффективность и вторая производная d2P/dt2 сигнала давления P(t) постоянны.
9. Способ по п.2, отличающийся тем, что упомянутая частота отсечки уменьшается по мере того, как возрастает вторая производная d2P/dt2 сигнала давления P(t), причем энергетическая эффективность и первая производная dP/dt сигнала давления P(t) постоянны.
10. Способ по п.2, отличающийся тем, что упомянутая частота отсечки изменяется в промежутке от 0,5 Гц до 10 Гц, предпочтительно в промежутке от 2 Гц до 80 Гц, более предпочтительно в промежутке от 3 Гц до 60 Гц.
11. Способ по п.2, отличающийся тем, что на этапе F сигнал давления P(t) отображают на дисплее.
12. Способ по п.2, отличающийся тем, что на этапе С RES отображают на дисплее.
13. Автоматическое устройство для обработки сигнала кровяного давления, отличающееся тем, что оно содержит средство обработки для выполнения этапов способа по любому из пп.1-12.
14. Запоминающий носитель, считываемый компьютерными средствами, на котором сохранена компьютерная программа, содержащая кодовые средства для выполнения этапов способа по любому из пп.1-12 при работе в соединении со средством обработки устройства.
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
ITRM2010A000468 | 2010-09-06 | ||
ITRM2010A000468A IT1402427B1 (it) | 2010-09-06 | 2010-09-06 | Metodo automatico di misura ed elaborazione della pressione sanguigna. |
PCT/IB2011/002025 WO2012032386A1 (en) | 2010-09-06 | 2011-09-02 | Method, apparatus and program for the automatic processing of blood pressure signals |
Publications (2)
Publication Number | Publication Date |
---|---|
RU2013115110A true RU2013115110A (ru) | 2014-10-20 |
RU2552685C2 RU2552685C2 (ru) | 2015-06-10 |
Family
ID=43739087
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
RU2013115110/14A RU2552685C2 (ru) | 2010-09-06 | 2011-09-02 | Способ, устройство и программа для автоматической обработки сигналов кровяного давления |
Country Status (12)
Country | Link |
---|---|
US (3) | US20130172761A1 (ru) |
EP (2) | EP2613691B1 (ru) |
CN (2) | CN103108585B (ru) |
AU (1) | AU2011300338B2 (ru) |
BR (2) | BR112013005065B1 (ru) |
CA (1) | CA2809930C (ru) |
ES (2) | ES2578998T3 (ru) |
HK (1) | HK1187225A1 (ru) |
IT (1) | IT1402427B1 (ru) |
PL (2) | PL2613691T3 (ru) |
RU (1) | RU2552685C2 (ru) |
WO (2) | WO2012032386A1 (ru) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
IT1402427B1 (it) | 2010-09-06 | 2013-09-04 | Romano | Metodo automatico di misura ed elaborazione della pressione sanguigna. |
ES2810801T3 (es) * | 2015-02-27 | 2021-03-09 | Preventicus Gmbh | Aparato y método para determinar la presión sanguínea |
WO2017013020A1 (en) * | 2015-07-17 | 2017-01-26 | Koninklijke Philips N.V. | Devices, systems, and methods for assessing a vessel |
KR20180010847A (ko) * | 2016-07-22 | 2018-01-31 | 엘지전자 주식회사 | 와치 타입의 이동 단말기 및 그 제어방법 |
US20220400959A1 (en) * | 2021-06-22 | 2022-12-22 | Apple Inc. | Non-invasive blood pressure measurement techniques based on wave shape change during an external pressure cycle |
Family Cites Families (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5810734A (en) * | 1994-04-15 | 1998-09-22 | Vital Insite, Inc. | Apparatus and method for measuring an induced perturbation to determine a physiological parameter |
US5743267A (en) * | 1995-10-19 | 1998-04-28 | Telecom Medical, Inc. | System and method to monitor the heart of a patient |
IT1315206B1 (it) * | 1999-04-27 | 2003-02-03 | Salvatore Romano | Metodo e apparato per la misura della portata cardiaca. |
US7025718B2 (en) * | 2002-11-19 | 2006-04-11 | Jonathan Williams | Method and device for correcting in-vivo sensor drift |
ITRM20030117A1 (it) | 2003-03-17 | 2004-09-18 | Matteo Bonan | Metodo automatizzato di discriminazione del battito cardiaco. |
RU2241373C1 (ru) * | 2003-12-29 | 2004-12-10 | Зао "Вниимп-Вита" | Автоматический измеритель артериального давления |
US20060064021A1 (en) * | 2004-09-20 | 2006-03-23 | David Hefele | Detection and correction of catheter line distortion in blood pressure measurements |
RU2378984C2 (ru) * | 2008-03-31 | 2010-01-20 | Закрытое Акционерное Общество "Поликониус-Центр" | Устройство контроля и оценки физиологических процессов |
IT1402427B1 (it) * | 2010-09-06 | 2013-09-04 | Romano | Metodo automatico di misura ed elaborazione della pressione sanguigna. |
-
2010
- 2010-09-06 IT ITRM2010A000468A patent/IT1402427B1/it active
-
2011
- 2011-09-02 ES ES11773535.7T patent/ES2578998T3/es active Active
- 2011-09-02 CN CN201180042812.7A patent/CN103108585B/zh active Active
- 2011-09-02 RU RU2013115110/14A patent/RU2552685C2/ru active
- 2011-09-02 PL PL11773535.7T patent/PL2613691T3/pl unknown
- 2011-09-02 US US13/820,853 patent/US20130172761A1/en active Pending
- 2011-09-02 EP EP11773535.7A patent/EP2613691B1/en active Active
- 2011-09-02 WO PCT/IB2011/002025 patent/WO2012032386A1/en active Application Filing
- 2011-09-02 BR BR112013005065-9A patent/BR112013005065B1/pt active IP Right Grant
- 2011-09-05 CN CN201180053466.2A patent/CN103200864B/zh active Active
- 2011-09-05 BR BR112013005200A patent/BR112013005200A2/pt not_active Application Discontinuation
- 2011-09-05 WO PCT/IT2011/000308 patent/WO2012032553A1/en active Application Filing
- 2011-09-05 US US13/818,985 patent/US20130150736A1/en not_active Abandoned
- 2011-09-05 PL PL11767797T patent/PL2613690T3/pl unknown
- 2011-09-05 EP EP11767797.1A patent/EP2613690B9/en active Active
- 2011-09-05 AU AU2011300338A patent/AU2011300338B2/en active Active
- 2011-09-05 ES ES11767797.1T patent/ES2503572T3/es active Active
- 2011-09-05 CA CA2809930A patent/CA2809930C/en active Active
-
2014
- 2014-01-09 HK HK14100242.0A patent/HK1187225A1/xx unknown
-
2016
- 2016-11-30 US US15/365,804 patent/US10271742B2/en active Active
Also Published As
Publication number | Publication date |
---|---|
EP2613691A1 (en) | 2013-07-17 |
CN103200864B (zh) | 2015-03-25 |
BR112013005200A2 (pt) | 2016-05-03 |
WO2012032386A1 (en) | 2012-03-15 |
PL2613690T3 (pl) | 2015-03-31 |
EP2613690A1 (en) | 2013-07-17 |
PL2613691T3 (pl) | 2016-11-30 |
US20130150736A1 (en) | 2013-06-13 |
CA2809930C (en) | 2021-05-04 |
WO2012032553A1 (en) | 2012-03-15 |
AU2011300338A1 (en) | 2013-03-14 |
CA2809930A1 (en) | 2012-03-15 |
AU2011300338B2 (en) | 2014-08-21 |
US20170086685A1 (en) | 2017-03-30 |
EP2613690B9 (en) | 2014-09-24 |
HK1187225A1 (en) | 2014-04-04 |
ITRM20100468A1 (it) | 2012-03-07 |
ES2503572T3 (es) | 2014-10-07 |
US20130172761A1 (en) | 2013-07-04 |
CN103200864A (zh) | 2013-07-10 |
RU2013115364A (ru) | 2014-10-20 |
US10271742B2 (en) | 2019-04-30 |
EP2613690B1 (en) | 2014-07-23 |
ES2578998T3 (es) | 2016-08-03 |
RU2552685C2 (ru) | 2015-06-10 |
IT1402427B1 (it) | 2013-09-04 |
CN103108585A (zh) | 2013-05-15 |
BR112013005065B1 (pt) | 2022-12-06 |
BR112013005065A2 (pt) | 2016-06-07 |
CN103108585B (zh) | 2015-04-29 |
EP2613691B1 (en) | 2016-05-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
RU2013115110A (ru) | Способ, устройство и программа для автоматической обработки сигналов кровяного давления | |
CN103654770B (zh) | 移动心电信号qrs波实时波检测方法及装置 | |
CN110236520B (zh) | 基于双卷积神经网络的心电图类型识别装置 | |
CN107432736A (zh) | 一种识别脉搏波形信号的方法 | |
CN104586381A (zh) | 一种基于物联网的心电监护系统 | |
CN105640545A (zh) | 一种胎儿心电信号提取方法及装置 | |
CN102835954A (zh) | 一种心拍波形模板生成、室性早搏检测处理方法及装置 | |
CN103584854A (zh) | 心电信号r波的提取方法 | |
RU2011128710A (ru) | Способ и устройство анализа баллистокардиографических сигналов | |
CN108143407A (zh) | 一种自动提取心音包络特征的心音分段方法 | |
CN104905785A (zh) | 室性/室上性早搏检测方法和系统 | |
CN204520670U (zh) | 一种基于物联网的心电监护系统 | |
CN102334986A (zh) | 动脉压信号中重搏切迹点识别方法 | |
RU2014126422A (ru) | Система и способ идентификации беременностей высокого риска | |
CN104983415A (zh) | 心拍识别的方法及装置 | |
CN109199361B (zh) | 心电信号数据中r峰时刻的提取方法、存储介质和装置 | |
US10750969B2 (en) | Heartbeat detection method and heartbeat detection device | |
RU2010136818A (ru) | Способ лечения фибрилляции предсердий у пациентов с артериальной гипертензией | |
Kang et al. | Heart sound segmentation toward automated heart murmur classification in pediatric patents | |
CN1759387A (zh) | 用于辨别心跳的自动化方法 | |
KR20120133793A (ko) | 특징점 추출을 통한 심실부정맥 검출방법 | |
Atbi et al. | Heart sounds and heart murmurs sepataion | |
RU2440023C1 (ru) | Способ выявления периодических составляющих в ритме сердца | |
RU2015131831A (ru) | Система и способ для уменьшения артефактов движения в экг-сигналах | |
RU2575308C2 (ru) | Автоматический способ измерения и обработки кровяного давления |