INTRODUCCIÓN
La monitorización de señales cardiorrespiratorias en diferentes condiciones (reposo, ambulatoria, ejercicio) y entornos (hospitalarios o no) facilita el diagnóstico, pero enfrenta la pérdida de información debido a perturbaciones que afectan el procesamiento y análisis, y a la carencia de métodos robustos a dichas perturbaciones.
La monitorización incluye las etapas de adquisición, procesamiento y análisis, en las cuales se realizan las diferentes funciones: acondicionamiento, preprocesamiento, reconocimiento de latidos, delineación de puntos característicos, análisis/cálculo de indicadores diagnósticos. Las 3 señales cardiorrespiratorias son: electrocardiograma (ECG), fonocardiograma (PCG) y seismocardiograma (SCG). El SCG y el PCG son el registro no invasivo de las vibraciones mecánicas de baja frecuencia (menor que 30 Hz) del cuerpo y de los sonidos cardíacos, respectivamente, ambos inducidos por los movimientos cardiacos. Por razones de fiabilidad, el interés se centra en analizar series temporales derivadas del: ECG (RR, QT, ritmo respiratorio: Rr), SCG (IMIM) o PCG (S1S1). El objetivo es mejorar los métodos y herramientas para la monitorización del ECG, el PCG y el SCG, enfocado al diagnóstico. Los resultados presentados han sido obtenidos en los últimos 15 años.
MÉTODOS
En la descripción de los métodos se usan diversos enfoques: uniderivacionales (SL) o multiderivacionales (ML), no supervisados o supervisados, aprendizaje automático (AA) o profundo (DL), robustez a derivas de líneas base (DLB), morfologías múltiples, respiración y ruido en la señal así como artefactos en las series temporales (falsos positivos: FP y falsos negativos: FN). Las bases de datos (BD) usadas depende de cada estudio. Así, para evaluar métodos ML, se usan BD con derivaciones ortogonales, con la señal respiratoria si se desea estimarla (PTB, CDB) o sin ella (CSE). Para métodos SL, se usan otras BD (MITDB, NSTD, QTDB, etc.). Para el ECG-SCG se usa la BD CEBS. También se usan BD propias de sujetos registrados con el Holter Excorde3C o con los sistemas desarrollados (ECG-PCG, ECG-WiFi), previo cumplimiento de: las normas de seguridad eléctrica, el consentimiento informado de los voluntarios, y la aprobación de los comités de ética correspondientes. Inicialmente se estudian las relaciones causa-efecto de algunas perturbaciones en el ECG o en sus series temporales, usando métodos existentes. Luego se proponen los nuevos métodos de procesamiento o análisis. Finalmente se describen las herramientas que integran los métodos anteriores.
Efecto del filtrado paso-alto (FPA) en el final de la onda T (Te): El FPA es el método más usado para reducir las DLB. En este estudio se analiza cómo influye la frecuencia de corte (fc) del FPA en la estimación en el Te (tema no publicado antes). 1
Caracterización espectral de las DLB en pruebas de esfuerzo: En este estudio se realiza lo anterior, usando señales con DLB debidas a la respiración sin y con artefactos en los electrodos pertenecientes a la BD NSTDB (tema no publicado antes). 2
Caracterización de los métodos para eliminar las derivas de línea base: Existen disimiles métodos para eliminar las DLB pero han sido validados (o comparados) con diferentes BD y métricas. En este estudio se comparan los 9 métodos más citados (IEEEXplore, WoS, Google-Scholar) en función de la similitud entre señales de entrada y salida. Se usan las mismas señales (QTDB y artificial) y métricas. 3
Caracterización de arquitecturas para reconocer artefactos en el ECG: Para analizar la variabilidad de la frecuencia cardiaca (VFC), cada QRS detectado debe corresponder a un latido sinusal sin bloqueo aurículo-ventricular, sino se considera un latido irregular (o artefacto en el ECG) y deberá excluirse. Hasta 2011 no se reporta 4 ninguna publicación con este enfoque (hay muchas para clasificar arritmias). Este estudio caracteriza el desempeño de 12 arquitecturas que combinan 3 métodos de extracción de características y 4 de agrupamiento, para investigar la combinación óptima con un enfoque SL. Se usa la BD “MIT-BIH Arrhythmia Database”. 4,5
Efecto de la respiración en las estrategias de delineación de la onda T: Se desconoce cuánto influye el efecto mecánico de la respiración (EMR) sobre la delineación SL o ML del Te, basadas en la transformada Wavelet discreta (TW). Para investigar lo anterior se usaron señales simuladas y reales (BD PTB) con 3 derivaciones ortogonales obtenidas mediante la transformación inversa de Dower. 6
Caracterización de la robustez a FP/FN de la Entropía Aproximada (ApEn): La aparición de FP/FN (artefactos) en series temporales (por fallos de los algoritmos de detección) causa errores diagnósticos. La solución tradicional es corregirlos mediante métodos automáticos (no infalibles) o visuales (tediosos). Una solución innovadora es usar indicadores robustos a artefactos. 7 Este estudio evalúa la robustez de la ApEn a FP/FN. Se usan series RR artificiales (no lineales) y reales (BD MITDB) a las que se les introduce (o contienen) FP, FN o ambos. 8 No se conoce un estudio similar.
MÉTODOS
Métodos propuestos: A continuación, se resume el principio de cada método propuesto.
Detección y delineación de QRS: El primer método (SL) estima el pico (R, QRSp), inicio (QRSon) y final del QRS (QRSoff) mediante la primera derivada y las transformadas de Hilbert (TH) y Wavelet Continua (CWT) usando un entrenamiento de 10 s para optimizar los parámetros de procesamiento (umbrales, escalas, etc.). 9 El segundo (ML) proyecta el bucle espacial de la escala 2 de la TW de 3 derivaciones ortogonales, sobre un vector director denominado, dirección de proyección. 10
Estimación del Te (delineación de onda T): Requiere alta exactitud temporal debido a la pequeña variabilidad del QT (menor que 20 ms) y a la lenta variación del ECG alrededor del Te, probablemente contaminado por ruido. El primer método es el único ML y usa el mismo principio del QRS, aunque la delineación de onda T (inicio To, pico Tp y final Te) se realiza en el bucle de las escalas 4 o 5 de la TW. 11 El segundo 12 calcula las áreas sucesivas de un trapecio rectangular con 3 vértices fijos y uno móvil que se desplaza por la señal; Te corresponde a la posición donde el área es máxima. El tercero es un delineador que combina un modelo matemático de función gaussiana asimétrica y el método anterior. 13 El cuarto es un enfoque inédito que consiste en buscar una función de regresión para estimar a Te usando redes neuronales artificiales simples 14) o combinadas con máquinas de vectores soporte (SVM). 15
Detección de latidos irregulares: Se proponen 4 métodos supervisados: 2 con enfoques SL y 2 con ML que usan tensores. El primero (SL) corresponde a la arquitectura de mejor desempeño que resulte de la caracterización descrita antes. 4 El segundo usa redes convolucionales (CNN), post-procesamiento, y sub-clasifica a latidos de origen ventricular. 16 El primero ML combina tensores y SVM, 17) y el segundo emplea descomposición de valores singulares secuencialmente truncada y SVM por mínimos cuadrados. 18
Obtención del ritmo respiratorio (Rr): Se estima el Rr usando las mismas direcciones finales de proyección del delineador ML con TW y el análisis espectral de los lazos espaciales de los puntos QRSp, To, Tp y Te. 10,11) Se compara con otros métodos SL que usan series RR, amplitudes del QRS y con el ML propuesto en Bailón R et al. 19
Obtención de indicadores robustos a artefactos: Se proponen 3 nuevos métodos/indicadores para analizar la VFC y uno de adaptación QT-RR. Los de VFC eliminan las colas inferiores (debajo de XL) y superior (encima de XH) del histograma RR mediante: (a) triangulación de áreas 20,21 (figura 1A), (b) umbrales iguales a la media(RR) ± 2 desviación-estándar (RR) 22) y (c) 2 simples click del especialista. 23 Luego, se calcula, en la nueva serie, los indicadores SDNN o SDSD. 24 El indicador de adaptación QT-RR 23 calcula la pendiente de regresión promediada ponderada, dando mayor peso a los puntos (RR, QT) cercanos a la nube central (figura 1B).
Detección de puntos AO/IM y análisis de la VFC con el SCG: Se realiza una detección gruesa mediante la CWT (Daubechies-4) y una regla de decisión, seguida de una búsqueda fina de máximos (AO) y mínimos (IM) locales en una ventana. 25 Para analizar la VFC, 26 se obtiene la serie IM-IM, se excluyen sus artefactos según 22) y se calculan SDNN, pNN50 y RMSDD 24 para las series RR e IM-IM.
Obtención de la frecuencia cardíaca instantánea (FCI): El método propuesto obtiene la serie FCI mediante un análisis espectro-temporal de la señal de SCG y se evalúa su similitud con la obtenida a partir de la señal de ECG. 27
Detección de puntos en el sonido S1 y análisis de la VFC con el PCG: Los nuevos métodos estiman la posición y duración de S1 mediante: (a) comparación de la envolvente del PCG con una plantilla gaussiana 28 o (b) la triangulación de áreas de la envolvente, 29 de forma similar al descrito en la figura 1. Luego, se evalúa la robustez al ruido. 28) Para analizar la VFC, se evalúa la similitud de las series S1S1 y RR, y la de los indicadores SDNN y SDSD calculados en ambas series. 30
Herramientas desarrolladas
Sistema inalámbrico de monitorización de ECG para dispositivos Android: La monitorización actual emplea dispositivos móviles. Para facilitar esto, se desarrolló un sistema de monitorización Wi-Fi de ECG (un canal), que usa un dispositivo portátil ubicado en el paciente (servidor) y una tableta con sistema operativo Android (cliente), comunicados entre sí. La tableta visualiza (en tiempo real), procesa (off-line) y muestra en pantalla la frecuencia cardíaca media y la desviación RR máxima. 31
Herramientas de adquisición y análisis ECG/PCG: La primera consta de un dispositivo de 2 canales (ECG/PCG) con aislamiento doble (señal, alimentación, cubierta plástica), protección vs desfibrilador, una tarjeta de adquisición (LabPC-1200) insertada en una PC con conductor de protección a tierra (Clase-I/IEC 60601-1) y una herramienta software (LabWindows-CVI) para adquirir/visualizar y crear la BD ECG-PCG. 29) La segunda, denominada CARDIOWAVE2.0, consta de 3 módulos: ECG, PCG, y ECG&PCG. 32 Cada módulo permite: leer registros (varios formatos), filtrar señales, detectar puntos característicos, analizar series temporales simples o cruzadas, corregir artefactos y obtener indicadores temporales y espectrales de VFC y de adaptación QT/RR.
Herramienta para análisis del intervalo QT (PyECG): La interfaz de usuario fue diseñada con la biblioteca PyQt. 33 Sus principales funcionalidades son: indica el canal con mejor calidad de señal (entre 2), filtrados diversos, detección de puntos (Q, R, Te), detección/eliminación de latidos irregulares, utilidades para corregir errores en los puntos detectados, cálculo/visualización de indicadores QTVI y de adaptación del QT.
RESULTADOS Y DISCUSIÓN
Hallazgos obtenidos: a) El FPA a fc mayor que 0,6 Hz distorsiona considerablemente a las 4 morfologías de onda T conllevando a errores inadmisibles en la estimación del Te y QT. Esto revela una inconsistencia de la recomendación AHA 1 que sugiere fc=0,67 Hz; b) el contenido espectral de las DLB durante ejercicio intenso se debe principalmente, al movimiento de los electrodos y, en menor medida, a la respiración; 2 c) la evaluación, en igualdad de condiciones, de los métodos más citados para eliminar las DLB, evidenció (para la morfología tradicional: T+) que el FPA FIR de Kaiser minimiza la distorsión a nivel de latido 3) (aunque no necesariamente en segmentos específicos del mismo), mejorando el reconocimiento de latidos de forma simple y con bajo costo computacional; d) el compromiso entre los índices de desempeño sugiere usar la arquitectura Análisis-componentes-principales-no-lineal/Perceptrón-Multicapa para reconocer latidos irregulares con enfoque SL; (4,5 e) el delineador ML basado en TW minimiza el error de estimación del Te debido al EMR garantizando un diagnóstico más robusto del QT 6 y f) el indicador ApEn de VFC es robusto a la presencia de artefactos en series RR. 8
Resultados obtenidos con los métodos y herramientas propuestos
Detección y delineación de QRS: El procesamiento usando la TW o la CWT minimiza la distorsión e incrementa las tasas de detección. El método de TH/CWT obtuvo altas tasas debido al entrenamiento inicial y al filtrado óptimo y recurrente que enfatiza los puntos característicos y favorece su detección de forma adaptativa, precisa, eficiente y robusta a ruido. 9 El método de delineación ML 10,11 superó en exactitud y precisión al único similar ML 34) que usa la mediana de los puntos R de las 3 derivaciones ortogonales (al ser pocas, no es un buen estimador estadístico).
Detección del punto Te y delineación de la onda T: El método ML superó al que usa la mediana 33) por las razones descritas y se sugiere cuando se dispone del vector-cardiograma. El método de áreas de trapecios (TR) 12) es exacto, no depende de ningún umbral empírico y es muy robusto al ruido blanco gaussiano, siendo idóneo en condiciones ambulatorias o ejercicio. El método de la función-gaussiana-asimétrica-TR logró un alto desempeño con un modelo simple, se adapta a diversas morfologías de onda T (con o sin marcapasos), 13) es robusto a ruido y su modelo no depende de ningún factor empírico, siendo generalizable a varios entornos clínicos. El uso de los métodos supervisados para estimar al Te, usando conjuntos de entrenamiento reducidos abre una alternativa a los cardiólogos que prefieren anotar en vez de corregir. 15 Todos los métodos obtuvieron desempeños similares o superiores al estado del arte.
Detección de latidos irregulares en el ECG: La arquitectura de mejor desempeño con enfoque SL usando AA se destacó en el hallazgo (d). 4 El método con CNN obtuvo un desempeño alto para detectar QRS y otro aceptable como clasificador (segundo lugar en sensibilidad) aunque aun admite mejoras para elevar su valor-predictivo-positivo. 16) Sin embargo, el resto de clasificadores localizan al QRS aplicando sofisticadas técnicas de extracción de características, mientras aquí no es necesario. Hasta nuestro conocimiento, este es el único método que no incorpora: filtrado, detección de puntos característicos, segmentación y extracción o selección de características. Además, su evaluación no se sustenta en segmentos de la señal, sino en todas las muestras de todos los registros. Los resultados sugieren a las CNN como una opción razonable para detectar/clasificar latidos cardiacos. Los métodos con tensores 17,18) obtuvieron altos desempeños con conjuntos de entrenamiento reducidos, siendo una alternativa viable para clasificar registros multicanales con un sólo clasificador sin usar reglas de decisión.
Obtención del ritmo respiratorio (Rr): El método propuesto supera integralmente (en precisión y costo computacional) al resto de los métodos y elimina la incertidumbre de las estimaciones SL. Facilita el diagnóstico cardiorrespiratorio, usando sólo al ECG. 10
Obtención de indicadores robustos a artefactos: La evaluación de los 4 métodos e indicadores propuestos evidenció su gran similitud con los indicadores estándares y su robustez a artefactos en series RR/QT. 20-23 El uso de estos indicadores fiables humaniza y potencia el trabajo de los cardiólogos.
Detección de puntos AO/IM y análisis de la VFC con el SCG: La evaluación del detector propuesto exhibió un alto desempeño y robustez al ruido. 25 Los indicadores SDNN, RMSDD y pNN50 calculados con el SCG y el ECG fueron semejantes. 26 El análisis de la VFC con el SCG reduce costos, contaminación ambiental (sensor económico reusable, sin electrodos), garantiza alta seguridad eléctrica y rápida generalización.
Obtención de la frecuencia cardíaca instantánea (FCI): El método propuesto estimó la FCI con elevada similitud (diferencia media cuadrática menor del 5 %) respecto al obtenido con ECG 27) sin necesidad de detectar latido a latido.
Detección del punto en el sonido S1 y análisis de la VFC con el PCG: El método de la plantilla gaussiana exhibió alta exactitud y robustez al ruido a los efectos deseados (monitoreo en reposo) 28) siendo superior al de triangulación 29) para estimar la posición de S1. El análisis de la VFC con el PCG y el ECG fue similar, 29 garantizando los mismos impactos atribuibles al SCG. Sistema inalámbrico de monitorización de ECG para dispositivos android: La figura 2 muestra la aplicación residente en el cliente y sus fases de ejecución, que culminan visualizando la serie RR y 2 indicadores de VFC. 31
El dispositivo hardware se alimenta de baterías, sin contacto eléctrico directo con el corazón por lo que, cumple las normas de seguridad eléctrica según la norma NTC-IEC 60601-1, siendo clasificado como un equipo de Clase III, tipo BF. El sistema desarrollado abre una alternativa para ampliar los servicios de salud en Cuba (incluyendo el monitoreo usando internet de las cosas) y aumentar la calidad de vida de personas de edad avanzada o con problemas de locomoción.
Herramientas de adquisición y análisis ECG/PCG: Las figuras 3 y 4 muestran las herramientas de adquisición y análisis de ambas señales. 30 Estas integran los métodos/indicadores propuestos. Se verificó la coincidencia numérica de los indicadores calculados con la herramienta de análisis y con programas elaborados en MATLAB. Ambas permitieron realizar la investigación de forma confiable, amistosa y cómoda.
Herramienta para análisis del intervalo QT (PyECG): La evaluación preliminar con registros del Holter cubano Excorde3C y de la BD QTDB fue exitosa. PyCG es la primera herramienta en software libre, 33) que brinda varios indicadores de variabilidad y dinamicidad del intervalo QT.
Conclusiones
Los novedosos métodos desarrollados han alcanzado desempeños similares o superiores a los mejores del estado del arte y son aplicables a diversos escenarios de la monitorización cardiorrespiratoria. Los hallazgos aportan nuevos conocimientos y sugieren procederes para un procesamiento y análisis más eficiente del ECG. Los hallazgos, métodos y herramientas obtenidos contribuyen a un diagnóstico más exacto, rápido e inclusivo usando dispositivos portables, no invasivos, seguros, cómodos, reusables, no contaminantes, de bajo costo y fácil generalización. Esto permite mejorar la calidad de vida de los pacientes y potenciar los servicios cardiológicos en Cuba, ahorrando considerables recursos para nuestro país.