Proyecto Transformadas
Proyecto Transformadas
UNIVERSIDAD EL BOSQUE
2018-1
II. INTRODUCCIÓN
I. RESUMEN El ECG es un registro de la actividad eléctrica del
En la actualidad las enfermedades cardiacas generan la corazón. Esta señal se registra desde la superficie
mayor cantidad de muertes a nivel mundial, por esta corporal del paciente y se dibuja en un papel mediante
razón su detección es de vital importancia y se ha una representación gráfica (figura 1), donde se observan
convertido en un prioridad para médicos e ingenieros diferentes ondas que representan los estímulos eléctricos
generar métodos más eficaces para la detección de las de las aurículas y los ventrículos [1].
mismas. Algunas de estas enfermedades se pueden
detectar por medio del análisis del complejo QRS, la
caracterización de sus segmentos, la frecuencia y la
amplitud de este complejo provisiona información
valiosa para la identificación de posibles anomalías
cardiacas. El presente documento muestra el diseño e
implementación de un prototipo capaz de adquirir
señales de electrocardiografía (ECG) realizando una Figura 1. Electrocardiograma.
detección en tiempo real del complejo QRS. [2]
abstract- currently, heart diseases generate the La actividad eléctrica del corazón es la encargada de la
highest number of deaths worldwide, for this reason contracción de la cavidades del corazón, para de esta
their detection is of vital importance and it has forma poder bombear toda la sangre al resto del cuerpo.
become a priority for doctors and engineers to Este impulso eléctrico nace desde el nodo sinusal hasta
generate more effective methods for detecting them. el nódulo atrioventricular, posteriormente viaja hasta el
Some of these diseases can be detected by means of haz de His, el cual se divide en dos ramas y de esta
the QRS complex analysis, the characterization of forma llevar el estímulo hacia ambos ventrìculos [3]. En
their segments, the frequency and the amplitude of la figura 2 se muestra como la secuencia de
this complex provides valuable information for the contracciones y su ubicación generan cada uno de los
identification of possible cardiac anomalies. This componentes que conforman una señal ECG.
document shows the design and implementation of a
prototype capable of acquiring electrocardiographic
(ECG) signals, real-time detection of the QRS
complex.
2
Figura 2. Generación de gráfica ECG según el ciclo cardiaco. Figura 4. Señal ECG con cada uno de sus componentes.
Existen diferentes formas de medir estos potenciales de Cada una de estas ondas responde a un impulso eléctrico
acción, estas formas se conocen como derivaciones ocurrido en el corazón y cualquier anomalía en alguno
cardiacas las cuales miden el diferencial eléctrico entre de los componentes puede ser causa de algún tipo de
un punto y otro. Estas variaciones están dadas según la patología cardiaca, por lo que es importante conocer el
ubicación de los electrodos al momento de la lectura y se origen de cada una de las ondas:
basan en el triángulo de Einthoven (Figura 4) [4].
● ONDA P: Se caracteriza por ser pequeña y
A partir de este triángulo se obtienen las derivaciones semicircular es la primera onda,y está situada
clásicas del ECG cada una de éstas presenta una por encima del cero. Representa la activación
inclinación diferente, por ende la forma y el auricular.
comportamiento de cada una de las ondas cambiará ● ONDA Q: Representa la primera curva negativa
según la derivada seleccionada. después de la onda P y el final del intervalo QT,
suele ser pequeña, no es profunda ni extensa.
Representa la activación ventricular.
● ONDA R: Es la primera onda positiva después
de la onda Q o en caso que la última no se
presente, de la onda P, es una onda corta pero la
de mayor amplitud de toda la señal. Representa
la activación ventricular.
● ONDA S:Es la primera onda negativa después
de la onda R, suele ser pequeña, al igual que la
onda Q. Representa la activación ventricular.
● ONDA T: Es una onda amplia, gruesa y de
forma semicircular. Representa la repolarización
de los ventrículos, marca el final de la acción
Figura 3. Representación triángulo de Einthoven. eléctrica del corazón.
● ONDA U: Esta onda no siempre está presente en
La señal ECG contiene diferentes componentes, cada la señal, es una onda muy pequeña y de forma
uno de estos presenta unas características específicas. En semicircular. Representa una oscilación al
la figura 5 se presenta el nombre con el cual se momento de repolarización de los ventrículos
diferencian cada uno de los picos que conforman una [5].
señal ECG.
El complejo QRS es la representación de la
despolarización del haz de his y de las fibras de purkinje.
La correcta detección de este complejo permitirá detectar
el eje eléctrico del corazón el cual está relacionado
3
directamente con este complejo QRS, en caso que exista electrodos. Por ejemplo, un exceso de electrodos en
alguna inclinación en este eje se puede reflejar algún tipo zonas precordiales podría llegar a simular patrones
de hipertrofia ventricular derecha o izquierda generados por un infarto de miocardio o retrasos
dependiendo la dirección de la inclinación. La duración intraventriculares. Por otra parte, la falta de electrodos
del QRS se desarrolla entre los 0.06 seg y los 0.10 seg, puede generar un pérdida de información de los
alguna anomalía en la amplitud de este complejo puede potenciales o la ausencia de alguna de las ondas que
estar relacionada con bloqueos ventriculares. Por último, componen el ECG [8].
el voltaje de la QRS presenta valores normales entre lo
Por último, una detección en tiempo real puede agilizar
25mm y -25 mm, los valores fuera de estos rangos están
el proceso de tratamiento al paciente, sin embargo, se
relacionados con alguna dificultad en la activación
presentan problemas debido a la cantidad de datos de la
ventricular. Todas estas posibles patologías presentes en
señal, por ejemplo, si se tiene un número de muestras de
el complejo QRS demuestran la importancia de la
la señal (M) y al hacerse el filtrado se realiza una
correcta detección del mismo, de forma que se pueda
convolución de la señal siendo dada por una operación
brindar un oportuno tratamiento de estas. [6]
matemática definida por un kernel (N), el vector
Ruidos presentes en un Señal de ECG resultante de esta operación es un vector de tamaño
N+M+1 muestras, esto evita que la señal tenga
● Interferencia señal de potencia: está
transitoriedad, cuando esto sucede la señal de ECG no se
compuesta por una señal sinusoidal de 60 Hz, se
puede filtrar ni tampoco se puede visualizar el complejo
considera una perturbación cuando su amplitud
QRS ya que los vectores no tienen las mismas
es superior al 50% de la amplitud pico a pico de
dimensiones por tanto no se puede visualizar debido a
la señal ECG.
estas diferencias. Al hacer uso de métodos
● Contacto piel-electrodo: es un ruido generado
computacionales automáticos se busca realizar la
por la pérdida de contacto entre el electrodo y la
optimización de procesos con base a la velocidad de
piel debido al movimiento del paciente suele
procesamiento y adquisición de la señales se vea
presentar una frecuencia de 60 Hz.
reducida y los recursos del ordenador se vean bien
● Desplazamiento de Línea base y Modulación
usados y por ende la velocidad de procesamiento
de amplitud de ECG con respiración: este
aumente.
desplazamiento por respiración puede ser
representado por una un sinusoidal en la
frecuencia respiratoria la cual se adhiere a la
señal ECG. Esta variable puede afectar la III. MARCO TEÓRICO
Anexo 1
amplitud de la onda en un 15% y a la frecuencia
0.15 Hz a 0.3 Hz.
● Ruido electro-quirúrgico: este tipo de ruido IV. METODOLOGÍA
puede destruir por completo la señal ECG
presentando frecuencias desde los 100 Hz hasta
los 1000 Hz, afectando la amplitud de la misma
hasta en un 200%, con una duración en los
1-10s.
[7]
Otro factor es la localización de los electrodos, este
último es de gran importancia, ya que si no se realiza de
forma adecuada se pueden arrojar falsas patologías, si se
invierte la posición de estos las ondas P y QRS pueden
aparecer invertidas en la señal, además, se pueden
presentar errores si no se coloca el número apropiado de
4
● Adquisición de la señal:
Las señales del ECG son obtenidas por medio de
electrodos desechables de materiales hipo
alergénicos, que son pequeños conductores eléctricos
que se adhiere a la piel para medir o causar actividad
eléctrica en el tejido debajo de él, convirtiendo así el
flujo iónico en corriente eléctrica, esto por medio de
una solución conductora que permite que haya un
flujo de electrones en la superficie. [9]
Figura 11. Filtro pasa bajas sallen key. de capacitancia, y son propensos a producir
fallos en las salidas (Al poner en las salidas un
código fuera de secuencia). Hay muchos de
estos utilizados para señales de video u otras
rápidas.
Una de las características más importantes que debe muestreo de 500 Hz y una cantidad mínima de 12 bits
cumplir el amplificador de instrumentación utilizado es para su procesamiento. Por esta razón se decidió trabajar
el valor de CMRR o el factor de rechazo de modo con la daq 6211 ya que cumple con los requerimientos
común. Durante la adquisición de la señal se presenta para la obtención de la señal ECG. [16]
una entrada común, a través de dicha entrada ingresa
ruido blanco proveniente del ambiente que rodea el Selección Amplificadores Operacionales para montaje
sistema, el objetivo del valor de CMRR es rechazar la del filtrado
mayor cantidad de ruido. Para el diseño de un ECG se
recomienda el uso de un amplificador cuyo valor de Por otra parte para el diseño de los filtros necesarios para
CMRR sea superior a los 100 dB. Por esta razón se el ECG se busca implementar amplificadores
decidió trabajar con el amplificador INA 128 el cual operaciones con un valor de CMRR alto y rango de
presenta el valor de CMRR màs alto, además, los alimentación que concuerde con la alimentación del
parámetros mencionados anteriormente como el voltaje amplificador de instrumentación.
de offset bajo, conceden una amplificación de precisión Referencia CMRR Offset V. Corriente de
de los potenciales obtenidos del cuerpo, al utilizar este alimentación polarización
amplificador nos permite obtener una configuración de entrada
diferencial lo que proporciona que haya una gran
precisión en la amplificación de la señal. [7] LM741 90 dB 1 mV ± 3.5 V - ± 200pA
18V
resolución. Para este último parámetro se utilizó la Tabla [Link]ón de amplificadores para filtro.
ecuación de resolución que se muestra a continuación:
Debido a que son los componentes que presentan el
menor valor de offset se trabajará con el LM741 y el
LM358. Además, presentan valores de CMRR dentro del
En donde Vf es el voltaje que ingresa a la tarjeta, en este rango o muy cercano al valor recomendado encontrado
caso para conocer la resolución de cada tarjeta se definió en la teoría.
este valor en 1V pico a pico y N es la cantidad de bits
que posee la tarjeta. ● Caracterización de datos:
[30] [Link]:
Para esta fase se adquirió una muestra proveniente de un
Referencia No. De Bits Frecuencia de Resolución de generador de señal donde se tuvo en cuenta amplitud
muestreo convertidor ADC (500mV), frecuencia(329 Hz) y valor medio de la señal
(Vin= 1Vpp) como se muestra en la figura 14.
Daq 6211 16 250 kHz 350mV
Para que este dispositivo realizara su adquisición de Figura 17. Adquisición de datos después de un pulso.
forma correcta fue necesario alimentarlo con una fuente
de voltaje dual (±2.25V). Así mismo se utilizó un buffer Según la siguiente ecuación se obtiene la resolución del
para lograr que la fuente original pudiera tener esta convertidor:
configuración.
Rango
Ancho de código = Ganancia×2Bits de resolución. Ec(2)
Cálculos
10
V out
Vs
= 1000
V out
Vs
= 1 + R1
R2
R1 =3.3k
R2 = (1 + 3.3k) − 1000
R2= 3.3k
Ec(3)
[Link]:
Para esta fase fue necesario tener en cuenta los
siguientes filtros con sus respectivas frecuencias de
corte: Pasabajas (10Hz), pasaltas (200Hz) y notch
(60Hz) que cuentan con una arquitectura Sallen-key de
segundo orden, estos realizan la eliminación del ruido Figura 25. Filtración de señal.
proveniente de fuentes externas como lo son artefactos
musculares, movimientos involuntario del usuario y la
alimentación.
De acuerdo a las frecuencias de corte varía el orden de
cada uno de los filtros para mejorar la precisión de estos
y tener una señal con el menor ruido posible.
[Link] altas:
C 1 = C 2 = 4.7μf
ω = 2π * 0.5Hz
1
R1 = R2 = 1/2π*0.5Hz *4.7μf
R1 = R2 = 3.3kΩ
[Link] bajas:
C 1 = C 2 = 4.7μf
ω = 2π * 200Hz
1
R1 = R2 = 1/2π*200Hz *4.7μf
R1 = R2 = 180Ω
Tablas
[Link] altas:
Figura 29. Gráfica pasabajos, diagrama de bode. Prueba pasa Vin Vout Ganancia
altas
1 2,20 760 1
2 2,24 670 1
12
transmitir un 1 y después un 0, hay que pasar de 5V a 0V afectan la detección del complejo QRS.
inmediatamente. ● Al identificar y caracterizar los parámetros de
Una manera de solucionar esto es codificando los bits de funcionamiento de los dispositivos de sensado
otra manera; por ejemplo, asociando cada bit a una señal de la señal podemos identificar la mejor forma
que el medio sí admita, como por ejemplo, senos y de establecer el cómo la señal de ECG y así
cosenos; si el medio limita a señales que se encuentren realizar al realizar análisis de la señal identificar
en el rango de 1 kHz y 4 kHz, podemos transmitir una los complejos QRS.
señal sinusoidal de 2 kHz para expresar un 1 y otra de 3
kHz para expresar un 0, lo que sería una manera Referencias:
primitiva de modulación FSK (codificación por [1]Ortega, D. (2016). Fundación del Corazón. [online]
[Link]. Available at:
desplazamiento en frecuencia). Estas señales tienen un
[Link]
tiempo de duración comúnmente llamado tiempo de ntes/metodos-diagnosticos/[Link].
símbolo T, de modo que cada T segundos se transmite [Accessed 5 Mar. 2018].
una de las dos señales. Como cada señal codifica 1 bit,
cada T segundos se transmite 1 bit, luego la tasa de bits [2] N. Dugarte, A. Alvarez, E. Dugarte and G. Álvarez,
es 1/T bps (bits por segundo), que en este caso coincide "Revista del Instituto Nacional de Higiene Rafael Rangel",
con los baudios. [Link], 2015. [Online]. Available:
la cantidad de baudios influyen debido a que la cantidad [Link]
0798-04772015000100007. [Accessed: 07- Mar- 2018].
de bits transmitidos depende de la velocidad o cantidad
de baudios transmitidos en el puerto donde se van a [3] “NI USB-621x User Manual” (2009) [online]
transmitir los datos obtenidos de la digitalización de la [Link]
señal.
[4]”Derivaciones cardiacas del electrocardiograma”. 2013
[Link]. [Online]. Available: [Link]/3Tngi5
[5] [Link] GmbH, G. (2018). Electrocardiograma
● Al observar las aplicaciones de las (ECG/EKG) Interpretación del ECG - [Link]. [online]
[Link]. Available at:
transformadas podemos identificar cual es la
[Link]
más certera a la hora de analizar el complejo [Link] [Accessed 6 Mar. 2018].
QRS.
● Se logra evidenciar como los filtros y [6][Link], 2018. [Online]. Available:
frecuencias de corte junto con la respuesta en [Link]
frecuencia de los dispositivos de adquisición y [Link]?sequence=8. [Accessed: 06- Mar- 2018].
conversión permiten visualizar la señal de
electrocardiograma en las aplicaciones [7] “ANÁLISIS DE ESQUEMAS DE FILTRADO
ANÁLOGO PARA SEÑALES ECG”. Scientia et Technica
observadas en el marco teórico.
Año XIII, No 37, Diciembre de 2007. Universidad
● Se logró determinar que con base a las Tecnológica de Pereira.
frecuencias de corte de la señal ECG en las
aplicaciones el uso de algunas transformadas no [8] “De las ondas de ECG a la patología” 2014 [Online].
permite la deteccion de patologias en el Available:
complejo QRS por ende se descarta. [Link]
● La determinación de frecuencias de corte para
los filtros permite establecer unos parámetros de [9] A. Escarcha, "[Solucionado] Explicación paso a paso de
cómo seguidor de | amplificador-operacional", [Link],
eliminación de ruido provenientes de otras
2015. [Online]. Available:
fuentes como artefactos y movimientos [Link]
involuntarios permite identificar de forma más a-paso-de-como-seguidor-de-voltaje-alcanza-el-estado-estacio
precisa la señal de electrocardiografía de lo que nario-mediante-retroalimentacion-negativa. [Accessed: 08-
no es y así eliminar estas interferencias que Mar- 2018].
14
[10] Carlos A. Alva, Wilfredo Reaño, Joel O. Castillo. [21] Rabanal, F. (2013). Universidad politécnica de Madrid.
“Diseño y Construcción de un Electrocardiógrafo de bajo [online] [Link]. Available at:
costo”. Facultad de Ingeniería, Universidad Ricardo Palma, [Link]
Lima-PERÚ (2011) [Online]. Available: [Link] [Accessed 5 Mar. 2018].
[Link]/pdf/ingenieria/.../CIR-11_Electrocardiografo
_de_Bajo_Costo.pdf [22]A. Sanagustín, "ECG normal (15): Complejo QRS",
[Link], 2018. [Online]. Available:
[11] "Electrodo - EcuRed", [Link], 2015. [Online]. [Link]
Available: [Link] [Accessed: 06- [Link]. [Accessed: 07- Mar- 2018].
Mar- 2018].
[23]S. Andrade, "cienciagora", [Link], 2006.
[12]W. Urbina Rojas and F. Martínez Santa, "Implementing [Online]. Available:
electromyographic equipment through USB interfaces", [Link]
[Link], 2012. [Online]. Available: ngenieria_Electronica_UPB.pdf. [Accessed: 06- Mar- 2018].
[Link]
0123-921X2012000300009. [Accessed: 07- Mar- 2018]. [23] "Texas instrument", [Link], 2018. [Online]. Available:
[Link] [Accessed: 06-
[13] "Penguat Operasional (Op-Amp)", Abi Sabrina, 2010. Mar- 2018].
[Online]. Available:
[Link] [24]S. Ramon, "Analisis y Reduccion del ruido en el ecg",
onal-op-amp/. [Accessed: 08- Mar- 2018]. [Link], 2018. [Online]. Available:
[Link]
[14] Imágenes tomadas de: J. Cabrera, "Filtros Activos", [Link]?sequence=8. [Accessed: 06- Mar- 2018] .
[Link], 2018. [Online]. Available:
[Link] [25] "Texas Instrument", [Link], 2013.
[Link]. [Accessed: 06- Mar- 2018]. [Online]. Available:
[Link]
[15] J. García Álvarez, "AS� FUNCIONA LA A-LM258-LM258A-LM358-LM358A-LM2904-LM2904V_T
CONVERSI�N ANAL�GICO-DIGITAL", [Link], [Link]. [Accessed: 07- Mar- 2018].
2013. [Online]. Available:
[Link] [26] "Electronica", electronicoscaldas, 2015. [Online].
[Link]. [Accessed: 09- Mar- 2018]. Available:
[Link]
[16] A. Galvez, "UNIDAD 4 - ANATOMIA Y FISIOLOGIA les/[Link]
CARDIACA", [Link], 2015. [Online]. Available: [Link]
[Link] es-y-sistemas/temas/tema-2-transformada-de-laplace.
[Accessed: 06- Mar- 2018]. [Accessed: 06- Mar- 2018].
[18]uliana M. Fernández. Rubén C. Acevedo. Carolina b. [28] "Conversión Analógico-Digital (ADC)", [Link],
Tabernig “Electrocardiografía”. Escuela de ingeniería de 2015. [Online]. Available:
Antioquia. Revista EIA, ISSN 1794-1237 Número 7, p. [Link]
111-119. Junio 2007. [Link]. [Accessed: 05- Mar- 2018].
[19]Lavry, Dan (2 de enero de 2012). “Sampling Theory for [29]Univesity of Virgina (2012), [online]
Digital Audio” Lavry Engineering, Inc. [Link]
iac/[Link]
[20] Teorema de muestreo de Nyquist-Shannon. Consultado el
19 de Abril, 2018 [Online] [Link]/WzqszV.
15
ANEXOS 1
ANÁLISIS DE SEÑALES Transformada de puede identificar no permite identificar GONZÁLEZ VILLADA, M. and
ELECTROCARDIOGRÁ fourier por medio de anomalías en la MARTINEZ OSORIO, X. (2009).
FICAS (ECG) CON análisis de frecuencia de las ANÁLISIS DE SEÑALES
ELECTROCARDIOGRÁFICAS
ISQUEMIA CARDIACA densidad espectral señal. (ECG) CON ISQUEMIA CARDIACA
USANDO TÉCNICAS DE fallas o anomalías USANDO TÉCNICAS DE
PROCESAMIENTO en el complejo PROCESAMIENTO DIGITAL DE
DIGITAL DE SEÑALES QRS. SEÑALES. [online]
[Link]. Available at:
[Link]
mlui/bitstream/handle/6789/156/180_G
onzalez_Villada_Marcela_2009.pdf?se
quence=1 [Accessed 6 Mar. 2018].