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

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

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

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

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

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

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

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

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

Reports on the activities of the Signal Processing Theory and Methods Technical Committee.@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

Diagnosis and treatment of atrial fibrillation can benefit from various signal processing approaches employed on atrial electrograms. However, the performance and interpretation of these approaches get highly degraded by far-field ventricular activities (VAs) that distort the morphology of the pure atrial activities (AAs). In this study, we aim to remove VAs from the recorded unipolar electrogram while preserving the AA components. To do so, we have developed a framework which first removes the VA-containing segments and interpolates the remaining samples. This will also partly remove the atrial components that overlap with VA signals, e.g., during atrial fibrillation. To reconstruct the AA components, we estimate them from the removed VA-containing segments based on a low-rank and sparse matrix decomposition and add them back to the electrograms. The presented framework is of rather low complexity, preserves AA components, and requires only a single EGM recording. Instrumental comparison to template matching and subtraction and independent component analysis shows that the proposed approach leads to smoother results with better similarity to the true atrial signal.

@en

Radio astronomy is known for its very large telescope dishes but is currently making a transition towards the use of a large number of small antennas. For example, the Low Frequency Array, commissioned in 2010, uses about 50 stations each consisting of 96 low band antennas and 768 or 1536 high band antennas. The low-frequency receiving system for the future Square Kilometre Array is envisaged to initially consist of over 131,000 receiving elements and to be expanded later. These instruments pose interesting array signal processing challenges. To present some aspects, we start by describing how the measured correlation data is traditionally converted into an image, and translate this into an array signal processing framework. This paves the way to describe self-calibration and image reconstruction as estimation problems. Self-calibration of the instrument is required to handle instrumental effects such as the unknown, possibly direction dependent, response of the receiving elements, as well a unknown propagation conditions through the Earth’s troposphere and ionosphere. Array signal processing techniques seem well suited to handle these challenges. Interestingly, image reconstruction, calibration and interference mitigation are often intertwined in radio astronomy, turning this into an area with very challenging signal processing problems.@en

Achieving high resolution time-of-arrival (TOA) estimation in multipath propagation scenarios from bandlimited observations of communication signals is challenging because the multipath channel impulse response (CIR) is not bandlimited. Modeling the CIR as a sparse sequence of Diracs, TOA estimation becomes a problem of parametric spectral inference from observed bandlimited signals. To increase resolution without arriving at unrealistic sampling rates, we consider multiband sampling approach, and propose a practical multibranch receiver for the acquisition. The resulting data model exhibits multiple shift invariance structures, and we propose a corresponding multiresolution TOA estimation algorithm based on the ESPRIT algorithm. The performance of the algorithm is compared against the derived Cramér Rao Lower Bound, using simulations with standardized ultra-wideband (UWB) channel models. We show that the proposed approach provides high resolution estimates while reducing spectral occupancy and sampling costs compared to traditional UWB approaches.

@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