RU2013115110A - Способ, устройство и программа для автоматической обработки сигналов кровяного давления - Google Patents

Способ, устройство и программа для автоматической обработки сигналов кровяного давления Download PDF

Info

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
Application number
RU2013115110/14A
Other languages
English (en)
Other versions
RU2552685C2 (ru
Inventor
Сальваторе РОМАНО
Original Assignee
Сальваторе РОМАНО
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 Сальваторе РОМАНО filed Critical Сальваторе РОМАНО
Publication of RU2013115110A publication Critical patent/RU2013115110A/ru
Application granted granted Critical
Publication of RU2552685C2 publication Critical patent/RU2552685C2/ru

Links

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
    • A61B5/021Measuring pressure in heart or blood vessels
    • A61B5/02108Measuring pressure in heart or blood vessels from analysis of pulse wave characteristics
    • A61B5/02125Measuring pressure in heart or blood vessels from analysis of pulse wave characteristics of pulse wave propagation time
    • 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
    • A61B5/021Measuring pressure in heart or blood vessels
    • A61B5/02108Measuring pressure in heart or blood vessels from analysis of pulse wave characteristics
    • 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
    • A61B5/021Measuring pressure in heart or blood vessels
    • A61B5/0215Measuring pressure in heart or blood vessels by means inserted into the body
    • 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/7235Details of waveform analysis
    • A61B5/7239Details of waveform analysis using differentiation including higher order derivatives
    • 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/7235Details of waveform analysis
    • A61B5/725Details of waveform analysis using specific filters therefor, e.g. Kalman or adaptive filters
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/74Details of notification to user or communication with user or patient ; user input means
    • A61B5/742Details 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 при работе в соединении со средством обработки устройства.
RU2013115110/14A 2010-09-06 2011-09-02 Способ, устройство и программа для автоматической обработки сигналов кровяного давления RU2552685C2 (ru)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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.

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) Автоматический способ измерения и обработки кровяного давления