ES2607721B2 - Método y aparato para estimar el tiempo de tránsito del pulso aórtico a partir de intervalos temporales medidos entre puntos fiduciales del balistocardiograma - Google Patents
Método y aparato para estimar el tiempo de tránsito del pulso aórtico a partir de intervalos temporales medidos entre puntos fiduciales del balistocardiograma Download PDFInfo
- Publication number
- ES2607721B2 ES2607721B2 ES201531414A ES201531414A ES2607721B2 ES 2607721 B2 ES2607721 B2 ES 2607721B2 ES 201531414 A ES201531414 A ES 201531414A ES 201531414 A ES201531414 A ES 201531414A ES 2607721 B2 ES2607721 B2 ES 2607721B2
- Authority
- ES
- Spain
- Prior art keywords
- bcg
- fiducial points
- interval
- ptt
- aortic
- 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
- 238000000034 method Methods 0.000 title claims description 37
- 238000012545 processing Methods 0.000 claims description 7
- 238000004364 calculation method Methods 0.000 claims description 5
- 230000002526 effect on cardiovascular system Effects 0.000 claims description 5
- 238000004891 communication Methods 0.000 claims description 4
- 238000012417 linear regression Methods 0.000 claims description 4
- 238000005259 measurement Methods 0.000 description 19
- 210000001367 artery Anatomy 0.000 description 14
- 230000029058 respiratory gaseous exchange Effects 0.000 description 8
- 210000000709 aorta Anatomy 0.000 description 7
- 230000036772 blood pressure Effects 0.000 description 6
- 230000000747 cardiac effect Effects 0.000 description 6
- 230000001020 rhythmical effect Effects 0.000 description 6
- 239000008280 blood Substances 0.000 description 5
- 210000004369 blood Anatomy 0.000 description 5
- 230000002123 temporal effect Effects 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 230000007211 cardiovascular event Effects 0.000 description 3
- 238000010989 Bland-Altman Methods 0.000 description 2
- 230000004872 arterial blood pressure Effects 0.000 description 2
- 230000001419 dependent effect Effects 0.000 description 2
- 230000005484 gravity Effects 0.000 description 2
- 210000003205 muscle Anatomy 0.000 description 2
- 238000002360 preparation method Methods 0.000 description 2
- 238000012552 review Methods 0.000 description 2
- 208000024172 Cardiovascular disease Diseases 0.000 description 1
- 206010020772 Hypertension Diseases 0.000 description 1
- 230000001133 acceleration Effects 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 238000009610 ballistocardiography Methods 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000037396 body weight Effects 0.000 description 1
- 210000000748 cardiovascular system Anatomy 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000004140 cleaning Methods 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000002354 daily effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 229940079593 drug Drugs 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 230000003203 everyday effect Effects 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000010197 meta-analysis Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000002107 myocardial effect Effects 0.000 description 1
- 208000010125 myocardial infarction Diseases 0.000 description 1
- 230000001991 pathophysiological effect Effects 0.000 description 1
- 238000000053 physical method Methods 0.000 description 1
- 230000002035 prolonged effect Effects 0.000 description 1
- 230000000250 revascularization Effects 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 208000011580 syndromic disease Diseases 0.000 description 1
- 230000009897 systematic effect Effects 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/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
- A61B5/1102—Ballistocardiography
-
- 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
-
- 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/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/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/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/352—Detecting R peaks, e.g. for synchronising diagnostic apparatus; Estimating R-R interval
-
- 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/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]
-
- 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/7246—Details of waveform analysis using correlation, e.g. template matching or determination of similarity
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Cardiology (AREA)
- Biomedical Technology (AREA)
- Medical Informatics (AREA)
- Physics & Mathematics (AREA)
- Veterinary Medicine (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Engineering & Computer Science (AREA)
- Public Health (AREA)
- Heart & Thoracic Surgery (AREA)
- General Health & Medical Sciences (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Physiology (AREA)
- Oral & Maxillofacial Surgery (AREA)
- Dentistry (AREA)
- Vascular Medicine (AREA)
- Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
- Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Description
DESCRIPCIÓN
MÉTODO Y APARATO PARA ESTIMAR EL TIEMPO DE TRÁNSITO DEL PULSO AÓRTICO A PARTIR DE INTERVALOS TEMPORALES MEDIDOS ENTRE PUNTOS
FIDUCIALES DEL BALISTOCARDIOGRAMA
SECTOR DE LA TÉCNICA
La presente invención concierne en general a los sistemas de medición de parámetros fisiológicos por métodos físicos y, en particular, a un método y aparato para estimar el tiempo de tránsito del pulso aórtico a partir de intervalos temporales medidos exclusivamente entre puntos fiduciales del balistocardiograma (BCG).
ESTADO DE LA TÉCNICA
El tiempo de tránsito del pulso (Pulse Transit Time o PTT) generado por la eyección de sangre del corazón al sistema arterial es un parámetro muy importante para diagnosticar el estado del sistema cardiovascular. Se define como el intervalo temporal entre la llegada de la onda de pulso a un punto proximal respecto al corazón y la llegada a otro punto distal, y a partir de él se puede evaluar, por ejemplo, la elasticidad arterial, que es un indicador cada vez más aceptado para predecir el riesgo de padecer enfermedades cardiovasculares. La elasticidad arterial ha sido asociada a la presencia de factores de riesgo cardiovascular y enfermedades arterioescleróticas, y su capacidad para predecir el riesgo de futuros eventos cardiovasculares tales como infarto de miocardio, derrame cerebral, revascularización o síndromes aórticos, entre otros, ha sido ampliamente corroborada, tal como se describe en el documento de C. Vlachopoulos, K. Aznaouridis, and C. Stefanadis, "Prediction of Cardiovascular Events and All-cause Mortality With Arterial Stiffness: a Systematic Review and Metaanalysis,” Journal American College Cardiology, vol. 55, no. 13, pp. 1318-27, Mar.
2010.
El grado de elasticidad de una arteria se evalúa normalmente a partir de la velocidad de la onda de pulso (Pulse Wave Velocity o PWV), de acuerdo con la fórmula descrita por Moens-Korteweg,
Eh
PW V =
------ ,
yj
2
rP
donde E es el módulo de elasticidad de la arteria, h es el grosor de la pared arterial, r es el radio de la arteria y p es la densidad de la sangre.
La medida de PWV en la aorta es la de mayor relevancia clínica porque la aorta y sus principales ramificaciones son responsables de la mayor parte de los efectos patofisiológicos derivados de la rigidez arterial, de modo que la PWV aórtica es un buen indicador del estado de rigidez de las arterias del sujeto. La alta predictividad de la PWV aórtica respecto a eventos cardiovasculares ha sido demostrada en distintos estudios epidemiológicos, tal como se describe en el documento de L. M. Van Bortel, S. Laurent, P. Boutouyrie, P. Chowienczyk, J. K. Cruickshank, et al., "Expert Consensus Document on the Measurement of Aortic Stiffness in Daily Practice Using Carotid-femoral Pulse Wave Velocity,” Journal Hypertension, vol. 30, no. 3, pp. 445 448, Mar. 2012.
La medida de la PWV en una arteria se realiza generalmente de forma no invasiva a partir del PTT en dicha arteria, de acuerdo con
D
PWV = —- ,
PTT
donde D es la distancia entre el punto proximal y el punto distal considerados. En el caso de la PWV aórtica, la medida se realiza generalmente entre el punto carotideo, localizado en la zona medial del borde anterior del músculo esternocleidomastoideo, y el punto femoral localizado en la zona medial del pliegue inguinal. Las arterias en dichos puntos son superficiales, por lo que son fácilmente accesibles para un sensor que establezca un buen contacto con la piel, y el tiempo de tránsito de la onda de pulso entre ellas refleja adecuadamente el PTT aórtico ya que incluye gran parte de la propagación aórtica y aórtico-ilíaca.
Otro parámetro que se puede medir a partir de la elasticidad de una arteria es la presión arterial, ya que el módulo de elasticidad se relaciona con los cambios en la presión arterial media P de acuerdo con
E = E0ekp ,
donde E0 es el módulo de elasticidad de la arteria a una presión arterial media de referencia y k es una constante dependiente de la arteria y cuyo valor está entre
0.016 mmHg"1 y 0.018 mmHg"1. Los cambios en el PTT inducidos por cambios en el módulo de elasticidad de la aorta, y en otras arterias, pueden ser usados para estimar cambios en la presión arterial, y también valores de presión absolutos mediante distintos métodos de calibración, tal como se describe por ejemplo en el documento de D. Buxi, J. M. Redouté, and M. R. Yuce, “A Survey on Signals and Systems in Ambulatory Blood Pressure Monitoring Using Pulse Transit Time,” Physiological Measurements, DOI 10.1088/0967-3334/36/3/R1.
El procedimiento habitual para medir el PTT aórtico de forma no invasiva conlleva la preparación (exposición, limpieza, colocación del sensor y conexión de cables) de los puntos carotideo y femoral para detectar en cada uno de ellos la llegada de la onda de pulso, por ejemplo mediante un fotopletismógrafo (PPG) o un pletismógrafo de impedancia (IPG) que detectan el cambio de volumen local debido a la llegada del pulso arterial, o mediante un tonómetro arterial que mide la presión que ejerce una arteria superficial sobre un sensor de fuerza en estrecho contacto con ella. Éstos y otros sensores capaces de detectar la llegada de la onda de pulso arterial exigen pericia en su colocación, conllevan lentitud en la medición y son incómodos para el sujeto. Además, la aplicación prolongada del sensor puede ocasionar molestias en el sujeto, lo que desaconseja la realización de la medida durante períodos de tiempo largos por las posibles repercusiones fisiológicas de la acción de medida.
Una forma alternativa de obtener información de la actividad mecánica derivada de la eyección cardíaca en la aorta, causante de la onda de pulso, y que requiere una menor preparación del sujeto, es determinar la posición temporal de puntos fiduciales del balistocardiograma (BCG), que refleja las variaciones que experimenta el centro de gravedad del cuerpo humano, ya sea en forma de desplazamiento, velocidad o aceleración, como resultado de la eyección de sangre en cada latido y la posterior propagación del pulso de presión a través de la red arterial. El BCG se puede obtener de múltiples formas, algunas de las cuales pueden implementarse con sensores incorporados en objetos de uso cotidiano, por ejemplo básculas de baño, sillas o camas, tal como se detalla en el documento de O. T. Inan, P. F. Migeotte, K.-S. Park, M. Etemadi, K. Tavakolian, et al., “Ballistocardiography and Seismocardiography: a Review of Recent Advances,” IEEE Journal of Biomedical Health and Informatics, DOI 10.1109/JBHI.2014.2361732, o incorporados en prendas de vestir, tales como calzado o calcetines. Estas soluciones permiten obtener el BCG de forma rápida y cómoda, y
en algún caso por períodos de tiempo largos, por cuanto no se colocan sensores sobre los puntos del cuerpo donde se desea detectar la llegada de la onda de pulso arterial sino que es el cuerpo el que entra en contacto de forma natural con un elemento (plataforma, báscula, silla, cama, prenda de vestir) donde están incorporados los sensores de BCG.
Hasta ahora, los puntos fiduciales del BCG han sido empleados como referencias temporales proximales en medidas de intervalos de tiempo debido a la relación que existe entre dicha señal y el momento de eyección de la sangre en cada latido. Ese es el caso de la patente US 20130310700 A1, en la que se propone usar puntos fiduciales del BCG obtenido a partir de un sistema integrado en una báscula doméstica como referencia temporal proximal para medir el PTT aórtico. El método descrito en dicha patente presenta el inconveniente de que necesita un sensor adicional para detectar la llegada de la onda de pulso arterial al punto distal considerado.
La obtención de información temporal proximal y distal en la misma señal del BCG permitiría medir el PTT aórtico de forma más rápida y cómoda incluso durante períodos de tiempo largos, lo que sería de gran utilidad para evaluar la elasticidad arterial y sus parámetros derivados. También serviría para calcular otros indicadores donde el PTT aórtico interviene junto con otros parámetros cardiovasculares, por ejemplo en la determinación de la contractilidad del miocardio a partir del periodo de pre-eyección (Pre-Ejection Period o PEP), obtenido calculando la diferencia entre el tiempo de llegada de la onda pulso (Pulse Arrival Time o PAT) y el PTT.
BREVE DESCRIPCIÓN DE LA INVENCIÓN
La invención consiste en un método y aparato para estimar el tiempo de tránsito del pulso (PTT) aórtico según se define en las reivindicaciones independientes. En las reivindicaciones dependientes se definen unas formas preferentes de realización de la invención. En la presente descripción y reivindicaciones se ha considerado que el PTT aórtico es el PTT entre el punto carotideo, localizado en la zona medial del borde anterior del músculo esternocleidomastoideo, y el punto femoral localizado en la zona medial del pliegue inguinal.
La solución innovadora propuesta por la presente invención la constituye la estimación
del PTT aórtico a partir de intervalos temporales medidos exclusivamente entre puntos fiduciales del BCG. Como esta señal se obtiene habitualmente mediante sensores integrados en un único elemento con el que entra en contacto el cuerpo del sujeto, el uso del BCG evita la necesidad de sensores de onda de pulso adicionales y la incomodidad de tener que colocar dichos sensores en las zonas específicas donde se desea detectar la llegada de la onda de pulso arterial.
Esta solución innovadora se fundamenta en que las ondas del BCG reflejan cambios en el centro de gravedad del cuerpo humano derivados de la superposición de los efectos de la eyección cardíaca y de la propagación de la onda de pulso arterial. Por ello, se espera que los puntos más iniciales del BCG respecto a la sístole cardíaca estén relacionados mayoritariamente con eventos ligados a la eyección cardíaca, mientras que es de esperar que los puntos fiduciales más alejados del inicio de la señal respecto a la sístole cardíaca estén más influidos por eventos relacionados con la llegada de la onda de pulso a zonas distales. Dado que la aorta es comparativamente la arteria con mayor volumen de sangre y su orientación es longitudinal (paralela al eje cabeza-pies), se espera que las ondas del BCG longitudinal estén especialmente influidas por la actividad mecánica derivada de la propagación de la onda de pulso que se produce en esa arteria principal.
En consecuencia de todo ello, se propone un método para estimar el PTT aórtico consistente, primero, en detectar dos puntos fiduciales de un BCG: un primer punto plausiblemente relacionado con la llegada de la onda de pulso a zonas más proximales respecto al corazón y un segundo punto posterior en el tiempo plausiblemente relacionado con la llegada de la onda de pulso a zonas más distales. A continuación, se mide el intervalo temporal existente entre dichos dos puntos fiduciales. Dicho intervalo se corresponde, en una primera manera de obtenerlo, directamente con el PTT aórtico. Una segunda manera alternativa de obtener dicho tiempo de tránsito a partir del intervalo de tiempo medido entre los dos puntos fiduciales es calibrar el intervalo obtenido del BCG usando como referencia el PTT aórtico obtenido simultáneamente con uno de los métodos ya conocidos pertenecientes al estado de la técnica. Usando la relación obtenida en la calibración, en mediciones posteriores se puede calcular el PTT aórtico a partir del intervalo obtenido exclusivamente del BCG, consiguiendo así una mayor exactitud que en la primera manera propuesta, aunque con un procedimiento inicial más lento y complejo.
Aplicando el método propuesto, los inventores han hallado que, concretamente, las ondas I y J del BCG son sistemáticamente coincidentes con la llegada de la onda de pulso a los puntos carotideo y femoral, respectivamente, por lo que el intervalo IJ resulta especialmente adecuado para obtener directamente de él el PTT aórtico, de acuerdo con la primera manera de obtención propuesta. El uso del segundo método, basado en la calibración del intervalo IJ, también es adecuado si se desea obtener una medida más exacta del PTT aórtico.
Por otro lado, se espera que los intervalos medidos entre otros puntos fiduciales del BCG longitudinal elegidos arbitrariamente sean igualmente sensibles a cambios en el PTT en la aorta, como por ejemplo el intervalo entre las ondas I y K, o el intervalo entre las ondas J y K, aunque la diferente duración de dichos intervalos respecto al PTT aórtico conllevará el uso del segundo método de obtención del PTT aórtico propuesto, basado en la calibración previa de la relación entre el intervalo considerado y el PTT aórtico medido con alguno de los métodos convencionales.
Aunque el método propuesto podría ser implementado por un experto que conociera las relaciones temporales entre las ondas del BCG y la llegada de la onda de pulso a distintos puntos de la red arterial descritas en esta patente de invención, y que identificara visualmente los puntos fiduciales, dentro de un mismo latido, en un registro del BCG y midiera manualmente el intervalo temporal entre ellos, una implementación óptima es mediante un aparato que contenga los sistemas de procesado de señal necesarios para detectar de forma automática los dos puntos fiduciales primero y segundo en una señal de BGC, los sistemas de cálculo necesarios para calcular el intervalo temporal entre dichos puntos fiduciales y para obtener el PTT aórtico a partir de dicho intervalo, y que contuviera, finalmente, un sistema de comunicación del PTT obtenido que se encargue de su representación en un elemento de visualización o de la comunicación del valor a otro dispositivo. Un algoritmo capaz de detectar y medir el intervalo temporal entre las ondas I y J a partir de únicamente la señal del BCG puede ser, por ejemplo, el propuesto en el documento de A. Akhbardeh, B. Kaminska y K. Tavakolian, "BSeg++: A modified Blind Segmentation Method for Ballistocardiogram Cycle Extraction," Annual Intemational Conference of the IEEE Engineering in Medicine and Biology Society (EMBC), 2007, pp.1896-1899. Otros algoritmos pertenecientes al estado de la técnica usan una señal cardiovascular adicional como
origen temporal para detectar las ondas I y J en vez de usar solamente el propio BCG. Por ejemplo, en el documento de Inan et al. antes citado (DOI 10.1109/JBHI.2014.2361732) se propone la detección de la onda J a partir de la medida del máximo local en un intervalo del BCG posterior al pico R del electrocardiograma (ECG). Este método es fácilmente replicable a partir de otras señales cardiovasculares como, por ejemplo, el PPG, el IPG medido localmente, o el IPG medido entre extremidades, señales que pueden ser obtenidas en zonas distales del cuerpo para mayor comodidad.
La invención aquí descrita tiene la ventaja principal de que permite obtener el PTT aórtico utilizando sólo puntos fiduciales del balistocardiograma, lo que facilita la medida del PTT aórtico de forma más sencilla, rápida y cómoda incluso durante períodos de tiempo largos respecto a aquellos sistemas que requieren señales distintas del BCG para obtener al menos uno de los dos puntos fiduciales entre los que se mide un intervalo de tiempo relacionado con el PTT aórtico, o que detectan la llegada de la onda de pulso arterial mediante la colocación de uno o más sensores en las zonas entre las que se desea medir el PTT.
DESCRIPCIÓN DE LOS DIBUJOS
Para complementar la descripción que se está realizando y con objeto de ayudar a una mejor comprensión de las características de la invención, se acompaña como parte integrante de esta descripción un juego de dibujos en donde con carácter ilustrativo y no limitativo, se ha representado lo siguiente:
Figura 1 - Muestra un diagrama de una báscula de baño capaz de obtener el BCG y que constituye el elemento con el que entra en contacto el sujeto en una de las realizaciones de la presente invención.
Figura 2 - Muestra la forma de onda típica del BCG, medida en un sujeto que está de pie y donde se detallan las principales ondas I, J, K, L, M y N que componen el BCG en cada latido y el intervalo IJ. El origen del eje de abscisas coincide con el pico R del electrocardiograma (ECG), aunque esta señal no es imprescindible para medir el intervalo IJ.
Figura 3 - Muestra, de arriba abajo, un registro del ECG (representado para facilitar la comprensión de la invención pero que no es imprescindible para medir intervalos entre puntos del BCG), el BCG obtenido con una báscula, el fotopletismograma (PPG) en el punto carotideo y el PPG en el punto femoral, medidos todos simultáneamente en un mismo sujeto.
Figura 4 - Muestra una traza del intervalo IJ del BCG y el PTT carotideo-femoral obtenidos simultáneamente en un mismo sujeto durante una maniobra de respiración rítmica para inducir cambios en la elasticidad arterial debidos a variaciones de la presión sanguínea asociadas a la respiración.
Figura 5 - Muestra el análisis de la regresión lineal y el análisis Bland-Altman de 407 pares de medidas simultáneas del intervalo IJ y el PTT carotideo-femoral.
Figura 6 - Muestra una traza del intervalo IK del BCG y el PTT carotideo-femoral obtenidos simultáneamente en un mismo sujeto durante una maniobra de respiración rítmica para inducir cambios en la elasticidad arterial debidos a variaciones de la presión sanguínea asociadas a la respiración.
MODOS DE REALIZACIÓN DE LA INVENCIÓN
En una realización preferente de la invención, un sistema integrado en una báscula de baño (1) obtiene un BCG longitudinal indicativo de la actividad mecánica derivada de la eyección cardíaca en la aorta, a partir de un sensor (2) formado por las propias galgas extensiométricas usadas por la báscula para medir el peso corporal y un circuito de procesado de señal analógico (3), según se muestra en la figura 1.
A partir del BCG obtenido en la salida del sistema descrito, el método para estimar el tiempo de tránsito del pulso aórtico consiste, primero, en detectar mediante procesado digital de la señal dos puntos fiduciales en dicho BCG: un primer punto relacionado con la llegada de la onda de pulso arterial a zonas más proximales, que en esta realización correspondería al mínimo de la onda I, y un segundo punto, relacionado con la llegada de la onda de pulso arterial a zonas más distales, que en esta realización correspondería al máximo de la onda J. A continuación, el sistema de
procesado digital de señal (4) que se está usando en esta realización preferente para detectar dichos puntos fiduciales mide el intervalo temporal existente entre ellos, que en esta realización preferente es el intervalo temporal entre el mínimo de la onda I y el máximo de la onda J, denominado intervalo IJ, en cada latido. Dicho intervalo IJ se correspondería, en una primera manera de obtenerlo, con el PTT aórtico. Finalmente, el módulo de comunicación (5) se encarga de comunicar el valor estimado del PTT aórtico del sujeto a través de un monitor LCD.
En la figura 2 se muestra un ejemplo del BCG de uno de dichos latidos, donde se indican las ondas I y J y el intervalo IJ medido es este caso. La figura 3 muestra simultáneamente en un mismo sujeto el ECG, el BCG, y los PPG carotideo y femoral que permiten obtener el PPT aórtico por el método tradicional que implica la necesidad de colocar los sensores de pulso (PPG) en las respectivas zonas del cuerpo donde se desea detectar la llegada de la onda de pulso arterial. En dicha figura se ilustra la correspondencia entre el mínimo de la onda I y el pie de la onda de pulso arterial en el punto carotideo (6), y entre el máximo de la onda J y el pie de la onda de pulso arterial en el punto femoral (7), y cómo, en consecuencia, es posible obtener el PTT aórtico a partir sólo de dos puntos fiduciales del BCG si se sigue el método propuesto en esta realización preferente de la invención.
Para ilustrar la correspondencia entre el intervalo IJ y el PTT aórtico en esta realización, se miden simultáneamente el PTT carotideo-femoral (obtenido de los PPG carotideo y femoral) y el intervalo IJ en un mismo sujeto durante una maniobra de respiración rítmica, que induce cambios en la elasticidad arterial debido a variaciones de presión sanguínea, tal como se muestra en la figura 4. La figura 5 muestra el análisis de la regresión lineal y el análisis Bland-Altman de 407 pares de medidas simultáneas del intervalo IJ y el PTT carotideo-femoral, obtenidas en diversos sujetos bajo respiración rítmica, que ilustran la buena concordancia entre ambos parámetros. Dado que la duración del intervalo IJ es similar a la duración del PTT aórtico y las variaciones en ambos parámetros introducidas por la respiración rítmica tienen una magnitud equivalente, en esta realización preferente se utiliza directamente el intervalo IJ como sustituto del PTT aórtico y la diferencia entre intervalos (media -5,2 ms y desviación estándar 13,2 ms en el grupo analizado en la figura 5) es considerada parte de la incertidumbre intrínseca de la medida.
1
Para mejorar la exactitud en la obtención del PTT aórtico, se propone una segunda manera de determinarlo en la que se obtiene mediante calibración previa la relación entre el intervalo IJ y el PTT aórtico. En esta realización preferente se calcula una regresión lineal entre el intervalo IJ y el PTT aórtico, obtenida a partir de la medida simultánea de ambos intervalos en un grupo de interés o una parte representativa de éste, que permite estimar con mayor exactitud en mediciones posteriores, el PTT aórtico a partir sólo del intervalo IJ a través de la ecuación de la recta obtenida.
La figura 6 muestra los datos obtenidos en otra realización preferente, en la que el PTT se estima a partir del intervalo entre el mínimo de la onda I y el mínimo de la onda K, que es el denominado intervalo IK. El resultado se muestra junto con el PTT carotideo-femoral medido simultáneamente en un mismo sujeto durante una maniobra de respiración rítmica. Tal como sucede con el intervalo IJ en la figura 4, el intervalo IK refleja los cambios en el PTT carotideo-femoral inducidos por la maniobra pero ahora están aumentados debido a que su duración es mayor que la del intervalo IJ, por lo que, en esta realización, el PTT carotideo-femoral se obtiene necesariamente a partir de la calibración del intervalo IK respecto al PTT carotideo-femoral medido con otros medios.
Una vez descrita suficientemente la invención, así como tres realizaciones preferentes, sólo debe añadirse que es posible realizar modificaciones en su constitución, materiales empleados, y en la elección de los sensores empleados para obtener el BCG y de los métodos para identificar los puntos fiduciales de este BCG, sin apartarse del alcance de la invención, definido en las siguientes reivindicaciones.
Claims (1)
- REIVINDICACIONES1 - Un método para estimar el tiempo de tránsito del pulso (PTT) aórtico a partir del intervalo entre puntos fiduciales de un balistocardiograma (BCG), caracterizado porque el método comprende:a) detectar, por un sistema de procesado digital de señal, un primer punto fiducial en un BCG;b) detectar, por dicho sistema de procesado digital de señal, un segundo punto fiducial en el BCG, posterior en el tiempo al primero, pero dentro de un mismo latido; c) medir, por un sistema de cálculo, un intervalo temporal existente entre el primer punto y el segundo punto fiducial; yd) estimar, por el sistema de cálculo, un PTT aórtico a partir del intervalo temporal medido entre los dos puntos fiduciales elegidos,en donde el PTT aórtico estimado se corresponde directamente con el intervalo obtenido entre los puntos fiduciales del BCG detectados, oen donde el PTT aórtico estimado se obtiene mediante la relación entre dicho PTT aórtico y el intervalo de tiempo entre puntos fiduciales del BCG, dicha relación obteniéndose de la calibración de dicho intervalo respecto de otro método de obtención del PTT aórtico.2 - El método según la reivindicación 1, caracterizado porque el primer y el segundo punto fiduciales del BCG se detectan empleando exclusivamente la señal BCG.3 - El método según la reivindicación 1, caracterizado porque los dos puntos fiduciales del BCG se identifican a partir de una señal cardiovascular auxiliar.4 - El método según la reivindicación 3, caracterizado porque se detecta un latido a partir del pico R del electrocardiograma y los dos puntos fiduciales del BCG se detectan en intervalos definidos respecto al pico R del electrocardiograma, correspondientes al mismo latido.5 - El método según la reivindicación 3, caracterizado porque se detecta un latido a partir del pie de la onda de pulso de un fotopletismograma o un pletismograma de impedancia y los dos puntos fiduciales del BCG se detectan en intervalos definidos respecto al pie de dicha onda de pulso correspondientes al mismo latido.6 - El método según una de las reivindicaciones 4 ó 5, caracterizado porque los dos puntos fiduciales del BCG detectados pertenecen a las ondas I y J.7 - El método según una de las reivindicaciones 4 ó 5, caracterizado porque los dos puntos fiduciales del BCG detectados pertenecen a las ondas I y K.8 - El método según una de las reivindicaciones 4 ó 5, caracterizado porque los dos puntos fiduciales del BCG detectados pertenecen a las ondas J y K.9 - El método según cualquiera de las reivindicaciones 1 a 6, caracterizado porque el intervalo entre puntos fiduciales del BCG se mide entre un mínimo de la onda I y un máximo de la onda J.10 - El método según cualquiera de las reivindicaciones 1 a 5 ó 7, caracterizado porque el intervalo entre puntos fiduciales del BCG se mide entre un mínimo de la onda I y un mínimo de la onda K.11 - El método según cualquiera de las reivindicaciones 1 a 5 u 8, caracterizado porque el intervalo entre puntos fiduciales del BCG se mide entre un máximo de la onda J y un mínimo de la onda K.12 - El método según la reivindicación 1, caracterizado porque dicha calibración comprende la realización de una regresión lineal entre el intervalo de tiempo entre puntos fiduciales del BCG y el PTT aórtico obtenido por dicho otro método, obtenido a partir de un fotopletismograma.13 - Un aparato para estimar de manera automática el tiempo de tránsito del pulso (PTT) aórtico a partir del intervalo entre puntos fiduciales de un balistocardiograma (BCG), que comprende:a) un sistema de procesado digital de señal (4) configurado para detectar de forma automática en una señal de BGC dos puntos fiduciales, dentro de un mismo latido, separados en el tiempo;1b) un sistema de cálculo (3) configurado para medir un intervalo temporal entre los dos puntos fiduciales y para estimar un PTT aórtico a partir del intervalo temporal medido entre los dos puntos fiduciales elegidos, en donde el PTT aórtico estimado se corresponde directamente con el intervalo obtenido entre los puntos fiduciales del BCG detectados, o en donde el PTT aórtico estimado se obtiene mediante la relación entre dicho PTT aórtico y el intervalo de tiempo entre puntos fiduciales del BCG, dicha relación obteniéndose de la calibración de dicho intervalo respecto de otro método de obtención del PTT aórtico; yc) un sistema de comunicación (5) configurado para comunicar el PTT aórtico estimado a un usuario o a otro aparato.14 - El aparato según la reivindicación 13, caracterizado porque comprende, adicionalmente, un segundo sistema de cálculo configurado para obtener el PTT aórtico a partir de dicho intervalo temporal.15 - El aparato según una de las reivindicaciones 13 ó 14, caracterizado porque comprende unos medios de entrada de una señal auxiliar y unos medios de determinación de un valor extremo en la señal de BCG posterior a un punto de referencia en dicha señal auxiliar.16 - El aparato según la reivindicación 15, caracterizado porque dicha señal auxiliar es un pletismograma de impedancia, un fotopletismograma o un electrocardiograma.
Priority Applications (7)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
ES201531414A ES2607721B2 (es) | 2015-10-02 | 2015-10-02 | Método y aparato para estimar el tiempo de tránsito del pulso aórtico a partir de intervalos temporales medidos entre puntos fiduciales del balistocardiograma |
KR1020187012248A KR102206785B1 (ko) | 2015-10-02 | 2016-09-30 | 심장탄도의 기준점 사이에서 계측된 시간 간격으로부터 대동맥 맥파 전달 시간을 추정하기 위한 방법 및 기구 |
PCT/ES2016/070692 WO2017055670A1 (es) | 2015-10-02 | 2016-09-30 | Método y aparato para estimar el tiempo de tránsito del pulso aórtico a partir de intervalos temporales medidos entre puntos fiduciales del balistocardiograma |
JP2018516114A JP6670376B2 (ja) | 2015-10-02 | 2016-09-30 | 心弾動図(bcg波)の複数の基準点の間で測定された時間差から大動脈の脈拍の伝搬時間(大動脈ptt)を見積もる方法と装置。 |
US15/765,293 US10925516B2 (en) | 2015-10-02 | 2016-09-30 | Method and apparatus for estimating the aortic pulse transit time from time intervals measured between fiducial points of the ballistocardiogram |
CN201680070423.8A CN108366743A (zh) | 2015-10-02 | 2016-09-30 | 根据在心冲击描记图的基准点之间测量的时间间隔来估算主动脉脉搏传导时间的方法和设备 |
EP16850429.8A EP3357418A4 (en) | 2015-10-02 | 2016-09-30 | Method and apparatus for estimating the transit time of the aortic pulse from time intervals measured between fiducial points of a ballistocardiogram |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
ES201531414A ES2607721B2 (es) | 2015-10-02 | 2015-10-02 | Método y aparato para estimar el tiempo de tránsito del pulso aórtico a partir de intervalos temporales medidos entre puntos fiduciales del balistocardiograma |
Publications (2)
Publication Number | Publication Date |
---|---|
ES2607721A1 ES2607721A1 (es) | 2017-04-03 |
ES2607721B2 true ES2607721B2 (es) | 2019-07-04 |
Family
ID=58422705
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
ES201531414A Active ES2607721B2 (es) | 2015-10-02 | 2015-10-02 | Método y aparato para estimar el tiempo de tránsito del pulso aórtico a partir de intervalos temporales medidos entre puntos fiduciales del balistocardiograma |
Country Status (7)
Country | Link |
---|---|
US (1) | US10925516B2 (es) |
EP (1) | EP3357418A4 (es) |
JP (1) | JP6670376B2 (es) |
KR (1) | KR102206785B1 (es) |
CN (1) | CN108366743A (es) |
ES (1) | ES2607721B2 (es) |
WO (1) | WO2017055670A1 (es) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20200359929A1 (en) * | 2018-01-31 | 2020-11-19 | Taiyo Yuden Co., Ltd. | Pulse wave propagation velocity measurement device and method for same |
CN108852322A (zh) * | 2018-02-26 | 2018-11-23 | 南京邮电大学 | 一种基于pwv的无约束血压测量系统 |
CN110420019B (zh) * | 2019-07-29 | 2021-04-20 | 西安电子科技大学 | 一种心冲击图信号的深度回归心率估计方法 |
KR102421435B1 (ko) * | 2019-09-27 | 2022-07-15 | 한양대학교 산학협력단 | 광용적맥파 신호를 보정하는 방법 및 웨어러블 디바이스 |
CN113827197B (zh) * | 2020-06-08 | 2023-05-05 | 华为技术有限公司 | 脉搏检测方法、终端设备和智能鞋 |
Family Cites Families (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4425920A (en) * | 1980-10-24 | 1984-01-17 | Purdue Research Foundation | Apparatus and method for measurement and control of blood pressure |
CN102197996B (zh) * | 2010-11-23 | 2013-04-17 | 中国科学院合肥物质科学研究院 | 利用桡动脉脉搏波获取主动脉脉搏波速度的方法 |
CA2825405A1 (en) * | 2011-01-27 | 2012-08-02 | The Board Of Trustees Of The Leland Stanford Junior University | Systems and methods for monitoring the circulatory system |
CN103596492B (zh) * | 2011-02-09 | 2017-02-22 | 麻省理工学院 | 耳戴式生命体征监视器 |
SG11201404113TA (en) * | 2012-01-16 | 2014-10-30 | Agency Science Tech & Res | Method and system for optical blood pressure monitoring |
US20160007870A1 (en) * | 2012-03-01 | 2016-01-14 | Koninklijke Philips N.V. | A method of processing a signal representing a physiological rhythm |
US10478082B2 (en) * | 2013-03-24 | 2019-11-19 | Seoul National University R&Db Foundation | Film-type biomedical signal measuring apparatus, blood pressure measuring apparatus using the same, cardiopulmonary fitness estimating apparatus, and personal authentication apparatus |
DE102013213481A1 (de) | 2013-07-10 | 2015-01-15 | Continental Automotive Gmbh | Drucksensorbaugruppe mit gerade verlaufenden Kontakten und senkrecht zur Wandung angeordnetem Schaltungsträger |
FI126008B (en) * | 2013-09-13 | 2016-05-31 | Murata Manufacturing Co | cardiac monitoring system |
CN104545863B (zh) * | 2013-10-10 | 2017-03-29 | 上海宽带技术及应用工程研究中心 | 基于模糊模式识别的bcg心率提取方法及系统 |
US20150112208A1 (en) * | 2013-10-23 | 2015-04-23 | Quanttus, Inc. | Medication management |
US20160081563A1 (en) * | 2014-09-23 | 2016-03-24 | PhysioWave, Inc. | Systems and methods to estimate or measure hemodynamic output and/or related cardiac output |
-
2015
- 2015-10-02 ES ES201531414A patent/ES2607721B2/es active Active
-
2016
- 2016-09-30 EP EP16850429.8A patent/EP3357418A4/en not_active Withdrawn
- 2016-09-30 JP JP2018516114A patent/JP6670376B2/ja active Active
- 2016-09-30 US US15/765,293 patent/US10925516B2/en active Active
- 2016-09-30 CN CN201680070423.8A patent/CN108366743A/zh active Pending
- 2016-09-30 WO PCT/ES2016/070692 patent/WO2017055670A1/es active Application Filing
- 2016-09-30 KR KR1020187012248A patent/KR102206785B1/ko active IP Right Grant
Also Published As
Publication number | Publication date |
---|---|
JP6670376B2 (ja) | 2020-03-18 |
KR102206785B1 (ko) | 2021-01-25 |
KR20180059911A (ko) | 2018-06-05 |
EP3357418A4 (en) | 2018-10-10 |
EP3357418A1 (en) | 2018-08-08 |
ES2607721A1 (es) | 2017-04-03 |
US20180279917A1 (en) | 2018-10-04 |
JP2018534023A (ja) | 2018-11-22 |
US10925516B2 (en) | 2021-02-23 |
WO2017055670A1 (es) | 2017-04-06 |
CN108366743A (zh) | 2018-08-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
ES2607721B2 (es) | Método y aparato para estimar el tiempo de tránsito del pulso aórtico a partir de intervalos temporales medidos entre puntos fiduciales del balistocardiograma | |
US9833151B2 (en) | Systems and methods for monitoring the circulatory system | |
Shriram et al. | Continuous cuffless blood pressure monitoring based on PTT | |
US20200000349A1 (en) | Pulse detection module and use-as-you-need blood pressure measurement device comprising the same | |
ES2525522T3 (es) | Aparato y método para medir la presión sanguínea | |
ES2946534T3 (es) | Dispositivo y método para medir la presión intracraneal | |
ES2398542A2 (es) | Método y aparato para obtener información cardiovascular en los pies | |
WO2013017718A2 (es) | Método y aparato para obtener información cardiovascular midiendo entre dos extremidades | |
US20190343407A1 (en) | An apparatus and method for determining a calibration parameter for a blood pressure measurement device | |
ES2616740B1 (es) | Método y aparato para estimar el tiempo de tránsito del pulso arterial a partir de medidas obtenidas en zonas distales de las extremidades | |
US20180092554A1 (en) | Wearable and unobtrusive multi-sensor array and method for pulse wave velocity imaging | |
Carek et al. | Robust sensing of distal pulse waveforms on a modified weighing scale for ubiquitous pulse transit time measurement | |
Rasool et al. | Continuous and noninvasive blood pressure estimation by two-sensor measurement of pulse transit time | |
Liu et al. | Comparison of simultaneous invasive and non-invasive measurements of blood pressure based upon MIMIC II database | |
ES2656765B1 (es) | Método y aparato para detectar eventos sistólicos mecánicos a partir del balistocardiograma | |
TWI615127B (zh) | 計算血壓的方法及其裝置 | |
Baek et al. | Validation of cuffless blood pressure monitoring using wearable device | |
Muehlsteff et al. | Experiences with Pulse Arrival Time as Surrogate for Systolic Blood Pressure | |
Hruskova et al. | Calculation of the pulse wave velocity from the waveform of the central aortic pressure pulse in young adults | |
Heinz et al. | ENHANCING NOCTURNAL BLOOD PRESSURE ASSESSMENT THROUGH HYDROSTATIC PRESSURE CORRECTION: A CLINICAL STUDY | |
KR20230099290A (ko) | 생체신호를 이용한 혈압 측정 장치 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
FC2A | Grant refused |
Effective date: 20181025 |
|
FG2A | Definitive protection |
Ref document number: 2607721 Country of ref document: ES Kind code of ref document: B2 Effective date: 20190704 |