Heart Artifact Removal in Electrohysterography Measurements Using Algebraic Differentiators

A continuación presento un resumen técnico detallado del artículo "Eliminación de artefactos cardíacos en mediciones de electrohisterografía utilizando diferenciadores algebraicos" en español.

1. Problema y Contexto

La Electrohisterografía (EHG) es un método no invasivo para registrar la actividad eléctrica uterina, crucial para el monitoreo del embarazo, la detección de parto prematuro y el análisis de contracciones uterinas durante el ciclo menstrual. Sin embargo, las grabaciones superficiales de la actividad eléctrica abdominal están severamente contaminadas por múltiples fuentes fisiológicas, siendo la más problemática el electrocardiograma (ECG).

• Desafío: Las señales cardíacas tienen una amplitud potencialmente grande que enmascara la actividad uterina.
• Limitaciones de métodos existentes:
  • Interpolación QRS: Riesgo de distorsionar la dinámica de la señal neural subyacente.
  • Filtrado adaptativo: Requiere canales de referencia de ECG auxiliares (no siempre disponibles).
  • Sustracción de plantillas: Asume morfología estacionaria del ECG; falla si la forma de onda varía.
  • Métodos basados en wavelets: Requieren recursos computacionales elevados o eliminan componentes espectrales de baja frecuencia importantes.
  • Aprendizaje automático: Necesitan grandes conjuntos de datos de entrenamiento y carecen de interpretabilidad.

El objetivo es extraer componentes específicos del útero eliminando los artefactos cardíacos sin canales de referencia externos, manteniendo la integridad de la señal y permitiendo procesamiento en tiempo real.

2. Metodología Propuesta

El trabajo propone un método basado en diferenciadores algebraicos, una herramienta de teoría de control para la estimación de derivadas sin modelo, adaptada para la detección y eliminación de impulsos esparsos (latidos cardíacos).

A. Modelo Matemático de la Señal

La señal observada en cada canal se modela como una superposición de:

1. Componente cardíaca: Trenes de impulsos (latidos) que se propagan a través del tejido.
2. Componente uterina y intestinal: Eventos discretos de menor frecuencia o diferentes características.
3. Línea base: Componentes lentos (respiración, deriva postural).
4. Interferencia de red eléctrica: Ruido de 50 Hz y sus armónicos.
5. Ruido de medición.

Se asume que los eventos impulsivos (latidos) son instantáneos en comparación con la escala temporal de la línea base, lo que permite representarlos mediante distribuciones delta de Dirac.

B. Diferenciación Algebraica

El núcleo del método es la aplicación de un diferenciador algebraico de orden $n$ ($n$-ésima derivada) a la señal.

• Operador: Se utiliza un núcleo causal de soporte compacto $g(t)$ construido a partir de polinomios de Jacobi ponderados.
• Efecto: La diferenciación suprime la línea base suave (cuya $n$-ésima derivada es cercana a cero) y atenúa la interferencia de red eléctrica si los parámetros se eligen correctamente, mientras que los impulsos cardíacos se transforman en firmas características dominantes en la señal derivada.
• Parámetros de diseño:
  • El orden de diferenciación $n$ y los parámetros del núcleo ($\alpha, \beta$) se seleccionan para maximizar la relación señal-ruido (SNR) y asegurar la atenuación de armónicos de la red eléctrica (ajustando la longitud de la ventana $T$ para anular la frecuencia fundamental de 50 Hz).

C. Detección y Agrupamiento (Clustering)

1. Detección de Umbrales: Se identifican los picos en la señal derivada que superan un umbral basado en percentiles. Se consolidan las detecciones cercanas para evitar múltiples conteos de un solo latido.
2. Clustering Multicanal: Se aprovecha la coherencia espacial. Los latidos cardíacos se detectan en la mayoría de los canales con pequeños retrasos de propagación, mientras que los artefactos localizados (movimiento, intestino) afectan a pocos canales.
   • Se agrupan las detecciones temporales de todos los canales en ventanas cortas ($\Delta t_{beat} \approx 20-50$ ms).
   • Un grupo se clasifica como cardíaco si está presente en un porcentaje alto de canales ($\rho \in [0.5, 0.7]$).
3. Estimación de Amplitud y Reconstrucción: Una vez identificados los latidos cardíacos, se estima su amplitud en cada canal mediante mínimos cuadrados. Finalmente, se reconstruye el componente cardíaco y se resta de la señal original filtrada, preservando la actividad uterina.

3. Contribuciones Clave

1. Modelo de señal distribucional: Representación formal de eventos fisiológicos impulsivos frente a una línea base suave.
2. Método sistemático de eliminación: Uso de diferenciadores algebraicos con diseño explícito de parámetros para suprimir simultáneamente interferencia de red y deriva de línea base.
3. Algoritmo de agrupamiento cruzado: Identificación de latidos cardíacos sin necesidad de un canal de referencia de ECG, basándose en la coherencia espacial entre canales.
4. Validación experimental: Comparación contra el método de sustracción de plantillas en grabaciones reales de sujetos humanos (hombre y mujer), demostrando baja complejidad computacional y ausencia de requisitos de entrenamiento.

4. Resultados Experimentales

El método se validó en grabaciones multicanal de EHG (8 electrodos, 5 kHz) de una mujer sana y un hombre sano.

• Rendimiento en sujetos masculinos (lejos del útero):
  • El método eliminó eficazmente los artefactos cardíacos (frecuencia cardíaca ~60 lpm) y el ruido de red.
  • Preservó la estructura de la señal subyacente.
  • Relación Señal-Ruido de salida ($SNR_{out}$): 5.14 dB y 7.42 dB en los canales mostrados, superando o igualando a la sustracción de plantillas sin distorsionar la señal.
• Rendimiento en sujetos femeninos (electrodos cerca del útero):
  • En electrodos muy cercanos al útero, la actividad uterina lenta viola la suposición de derivada base cero, lo que genera falsos positivos (la actividad uterina se confunde con latidos).
  • $SNR_{out}$ negativo (-2.15 dB) en este caso específico, indicando que el método es menos efectivo cuando la señal de interés tiene componentes de baja frecuencia similares a la línea base.
• Propagación Espacial:
  • Se estimaron retrasos intercanales de los pulsos cardíacos (rango de -11.6 ms a 3.0 ms).
  • La coherencia de estos retrasos (rango intercuartílico de 1-2 ms) valida la suposición de que los latidos son eventos globales y coherentes, diferenciándolos de artefactos locales.

5. Significado y Conclusiones

Este trabajo presenta una solución causal, de bajo costo computacional y sin necesidad de aprendizaje para el problema de los artefactos cardíacos en EHG.

• Ventajas principales: No requiere canales de ECG auxiliares, funciona en tiempo real (filtro FIR causal) y es robusto frente a variaciones de morfología del ECG.
• Limitaciones: El método asume que la derivada de la línea base es cero. Esto funciona bien para electrodos alejados del útero, pero falla cerca del útero donde la actividad uterina lenta genera derivadas no nulas que el algoritmo interpreta erróneamente como latidos.
• Futuro: Se sugiere extender el modelo de señal para acomodar derivadas de línea base no nulas (posiblemente mediante identificación de sistemas) y realizar validaciones en cohortes más grandes de pacientes.

En resumen, la propuesta ofrece una herramienta teóricamente sólida y práctica para limpiar señales EHG, facilitando un monitoreo más preciso de la actividad uterina sin la complejidad de los métodos actuales basados en referencias externas o aprendizaje profundo.