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.

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.

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.

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.

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.

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.

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.

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.

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

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.

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.

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.

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.

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. In the proposed method, a LiDAR sensor is added to a loudspeaker to improve wall detection accuracy and robustness. This is done in two ways. First, the model mismatch introduced by horizontal reflectors can be resolved by detecting reflectors with the LiDAR sensor to enable elimination of their detrimental influence from the 2D problem in pre-processing. Second, a LiDAR-based method is proposed to compensate for the challenging directions where the directive loudspeaker emits little energy. We show via simulations that this multi-modal approach, i.e., combining microphone and LiDAR sensors, improves the robustness and accuracy of wall detection.

Mathematical models of the electrophysiology of cardiac tissue play an important role when studying heart rhythm disorders like atrial fibrillation. Model parameters such as conductivity, activation time, and anisotropy ratio are useful parameters to determine the arrhythmogenic substrate that causes abnormalities in the atrial tissue. Existing methods often estimate the model parameters separately and assume some of the parameters to be known as a priori knowledge. In this work, we propose an efficient method to jointly estimate the parameters of interest from the cross power spectral density matrix (CPSDM) model of the electrograms. By applying confirmatory factor analysis (CFA) to the CPSDMs of multi-electrode electrograms, we can make use of the spatial information of the data and analyze the relationship between the desired resolution and the required amount of data. With the reasonable assumptions that the conductivity parameters and the anisotropy parameters are constant across different frequencies and heart beats, we estimate these parameters using multiple frequencies and multiple heart beats simultaneously to easier satisfy the identifiability conditions in the CFA problem. Results on the simulated data show that using multiple heart beats decreases the estimation errors of the conductivity and the estimated activation time parameters. The experimental results on clinical data show that using multiple heart beats for parameter estimation can reduce the reconstruction errors of the clinical electrograms, which further demonstrates the robustness of the proposed method.

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

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.

Many multi-microphone algorithms depend on knowing the relative acoustic transfer functions (RTFs) of the individual sound sources in the acoustic scene. However, accurate joint RTF estimation for multiple sources is a challenging problem. Existing methods to jointly estimate the RTF for multiple sources have either no satisfying performance, or, suffer from a very large computational complexity. In this paper, we propose a method for robust estimation of the individual RTFs in a multi-source acoustic scenario. The presented algorithm is based on linear algebraic concepts and therefore of lower computational complexity compared to a recently presented state-of-the-art algorithm, while having a similar performance. Experimental results are presented to demonstrate the RTF estimation performance as well as the noise reduction performance when combining the estimated RTFs with a beamformer.

The sound field in a room can be represented by a weighted sum of room modes. To estimate the room modes, current literature uses on-the-grid, sparse reconstruction methods. However, these on-the-grid methods are known to suffer from basis mismatch. In this work, we investigate the use of a gridless framework for estimating room modes using atomic norm minimization, a gridless method. The advantage of this approach would be that it does not suffer from this basis mismatch problem. We derive a bound for the sound field reconstruction problem in a one-dimensional room with rigid walls and relate this to the frequency separation that is required by the atomic norm. We conclude that for perfect reconstruction based on the investigated gridless approach, additional prior knowledge about the signal model is required. We show how recovery is possible in a one-dimensional setting by exploiting both the structure of the sound field and the acquisition method.

Background: An increasing number of wearables are capable of measuring electrocardiograms (ECGs), which may help in early detection of atrial fibrillation (AF). Therefore, many studies focus on automated detection of AF in ECGs. A major obstacle is the required amount of manually labelled data. This study aimed to provide an efficient and reliable method to train a classifier for AF detection using large datasets of real-life ECGs. Method: Human-controlled semi-supervised learning was applied, consisting of two phases: the pre-training phase and the semi-automated training phase. During pre-training, an initial classifier was trained, which was used to predict the classes of new ECG segments in the semi-automated training phase. Based on the degree of certainty, segments were added to the training dataset automatically or after human validation. Thereafter, the classifier was retrained and this procedure was repeated. To test the model performance, a real-life telemetry dataset containing 3,846,564 30-s ECG segments of hospitalized patients (n = 476) and the CinC Challenge 2017 database were used. Results: After pre-training, the average F1-score on a hidden testing dataset was 89.0%. Furthermore, after the pre-training phase 68.0% of all segments in the hidden test set could be classified with an estimated probability of successful classification of 99%, providing an F1-score of 97.9% for these segments. During the semi-automated training phase, this F1-score showed little variation (97.3%–97.9% in the hidden test set), whilst the number of segments which could be automatically classified increased from 68.0% to 75.8% due to the enhanced training dataset. At the same time, the overall F1-score increased from 89.0% to 91.4%. Conclusions: Human-validated semi-supervised learning makes training a classifier more time efficient without compromising on accuracy, hence this method might be valuable in the automated detection of AF in real-life ECGs.