CutFEM-based MEG forward modeling improves source separability and sensitivity to quasi-radial sources: A somatosensory group study.

Source analysis of magnetoencephalography (MEG) data requires the computation of the magnetic fields induced by current sources in the brain. This so-called MEG forward problem includes an accurate estimation of the volume conduction effects in the human head. Here, we introduce the Cut finite element method (CutFEM) for the MEG forward problem. CutFEM's meshing process imposes fewer restrictions on tissue anatomy than tetrahedral meshes while being able to mesh curved geometries contrary to hexahedral meshing. To evaluate the new approach, we compare CutFEM with a boundary element method (BEM) that distinguishes three tissue compartments and a 6-compartment hexahedral FEM in an n = 19 group study of somatosensory evoked fields (SEF). The neural generators of the 20 ms post-stimulus SEF components (M20) are reconstructed using both an unregularized and a regularized inversion approach. Changing the forward model resulted in reconstruction differences of about 1 centimeter in location and considerable differences in orientation. The tested 6-compartment FEM approaches significantly increase the goodness of fit to the measured data compared with the 3-compartment BEM. They also demonstrate higher quasi-radial contributions for sources below the gyral crowns. Furthermore, CutFEM improves source separability compared with both other approaches. We conclude that head models with 6 compartments rather than 3 and the new CutFEM approach are valuable additions to MEG source reconstruction, in particular for sources that are predominantly radial.

[Developments of ex vivo cardiac electrical mapping and intelligent labeling of atrial fibrillation substrates].

Cardiac three-dimensional electrophysiological labeling technology is the prerequisite and foundation of atrial fibrillation (AF) ablation surgery, and invasive labeling is the current clinical method, but there are many shortcomings such as large trauma, long procedure duration, and low success rate. In recent years, because of its non-invasive and convenient characteristics, ex vivo labeling has become a new direction for the development of electrophysiological labeling technology. With the rapid development of computer hardware and software as well as the accumulation of clinical database, the application of deep learning technology in electrocardiogram (ECG) data is becoming more extensive and has made great progress, which provides new ideas for the research of ex vivo cardiac mapping and intelligent labeling of AF substrates. This paper reviewed the research progress in the fields of ECG forward problem, ECG inverse problem, and the application of deep learning in AF labeling, discussed the problems of ex vivo intelligent labeling of AF substrates and the possible approaches to solve them, prospected the challenges and future directions for ex vivo cardiac electrophysiology labeling.

Non-Contact Monitoring of ECG in the Home Environment-Selecting Optimal Electrode Configuration.

Capacitive electrocardiography (cECG) is most often used in wearable or embedded measurement systems. The latter is considered in the paper. An optimal electrocardiographic lead, as an individual feature, was determined based on model studies. It was defined as the possibly highest value of the R-wave amplitude measured on the back of the examined person. The lead configuration was also analyzed in terms of minimizing its susceptibility to creating motion artifacts. It was found that the direction of the optimal lead coincides with the electrical axis of the heart. Moreover, the electrodes should be placed in the areas preserving the greatest voltage and at the same time characterized by the lowest gradient of the potential. Experimental studies were conducted using the developed measurement system on a group of 14 people. The ratio of the R-wave amplitude (as measured on the back and chest, using optimal leads) was less than 1 while the SNR reached at least 20 dB. These parameters allowed for high-quality QRS complex detection with a PPV of 97%. For the "worst" configurations of the leads, the signals measured were practically uninterpretable.

A Finite-Difference Solution for the EEG Forward Problem in Inhomogeneous Anisotropic Media.

Accurate source localization of electroencephalographic (EEG) signals requires detailed information about the geometry and physical properties of head tissues. Indeed, these strongly influence the propagation of neural activity from the brain to the sensors. Finite difference methods (FDMs) are head modelling approaches relying on volumetric data information, which can be directly obtained using magnetic resonance (MR) imaging. The specific goal of this study is to develop a computationally efficient FDM solution that can flexibly integrate voxel-wise conductivity and anisotropy information. Given the high computational complexity of FDMs, we pay particular attention to attain a very low numerical error, as evaluated using exact analytical solutions for spherical volume conductor models. We then demonstrate the computational efficiency of our FDM numerical solver, by comparing it with alternative solutions. Finally, we apply the developed head modelling tool to high-resolution MR images from a real experimental subject, to demonstrate the potential added value of incorporating detailed voxel-wise conductivity and anisotropy information. Our results clearly show that the developed FDM can contribute to a more precise head modelling, and therefore to a more reliable use of EEG as a brain imaging tool.

New Strategy for Finite Element Mesh Generation for Accurate Solutions of Electroencephalography Forward Problems.

The finite element method (FEM) is a numerical method that is often used for solving electroencephalography (EEG) forward problems involving realistic head models. In this study, FEM solutions obtained using three different mesh structures, namely coarse, densely refined, and adaptively refined meshes, are compared. The simulation results showed that the accuracy of FEM solutions could be significantly enhanced by adding a small number of elements around regions with large estimated errors. Moreover, it was demonstrated that the adaptively refined regions were always near the current dipole sources, suggesting that selectively generating additional elements around the cortical surface might be a new promising strategy for more efficient FEM-based EEG forward analysis.

Using reciprocity for relating the simulation of transcranial current stimulation to the EEG forward problem.

To explore the relationship between transcranial current stimulation (tCS) and the electroencephalography (EEG) forward problem, we investigate and compare accuracy and efficiency of a reciprocal and a direct EEG forward approach for dipolar primary current sources both based on the finite element method (FEM), namely the adjoint approach (AA) and the partial integration approach in conjunction with a transfer matrix concept (PI). By analyzing numerical results, comparing to analytically derived EEG forward potentials and estimating computational complexity in spherical shell models, AA turns out to be essentially identical to PI. It is then proven that AA and PI are also algebraically identical even for general head models. This relation offers a direct link between the EEG forward problem and tCS. We then demonstrate how the quasi-analytical EEG forward solutions in sphere models can be used to validate the numerical accuracies of FEM-based tCS simulation approaches. These approaches differ with respect to the ease with which they can be employed for realistic head modeling based on MRI-derived segmentations. We show that while the accuracy of the most easy to realize approach based on regular hexahedral elements is already quite high, it can be significantly improved if a geometry-adaptation of the elements is employed in conjunction with an isoparametric FEM approach. While the latter approach does not involve any additional difficulties for the user, it reaches the high accuracies of surface-segmentation based tetrahedral FEM, which is considerably more difficult to implement and topologically less flexible in practice. Finally, in a highly realistic head volume conductor model and when compared to the regular alternative, the geometry-adapted hexahedral FEM is shown to result in significant changes in tCS current flow orientation and magnitude up to 45° and a factor of 1.66, respectively.

The role of blood vessels in high-resolution volume conductor head modeling of EEG.

Reconstruction of the electrical sources of human EEG activity at high spatio-temporal accuracy is an important aim in neuroscience and neurological diagnostics. Over the last decades, numerous studies have demonstrated that realistic modeling of head anatomy improves the accuracy of source reconstruction of EEG signals. For example, including a cerebro-spinal fluid compartment and the anisotropy of white matter electrical conductivity were both shown to significantly reduce modeling errors. Here, we for the first time quantify the role of detailed reconstructions of the cerebral blood vessels in volume conductor head modeling for EEG. To study the role of the highly arborized cerebral blood vessels, we created a submillimeter head model based on ultra-high-field-strength (7T) structural MRI datasets. Blood vessels (arteries and emissary/intraosseous veins) were segmented using Frangi multi-scale vesselness filtering. The final head model consisted of a geometry-adapted cubic mesh with over 17×10(6) nodes. We solved the forward model using a finite-element-method (FEM) transfer matrix approach, which allowed reducing computation times substantially and quantified the importance of the blood vessel compartment by computing forward and inverse errors resulting from ignoring the blood vessels. Our results show that ignoring emissary veins piercing the skull leads to focal localization errors of approx. 5 to 15mm. Large errors (>2cm) were observed due to the carotid arteries and the dense arterial vasculature in areas such as in the insula or in the medial temporal lobe. Thus, in such predisposed areas, errors caused by neglecting blood vessels can reach similar magnitudes as those previously reported for neglecting white matter anisotropy, the CSF or the dura - structures which are generally considered important components of realistic EEG head models. Our findings thus imply that including a realistic blood vessel compartment in EEG head models will be helpful to improve the accuracy of EEG source analyses particularly when high accuracies in brain areas with dense vasculature are required.

The Role of Skull Modeling in EEG Source Imaging for Patients with Refractory Temporal Lobe Epilepsy.

We investigated the influence of different skull modeling approaches on EEG source imaging (ESI), using data of six patients with refractory temporal lobe epilepsy who later underwent successful epilepsy surgery. Four realistic head models with different skull compartments, based on finite difference methods, were constructed for each patient: (i) Three models had skulls with compact and spongy bone compartments as well as air-filled cavities, segmented from either computed tomography (CT), magnetic resonance imaging (MRI) or a CT-template and (ii) one model included a MRI-based skull with a single compact bone compartment. In all patients we performed ESI of single and averaged spikes marked in the clinical 27-channel EEG by the epileptologist. To analyze at which time point the dipole estimations were closer to the resected zone, ESI was performed at two time instants: the half-rising phase and peak of the spike. The estimated sources for each model were validated against the resected area, as indicated by the postoperative MRI. Our results showed that single spike analysis was highly influenced by the signal-to-noise ratio (SNR), yielding estimations with smaller distances to the resected volume at the peak of the spike. Although averaging reduced the SNR effects, it did not always result in dipole estimations lying closer to the resection. The proposed skull modeling approaches did not lead to significant differences in the localization of the irritative zone from clinical EEG data with low spatial sampling density. Furthermore, we showed that a simple skull model (MRI-based) resulted in similar accuracy in dipole estimation compared to more complex head models (based on CT- or CT-template). Therefore, all the considered head models can be used in the presurgical evaluation of patients with temporal lobe epilepsy to localize the irritative zone from low-density clinical EEG recordings.

Influence of the head model on EEG and MEG source connectivity analyses.

The results of brain connectivity analysis using reconstructed source time courses derived from EEG and MEG data depend on a number of algorithmic choices. While previous studies have investigated the influence of the choice of source estimation method or connectivity measure, the effects of the head modeling errors or simplifications have not been studied sufficiently. In the present simulation study, we investigated the influence of particular properties of the head model on the reconstructed source time courses as well as on source connectivity analysis in EEG and MEG. Therefore, we constructed a realistic head model and applied the finite element method to solve the EEG and MEG forward problems. We considered the distinction between white and gray matter, the distinction between compact and spongy bone, the inclusion of a cerebrospinal fluid (CSF) compartment, and the reduction to a simple 3-layer model comprising only the skin, skull, and brain. Source time courses were reconstructed using a beamforming approach and the source connectivity was estimated by the imaginary coherence (ICoh) and the generalized partial directed coherence (GPDC). Our results show that in both EEG and MEG, neglecting the white and gray matter distinction or the CSF causes considerable errors in reconstructed source time courses and connectivity analysis, while the distinction between spongy and compact bone is just of minor relevance, provided that an adequate skull conductivity value is used. Large inverse and connectivity errors are found in the same regions that show large topography errors in the forward solution. Moreover, we demonstrate that the very conservative ICoh is relatively safe from the crosstalk effects caused by imperfect head models, as opposed to the GPDC.

A guideline for head volume conductor modeling in EEG and MEG.

For accurate EEG/MEG source analysis it is necessary to model the head volume conductor as realistic as possible. This includes the distinction of the different conductive compartments in the human head. In this study, we investigated the influence of modeling/not modeling the conductive compartments skull spongiosa, skull compacta, cerebrospinal fluid (CSF), gray matter, and white matter and of the inclusion of white matter anisotropy on the EEG/MEG forward solution. Therefore, we created a highly realistic 6-compartment head model with white matter anisotropy and used a state-of-the-art finite element approach. Starting from a 3-compartment scenario (skin, skull, and brain), we subsequently refined our head model by distinguishing one further of the above-mentioned compartments. For each of the generated five head models, we measured the effect on the signal topography and signal magnitude both in relation to a highly resolved reference model and to the model generated in the previous refinement step. We evaluated the results of these simulations using a variety of visualization methods, allowing us to gain a general overview of effect strength, of the most important source parameters triggering these effects, and of the most affected brain regions. Thereby, starting from the 3-compartment approach, we identified the most important additional refinement steps in head volume conductor modeling. We were able to show that the inclusion of the highly conductive CSF compartment, whose conductivity value is well known, has the strongest influence on both signal topography and magnitude in both modalities. We found the effect of gray/white matter distinction to be nearly as big as that of the CSF inclusion, and for both of these steps we identified a clear pattern in the spatial distribution of effects. In comparison to these two steps, the introduction of white matter anisotropy led to a clearly weaker, but still strong, effect. Finally, the distinction between skull spongiosa and compacta caused the weakest effects in both modalities when using an optimized conductivity value for the homogenized compartment. We conclude that it is highly recommendable to include the CSF and distinguish between gray and white matter in head volume conductor modeling. Especially for the MEG, the modeling of skull spongiosa and compacta might be neglected due to the weak effects; the simplification of not modeling white matter anisotropy is admissible considering the complexity and current limitations of the underlying modeling approach.

Forward problem of electrocardiography based on cardiac source vector orientations.

To localize the unusual cardiac activities non-invasively, one has to build a prior forward model that relates the heart, torso, and detectors. This model has to be constructed to mathematically relate the geometrical and functional activities of the heart. Several methods are available to model the prior sources in the forward problem, which results in the lead field matrix generation. In the conventional technique, the lead field assumed the fixed prior sources, and the source vector orientations were presumed to be parallel to the detector plane with the unit strength in all directions. However, the anomalies cannot always be expected to occur in the same location and orientation, leading to misinterpretation and misdiagnosis. To overcome this, the work proposes a new forward model constructed using the VCG signals of the same subject. Furthermore, three transformation methods were used to extract VCG in constructing the time-varying lead fields to steer to the orientation of the source rather than just reconstructing its activities in the inverse problem. In addition, the unit VCG loop of the acute ischemia patient was extracted to observe the changes compared to the normal subject. The abnormality condition was achieved by delaying the depolarization time by 15ms. The results involving the unit vectors of VCG demonstrated the anisotropic nature of cardiac source orientations, providing information about the heart's electrical activity.

Fast and Accurate EEG/MEG BEM-Based Forward Problem Solution for High-Resolution Head Models.

A BEM (boundary element method) based approach is developed to accurately solve an EEG/MEG forward problem for a modern high-resolution head model in approximately 60 seconds using a common workstation. The method utilizes a charge-based BEM with fast multipole acceleration (BEM-FMM) and a "smart" mesh pre-refinement (called b-refinement) close to the singular source(s). No costly matrix-filling or direct solution steps typical for the standard BEM are required; the method generates on-skin voltages as well as MEG magnetic fields for high-resolution head models in approximately 60 seconds after initial model assembly. The method is verified both theoretically and experimentally.

A hybrid boundary element-finite element approach for solving the EEG forward problem in brain modeling.

This article introduces a hybrid BE-FE method for solving the EEG forward problem, leveraging the strengths of both the Boundary Element Method (BEM) and Finite Element Method (FEM). FEM accurately models complex and anisotropic tissue properties for realistic head geometries, while BEM excels in handling isotropic tissue regions and dipolar sources efficiently. The proposed hybrid method divides regions into homogeneous boundary element (BE) regions that include sources and heterogeneous anisotropic finite element (FE) regions. So, BEM models the brain, including dipole sources, and FEM models other head layers. Validation includes inhomogeneous isotropic/anisotropic three- and four-layer spherical head models, and a four-layer MRI-based realistic head model. Results for six dipole eccentricities and two orientations are computed using BEM, FEM, and hybrid BE-FE method. Statistical analysis, comparing error criteria of RDM and MAG, reveals notable improvements using the hybrid FE-BE method. In the spherical head model, the hybrid BE-FE method compared with FEM demonstrates enhancements of at least 1.05 and 38.31% in RDM and MAG criteria, respectively. Notably, in the anisotropic four-layer head model, improvements reach a maximum of 88.3% for RDM and 93.27% for MAG over FEM. Moreover, in the anisotropic four-layer realistic head model, the proposed hybrid method exhibits 55.4% improvement in RDM and 89.3% improvement in MAG compared to FEM. These findings underscore the proposed method is a promising approach for solving the realistic EEG forward problems, advancing neuroimaging techniques and enhancing understanding of brain function.

CutFEM forward modeling for EEG source analysis.

Introduction: Source analysis of Electroencephalography (EEG) data requires the computation of the scalp potential induced by current sources in the brain. This so-called EEG forward problem is based on an accurate estimation of the volume conduction effects in the human head, represented by a partial differential equation which can be solved using the finite element method (FEM). FEM offers flexibility when modeling anisotropic tissue conductivities but requires a volumetric discretization, a mesh, of the head domain. Structured hexahedral meshes are easy to create in an automatic fashion, while tetrahedral meshes are better suited to model curved geometries. Tetrahedral meshes, thus, offer better accuracy but are more difficult to create. Methods: We introduce CutFEM for EEG forward simulations to integrate the strengths of hexahedra and tetrahedra. It belongs to the family of unfitted finite element methods, decoupling mesh and geometry representation. Following a description of the method, we will employ CutFEM in both controlled spherical scenarios and the reconstruction of somatosensory-evoked potentials. Results: CutFEM outperforms competing FEM approaches with regard to numerical accuracy, memory consumption, and computational speed while being able to mesh arbitrarily touching compartments. Discussion: CutFEM balances numerical accuracy, computational efficiency, and a smooth approximation of complex geometries that has previously not been available in FEM-based EEG forward modeling.

EEG-based source localization with enhanced virtual aperture using second order statistics.

For the past few decades source localization, based on EEG modality, has been a very active area of research. EEG signal provides temporal resolution in millisecond range that can capture rapidly changing patterns of brain activity but it has a low spatial resolution as compared to techniques like fMRI, PET, CT scan, etc. So, one of the motives of this research is to improve the spatial resolution of the EEG signal. Many successful attempts have been made to localise the active neural sources using EEG signals with the introduction of techniques like MNE, LORETA, sLORETA, FOCUSS, etc. But these techniques require a large number of electrodes for correct localization of a few sources. This paper aims at providing a new method for the localization of EEG sources with a fewer electrode. This is achieved by exploiting the second-order statistics to enhance the aperture and solve the EEG localization problem. The comparison of the proposed method with the state-of-the-art methods is done by observing the localization error with variation in SNR, number of snapshots (time samples), number of active sources, and number of electrodes. The results show that the proposed method can detect a greater number of sources with fewer electrodes and with higher accuracy as compared to methods available in the literature. Real -time EEG signal during an arithmetic task is considered and the proposed algorithm clearly shows a sparse activity in the frontal region.

Impact of skull sutures, spongiform bone distribution, and aging skull conductivities on the EEG forward and inverse problems.

Objective. Source imaging is a principal objective for electroencephalography (EEG), the solutions of which require forward problem (FP) computations characterising the electric potential distribution on the scalp due to known sources. Additionally, the EEG-FP is dependent upon realistic, anatomically correct volume conductors and accurate tissue conductivities, where the skull is particularly important. Skull conductivity, however, deviates according to bone composition and the presence of adult sutures. The presented study therefore analyses the effect the presence of adult sutures and differing bone composition have on the EEG-FP and inverse problem (IP) solutions.Approach. Utilising a well-established head atlas, detailed head models were generated including compact and spongiform bone and adult sutures. The true skull conductivity was considered as inhomogeneous according to spongiform bone proportion and sutures. The EEG-FP and EEG-IP were solved and compared to results employing homogeneous skull models, with varying conductivities and omitting sutures, as well as using a hypothesised aging skull conductivity model.Main results. Significant localised FP errors, with relative error up to 85%, were revealed, particularly evident along suture lines and directly related to the proportion of spongiform bone. This remained evident at various ages. Similar EEG-IP inaccuracies were found, with the largest (maximum 4.14 cm) across suture lines.Significance. It is concluded that modelling the skull as an inhomogeneous layer that varies according to spongiform bone proportion and includes differing suture conductivity is imperative for accurate EEG-FP and source localisation calculations. Their omission can result in significant errors, relevant for EEG research and clinical diagnosis.

Improvement of performance and sensitivity of 2D and 3D image reconstruction in EIT using EFG forward model.

This paper presents a pure element-free Galerkin method (EFGM) forward model for image reconstruction in 2D and 3D electrical impedance tomography (EIT) using an adaptive current injection method. In EIT systems with the adapting current injection method, both static and dynamic images can be reconstructed; however, determination of electrode contact impedances in the complete electrode model is difficult and the Gap model is used. In this paper, in the EIT forward problem a weak form functional based on the Gap model and a pure EFGM approach are developed, and in the EIT inverse problem, Jacobian matrix is computed by the EFGM, and a fast integration technique is introduced to calculate the entries of the Jacobian matrix within an adequate computation time. The influence of increasing the density of nodes at and near the electrodes with steep electric potential gradients on the accuracy of FEM and EFGM forward solutions is investigated, and the performance of the image reconstruction algorithm with the proposed fast integration technique is examined. The numerical results reveal that the proposed EFGM forward model with the fast integration technique has an efficient performance both in terms of mean relative imaging errors and computational time.

The electrocardiographic forward problem: A benchmark study.

BACKGROUND: Electrocardiographic forward problems are crucial components for noninvasive electrocardiographic imaging (ECGI) that compute torso potentials from cardiac source measurements. Forward problems have few sources of error as they are physically well posed and supported by mature numerical and computational techniques. However, the residual errors reported from experimental validation studies between forward computed and measured torso signals remain surprisingly high. OBJECTIVE: To test the hypothesis that incomplete cardiac source sampling, especially above the atrioventricular (AV) plane is a major contributor to forward solution errors. METHODS: We used a modified Langendorff preparation suspended in a human-shaped electrolytic torso-tank and a novel pericardiac-cage recording array to thoroughly sample the cardiac potentials. With this carefully controlled experimental preparation, we minimized possible sources of error, including geometric error and torso inhomogeneities. We progressively removed recorded signals from above the atrioventricular plane to determine how the forward-computed torso-tank potentials were affected by incomplete source sampling. RESULTS: We studied 240 beats total recorded from three different activation sequence types (sinus, and posterior and anterior left-ventricular free-wall pacing) in each of two experiments. With complete sampling by the cage electrodes, all correlation metrics between computed and measured torso-tank potentials were above 0.93 (maximum 0.99). The mean root-mean-squared error across all beat types was also low, less than or equal to 0.10 mV. A precipitous drop in forward solution accuracy was observed when we included only cage measurements below the AV plane. CONCLUSION: First, our forward computed potentials using complete cardiac source measurements set a benchmark for similar studies. Second, this study validates the importance of complete cardiac source sampling above the AV plane to produce accurate forward computed torso potentials. Testing ECGI systems and techniques with these more complete and highly accurate datasets will improve inverse techniques and noninvasive detection of cardiac electrical abnormalities.

A New Method of Haemorrhagic Stroke Detection Via Deep Magnetic Induction Tomography.

Hemorrhage imaging is one of the most common applications of magnetic induction tomography (MIT). Depth and the mass of stroke stimulated (MSS) are the most important issues that need to be solved for this application. Transcranial magnetic stimulation (TMS) is a technique belonging to the deep brain stimulation (DBS) field, which aims at overcoming human diseases such as depression. TMS coils, namely, circular, figure-8, and H-coils, play an important role in TMS. Among these, H-coils individually focus on the issues of achieving effective stimulation of deep region. MIT and TMS mechanisms are similar. Herein, for the first time, improved TMS coils, including figure-8 and H-coils, are applied as MIT excitation coils to study the possibility of achieving the mass of stroke stimulated and deep detection through MIT. In addition, the configurations of the detection coils are varied to analyze their influence and determine the optimal coils array. Finally, MIT is used to detect haemorrhagic stroke occurring in humans, and the application of deep MIT to the haemorrhagic stroke problem is computationally explored. Results show that among the various coils, the improved H-coils have MSS and depth characteristics that enable the detection of deep strokes through MIT. Although the detecting depth of the figure-8 coil is weaker, its surface signal is good. The deep MIT technique can be applied to haemorrhagic detection, providing a critical base for deeper research.

Computational imaging of the cardiac activities using magnetocardiography.

The electrical impulses of the heart will generate a tiny magnetic field outside the thorax that is measured as Magnetocardiographic signals. The challenging study is to estimate the cardiac activities in terms of depolarisation and repolarization maps from the measured signals called as inverse problem. This is computed only if one has solved generic or subject- specific prior models using the anatomical structures of the myocardium, the torso and the detectors called as forward problem. In this study, the Discretised heart is priorily assumed as the dipolar sources forming a double layer. The thorax structure modelled with finite element meshes is considered in the forward study. The magnetocardiographic data are simulated using uniform double layer model representing transmembrane distribution on the epicardium and endocardium. Using this data, the activation maps are non-invasively imaged on the heart surface using Tikhonov's regularisation technique. The inverse study is extended to reconstruct the depolarisation sequences of the abnormal cases.