Objective: The severity of atrial fibrillation (AF) can be assessed from intra-operative epicardial measurements (high-resolution electrograms), using metrics such as conduction block (CB) and continuous conduction delay and block (cCDCB). These features capture differences in conduction velocity and wavefront propagation, but ignore complementary properties such as the morphology of the action potentials. In this work, we focus on such complementary properties, and derive features to detect variations in the atrial potential waveforms. Methods: We show that the spatial variation of atrial potential morphology during a single beat may be described by changes in the singular values of the epicardial measurement matrix. The method is non-parametric and requires little preprocessing. A corresponding singular value map points at areas subject to fractionation and block. Further, we developed an experiment where we simultaneously measure electrograms (EGMs) and a multi-lead ECG. Results: The captured data showed that the normalized singular values of the heartbeats during AF are higher than during SR, and that this difference is more pronounced for the (non-invasive) ECG data than for the EGM data, if the electrodes are positioned at favorable locations. Conclusion: Overall, the singular value-based features are a useful indicator to detect and evaluate AF. Significance: The proposed method might be beneficial for identifying electropathological regions in the tissue without estimating the local activation time.

@en

In wireless networks, an essential step for precise range-based localization is the high-resolution estimation of multipath channel delays. The resolution of traditional delay estimation algorithms is inversely proportional to the bandwidth of the training signals used for channel probing. Considering that typical training signals have limited bandwidth, delay estimation using these algorithms often leads to poor localization performance. To mitigate these constraints, we exploit the multiband and carrier frequency switching capabilities of wireless transceivers and propose to acquire channel state information (CSI) in multiple bands spread over a large frequency aperture. The data model of the acquired measurements has a multiple shift-invariance structure, and we use this property to develop a high-resolution delay estimation algorithm. We derive the Cramér-Rao Bound (CRB) for the data model and perform numerical simulations of the algorithm using system parameters of the emerging IEEE 802.11be standard. Simulations show that the algorithm is asymptotically efficient and converges to the CRB. To validate modeling assumptions, we test the algorithm using channel measurements acquired in real indoor scenarios. From these results, it is seen that delays (ranges) estimated from multiband CSI with a total bandwidth of 320 MHz show an average RMSE of less than 0.3 ns (10 cm) in 90% of the cases.

@en

We aim to provide a critical appraisal of basic concepts underlying signal recording and processing technologies applied for (i) atrial fibrillation (AF) mapping to unravel AF mechanisms and/or identifying target sites for AF therapy and (ii) AF detection, to optimize usage of technologies, stimulate research aimed at closing knowledge gaps, and developing ideal AF recording and processing technologies. Recording and processing techniques for assessment of electrical activity during AF essential for diagnosis and guiding ablative therapy including body surface electrocardiograms (ECG) and endo- or epicardial electrograms (EGM) are evaluated. Discussion of (i) differences in uni-, bi-, and multi-polar (omnipolar/Laplacian) recording modes, (ii) impact of recording technologies on EGM morphology, (iii) global or local mapping using various types of EGM involving signal processing techniques including isochronal-, voltage- fractionation-, dipole density-, and rotor mapping, enabling derivation of parameters like atrial rate, entropy, conduction velocity/direction, (iv) value of epicardial and optical mapping, (v) AF detection by cardiac implantable electronic devices containing various detection algorithms applicable to stored EGMs, (vi) contribution of machine learning (ML) to further improvement of signals processing technologies. Recording and processing of EGM (or ECG) are the cornerstones of (body surface) mapping of AF. Currently available AF recording and processing technologies are mainly restricted to specific applications or have technological limitations. Improvements in AF mapping by obtaining highest fidelity source signals (e.g. catheter-electrode combinations) for signal processing (e.g. filtering, digitization, and noise elimination) is of utmost importance. Novel acquisition instruments (multi-polar catheters combined with improved physical modelling and ML techniques) will enable enhanced and automated interpretation of EGM recordings in the near future.

@en

Atrial fibrillation (AF) is the most common arrhythmia. Although the exact cause is unclear, electropathology of atrial tissue is one contributing factor. Electropathological characteristics derived from intra-operative epicardial measurements, such as conduction block (CB) and continues conduction delay and block (cCDCB), can be used to assess the severity of AF. In sinus rhythm, however, these parameters do not indicate significant difference between different development stages of AF, such as paroxysmal and persistent AF. Therefore, we propose a methodology to improve AF severity detection using intra-operative electrograms. We propose a model that describes the spatial diversity of atrial potential waveforms during a single beat on the multi-channel electrograms. Based on this model, we derive two novel features. During sinus rhythm, we used 293 beats from patients with a history of paroxysmal or persistent AF. Using a random forest classifier, we achieved 78.42% classification accuracy, while classification based on the CB and cCDCB leads to an accuracy of 58.34%.

@en

Atrial fibrillation (AF) is the most sustained arrhythmia in the heart and also the most common complication developed after cardiac surgery. Due to its progressive nature, timely detection of AF is important. Currently, physicians use a surface electrocardiogram (ECG) for AF diagnosis. However, when the patient develops AF, its various development stages are not distinguishable for cardiologists based on visual inspection of the surface ECG signals. Therefore, severity detection of AF could start from differentiating between short-lasting AF and long-lasting AF. Here, de novo post-operative AF (POAF) is a good model for short-lasting AF while long-lasting AF can be represented by persistent AF. Therefore, we address in this paper a binary severity detection of AF for two specific types of AF. We focus on the differentiation of these two types as de novo POAF is the first time that a patient develops AF. Hence, comparing its development to a more severe stage of AF (e.g., persistent AF) could be beneficial in unveiling the electrical changes in the atrium. To the best of our knowledge, this is the first paper that aims to differentiate these different AF stages. We propose a method that consists of three sets of discriminative features based on fundamentally different aspects of the multi-channel ECG data, namely based on the analysis of RR intervals, a greyscale image representation of the vectorcardiogram, and the frequency domain representation of the ECG. Due to the nature of AF, these features are able to capture both morphological and rhythmic changes in the ECGs. Our classification system consists of a random forest classifier, after a feature selection stage using the ReliefF method. The detection efficiency is tested on 151 patients using 5-fold cross-validation. We achieved 89.07% accuracy in the classification of de novo POAF and persistent AF. The results show that the features are discriminative to reveal the severity of AF. Moreover, inspection of the most important features sheds light on the different characteristics of de novo post-operative and persistent AF.

@en

Atrial Fibrillation (AF) is the most sustained arrhythmia in the heart. On the surface electrocardiogram (ECG), AF is characterised by the irregular RR intervals and by fibrillatory waves or the absence of a P wave. Since AF is a progressive disease, timely and correct detection is crucial for AF treatment. Detailed insight into the areas of arrhythmia-related electropathology can be obtained by analyzing high-resolution (inter-electrode distance 2mm) electrograms (EGMs). However, these measurements are rather invasive. By integration of high-resolution epicardial mapping data and surface ECG data, we hope to learn how different stages of AF represent themselves on the ECG. Eventually this can help to guide to identify areas of electropathology as target sites of ablation therapy on the less invasive ECG. A first step in this direction is to learn how to reconstruct the ECG based on EGM measurements. In practice, however, EGMs are only measured at few atrial locations, not covering the complete atria. An important question therefore is: How can we reconstruct ECG based on the observations from a limited part of the heart? To answer this question, we propose two methods. In the first method, we increase the number of observations from a part of the right atrium (RA) to the whole RA by synchronizing EGMs that are measured at different moments in time based on the local activation time (LAT). In the second method, under the assumption that atrial EGMs measured at different spatial areas are linearly related, the conductivity matrix is estimated for the whole atrium which enables us to reconstruct the ECGs from the limited observations. The second method brings twofold benefits. First, the conductivity matrix can be used as a novel diagnostic tool to detect AF as well as areas of electropathology. Second, it provides a practical solution to reconstruct epicardial potentials from ECGs, non- This research was funded in part by the Medical Delta Cardiac Arrhythmia Lab (CAL), the Netherlands. invasively. The results show that method one increases the reconstruction accuracy. Furthermore, the conductivity matrix reveals the structural differences between sinus rhythm (SR) and AF episodes which could be the first step to interpret the underlying electropathology of AF @en

The presence of rich scattering in indoor and urban radio propagation scenarios may cause a high arrival density of multipath components (MPCs). Often the MPCs arrive in clusters at the receiver, where MPCs within one cluster have similar angles and delays. The MPCs arriving within a single cluster are typically unresolvable in the delay domain. In this paper, we analyze the effects of unresolved MPCs on the bias of the delay estimation with a multiband subspace fitting algorithm. We treat the unresolved MPCs as a model error that results in perturbed subspace estimation. Starting from the first-order approximation of the perturbations, we derive the bias of the delay estimate of the line-of-sight (LOS) component. We show that it depends on the power and relative delay of the unresolved MPCs in the first cluster compared to the LOS component. Numerical experiments are included to show that the derived expression for the bias well describes the effects of unresolved MPCs on the delay estimation.@en

Objective: Over the last two decades, radar-based contactless monitoring of vital signs (heartbeat and respiration rate) has raised increasing interest as an emerging and added value to health care. However, until now, the flaws caused by indoor multipath propagation formed a fundamental hurdle for the adoption of such technology in practical healthcare applications where reliability and robustness are crucial. Multipath reflections, originated from one person, combine with the direct signals and multipaths of other people and stationary objects, thus jeopardizing individual vital signs extraction and localization. This work focuses on tackling indoor multipath propagation. Methods: We describe a methodology, based on accurate models of the indoor multipaths and of the radar signals, that enables separating the undesired multipaths from desired signals of multiple individuals, removing a key obstacle to real-world contactless vital signs monitoring and localization. Results: We also demonstrated it by accurately measure individual heart rates, respiration rates, and absolute distances (range information) of paired volunteers in a challenging real-world office setting. Conclusion: The approach, validated using a frequency-modulated continuous wave (FMCW) radar, was shown to function in an indoor environment where radar signals are severely affected by multipath reflections. Significance: Practical applications arise for health care, assisted living, geriatric and quarantine medicine, rescue and security purposes.

@en

In Principal Component Analysis (PCA), the dimension of the signal subspace is detected by counting the number of eigenvalues of a covariance matrix that are above a threshold. Random matrix theory provides accurate estimates for this threshold if the underlying data matrix has independent identically distributed columns. However, in time series analysis, the underlying data matrix has a Hankel or Toeplitz structure, and the columns are not independent. Using an empirical approach, we observe that the largest eigenvalue is fitted well by a Generalized Extreme Value (GEV) distribution, and we obtain accurate estimates for the thresholds to be used in a sequential rank detection test. In contrast to AIC or MDL, this provides a parameter that controls the probability of false alarm. Also a lower bound is presented for the rank detection rate of threshold-based detection for rank-1 problems. @en

Atrial fibrillation (AF) is the most common arrhythmia in the heart. Two main types of AF are defined as paroxysmal and persistent. In this paper, we present a method to discriminate between the characteristics of paroxysmal and persistent using tensor decompositions of a multi-channel electrocardiogram (ECG) signal. For this purpose, ECG signals are segmented by applying a Hilbert transform on the thresholded signal. Dynamic time warping is used to align the separated segments of each channel and then a tensor is constructed with three dimensions as time, heartbeats and channels. A Canonical polyadic decomposition with rank 2 is computed from this tensor and the resulting loading vectors describe the characteristics of paroxysmal and persistent AF in these three dimensions. The time loading vector reveals the pattern of a single P wave or abnormal AF patterns. The heartbeat loading vector shows whether the pattern is present or absent in a specific beat. The results can be used to distinguish between the patterns in paroxysmal AF and persistent AF. @en

The past two decades have witnessed a renewed interest in low frequency radio astronomy, with a particular focus on frequencies above 30 MHz e.g., LOFAR (LOw Frequency ARray) in the Netherlands and its European extension ILT, the International LOFAR Telescope. However, at frequencies below 30 MHz, Earth-based observations are limited due to a combination of severe ionospheric distortions, almost full reflection of radio waves below 10 MHz, solar eruptions and the radio frequency interference (RFI) of human-made signals. Moreover, there are interesting scientific processes which naturally occur at these low frequencies. A space or Lunar-based ultra-low-frequency (also referred to as ultra-long-wavelength, ULW) radio array would suffer significantly less from these limitations and hence would open up the last, virtually unexplored frequency domain in the electromagnetic spectrum. A roadmap has been initiated by astronomers and researchers in the Netherlands to explore the opportunity of building a swarm of satellites to observe at the frequency band below 30 MHz. This roadmap dubbed Orbiting Low Frequency Antennas for Radio Astronomy (OLFAR), a space-based ultra-low frequency radio telescope that will explore the Universe's so-called dark ages, map the interstellar medium, and study planetary and solar bursts in the solar system and search them in other planetary systems. Such a radio astronomy system will comprise of a swarm of hundreds to thousands of satellites, working together as a single aperture synthesis instrument deployed sufficiently far away from Earth to avoid terrestrial RFI. The OLFAR telescope is a novel and complex system, requiring yet to be proven engineering solutions. Therefore, a number of key technologies are still required to be developed and proven. The first step in this roadmap is the NCLE (Netherlands China Low Frequency Explorer) experiment, which was launched in May 2018 on the Chinese Chang'e 4 mission. The NCLE payload consists of a three monopole antenna system for low frequency observations, from which the first data stream is expected in the second half of 2019, which will provide important feedback for future science and technology opportunities. In this paper, the roadmap towards OLFAR, a brief overview of the science opportunities, and the technological and programmatic challenges of the mission are presented.

@en

In this paper, we focus on the problem of blind joint calibration of multiband transceivers and time-delay (TD) estimation of multipath channels. We show that this problem can be formulated as a particular case of covariance matching. Although this problem is severely ill-posed, prior information about radio-frequency chain distortions and multipath channel sparsity is used for regularization. This approach leads to a biconvex optimization problem, which is formulated as a rank-constrained linear system and solved by a simple group Lasso algorithm. Numerical experiments show that the proposed algorithm provides better calibration and higher resolution for TD estimation than current state-of-the-art methods.@en

Background: Local activation time (LAT) annotation in unipolar electrograms is complicated by interference from nonlocal atrial activities of neighboring tissue. This happens due to the spatial blurring that is inherent to electrogram recordings. In this study, we aim to exploit multi-electrode electrogram recordings to amplify the local activity in each electrogram and subsequently improve the annotation of LATs. Methods: An electrogram array can be modeled as a spatial convolution of per cell transmembrane currents with an appropriate distance kernel, which depends on the cells’ distances to the electrodes. By deconvolving the effect of the distance kernel from the electrogram array, we undo the blurring and estimate the underlying transmembrane currents as our desired local activities. However, deconvolution problems are typically highly ill-posed and result in unstable solutions. To overcome this issue, we propose to use a regularization term that exploits the sparsity of the first-order time derivative of the transmembrane currents. Results: We perform experiments on simulated two-dimensional tissues, as well as clinically recorded electrograms during paroxysmal atrial fibrillation. The results show that the proposed approach for deconvolution can improve the annotation of the true LAT in the electrograms. We also discuss, in summary, the required electrode array specifications for an appropriate recording and subsequent deconvolution. Conclusion: By ignoring small but local deflections, algorithms based on steepest descent are prone to generate smoother activation maps. However, by exploiting multi-electrode recordings, we can efficiently amplify small but local deflections and reveal new details in the activation maps that were previously missed.

@en

Reports on the activities of the Signal Processing Theory and Methods Technical Committee.@en

The multipath radio channel is considered to have a non-bandlimited channel impulse response. Therefore, it is challenging to achieve high resolution time-delay (TD) estimation of multipath components (MPCs) from bandlimited observations of communication signals. It this paper, we consider the problem of multiband channel sampling and TD estimation of MPCs. We assume that the nonideal multi-branch receiver is used for multiband sampling, where the noise is nonuniform across the receiver branches. The resulting data model of Hankel matrices formed from acquired samples has multiple shift-invariance structures, and we propose an algorithm for TD estimation using weighted subspace fitting. The subspace fitting is formulated as a separable nonlinear least squares (NLS) problem, and it is solved using a variable projection method. The proposed algorithm supports high resolution TD estimation from an arbitrary number of bands, and it allows for nonuniform noise across the bands. Numerical simulations show that the algorithm almost attains the Cramér Rao Lower Bound, and it outperforms previously proposed methods such as multiresolution TOA, MI-MUSIC, and ESPRIT.@en

Synchronization and ranging in internet of things (IoT) networks are challenging due to the narrowband nature of signals used for communication between IoT nodes. Recently, several estimators for range estimation using phase difference of arrival (PDoA) measurements of narrowband signals have been proposed. However, these estimators are based on data models which do not consider the impact of clock-skew on the range estimation. In this paper, clock-skew and range estimation are studied under a unified framework. We derive a novel and precise data model for PDoA measurements which incorporates the unknown clock-skew effects. We then formulate joint estimation of the clock-skew and range as a two-dimensional (2-D) frequency estimation problem of a single complex sinusoid. Furthermore, we propose: (i) a two-way communication protocol for collecting PDoA measurements and (ii) a weighted least squares (WLS) algorithm for joint estimation of clock-skew and range leveraging the shift invariance property of the measurement data. Finally, through numerical experiments, the performance of the proposed protocol and estimator is compared against the Cramér Rao lower bound demonstrating that the proposed estimator is asymptotically efficient.

@en

Finding the hidden parameters of the cardiac electrophysiological model would help to gain more insight on the mechanisms underlying atrial fibrillation, and subsequently, facilitate the diagnosis and treatment of the disease in later stages. In this work, we aim to estimate tissue conductivity from recorded electrograms as an indication of tissue (mal)functioning. To do so, we first develop a simple but effective forward model to replace the computationally intensive reaction-diffusion equations governing the electrical propagation in tissue. Using the simplified model, we present a compact matrix model for electrograms based on conductivity. Subsequently, we exploit the simplicity of the compact model to solve the ill-posed inverse problem of estimating tissue conductivity. The algorithm is demonstrated on simulated data as well as on clinically recorded data. The results show that the model allows to efficiently estimate the conductivity map. In addition, based on the estimated conductivity, realistic electrograms can be regenerated demonstrating the validity of the model.

@en

The estimation of the coordinates of nodes their proximity (or distance) measurements, is a principal challenge in numerous fields. Conventionally, when localizing a static network of immobile nodes, non-linear dimensionality reduction techniques are applied on the measured distances to obtain the relative coordinates up to a rotation and translation. In this article, we consider an anchorless network of mobile nodes, where the distance measurements between the mobile nodes are time-varying. In such an anchorless framework, where the absolute knowledge of any node position, motion or reference frame is absent, we aim to estimate the relative positions using the measured time-varying distances. To this end, we derive a data model which relates the time-varying distances to the time-varying relative positions of an anchorless network. Given this data model, we estimate the relative (position, velocity) and higher order derivatives, which are collectively termed as the relative kinematics of the anchorless network. The derived data model is inherently ill-posed, however under certain immobility constraints, we propose closed-form solutions to recursively estimate the relative kinematics. For the sake of completeness, we also estimate the absolute node kinematics, given reference anchors. Theoretical bounds are derived, and simulations are conducted to analyse the performance of proposed solutions.

@en