RH

R.C. Hendriks

info

Please Note

83 records found

This paper investigates the positioning of the pilot symbols, as well as the power distribution between the pilot and the communication symbols for the orthogonal time frequency space (OTFS) modulation scheme. We analyze the pilot placements that minimize the mean squared error (MSE) in estimating the channel taps. This allows us to identify two new pilot allocations for OTFS that save approximately 50% of the pilot overhead compared to existing allocations. In addition, we optimize the average channel capacity by adjusting the power distribution. We show that this leads to a significant increase in average capacity. The results provide valuable guidance for designing the OTFS parameters to achieve maximum capacity. Numerical simulations are performed to validate the findings. ...
Journal article (2025) - Giovanni Bologni, Richard C. Hendriks, Richard Heusdens
This article focuses on estimating relative transfer functions (RTFs) for beamforming applications. Traditional methods often assume that spectra are uncorrelated, an assumption that is often violated in practical scenarios due to factors such as time-domain windowing or the non-stationary nature of signals, as observed in speech. To overcome these limitations, we propose an RTF estimation technique that leverages spectral and spatial correlations through subspace analysis. Additionally, we derive Cramér–Rao bounds (CRBs) for the RTF estimation task, providing theoretical insights into the achievable estimation accuracy. These bounds reveal that channel estimation can be performed more accurately if the noise or the target signal exhibits spectral correlations. Experiments with both real and synthetic data show that our technique outperforms the narrowband maximum-likelihood estimator, known as covariance whitening (CW), when the target exhibits spectral correlations. Although the proposed algorithm generally achieves accuracy close to the theoretical bound, there is potential for further improvement, especially in scenarios with highly spectrally correlated noise. While channel estimation has various applications, we demonstrate the method using a minimum variance distortionless (MVDR) beamformer for multichannel speech enhancement. A free Python implementation is also provided. ...
Journal article (2025) - Lianne N. van Staveren, Yannick J.H.J. Taverne, Richard Hendriks, Natasja M.S. de Groot
Background: Short atrial fibrillation cycle lengths (AFCLs) and regular activation patterns are associated with drivers of atrial fibrillation, although the relation with underlying patterns of activation is incompletely understood. Previous studies used automated assessment of electrograms to determine fast and regular fibrillatory rates. Objective: We investigated the relation among AFCL, temporal variation in AFCL, and the occurrence of driver-like patterns of activation using high-density local activation time mapping. Methods: High-density epicardial mapping of the right atrium and left atrial ventricular groove including Bachmann's bundle was performed in 71 patients admitted for elective cardiac surgery. Recording sites with the shortest median AFCL or the smallest standard deviation of AFCL were identified. Patterns of activation included focal or rotational activation, smooth propagation, propagation with conduction block (CB), collision, and remnant activity. Results: There was a higher number of fibrillation waves with CB (81% [interquartile range (IQR) 76%–85%] vs 74% [68%–76%]; P < .001) and fractionated potentials (22% [12%–37%] vs 12% [9%–15%]; P < .001) at shortest median AFCL than at other recording sites. Smallest standard deviation sites harbored more smoothly propagating waves (33% [24%–54%] vs 17% [11%–25%]; P < .001) and a higher proportion of single potentials (76% [60%–89%] vs 59% [54%–65%]; P < .001). Both highly regular and fastest reactivated sites did not correspond to the origin of (repetitive) focal fibrillation waves. Conclusion: During extensive mapping, the fastest or most regularly activated areas are characterized by CB and smoothly propagating fibrillation waves instead of repetitive occurrence of focal or rotational activation patterns. This study rejects the concept of detecting drivers by identifying the fastest or most regularly activated recording sites. ...
Journal article (2024) - K. R. Stunnenberg, R. C. Hendriks, J. L. Vroegop, M. L. Adank, B. Hunyadi
The pursuit of sensitive and dependable biomarkers capable of capturing the neural processes associated with cognition is a prominent area of interest. Event-related potentials (ERPs) hold significant promise for assessing cognitive dysfunction in various neurological disorders. However, existing data analysis techniques often underutilize the available data and may benefit from potential enhancements. In this paper, we investigate biomarker extraction methods based on two ERP experiments. First, we derive average ERPs from the electroencephalography (EEG) recorded during each experiment and store them in third-order tensors with subjects, channels and time samples along the three modes. Then, we extract biomarkers from these datasets via tensor decompositions. We compare single tensor decompositions and joint tensor decompositions that fuse the data from the individual tensors. In a simulated ERP experiment we compare the benefits and limitations of different tensor-based data fusion methods. Finally, we investigate their performance on a real dataset obtained from schizophrenia patients. ...
The main focus of this paper is an active sensing application that involves selecting transmit and receive sensors to optimize the Cramér-Rao bound (CRB) on target parameters. Although the CRB is non-convex in the transmit and receive selection, we demonstrate that it is convex in the virtual array weight vector, which describes the multiplicity of the virtual array elements. Based on this finding, we propose a novel algorithm that optimizes the virtual array weight vector first and then finds a matching transceiver array. This greatly enhances the efficiency of the transmit and receive sensor selection problem. ...
Journal article (2024) - Niels De Koeijer, Martin Bo Moller, Jorge Martinez, Pablo Martinez-Nuevo, Richard C. Hendriks
Sound zone algorithms control the inputs to a loudspeaker array such that spatially distinct zones, each with separate audio content, are created. This work proposes a sound zone approach which includes a model of human auditory perception in the optimization problem designing the loudspeaker control filters. The control filters are therefore optimized directly for human experience, rather than by proxy through sound pressure, as is done in typical approaches. The proposed optimization problem features a perceptually weighted constraint on the bright zone reproduction error, which allows the user of the algorithm to specify the desired bright zone quality. The proposed method achieves 2 to 4 dB of additional acoustic contrast and is expected to yield less distracting dark-zone interference for the same perceived quality when compared to a traditional approach. ...
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. ...
Journal article (2024) - Johannes W. de Vries, Steven van de Par, Geert Leus, Richard Heusdens, Richard C. Hendriks
Hearing impairment is a prevalent problem with daily challenges like impaired speech intelligibility and sound localisation. One of the shortcomings of spatial filtering in hearing aids is that speech intelligibility is often not optimised directly, meaning that different auditory processes contributing to intelligibility are often not considered. One example is the perceptual phenomenon known as spatial release from masking (SRM). This paper develops a signal model that explicitly considers SRM in the beamforming design, achieved by transforming the binaural intelligibility prediction model (BSIM) into a signal processing framework. The resulting extended signal model is used to analyse the performance of reference beamformers and design a novel beamformer that more closely considers how the auditory system perceives binaural sound. It can be shown that the binaural minimum variance distortionless response (BMVDR) beamformer is also an optimal solution for the extended, perceived model, suggesting that SRM does not play a significant role in intelligibility enhancement after optimal beamforming. However, the optimal beamformer is no longer unique in the extended signal model. The additional secondary degrees of freedom can be used to preserve binaural cues of interfering sources while still achieving the same perceived performance of the BMVDR beamformer, though with a possible high sensitivity to intelligibility model mismatch errors. ...
Journal article (2023) - Lianne N. van Staveren, Richard C. Hendriks, Yannick J.H.J. Taverne, Natasja M.S. de Groot
Background: Dominant frequencies (DFs) or complex fractionated atrial electrograms (CFAEs), indicative of focal sources or rotational activation, are used to identify target sites for atrial fibrillation (AF) ablation in clinical studies, although the relationship among DF, CFAE, and activation patterns remains unclear. Objectives: This study sought to investigate the relationship between patterns of activation underlying DF and CFAE sites during AF. Methods: Epicardial high-resolution mapping of the right and left atrium including Bachmann's bundle was performed in 71 participants. We identified the highest dominant frequency (DF max) and highest degree of CFAE (CFAE max) with the use of existing clinical criteria and classified patterns of activation as focal or rotational activation and smooth propagation, conduction block (CB), collision and remnant activity, and fibrillation potentials as single, double, or fractionated potentials containing, respectively, 1, 2, or 3 or more negative deflections. Relationships among activation patterns, DF max, and potential types were investigated. Results: DF max were primarily located at the left atrioventricular groove and did not harbor focal activation (proportion focal waves: 0% [IQR: 0%-2%]). Compared with non-DF max sites, DF max were characterized by more frequent smooth propagation (22% [IQR: 7%-48%] vs 17% [IQR: 11%-24%]; P = 0.001), less frequent conduction block (69% [IQR: 51%-81%] vs 74% [IQR: 69%-78%]; P = 0.006), a higher proportion of single potentials (72% [IQR: 55%-84%] vs 6%1 [IQR: 55%-65%]; P = 0.003), and a lower proportion of fractionated potentials (4% [IQR: 1%-11%] vs 12% [IQR: 9%-15%]; P = 0.004). CFAE max were mainly found at the pulmonary veins area, and only 1% [IQR: 0%-2%] of all CFAE max contained focal activation. Compared with non-CFAE max sites, CFAE max sites were characterized by less frequent smooth propagation (1% [IQR: 0%-1%] vs 17% [IQR: 12%-24%]; P < 0.001) and more frequent remnant activity (20% [IQR: 12%-29%] vs 8% [IQR: 5%-10%]; P < 0.001), and harbored predominantly fractionated potentials (52% [IQR: 43%-66%] vs 12% [IQR: 9%-14%]; P < 0.001). Conclusions: Focal or rotational patterns of activation were not consistently detected at DF max domains and CFAE max sites. These findings do not support the concept of targeting DF max or CFAE max according to existing criteria for AF ablation. ...
Conference paper (2023) - Johannes W. de Vries, Miao Sun, Natasja M.S. de Groot, Richard C. Hendriks
Estimating tissue conductivity parameters from electrograms (EGMs) could be an important tool for diagnosing and treating heart rhythm disorders such as atrial fibrillation (AF). One of these parameters is the fibre direction, often assumed to be known in conductivity estimation methods. In this paper, a novel method to estimate the fibre direction from EGMs is presented. This method is based on local conduction slowness vectors of a propagating activation wave. These conduction slowness vectors follow an elliptical pattern that depends on the underlying conductivity parameters. The fibre direction and conductivity anisotropy ratio can therefore be estimated by fitting an ellipse to the conduction slowness vectors. Applying the presented method on simulated data shows that it can estimate the fibre direction more accurately than existing methods, and that its performance depends mostly on the range of wavefront directions present in the measurement area. The main advantage of the presented method is that it still functions relatively well in the presence of conduction blocks, as long as the surrounding tissue is approximately homogeneous. ...
Journal article (2023) - I. van der Werf, H. S. Dol, K. C. H. Blom, R. Heusdens, R. C. Hendriks, G. J. T. Leus
In this paper, we show the mathematical equivalence of two popular modulation schemes: OSDM and OTFS. The former is mainly used in underwater acoustic communications, while the latter scheme is a promising modulation technique in radio-frequency communications. Although literature suggests a link between the two modulation schemes by connecting them to related modulation schemes like V-OFDM and A-OFDM, to the best of the authors’ knowledge, a direct mathematical comparison between the schemes has not been presented yet. The main purpose of this paper is therefore to show the mathematical equivalence of the two schemes. In addition, by combining the knowledge of acoustic and radio-frequency communications, we give insight in the performance of OSDM/OTFS in terms of intersymbol interference (ISI) and intercarrier interference (ICI) by analyzing its signal structure. ...
In this paper, we present an improvement of an existing method to estimate the channel parameters of an underwater acoustic channel in a wideband scenario. For underwater acoustic communications, the channel is characterized by paths having both time delay and Doppler scale. In this paper we focus on using a linear frequency modulated (LFM) preamble to estimate the channel parameters. We propose an improvement to the state-of-the-art method, which utilizes the fractional Fourier transform (FrFT), and show that matched filtering and the FrFT method can both be viewed as integral transforms having the same kernel. We argue that, contrary to the prevailing belief in the literature, matched filtering is superior compared to the FrFT method for channel parameter estimation. We support our findings using numerical experiments. ...
Estimation of the acoustic-scene related parameters such as relative transfer functions (RTFs) from source to microphones, source power spectral densities (PSDs) and PSDs of the late reverberation is essential and also challenging. Existing maximum likelihood estimators typically consider only subsets of these parameters and use each time frame separately. In this paper we explicitly focus on the single source scenario and first propose a joint maximum likelihood estimator (MLE) to estimate all parameters jointly using a single time frame. Since the RTFs are typically invariant for a number of consecutive time frames we also propose a joint maximum likelihood estimator (MLE) using multiple time frames which has similar estimation performance compared to a recently proposed reference algorithm called simultaneously confirmatory factor analysis (SCFA), but at a much lower complexity. Moreover, we present experimental results which demonstrate that the estimation accuracy, together with the performance of noise reduction, speech quality and speech intelligibility, of our proposed joint MLE outperform those of existing MLE based approaches that use only a single time frame. ...
Conference paper (2023) - Changheng Li, Richard Hendriks
Spatial filtering techniques typically rely on estimates of the target relative transfer function (RTF). However, the target speech signal is typically corrupted by late reverberation and ambient noise, which complicates RTF estimation. Existing methods subtract the noise covariance matrix to obtain the target plus late reverberation covariance matrix, from where the RTF is estimated. However, the noise covariance matrix is typically unknown. More specifically, the noise power spectral density (PSD) is typically unknown, while the spatial coherence matrix can be assumed known as it might remain time-invariant for a longer time. Using the spatial coherence matrices we simplify the signal model such that the off-diagonal elements are not affected by the PSDs of the late reverberation and the ambient noise. Then we use these elements to estimate the target covariance matrix, from where the RTF can be obtained. Hence, the resulting estimate of the RTF is insensitive to the noise PSD. Experiments demonstrate the estimation performance of our proposed method. ...
Conference paper (2023) - Changheng Li, Richard C.Hendriks
Estimating the parameters that describe the acoustic scene is very important for many microphone array applications. For example, consider the power spectral densities (PSDs) or relative acoustic transfer functions (RTFs) that are required when estimating a particular sound source using multi-microphone noise reduction. State-of-the-art algorithms estimate the pa-rameters per segment, where each segment consists of a fixed number of time frames. These algorithms exploit the assumption that PSDs are constant per time frame, and RTFs are constant per segment. However, in practice, sound sources will move relative to the microphone array. Improved per-formance is therefore expected when the actual time frames that are used to form the segments are adapted such that time frames all share the same (unknown) RTF. In this paper, we therefore present an algorithm to obtain an optimal adaptive time segmentation and combine this with our previously pub-lished joint maximum likelihood estimator (JMLE) for jointly estimating the RTF, source PSD and late reverberation PSD of a single source in a reverberant environment. ...
Journal article (2023) - Changheng Li, Richard C. Hendriks
Acoustic-scene-related parameters such as relative transfer functions (RTFs) and power spectral densities (PSDs) of the target source, late reverberation and ambient noise are essential for microphone array signal processing and are challenging to estimate. Existing methods typically only estimate a subset of the parameters by assuming the other parameters are known. This can lead to unmatched scenarios and reduced estimation performance on the parameters of interest. Moreover, many methods process time frames independently, despite they share common information such as the same RTF. In this work, we consider a noisy scenario by modelling the noise component as a spatially homogeneous sound field with a time-invariant spatial coherence matrix and time-varying PSD. We first modify an existing alternating least squares (ALS) method to obtain more accurate estimates using a single time frame. Then, we extend the method to use multiple time frames that share the same RTF. Furthermore, we propose more robust constraints on the PSDs to avoid large estimation errors. We compare our proposed methods to the state-of-the-art simultaneously confirmatory factor analysis (SCFA) method, a joint maximum likelihood estimation (JMLE) method and an existing ALS-based method. The experimental results in terms of estimation accuracy, noise reduction performance, predicted speech quality, and predicted speech intelligibility demonstrate that our proposed methods achieve similar performance compared to the state-of-the-art SCFA method, which outperforms the existing ALS method in all scenarios and outperforms the JMLE method particularly in low SNR scenarios. Moreover, our proposed methods have significantly lower computational complexity than SCFA. ...
Journal article (2023) - S. Difrancesco, J.U. van Baardewijk, A.S. Cornelissen, C. Varon, R.C. Hendriks, A.M. Bouwer
Wearable sensors offer new opportunities for the early detection and identification of toxic chemicals in situations where medical evaluation is not immediately possible. We previously found that continuously recorded physiology in guinea pigs can be used for early detection of exposure to an opioid (fentanyl) or a nerve agent (VX), as well as for differentiating between the two. Here, we investigated how exposure to these different chemicals affects the interactions between ECG and respiration parameters as determined by Granger causality (GC). Features reflecting such interactions may provide additional information and improve models differentiating between chemical agents. Traditional respiration and ECG features, as well as GC features, were extracted from data of 120 guinea pigs exposed to VX (n = 61) or fentanyl (n = 59). Data were divided in a training set (n = 99) and a test set (n = 21). Minimum Redundancy Maximum Relevance (mRMR) and Support Vector Machine (SVM) algorithms were used to, respectively, perform feature selection and train a model to discriminate between the two chemicals. We found that ECG and respiration parameters are Granger-related under healthy conditions, and that exposure to fentanyl and VX affected these relationships in different ways. SVM models discriminated between chemicals with accuracy of 95% or higher on the test set. GC features did not improve the classification compared to traditional features. Respiration features (i.e., peak inspiratory and expiratory flow) were the most important to discriminate between different chemical’s exposure. Our results indicate that it may be feasible to discriminate between chemical exposure when using traditional physiological respiration features from wearable sensors. Future research will examine whether GC features can contribute to robust detection and differentiation between chemicals when considering other factors, such as generalizing results across species ...
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 ...

A Multi-Modal Approach to Acoustic Reflector Estimation

Conference paper (2022) - E.H.J. Riemens, Pablo Martinez-Nuevo, Jorge Martinez, Martin Bo Møller, R.C. Hendriks
Loudspeakers are usually placed in an environment unknown to the loudspeaker designers. Having knowledge on the room acoustic properties, e.g., the location of acoustic reflectors, allows to better reproduce the sound field as intended. Current state-of-the-art methods for room boundary detection using microphone measurements typically focus on a two-dimensional setting, causing a model mismatch when employed in real-life scenarios. Detection of arbitrary reflectors in three dimensions encounters practical limitations, e.g., the need for a spherical array and the increased computational complexity. Moreover, loudspeakers may not have an omnidirectional directivity pattern, as usually assumed in the literature, making the detection of acoustic reflectors in some directions more challenging. ...
Conference paper (2022) - Tobias Roest, Jorge Martinez, Han Oosterom, Richard Hendriks
Moving-coil electrodynamic loudspeakers and dynamic microphones use the same linear actuator technology at the core of their operation. Utilising this similarity, loudspeakers have a possible use as recording devices in cases where using dedicated microphones is not feasible. Such a use case exists in public address and voice alarm systems. This paper evaluates the feasibility of using the loudspeakers already in place in these systems as recording devices to provide information back to the system. A system using a single loudspeaker as both a playback and recording device simultaneously is analysed, modelled and simulated. The results show that using a current measuring set-up with an analogue-to-digital converter capable of detecting a range of roughly 120 dB, a speech signal incident at 46 dBSPL in a cone of 150° from a loudspeaker can be successfully estimated in an office room with an announcement playing at 88 dBSPL and background interference present at the same time. As the estimated signal is unknown to the system, the solution generalises to other signal types as well. ...