Vessela Krasteva

Real time detection of ventricular fibrillation and tachycardia

The automatic external defibrillator (AED) is a lifesaving device, which processes and analyses the electrocardiogram (ECG) and delivers a defibrillation shock to terminate ventricular fibrillation or tachycardia above 180 bpm. The built-in algorithm for ECG analysis has to discriminate between shockable and non-shockable rhythms and its accuracy, represented by sensitivity and specificity, is aimed at approaching the maximum values of 100%. An algorithm for VF/VT detection is proposed using a band-pass digital filter with integer coefficients, which is very simple to implement in real-time operation. A branch for wave detection is activated for heart rate measurement and an auxiliary parameter calculation. The method was tested with ECG records from the widely recognized databases of the American Heart Association (AHA) and the Massachusetts Institute of Technology (MIT). A sensitivity of 95.93% and a specificity of 94.38% were obtained.

Estimation of current density distribution under electrodes for external defibrillation.

BackgroundTransthoracic defibrillation is the most common life-saving technique for the restoration of the heart rhythm of cardiac arrest victims. The procedure requires adequate application of large electrodes on the patient chest, to ensure low-resistance electrical contact. The current density distribution under the electrodes is non-uniform, leading to muscle contraction and pain, or risks of burning. The recent introduction of automatic external defibrillators and even wearable defibrillators, presents new demanding requirements for the structure of electrodes.Method and ResultsUsing the pseudo-elliptic differential equation of Laplace type with appropriate boundary conditions and applying finite element method modeling, electrodes of various shapes and structure were studied. The non-uniformity of the current density distribution was shown to be moderately improved by adding a low resistivity layer between the metal and tissue and by a ring around the electrode perimeter. The inclusion of openings in long-term wearable electrodes additionally disturbs the current density profile. However, a number of small-size perforations may result in acceptable current density distribution.ConclusionThe current density distribution non-uniformity of circular electrodes is about 30% less than that of square-shaped electrodes. The use of an interface layer of intermediate resistivity, comparable to that of the underlying tissues, and a high-resistivity perimeter ring, can further improve the distribution. The inclusion of skin aeration openings disturbs the current paths, but an appropriate selection of number and size provides a reasonable compromise.

Peripheral nerve magnetic stimulation: influence of tissue non-homogeneity

BackgroundPeripheral nerves are situated in a highly non-homogeneous environment, including muscles, bones, blood vessels, etc. Time-varying magnetic field stimulation of the median and ulnar nerves in the carpal region is studied, with special consideration of the influence of non-homogeneities.MethodsA detailed three-dimensional finite element model (FEM) of the anatomy of the wrist region was built to assess the induced currents distribution by external magnetic stimulation. The electromagnetic field distribution in the non-homogeneous domain was defined as an internal Dirichlet problem using the finite element method. The boundary conditions were obtained by analysis of the vector potential field excited by external current-driven coils.ResultsThe results include evaluation and graphical representation of the induced current field distribution at various stimulation coil positions. Comparative study for the real non-homogeneous structure with anisotropic conductivities of the tissues and a mock homogeneous media is also presented. The possibility of achieving selective stimulation of either of the two nerves is assessed.ConclusionThe model developed could be useful in theoretical prediction of the current distribution in the nerves during diagnostic stimulation and therapeutic procedures involving electromagnetic excitation. The errors in applying homogeneous domain modeling rather than real non-homogeneous biological structures are demonstrated. The practical implications of the applied approach are valid for any arbitrary weakly conductive medium.

Real-time detection of pathological cardiac events in the electrocardiogram.

The development of accurate and fast methods for real-time electrocardiogram (ECG) analysis is mandatory in handheld fully automated monitoring devices for high-risk cardiac patients. The present work describes a simple software method for fast detection of pathological cardiac events. It implements real-time procedures for QRS detection, interbeat RR-intervals analysis, QRS waveform evaluation and a decision-tree beat classifier. Two QRS descriptors are defined to assess (i) the RR interval deviation from the mean RR interval and (ii) the QRS waveform deviation from the QRS pattern of the sustained rhythm. The calculation of the second parameter requires a specific technique, in order to satisfy the demand for straight signal processing with minimum iterations and small memory size. This technique includes fast and resource efficient estimation of a histogram matrix, which accumulates dynamically the amplitude-temporal distribution of the successive QRS pattern waveforms. The pilot version of the method is developed in Matlab and it is tested with internationally recognized ECG databases. The assessment of the online single lead QRS detector showed sensitivity and positive predictivity of above 99%. The classification rules for detection of pathological ventricular beats were defined empirically by statistical analysis. The attained specificity and sensitivity are about 99.5% and 95.7% for all databases and about 99.81% and 98.87% for the noise free dataset. The method is applicable in low computational cost systems for long-term ECG monitoring, such as intelligent holters, automatic event/alarm recorders or personal devices with intermittent wireless data transfer to a central terminal.

Assessment of ECG frequency and morphology parameters for automatic classification of life-threatening cardiac arrhythmias

The reliable recognition and adequate electrical shock therapy of life-threatening cardiac states depend on the electrocardiogram (ECG) descriptors which are used by the defibrillator-embedded automatic arrhythmia analysis algorithms. We propose a method for real-time ECG processing and parameter set extraction using band-pass digital filtration and ECG peak detection. Twelve parameters were derived: (i) seven parameters from the band-pass filter output-six threshold parameters and one peak counter; (ii) five parameters from the ECG peak detection branch, which assess the heart rate, the periodicity and the amplitude/slope symmetry of the ECG peaks. The statistical assessment for more than 36 h of cardiac arrhythmia episodes collected from the public AHA and MIT databases showed that some of the parameters achieved high specificity and sensitivity, but there was no parameter providing 100% separation between non-shockable and shockable rhythms. In order to estimate the influence of the wide variety of cardiac arrhythmias and the different artifacts in real recording conditions, we performed a more detailed study for eight non-shockable and four shockable rhythm categories. The combination of the six top-ranked parameters provided specificity: (i) more than 99% for rhythms with narrow supraventricular complexes, premature ventricular contractions, paced beats and bradycardias; (ii) almost 95% for supraventricular tachycardias; (iii) 91.5% for bundle branch blocks; (iv) 92.2% for slow ventricular tachycardias. The attained sensitivity was above 98% for coarse and fine ventricular fibrillations and 94% for the rapid ventricular tachycardias. The accuracy for the noise contaminated non-shockable and shockable signals exceeded 93%. The proposed parameter set guarantees an accuracy that meets the AHA performance goal for each rhythm category and could be a reliable facility for AED shock-advisory algorithms.

Transthoracic impedance study with large self-adhesive electrodes in two conventional positions for defibrillation

External defibrillation requires the application of high voltage electrical impulses via large external electrodes, placed on selected locations on the thorax surface. The position of the electrodes is one of the major determinants of the transthoracic impedance (TTI) which influences the intracardiac current flow during electric shock and defibrillation success. The variety of factors which influence TTI measurements raised our interest to investigate the range of TTI values and the temporal TTI variance during long-term application of defibrillation self-adhesive electrodes in two conventional positions on the patients chest--position 1 (sub-clavicular/sub-axillar position) and position 2 (antero-posterior position). The prospective study included 86 randomly selected volunteers (39 male and 49 female, 67 patients with normal skin, 13 patients with dry skin and 6 patients with greasy skin, 16 patients with chest pilosity and 70 patients without chest pilosity). The TTI was measured according to the interelectrode voltage drop obtained by passage of a low-amplitude high-frequency current (32 kHz) between the two self-adhesive electrodes (active area about 92 cm2). For each patient, the TTI values were measured within 10 s, 1 min and 5 min after sticking the electrodes to the skin surface, independently for the two tested electrode positions. We found that the expected TTI range is between 58 Omega and 152 Omega for position 1 and between 55 Omega and 149 Omega for position 2. Although the two TTI ranges are comparable, we measured significantly higher TTI mean of about (107.2 +/- 22.3) Omega for position 1 compared to (96.6 +/- 19.2) Omega for position 2 (p = 0.001). This fact suggested that the antero-posterior position of the electrodes is favourable for defibrillation. Within the investigated time interval of 5 min, we observed a significant TTI reduction with about 6.9% (7.4 Omega/107.2 Omega) for position 1 and about 5.3% (5.1 Omega/96.6 Omega) for position 2. We suppose that the long-term application of self-adhesive electrodes would lead to improvement of the physical conditions for conduction of the defibrillation current and to diminution of energy loss in the electrode-skin contact impedance. We found that gender is important when position 1 is used because women have significantly higher TTI (111 +/- 20.3) Omega compared to the TTI of men (102.6 +/- 24) Omega (p = 0.0442). Although we found some specifics of the electrode-skin contact layer, we can conclude that because of the insignificant differences in TTI, the operator of the defibrillator paddles does not need to take into consideration the skin type and pilosity of the patients. Analysis of the correlations between TTI and the individual patient characteristics (chest size, weight, height, age) showed that these patient characteristics are unreliable factors for prediction of the TTI values and optimal defibrillation pulse parameters and energy.

Bench study of the accuracy of a commercial AED arrhythmia analysis algorithm in the presence of electromagnetic interferences.

This paper presents a bench study on a commercial automated external defibrillator (AED). The objective was to evaluate the performance of the defibrillation advisory system and its robustness against electromagnetic interferences (EMI) with central frequencies of 16.7, 50 and 60 Hz. The shock advisory system uses two 50 and 60 Hz band-pass filters, an adaptive filter to identify and suppress 16.7 Hz interference, and a software technique for arrhythmia analysis based on morphology and frequency ECG parameters. The testing process includes noise-free ECG strips from the internationally recognized MIT-VFDB ECG database that were superimposed with simulated EMI artifacts and supplied to the shock advisory system embedded in a real AED. Measurements under special consideration of the allowed variation of EMI frequency (15.7-17.4, 47-52, 58-62 Hz) and amplitude (1 and 8 mV) were performed to optimize external validity. The accuracy was reported using the American Heart Association (AHA) recommendations for arrhythmia analysis performance. In the case of artifact-free signals, the AHA performance goals were exceeded for both sensitivity and specificity: 99% for ventricular fibrillation (VF), 98% for rapid ventricular tachycardia (VT), 90% for slow VT, 100% for normal sinus rhythm, 100% for asystole and 99% for other non-shockable rhythms. In the presence of EMI, the specificity for some non-shockable rhythms (NSR, N) may be affected in some specific cases of a low signal-to-noise ratio and extreme frequencies, leading to a drop in the specificity with no more than 7% point. The specificity for asystole and the sensitivity for VF and rapid VT in the presence of any kind of 16.7, 50 or 60 Hz EMI simulated artifact were shown to reach the equivalence of sensitivity required for non-noisy signals. In conclusion, we proved that the shock advisory system working in a real AED operates accurately according to the AHA recommendations without artifacts and in the presence of EMI. The results may be affected for specificity in the case of a low signal-to-noise ratio or in some extreme frequency setting.

Magnetic stimulation for non-homogeneous biological structures

BackgroundMagnetic stimulation has gained relatively wide application in studying nervous system structures. This technology has the advantage of reduced excitation of sensory nerve endings, and hence results in quasi-painless action. It has become clinically accepted modality for brain stimulation. However, theoretical and practical solutions for assessment of induced current distribution need more detailed and accurate consideration.Some possible analyses are proposed for distribution of the current induced from excitation current contours of different shape and disposition. Relatively non-difficult solutions are shown, applicable for two- and three-dimensional analysis.MethodsThe boundary conditions for field analysis by the internal Dirichlet problem are introduced, based on the vector potential field excited by external current coils. The feedback from the induced eddy currents is neglected. Finite element modeling is applied for obtaining the electromagnetic fields distribution in a non-homogeneous domain.ResultsThe distributions were obtained in a non-homogeneous structure comprised of homogeneous layers. A tendency was found of the induced currents to follow paths in lower resistivity layers, deviating from the expected theoretical course for a homogeneous domain. Current density concentrations occur at the boundary between layers, suggesting the possibility for focusing on, or predicting of, a zone of stimulation.ConclusionThe theoretical basis and simplified approach for generation of 3D FEM networks for magnetic stimulation analysis are presented, applicable in non-homogeneous and non-linear media. The inconveniences of introducing external excitation currents are avoided. Thus, the possibilities are improved for analysis of distributions induced by time-varying currents from contours of various geometry and position with respect to the medium.

Modelling transthoracic defibrillation waveforms

Recent investigations connected with implantable defibrillators yielded new data on heart electrophysiology, resulting in reassessment of existing and advancing of new types of electrical impulses. Different electrical equivalent circuits were proposed for modelling intracardiac and transthoracic defibrillation pulse waveforms, comprising generator, electrode interface and tissue resistances. We attempted modelling of the transmembrane voltage Vm time course, induced by different applied voltage Vs waveforms, taking into account only the shapes and the relative Vs and Vm amplitudes. The excitable cell membrane impedance Zm was modelled with higher resistance and lower capacitance, so that a shunting effect on the generator and tissue resistances was avoided. The result was a very simple equivalent circuit. We proposed criteria for efficient defibrillation pulse waveforms yielding a straightforward approach to model existing and new pulses and to assess their efficiency.Recent investigations connected with implantable defibrillators yielded new data on heart electrophysiology, resulting in reassessment of existing and advancing of new types of electrical impulses. Different electrical equivalent circuits were proposed for modelling intracardiac and transthoracic defibrillation pulse waveforms, comprising generator, electrode interface and tissue resistances. We attempted modelling of the transmembrane voltage Vm time course, induced by different applied voltage Vs waveforms, taking into account only the shapes and the relative Vs and Vm amplitudes. The excitable cell membrane impedance Z was modeled with higher resistance and lower capacitance, so that a shunting effect on the generator and tissue resistances was avoided. The result was a very simple equivalent circuit. We proposed criteria for efficient defibrillation pulse waveforms yielding a straightforward approach to model existing and new pulses and to assess their efficiency.

Automatic adjustment of biphasic pulse duration in transthoracic defibrillation.

Many studies have proven that biphasic defibrillation pulses are more efficient than the damped sinusoid monopolar waveform. Transthoracic resistance was shown to change during the two phases. On the other hand, it was proven that transthoracic resistance plays an important role in the defibrillation process, yielding the current for selected energy or voltage. Pre-shock measurement of the resistance may lead to improved selection. Stabilized current defibrillators are of low stored-to-delivered energy ratio. Therefore, automatic dynamic adjustment of some defibrillator parameters with respect to transthoracic resistance changes seems rational. An approach is known for modifying the pulse duration, in order to deliver a selected energy. A method is proposed here and an experimental defibrillator is developed for dynamic pulse duration adjustment with the purpose of obtaining a desired optimal time-course of the cardiac cell transmembrane potential.Many studies have proven that biphasic defibrillation pulses are more efficient than the damped sinusoid monopolar waveform. Transthoracic resistance was shown to change during the two phases. On the other hand, it was proven that transthoracic resistance plays an important role in the defibrillation process, yielding the current for selected energy or voltage. Pre-shock measurement of the resistance may lead to improved selection. Stabilized current defibrillators are of low stored-to-delivered energy ratio. Therefore, automatic dynamic adjustment of some defibrillator parameters with respect to transthoracic resistance changes seems rational. An approach is known for modifying the pulse duration, in order to deliver a selected energy. A method is proposed here and an experimental defibrillator is developed for dynamic pulse duration adjustment with the purpose of obtaining a desired optimal timecourse of the cardiac cell transmembrane potential.