The motion of and interaction between phase singularities that lie at the centers of spiral waves captures many qualitative and, in some cases, quantitative features of complex dynamics in excitable systems. Being able to accurately reconstruct their position is thus quite important, even if the data are noisy and sparse, as in electrophysiology studies of cardiac arrhythmias, for instance. A recently proposed global topological approach [Marcotte & Grigoriev, Chaos 27, 093936 (2017)] promises to meaningfully improve the quality of the reconstruction compared with traditional, local approaches. Indeed, we found that this approach is capable of handling noise levels exceeding the range of the signal with minimal loss of accuracy. Moreover, it also works successfully with data sampled on sparse grids with spacing comparable to the mean separation between the phase singularities for complex patterns featuring multiple interacting spiral waves.
Keywords: phase singularity, cardiac arrhythmia, spiral wave chaosCatheter ablation has recently emerged as a leading medical treatment for a range of cardiac arrhythmias, especially atrial fibrillation. The premise of the treatment is that certain localized regions of heart tissue can become sources of spiral excitation waves -or rotors -competing with the heart's natural pacemaker, i.e., the sinoatrial node for the atria or the atrio-ventricular node for the ventricles. The success of the ablation procedure then critically depends on the precision with which these sources are located based on electrograms obtained using intra-cardiac multielectrode catheters. This paper explains how the sources of excitation waves in a numerical model of atrial fibrillation can be reliably located with subgrid precision using sparse and noisy measurements of the transmembrane voltage. A similar approach could be used to improve the quality of rotor mapping in a clinical setting.