Introducción
La automatización de la medición de la Presión Sanguínea Arterial (PSA) sigue siendo un reto actual; debido a los métodos que utilizan las diferentes propuestas [8, 14], 15]. Entre estos se encuentra el método oscilométrico, cuyas debilidades son evidentes ante artefactos de movimiento, que alteran la forma de la envolvente de dicha señal [4]. La medición automatizada de la PSA en los equipos comerciales emplea, más usualmente, dicho método. Este método es no invasivo, basado en la colocación e inflado de un brazalete hasta una presión superior a la PSA sistólica del sujeto; y su posterior desinflado gradual a razón constante hasta por debajo de la PSA diastólica. Actualmente la comunidad científica internacional enfoca sus esfuerzos cada vez más, en métodos no invasivos que presenten robustez ante artefactos de movimiento [1]. Incursionando en el uso de diferentes marcas fiduciales que sostengan métodos fiables para la medición automatizada de presión arterial [1, 5, 7].
En un trabajo reciente, se presenta una propuesta para automatizar la medición de la PSA basada en el tiempo transcurrido entre el complejo QRS del electrocardiograma (ECG) y la máxima pendiente en la onda oscilométrica, denominándose tiempo de tránsito del pulso o PTT (por sus siglas en inglés, pulse transit time) [1]. Este método sugiere información que puede delimitar mejor los patrones para determinar los instantes sistólicos, diastólicos y de presión media, con el uso de la señal oscilométrica; y a su vez introduce el uso del ECG como referencia. Como consecuencia, se obtiene una envolvente en función del tiempo PTT.
Otro trabajo reciente para automatizar la medición de la PSA se basa en el tiempo de retraso transcurrido entre el pico R del complejo QRS del electrocardiograma (ECG) y el pulso de la onda oscilométrica [3], [2]. La envolvente en este método se obtiene a partir de cálculos de la media y la desviación estándar de la velocidad de propagación del pulso. Resaltando el valor de una referencia temporal al ECG, sin la necesidad de una medición continúa. No obstante, no cuenta con una referencia a la amplitud del pulso ni ofrece un criterio para pacientes no sanos [12].
Con la finalidad de lograr un sistema de medición robusto de la presión arterial, se desarrolló en el Centro de Estudios de Electrónica y Tecnologías de la Información (CEETI) el método de Contornos de la Imagen Tiempo-Latido Oscilométrico (CITLO) [10]. Este método comprende la adquisición de señales, la obtención de la Imagen Tiempo-Latido Oscilométrico (ITLO), el procesamiento bidimensional de la ITLO y la estimación de los instantes sistólico y diastólico, a partir de indicadores basados en los contornos de la ITLO.
Ambos trabajos [1, 10] aportan individualmente información que puede delimitar mejor los patrones para determinar los instantes sistólicos, diastólicos y de presión media, con el uso de la señal oscilométrica. En la presente investigación se presenta una comparación de las bondades y límites entre dichos métodos. Remarcando el empleo del PTT como referencia temporal, como sugieren otros trabajos investigativos actuales [2,4]; y el uso de la envolvente oscilométrica en la estimación de la presión arterial [8, 9, 10, 11, 12 y 13]. Con la finalidad de ofrecer los fundamentos para implementarlo en un algoritmo de medición automatizada y ambulatoria de PSA.
Materiales y métodos
El método Puntos de Identificación de Parámetros (PIP) será el punto de partida del análisis [5, 6, 7, 8, 9], debido a que básicamente es la línea que sigue el método en cuestión. Este se basa fundamentalmente en la obtención de la posición de cada latido o pulso y en función de ello, determinar la envolvente de la señal oscilométrica detectando los instantes en que ocurren los momentos sistólicos y diastólicos [5, 6, 7, 8, 9].
Descripción de la base de datos empleada
La base de señales utilizada contiene un total de 99 capturas a 33 sujetos, tomando tres mediciones de cada uno. Estas se extrajeron al azar de un conjunto original formado por 108 totales a la que se le agregó 66 como resultado de tomar cada una de las 33 y contaminarla con artefactos de movimiento en dos ocasiones [5]. Se parte de la consideración de que para un sujeto con frecuencia cardiaca de 1Hz (60 pulsaciones/s), ello representa alrededor de 0,5 s con frecuencia de muestreo 500 muestras/s, lo que resulta un valor representativo del pulso, incluso si la frecuencia cardiaca variara entre los valores máximos y mínimos estadísticamente reconocidos de 0,5 a 3,5 [9-10]. Esta base de datos está conformada por las siguientes señales: presión del brazalete, señal de latidos de oscilometría, sonidos de Korotkoff y señal electrocardiografía (derivación II). Además de contar con la medición de referencia, dada por dos expertos respecto a las marcas de los instantes sistólicos (ts1 y ts2) y diastólicos (td1 y td2) como muestra la Fig. 1.
Algoritmo de detección del complejo qrs
Para detectar los complejos QRS de la señal ECG se seleccionó el algoritmo de Pan y Tompkins [2]. Este incluye un filtro pasabanda (pasoalto y pasobajo en cascada), que aísla la energía predominante en el QRS, centrada en 10 Hz, y atenúa las otras componentes ECG, el ruido electromiográfico y las interferencias de la línea de potencia. Posteriormente realiza una diferenciación para encontrar las pendientes altas que normalmente distinguen los complejos QRS de otras ondas del ECG, y seguidamente una potenciación al cuadrado, punto a punto. Después, la señal resultante de los pasos anteriores pasa a través de un integrador de ventana móvil de 150 ms, lo cual abarca todo el complejo QRS sin superposición con la onda T. Esto permite el uso de umbrales, justo encima de los niveles de ruido máximos, por encima de los cuales se tendrán los picos R correspondientes al complejo QRS como el punto máximo de dichos picos, otorgando adicionalmente la posibilidad de localizar las ondas Q y S en complejos de este tipo que tengan duración normal.
Algoritmo de detección del pulso oscilométrico
Para detectar las marcas fiduciarias en la señal oscilométrica, se deriva esta señal y se buscan las localizaciones y amplitudes de los picos máximos, excluyendo los que estén por debajo de la media. Se comprueba que el ancho de los pulsos (a cierta amplitud umbral) sea mayor que 10 muestras, para discriminar ruido. Encontrando los valles en la señal oscilométrica se determina la máxima pendiente.
Conformación de envolvente oscilométrica y del ptt
La envolvente oscilométrica se obtiene buscando los máximos de la señal oscilométrica delimitados por los dos picos consecutivos del ECG (onda R) y los mínimos en dicha señal, pero en esta ocasión entre los dos picos consecutivos de la propia onda oscilométrica. La sustracción de estos será el valor de máxima pendiente que marcará cada pulso [5, 6, 7, 8, 9, 10].
Por su parte, la envolvente del PTT se conforma mediante la determinación del tiempo transcurrido desde el pico de la onda R (ECG) y el máximo valor de la derivada de la señal oscilométrica, acotada igualmente entre los dos picos consecutivos de la onda R.
Resultados y discusión
En la Fig. 2 se contrasta la variabilidad descrita por la envolvente conformada mediante las consideraciones tradicionales del método PIP (máximos oscilométricos) y las de máxima pendiente (empleada en el método PTT para evaluar resultados). Mientras en la Fig. 3 se representan de manera simultánea para ofrecer un primer contraste entre las envolventes obtenidas por cada método.
Para una mejor apreciación de este efecto se ilustran los resultados independientes por el método PIP bajos sus consideraciones tradicionales y empleando el PTT en las Fig. 4 y 5. El inconveniente fundamental del método PIP es que, en cualquiera de sus variantes, el cálculo de la presión sistólica y diastólica es dependiente de la envolvente de los pulsos oscilométricos, que puede variar en función de artefactos de movimiento del sujeto objeto de la medición.
Evaluando la efectividad de este último método en la determinación de los instantes sistólicos y diastólicos; se encuentra que en el primer caso no existen dificultades (significación asintótica es mayor que 0.005, Sig. 0.078), pero en la determinación de la presión diastólica, contrastado con el estándar dorado, existen diferencias significativas en cuanto a su valor medio entre ellos, (significación asintótica es menor que 0.005, Sig. 0.000) ver Tablas 1, 2, 3 y 4. Dicha evaluación se realiza mediante pruebas no paramétricas pues la serie de datos no sigue una distribución normal (observando el estadístico de Shapiro - Wilk, ver Tablas 5 y 6, pues las muestras tienen 99 observaciones); más específicamente, la prueba de rangos con signo de Wilcoxon.
N | Rango promedio | Suma de rangos | ||
---|---|---|---|---|
Ps PIP - Ps OBS | Rangos negativos | 40a | 49,25 | 1970,00 |
Rangos positivos | 59b | 50,51 | 2980,00 | |
Empates | 0c | |||
Total | 99 |
a. Ps PIP < Ps OBS b. Ps PIP > Ps OBS c. Ps PIP = Ps OBS
Ps PIP - Ps OBS | |
---|---|
Z | -1,763a |
Sig. asintót. (bilateral) | ,078 |
a. Basado en los rangos negativos. b. Prueba de los rangos con signo de Wilcoxon
Rangos | N | Rango promedio | Suma de rangos |
---|---|---|---|
Rangos negativos | 40a | 49,25 | 1970,00 |
Rangos positivos | 59b | 50,51 | 2980,00 |
Empates | 0c | ||
Total | 99 |
a. Pd PIP < Pd OBS b. Pd PIP > Pd OBS c. Pd PIP = Pd OBS
Estadísticos de contraste | Pd PIP - Pd OBS |
---|---|
Z | -5,858a |
Sig. asintót. (bilateral) | ,000 |
a. Basado en los rangos negativos. b. Prueba de los rangos con signo de Wilcoxon
Pruebas de normalidad | Kolmogorov-Smirnova | Shapiro-Wilk | ||||
---|---|---|---|---|---|---|
Estadístico | gl | Sig. | Estadístico | gl | Sig. | |
Ps OBS | ,113 | 99 | ,003 | ,907 | 99 | ,000 |
Ps PIP | ,082 | 99 | ,096 | ,961 | 99 | ,005 |
a. Corrección de la significación de Lilliefors
Pruebas de normalidad | ||||||
---|---|---|---|---|---|---|
Kolmogorov-Smirnova | Shapiro-Wilk | |||||
Estadístico | gl | Sig. | Estadístico | gl | Sig. | |
Pd OBS | ,067 | 99 | ,200* | ,981 | 99 | ,172 |
Pd PIP | ,102 | 99 | ,013 | ,970 | 99 | ,025 |
a. Corrección de la significación de Lilliefors
*. Este es un límite inferior de la significación verdadera.
En el caso del PTT, puede apreciarse en la Fig. 5, que no se logra un óptimo resultado (en líneas discontinuas se representa la envolvente real obtenida). Considérese que para la obtención del PTT se emplea la derivada del pulso oscilométrico lo cual acentuará el ruido de alta frecuencia introducido; afectando negativamente el resultado final.
Al graficarse simultáneamente, puede ser apreciado cómo se pierde información respecto a la amplitud del pulso (Fig. 6); además de que se distorsiona la envolvente oscilométrica (Tabla 7). Obsérvense los valores de varianza y medias.
En la envolvente obtenida mediante el método PTT la dispersión de las muestras es significativa. Por lo tanto, aunque dicho método ofrece una referencia temporal al ECG, conlleva un procesamiento posterior de la envolvente obtenida que degrada la fidelidad de la información respecto a la amplitud que esta aporta. Sin embargo, en la Fig. 7, puede observarse cómo la referencia a la señal ECG que ofrece el PTT, puede ser muy válida para análisis posteriores si se emplea en colaboración con la envolvente oscilométrica obtenida siguiendo las consideraciones del método PIP.
N | Mínimo | Máximo | Media | Desv. típ. | Varianza | |
---|---|---|---|---|---|---|
Envolvente PIP | 312 | 0,08133 | 31,449 | 14,9357 | 8,66781 | 75,131 |
Envolvente PTT | 312 | -82,50 | 222,50 | 52,7195 | 73,6153 | 5419,2 |
N válido (según lista) | 312 |
Conclusiones
Al evaluar la conformación de la envolvente a través del método PIP y el PTT, se hacen visibles las flaquezas de este último si es empleado exclusivamente en la medición automatizada de presión arterial.
Esta propuesta del empleo del método PTT ofrece una referencia a la señal ECG válida para posteriores análisis, pero muy variable en la conformación de la envolvente de la onda oscilométrica y desestima la información que la amplitud del pulso representa al utilizar procesamientos digitales de señales para la supresión de ruido y aproximar dicha envolvente. Además, encarece y resta movilidad al sistema de medida para la automatización de la medición de PSA.
Por lo tanto, se propone la combinación de la información que el PTT ofrece con otro método de medición de PSA que presente mayor robustez ante artefactos de movimiento como el método CITLO. Se debe evaluar e integrar los criterios ilustrados, para localizar los instantes sistólico y diastólico, así como continuar ampliando la base de datos. Ello permitirá corroborar la utilidad del PTT como marca fiduciaria en la medición de PSA.
En la presente investigación se consideraron las limitantes y bondades de dichos métodos. Además, se propuso un enfoque novedoso para su empleo, obteniéndose una referencia temporal al ECG sin prescindir de la información aportada por la amplitud del pulso oscilométrico.