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 DESCRIPCIÓN
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 Meta- analysis," 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,
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 EQ 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 Signáis 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.1 109/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 International 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.1 109/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-AItman 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-AItman 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.
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 I J 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.