Av

A.J. van der Veen

info

Please Note

49 records found

Multi-agent systems, such as fleets of robots or drones, are increasingly deployed in logistics, inspection, and surveillance. These systems must reach their targets while maintaining safe separation, even under uncertain dynamics. This is challenging because unmodeled effects, disturbances, and sensor noise can degrade tracking performance and compromise safety. Model Predictive Control (MPC) is well suited for multi-agent navigation since it optimizes trajectories over a prediction horizon while enforcing input and state constraints. However, its performance depends on accurate models, and centralized formulations suffer from poor scalability and a single point of failure. We propose a cooperative Gaussian Process–augmented MPC (GP-MPC) framework that combines learning, chance-constrained safety, and distributed optimization. Each agent uses a Gaussian Process to learn its residual dynamics and quantify local uncertainty, incorporates this uncertainty into a chance-constrained collision-avoidance scheme, and coordinates only with neighbors through an ADMM-based distributed optimization method. This integration provides robustness to model errors and scalability to larger teams. The framework enables collision avoidance using only local uncertainty estimates, removing the need to share covariance information. By restricting computation and communication to each agent’s neighborhood, it maintains scalability and efficiency. Simulations show that the approach yields smoother and more efficient trajectories, faster convergence to targets, and reliable probabilistic safety compared to nominal and nonlinear MPC baselines. Convergence analysis further confirms robust consensus across a range of tuning parameters. ...
Journal article (2025) - Metin Calis, Massimo Mischi, Alle-Jan van der Veen, Raj Thilak Rajan, Borbàla Hunyadi
Tensor decomposition methods for signal processing applications are an active area of research. Real data are often low-rank, noisy, and come in a higher-order format. As such, low-rank tensor approximation methods that account for the high-order structure of the data are often used for denoising. One way to represent a tensor in a low-rank form is to decompose the tensor into a set of orthonormal factor matrices and an all-orthogonal core tensor using a higher-order singular value decomposition. Under noisy measurements, the lower bound for recovering the factor matrices and the core tensor is unknown. In this paper, we exploit the well-studied constrained Cramér-Rao bound to calculate a lower bound on the mean squared error of the unbiased estimates of the components of the multilinear singular value decomposition under additive white Gaussian noise, and we validate our approach through simulations. ...
Journal article (2025) - Metin Calis, Massimo Mischi, Alle Jan van der Veen, Borbala Hunyadi
Dynamic contrast-enhanced ultrasound (DCEUS) is an imaging modality for assessing microvascular perfusion and dispersion kinetics. However, the presence of speckle noise may hamper the quantitative analysis of the contrast kinetics. Common speckle denoising techniques based on low-rank approximations typically model the speckle noise as white Gaussian noise (WGN) after the log transformation and apply matrix-based algorithms. We address the high dimensionality of the 4D DCEUS data and apply low-rank tensor decomposition techniques to denoise speckles. Although there are many tensor decompositions that can describe low rankness, we limit our research to multilinear rank and tubal rank. We introduce a gradient-based extension of the multilinear singular value decomposition to model low multilinear rankness, assuming that the log-transformed speckle noise follows a Fisher-tippet distribution. In addition, we apply an algorithm based on tensor singular value decomposition to model low tubal rankness, assuming that the log-transformed speckle noise is WGN with sparse outliers. The effectiveness of the methods is evaluated through simulations and phantom studies. Additionally, the tensor-based algorithms’ real-world performance is assessed using DCEUS prostate recordings. Comparative analyses with existing DCEUS denoising literature are conducted, and the algorithms’ capabilities are showcased in the context of prostate cancer classification. The addition of Fisher-tippet distribution did not improve the results of tr-MLSVD in the in vivo case. However, most cancer markers are better distinguishable when using a tensor denoising technique than state-of-the-art approaches. ...
Journal article (2024) - Hanie Moghaddasi, Richard C. Hendriks, Borbala Hunyadi, Paul Knops, Mathijs S. Van Schie, Natasja M.S. De Groot, Alle Jan Van Der Veen
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. ...

Celebrating 75 Years of Remarkable Achievements [From the Guest Editors]

Journal article (2023) - Rodrigo Capobianco Guido, Tulay Adali, Emil Bjornson, Laure Blanc-Feraud, Ulisses Braga-Neto, Behnaz Ghoraani, Christian Jutten, Alle Jan Van Der Veen, Hong Vicky Zhao, Xiaoxing Zhu

Celebrating 75 Years of Remarkable Achievements (Part 2) [From the Guest Editors]

Journal article (2023) - Rodrigo Capobianco Guido, Tulay Adali, Emil Bjornson, Laure Blanc-Feraud, Ulisses Braga-Neto, Behnaz Ghoraani, Christian Jutten, Alle Jan Van Der Veen, Hong Vicky Zhao, Xiaoxing Zhu
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%. ...
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 ...
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. ...
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. ...

A position paper from the European Heart Rhythm Association and European Society of Cardiology Working Group on eCardiology in collaboration with the Heart Rhythm Society, Asia Pacific Heart Rhythm Society, Latin American Heart Rhythm Society and Computing in Cardiology

Journal article (2022) - Natasja M.S. de Groot, Dipen Shah, Patrick M. Boyle, Elad Anter, Gari D. Clifford, Isabel Deisenhofer, Pascal van Dessel, Polychronis Dilaveris, Alle Jan van der Veen, More authors...
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. ...
Journal article (2021) - Marco Mercuri, Yiting Lu, Salvatore Polito, Fokko Wieringa, Alle-Jan van der Veen, Chris Van Hoof, Tom Torfs
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. ...
Conference paper (2021) - Tarik Kazaz, Jac Romme, Gerard J.M. Janssen, Alle-Jan van der Veen
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. ...
Conference paper (2020) - A.J. van der Veen, J.P.A. Romme, Ye Cui
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. ...
Conference paper (2020) - H. Moghaddasi, A.J. van der Veen, N.M.S. de Groot, B. Hunyadi
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. ...
Journal article (2020) - Jose C.M. Bermudez, Monica F. Bugallo, Alle-Jan Van Der Veen
Reports on the activities of the Signal Processing Theory and Methods Technical Committee. ...
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. ...
Journal article (2020) - M. J. Bentum, M. K. Verma, R. T. Rajan, A. J. Boonstra, C. J.M. Verhoeven, E. K.A. Gill, A. J. van der Veen, H. Falcke, L. I. Gurvits, More Authors...
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. ...
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. ...
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. ...