Ctrl+K

Generating Abnormal 12-leads ECG

Generating Abnormal 12-leads ECG#

This example can be referenced by citing the package.

import numpy as np
import neurokit2 as nk

Normal Multi-lead ECG#

We can generate a normal multi-lead ECG signal with the following code:

ecg12 = nk.ecg_simulate(duration=10, method="multileads")

# Visualize results
ecg12[0:10000].plot(subplots=True);
../../_images/141780adf0e2d009fa78a17ac3509f7b7a7f610a3f1ff1c0d398a04e88c06ebe.png

Abnormal Multi-lead ECG#

Abnormal ECG can be generated by modifying the simulation parameters. We will draw random parameters from a distribution.

# Normal parameters (used by default)
# ===================================
# t, the starting position along the circle of each interval in radius
ti = np.array((-70, -15, 0, 15, 100))
# a, the amplitude of each spike
ai = np.array((1.2, -5, 30, -7.5, 0.75))
# b, the width of each spike
bi = np.array((0.25, 0.1, 0.1, 0.1, 0.4))

# Add noise
# ===============
ti = np.random.normal(ti, np.ones(5) * 3)
ai = np.random.normal(ai, np.abs(ai / 5))
bi = np.random.normal(bi, np.abs(bi / 5))

ecg12 = nk.ecg_simulate(duration=10, method="multileads", ti=ti, ai=ai, bi=bi)

# Visualize results
ecg12[0:10000].plot(subplots=True);
../../_images/39ce66b67c4898d062918c620d0841244bb9c3de52e4595ba0c7ec555619078c.png

Explanation#

For a better understanding of the above parameters, please read the following.

We can see from the above table that each interval in the 3D trajectory can be fixed by 3 parameters:

  • ti: The starting position $\theta_i/t_i$ along the circle of each interval in radius

  • ai: The amplitude of each spike a

  • bi: The width of each wave b. By altering these 3 parameters we can change the shape of the 3D trajectory and thus change the waveform of the resulting ECG.

Prebuilt Synthetic Data Prebuilt synthetic data can be found at the following link”:

  • https://drive.google.com/drive/folders/1iqyAlyHAvNWdOvjEGn8Y6C8kItsdZ_GC?usp=sharing

Contents