Research software has become a central asset in academic research. It optimizes existing and enables new research methods, implements and embeds research knowledge, and constitutes an essential research product in itself. Research software must be sustainable in order to understand, replicate, reproduce, and build upon existing research or conduct new research effectively. In other words, software must be available, discoverable, usable, and adaptable to new needs, both now and in the future. Research software therefore requires an environment that supports sustainability. Hence, a change is needed in the way research software development and maintenance are currently motivated, incentivized, funded, structurally and infrastructurally supported, and legally treated. Failing to do so will threaten the quality and validity of research. In this paper, we identify challenges for research software sustainability in Germany and beyond, in terms of motivation, selection, research software engineering personnel, funding, infrastructure, and legal aspects. Besides researchers, we specifically address political and academic decision-makers to increase awareness of the importance and needs of sustainable research software practices. In particular, we recommend strategies and measures to create an environment for sustainable research software, with the ultimate goal to ensure that software-driven research is valid, reproducible and sustainable, and that software is recognized as a first class citizen in research. This paper is the outcome of two workshops run in Germany in 2019, at deRSE19 - the first International Conference of Research Software Engineers in Germany - and a dedicated DFG-supported follow-up workshop in Berlin.
OBJECTIVE: Unipolar intracardiac electrograms (uEGMs) measured inside the atria during electro-anatomic mapping contain diagnostic information about cardiac excitation and tissue properties. The ventricular far field (VFF) caused by ventricular depolarization compromises these signals. Current signal processing techniques require several seconds of local uEGMs to remove the VFF component and thus prolong the clinical mapping procedure. We developed an approach to remove the VFF component using data obtained during initial anatomy acquisition. METHODS: We developed two models which can approximate the spatio-temporal distribution of the VFF component based on acquired EGM data: Polynomial fit, and dipole fit. Both were benchmarked based on simulated cardiac excitation in two models of the human heart and applied to clinical data. RESULTS: VFF data acquired in one atrium were used to estimate model parameters. Under realistic noise conditions, a dipole model approximated the VFF with a median deviation of 0.029mV, yielding a median VFF attenuation of 142. In a different setup, only VFF data acquired at distances of more than 5mm to the atrial endocardium were used to estimate the model parameters. The VFF component was then extrapolated for a layer of 5mm thickness lining the endocardial tissue. A median deviation of 0.082mV (median VFF attenuation of 49x) was achieved under realistic noise conditions. CONCLUSION: It is feasible to model the VFF component in a personalized way and effectively remove it from uEGMs. SIGNIFICANCE: Application of our novel, simple and computationally inexpensive methods allows immediate diagnostic assessment of uEGM data without prolonging data acquisition.
Catheter ablation is a curative therapeutic approach for atrial fibrillation (AF). Ablation of rotational sources based on basket catheter measurements has been proposed as a promising approach in patients with persistent AF to complement pulmonary vein isolation. However, clinically reported success rates are equivocal calling for a mechanistic investigation under controlled conditions. We present a computational framework to benchmark ablation strategies considering the whole cycle from excitation propagation to electrogram acquisition and processing to virtual therapy. Fibrillation was induced in a patient-specific 3D volumetric model of the left atrium, which was homogeneously remodelled to sustain reentry. The resulting extracellular potential field was sampled using models of grid catheters as well as realistically deformed basket catheters considering the specific atrial anatomy. Virtual electrograms were processed to compute phase singularity density maps to target rotor tips with up to three circular ablations. Stable rotors were successfully induced in different regions of the homogeneously remodelled atrium showing that rotors are not constrained to unique anatomical structures or locations. Phase singularity density maps correctly identified and located the rotors (deviation < 10 mm) based on catheter recordings only for sufficient resolution (inter-electrode distance = 3 mm) and proximity to the wall (< 10 mm). Targeting rotor sites with ablation did not stop reentries in the homogeneously remodelled atria independent from lesion size (1-7 mm radius), from linearly connecting lesions with anatomical obstacles, and from the number of rotors targeted sequentially (up to 3). Our results show that phase maps derived from intracardiac electrograms can be a powerful tool to map atrial activation patterns, yet they can also be misleading due to inaccurate localization of rotor tips depending on electrode resolution and distance to the wall. This should be considered to avoid ablating regions that are in fact free of rotor sources of AF. In our experience, ablation of rotor sites was not successful to stop fibrillation. Our comprehensive simulation framework provides the means to holistically benchmark ablation strategies in silico under consideration of all steps invol
Conference Contributions (10)
L. A. Unger, A. Luik, A. Haas, and O. Dössel. Comparison of Morphology-Based and Delay-Based Measures for Reference Beat Classification during Atrial Tachycardia. In Computing in Cardiology Conference (CinC), 2020
Beat acceptance and rejection during atrial tachycardia are crucial for the compilation of meaningful electroanatomical maps during an electrophysiological study. State of the art methods compare the delays in activation time between two or more electrograms recorded with electrodes of a spatially stable reference catheter. This work introduces morphology-based measures for beat selection in the context of mapping atrial tachycardia. Active segments were extracted from bipolar reference electrograms with the help of the non-linear energy operator. After prealignment by means of maximum cross-correlation, the correlation coefficient as well as the normalized 1-norm distance yielded a similarity measure for each pair of prealigned active segments. The morphology-based measures were then compared to the delay-based measure. In an exemplary patient with 5163 recorded beats, the delay-based measures were strongly dependent on the accuracy of the local activation times as well as on the selection of reference leads. The morphology-based measures emphasized changes in the target tachycardia which were not detectable by the delay-based method. The correlation and the distance measure showed similar behavior but stressed different aspects of morphological changes. Ventricular components in active segments caused minor changes in morphology which were also reflected in the morphology-based measures. The morphology-based measures introduced in this work enhanced beat selection in the exemplary patient. A follow-up study with a representative patient cohort needs to quantify the improvement across patients and translate the measure to clinical practice. A combination of activation delays and morphological similarity is strongly expected to exploit the advantages of both methods for beat selection.
Intracardiac electrograms (EGMs) form the basis for the diagnosis of arrhythmia mechanisms. Bipolar EGMs dominate clinical practice despite major disadvantages over unipolar EGMs since noise strongly distorts the latter. In this study, we quantified and reduced the noise level of uni- and bipolar EGMs recorded with Rhythmia HDx and the Orion catheter. Distinct noise frequencies in the power spectral density (PSD) were detected with a sliding win- dow of variable width and subsequently removed by notch filtering. The absolute peak to peak voltage remaining in the inactive segments after baseline removal quantified the noise level of the system. An international, multi-center selection of 33 patients served as a broad sample cohort. The case-specific detection and removal of noise peaks reduced the noise level in unipolar EGMs by 30% to 0.076 mV compared to standard clinical filtering. With a bipolar noise level of 0.01 mV, we saw that Rhythmia HDx meets the low noise floor claimed in the system specifica- tions. Certain noise frequencies presented permanently in all cases whereas others showed up only intermittently or in individual cases. The suggested extension of filter settings lowers the noise level, enhances the detailed segmentation of low volt- age areas, and encourages to exploit the advantages of unipolar over bipolar EGMs in clinical practice.
Acquiring adequate mapping data in patients with atrial fibrillation is still one of the main obstacles in the treatment of this atrial arrhythmia. Due to the lack of catheters with both a panoramic field of view and sufficient electrode density for simultaneous mapping, electrophysiologists are forced to fall back on sequential mapping techniques. But, because activation patterns change rapidly during atrial fibrillation, they cannot be mapped sequentially. We propose that mapping tissue properties which are time independent, in contrast, allows a sequential approach. Here, we use the shortest measured electrogram cycle length to estimate the effective refractory period of the underlying tissue in a simulation study. Atrial fibrillation was simulated in a spherical model of the left atrium comprised of regions with varied refractory period. We found that the minimal measured electrogram cycle length correlates with the effective refractory period of the underlying tissue if the regions with distinct refractory properties are large enough and if the absolute difference in effective refractory periods is sufficient. This approach is capable of identifying regions of lowered effective refractory period without the need for cardioversion. Those regions are likely to harbor drivers of atrial fibrillation, which emphasizes the necessity of their localization.
L. A. Unger, M. Rottmann, G. Seemann, and O. Dössel. Detecting phase singularities and rotor center trajectories based on the Hilbert transform of intraatrial electrograms in an atrial voxel model. In Current Directions in Biomedical Engineering, vol. 1(1) , pp. 38-41, 2015
This work aimed at the detection of rotor centers within the atrial cavity during atrial fibrillation on the basis of phase singularities. A voxel based method was established which employs the Hilbert transform and the phase of unipolar electrograms. The method provides a 3D overview of phase singularities at the endocardial surface and within the blood volume. Mapping those phase singularities from the inside of the atria at the endocardium yielded rotor center trajectories.We discuss the results for an unstable and a more stable rotor. The side length of the areas covered by the trajectories varied from 1.5mm to 10 mm. These results are important for cardiologists who target rotors with RF ablation in order to cure atrial fibrillation.
C. Nagel, N. Pilia, L. Unger, and O. Dössel. Performance of Different Atrial Conduction Velocity Estimation Algorithms Improves with Knowledge about the Depolarization Pattern. In Current Directions in Biomedical Engineering, vol. 5(1) , pp. 101-104, 2019
Quantifying the atrial conduction velocity (CV) reveals important information for targeting critical arrhythmia sites that initiate and sustain abnormal electrical pathways, e.g. during atrial flutter. The knowledge about the local CV distribution on the atrial surface thus enhances clinical catheter ablation procedures by localizing pathological propagation paths to be eliminated during the intervention. Several algorithms have been proposed for estimating the CV. All of them are solely based on the local activation times calculated from electroanatomical mapping data. They deliver false values for the CV if applied to regions near scars or wave collisions. We propose an extension to all approaches by including a distinct preprocessing step. Thereby, we first identify scars and wave front collisions and provide this information for the CV estimation algorithm. In addition, we provide reliable CV values even in the presence of noise. We compared the performance of the Triangulation, the Polynomial Fit and the Radial Basis Functions approach with and without the inclusion of the aforementioned preprocessing step. The evaluation was based on different activation patterns simulated on a 2D synthetic triangular mesh with different levels of noise added. The results of this study demonstrate that the accuracy of the estimated CV does improve when knowledge about the depolarization pattern is included. Over all investigated test cases, the reduction of the mean velocity error quantified to at least 25 mm/s for the Radial Basis Functions, 14 mm/s for the Polynomial Fit and 14 mm/s for the Triangulation approach compared to their respective implementations without the preprocessing step. Given the present results, this novel approach can contribute to a more accurate and reliable CV estimation in a clinical setting and thus improve the success of radio-frequency ablation to treat cardiac arrhythmias.
The outcomes of ablation targeting either reentry activations or fractionated activity during persistent atrial fibrillation (AF) therapy remain suboptimal due to, among others, the intricate underlying AF dynamics. In the present work, we sought to investigate such AF dynamics in a heterogeneous simulation setup using recurrence quantification analysis (RQA). AF was simulated in a spherical model of the left atrium, from which 412 unipolar atrial electrograms (AEGs) were extracted (2 s duration; 5 mm spacing). The phase was calculated using the Hilbert transform, followed by the identification of points of singularity (PS). Three regions were defined according to the occurrence of PSs: 1) no rotors; 2) transient rotors and; 3) long-standing rotors. Bipolar AEGs (1114) were calculated from pairs of unipolar nodes and bandpass filtered (30-300 Hz). The CARTO criterion (Biosense Webster) was used for AEGs classification (normal vs. fractionated). RQA attributes were calculated from the filtered bipolar AEGs: determinism (DET); recurrence rate (RR); laminarity (LAM). Sample entropy (SampEn) and dominant frequency (DF) were also calculated from the AEGs. Regions with longstanding rotors have shown significantly lower RQA attributes and SampEn when compared to the other regions, suggesting a higher irregular behaviour (P≤0.01 for all cases). Normal and fractionated AEGs were found in all regions (respectively; Region 1: 387 vs. 15; Region 2: 221 vs. 13; Region 3: 415 vs. 63). Region 1 vs. Region 3 have shown significant differences in normal AEGs (P≤0.0001 for all RQA attributes and SampEn), and significant differences in fractionated AEGs for LAM, RR and SampEn (P=0.0071, P=0.0221 and P=0.0086, respectively). Our results suggest the co-existence of normal and fractionated AEGs within long-standing rotors. RQA has unveiled distinct dynamic patterns–irrespective of AEGs classification–related to regularity structures and their nonstationary behaviour in a rigorous deterministic context.
Cardiac arrhythmias such as atrial fibrillation occur frequently in industrialized countries. Radiofrequency ablation (RFA) is a standard treatment if drug therapy fails. This minimally invasive surgery aims at stabilizing the heart rhythm on a permanent basis. However, the procedure commonly needs to be repeated because of the high recurrence rate of arrhythmias. Non-transmural lesions as well as gaps within linear lesions are among the main problems during the RFA. The assessment of lesion formation is not adequate in state of the art procedures. Therefore, the aim of this study is to investigate the short-term reversibility of lesions using human electrograms recorded by a high-density mapping system during an electrophysiological study (EPS). A predefined measurement protocol was executed during the EPS in order to create three ablation points in the left atrium. Subsequently, after preprocessing the recorded signals, electrogram (EGM) paths were formed along the endocardial surface of the atrium. By analyzing changes of peak to peak amplitudes of unipolar EGMs before and after ablation, it was possible to distinguish lesion area and healthy myocardium. The peak to peak amplitudes of the EGMs decreased by 40-61% after 30 seconds of ablation. Furthermore, we analyzed the morphological changes of EGMs surrounding the lesion. High-density mapping data showed that not only the tissue, which had direct contact with the catheter tip during the RFA, but also the surrounding tissue was affected. This was demonstrated by low peak to peak amplitudes in large areas with a width of 14 mm around the center of the ablation lesion. After right pulmonary vein isolation, high-density mapping was repeated on the previous lesions. The outer region of RFA-treated tissue appears to recover as opposed to the central core of the ablation point. This observation suggests that the meaningfulness of an immediate remap after ablation during an EPS may lead the physician to false conclusions.
O. Dössel, T. Oesterlein, L. Unger, A. Loewe, C. Schmitt, and A. Luik. Spatio-temporal Analysis of Multichannel Atrial Electrograms Based on a Concept of Active Areas. In Annual International Conference of the IEEE Engineering in Medicine and Biology Society. IEEE Engineering in Medicine and Biology Society. Annual Conference, vol. 2018, pp. 490-493, 2018
Atrial tachycardia and atrial flutter are frequent arrhythmia that occur spontaneously and after ablation of atrial fibrillation. Depolarization waves that differ significantly from sinus rhythm propagate across the atria with high frequency (typically 140 to 220 beats per minute). A detailed and personalized analysis of the spread of depolarization is imperative for a successful ablation therapy. Thus, catheters with several electrodes are employed to measure multichannel electrograms inside the atria. Here we propose a new concept for spatio-temporal analysis of multichannel electrograms during atrial tachycardia and atrial flutter. It is based on the calculation of simultaneously active areas. The method allows to identify atrial tachycardia and to automatically distinguish between subtypes of focal activity, micro-reentry and macro-reentry.
Atrial fibrillation is a common irregular heart rhythm. Until today there is still a need for research to quantify typical signal characteristics of rotors, which can induce atrial fibrillation. In this work, signal characteristics of a stable and a more unstable rotor in a realistic heart model including fiber orientation were analyzed with the following methods: peak-to-peak amplitude, Hilbert phase, approximate entropy and RS-difference. In this simulation model the stable rotor rotated with a cycle length of 145 ms and stayed in an area of 1.5 mm x 3 mm. Another more unstable rotor with a cycle length of 190 ms moved in an area of 10 mm × 4 mm. In a distance of 2 mm to the rotor tip, the peak-to-peak amplitude decreased significantly, whereas the RS-difference and the approximate entropy were maximal. The rotor center trajectories were detected by phase singularity points determined by the Hilbert transform. We showed that more unstable rotors resulted in more amplitude changes over time and also the cycle length differed more. Furthermore, we presented typical activation time patterns of the Lasso catheter centered at the rotor tip and in different distances to the rotor tip. We suggest that cardiologists use a combination of the described methods to determine a rotor tip position in a more robust manner.
Student Theses (2)
L. A. Unger. Substrate Mapping During Atrial Fibrillation - A Combined in Silico and Clinical Test of Concept Study. Institute of Biomedical Engineering, Karlsruhe Institute of Technology (KIT). Masterarbeit. 2017
Acquiring adequate mapping data in patients with atrial fibrillation (AFib) is one of the main obstacles in the treatment of this arrhythmia. Due to the lack of catheters with both a panoramic field of view and sufficient electrode density for simultaneous mapping, electrophysiologists are forced to fall back on sequential mapping techniques to identify activation patterns. However, this approach is insufficient for rapidly changing patterns as they typically occur during AFib. In contrast to activation time mapping, substrate mapping avoids this drawback by analyzing time independent tissue properties. While most of the existing methods for substrate mapping do not reflect actual tissue properties but certain electrogram features, the results suffer from dependencies on parameters of data analysis or are limited to harmonic signals. Here, we investigate the potential and limitations of measured electrogram cycle lengths to derive information about the effective refractory period of the underlying tissue during sequences of AFib. Following theoretical considerations, areas with decreased effective refractory period are likely to harbor AFib drivers.In a first step, different parametrizations of the Courtemanche-Ramirez-Nattel model with varying ERP served as substrates for bidomain simulations of AFib in a spherical model of the left atrium. Circular regions of deviating effective refractory period with radii between 9 mm and 19 mm were imposed. We found that the minimal measured electrogram cycle length correlates with the effective refractory period of the underlying tissue if the regions with distinct refractory properties are large enough and if the absolute difference in effective refractory periods is significant. Rhythms of high complexity which cause many other mapping approaches to come up against limiting factors favor the here introduced method as statistics profit from the variability in observed events.In a second step, the clinical feasibility of using measured cycle length statistics to conclude on the underlying ERP of the tissue was investigated. Lacking ground truth data providing reliable in vivo information on the ERP of the tissue, final validations of our hypotheses remain due. The 25 % quantile of cycle lengths was used rather than the minimum in favor of improved robustness in clinical application. The cycle length analysis in patient data yielded physiologically reasonable results with locally low gradients but globally differing statistics. Both the in silico and the clinical test of concept study suggested that applying statistical measures such as the 25 % quantile to measured electrogram cycle lengths is capable of revealing information on the effective refractory period of the underlying tissue. This, in turn, is of particular clinical interest, as regions of lowered effective refractory period are likely to harbor AFib drivers.
L. A. Unger. Simulation based Estimation of Parameters for Reentries in Human Atria. Institute of Biomedical Engineering, Karlsruhe Institute of Technology (KIT). Bachelorarbeit. 2014