Integration of amplitude and phase statistics for complete artifact removal in independent components of neuromagnetic recordings

In magnetoencephalography (MEG) and electroencephalography (EEG), independent component analysis is widely applied to separate brain signals from artifact components. A number of different methods have been proposed for the automatic or semiautomatic identification of artifact components. Most of the proposed methods are based on amplitude statistics of the decomposed MEG/EEG signal. We present a fully automated approach based on amplitude and phase statistics of decomposed MEG signals for the isolation of biological artifacts such as ocular, muscle, and cardiac artifacts (CAs). The performance of different artifact identification measures was investigated. In particular, we show that phase statistics is a robust and highly sensitive measure to identify strong and weak components that can be attributed to cardiac activity, whereas a combination of different measures is needed for the identification of artifacts caused by ocular and muscle activity. With the introduction of a rejection performance parameter, we are able to quantify the rejection quality for eye blinks and CAs. We demonstrate in a set of MEG data the good performance of the fully automated procedure for the removal of cardiac, ocular, and muscle artifacts. The new approach allows routine application to clinical measurements with small effect on the brain signal.      

Imaging brain dynamics using independent component analysis

The analysis of electroencephalographic and magnetoencephalographic recordings is important both for basic brain research and for medical diagnosis and treatment. Independent component analysis (ICA) is an effective method for removing artifacts and separating sources of the brain signals from these recordings. A similar approach is proving useful for analyzing functional magnetic resonance brain imaging data. In this paper, we outline the assumptions underlying ICA and demonstrate its application to a variety of electrical and hemodynamic recordings from the human brain     

Motion artifact cancellation in NIR spectroscopy using Wiener filtering

We present a Wiener filtering based algorithm for the elimination of motion artifacts present in Near Infrared (NIR) spectroscopy measurements. Until now, adaptive filtering was the only technique used in the noise cancellation in NIR studies. The results in this preliminary study revealed that the proposed method gives better estimates than the classical adaptive filtering approach without the need for additional sensor measurements. Moreover, this novel technique has the potential to filter out motion artifacts in functional near infrared (fNIR) signals, too.     

CPR artifact removal from human ECG using optimal multichannel filtering

The purpose of this study was to assess whether the artifacts presented by precordial compressions during cardiopulmonary resuscitation could be removed from the human electrocardiogram (ECG) using a filtering approach. This would allow analysis and defibrillator charging during ongoing precordial compressions yielding a very important clinical improvement to the treatment of cardiac arrest patients. In this investigation we started with noise-free human ECGs with ventricular fibrillation (VF) and ventricular tachycardia (VT) records. To simulate a realistic resuscitation situation, we added a weighted artifact signal to the human ECG, where the weight factor was chosen to provide the desired signal-to-noise ratio (SNR) level. As artifact signals we used ECGs recorded from animals in asystole during precordial compressions at rates 60, 90, and 120 compressions/min. The compression depth and the thorax impedance was also recorded. In a real-life situation such reference signals are available and, using an adaptive multichannel Wiener filter, we construct an estimate of the artifact signal, which subsequently can be subtracted from the noisy human ECG signal. The success of the proposed method is demonstrated through graphic examples, SNR, and rhythm classification evaluations.     

Artifacts in Wiener kernels estimated using Gaussian white noise

Wiener's nonlinear system identification theory characterizes a system function with a set of kernels of integrals. One method of determining these Wiener kernels is the cross-correlation technique proposed by Lee and Schetzen, which uses Gaussian white noise as the input to the unknown system. Because a test stimulus is only an approximation of infinitely long Gaussian white noise, it is possible that artifacts are generated during the estimation of the kernels. To help identify and characterize these artifacts, Wiener kernel estimates for two simple nonlinear model systems were made using a pseudorandom Gaussian white noise sequence. The results showed that because of the approximation of a Gaussian distribution, artifacts appear in the estimated kernels due to a form of aliasing. These artifacts can be reduced by increasing the sequence length of the input noise.     

On Nonparametric Identification of Wiener Systems

In this paper, a new method for the identification of the Wiener nonlinear system is proposed. The system, being a cascade connection of a linear dynamic subsystem and a nonlinear memoryless element, is identified by a two-step semiparametric approach. The impulse response function of the linear part is identified via the nonlinear least-squares approach with the system nonlinearity estimated by a pilot nonparametric kernel regression estimate. The obtained estimate of the linear part is then used to form a nonparametric kernel estimate of the nonlinear element of the Wiener system. The proposed method permits recovery of a wide class of nonlinearities which need not be invertible. As a result, the proposed algorithm is computationally very efficient since it does not require a numerical procedure to calculate the inverse of the estimate. Furthermore, our approach allows non-Gaussian input signals and the presence of additive measurement noise. However, only linear systems with a finite memory are admissible. The conditions for the convergence of the proposed estimates are given. Computer simulations are included to verify the basic theory     

Second-order adaptive Volterra system identification based ondiscrete nonlinear Wiener model

The authors present the nonlinear LMS adaptive filtering algorithm based on the discrete nonlinear Wiener (1942) model for second-order Volterra system identification application. The main approach is to perform a complete orthogonalisation procedure on the truncated Volterra series. This allows the use of the LMS adaptive linear filtering algorithm for calculating all the coefficients with efficiency. This orthogonalisation method is based on the nonlinear discrete Wiener model. It contains three sections: a single-input multi-output linear with memory section, a multi-input, multi-output nonlinear no-memory section and a multi-input, single-output amplification and summary section. For a white Gaussian noise input signal, the autocorrelation matrix of the adaptive filter input vector can be diagonalised unlike when using the Volterra model. This dramatically reduces the eigenvalue spread and results in more rapid convergence. Also, the discrete nonlinear Wiener model adaptive system allows us to represent a complicated Volterra system with only few coefficient terms. In general, it can also identify the nonlinear system without over-parameterisation. A theoretical performance analysis of steady-state behaviour is presented. Computer simulations are also included to verify the theory     

Optimal reduction of MCG in fetal MEG recordings

Recording fetal magnetoencephalographic (fMEG) signals in-utero is a demanding task due to biological interference, especially maternal and fetal magnetocardiographic (MCG) signals. A method based on orthogonal projection of MCG signal space vectors (OP) was evaluated and compared with independent component analysis (ICA). The evaluation was based on MCG amplitude reduction and signal-to-noise ratio of fetal brain signals using exemplary datasets recorded during ongoing studies related to auditory evoked fields. The results indicate that the OP method is the preferable approach for attenuation of MCG and for preserving the fetal brain signals in fMEG recordings.     

Body movement activity recognition for ambulatory cardiac monitoring

Wearable electrocardiogram (W-ECG) recorders are increasingly in use by people suffering from cardiac abnormalities who also choose to lead an active lifestyle. The challenge presently is that the ECG signal is influenced by motion artifacts induced by body movement activity (BMA) of the wearer. The usual practice is to develop effective filtering algorithms which will eliminate artifacts. Instead, our goal is to detect the motion artifacts and classify the type of BMA from the ECG signal itself. We have recorded the ECG signals during specified BMAs, e.g., sitting still, walking, movements of arms and climbing stairs, etc. with a single-lead system. The collected ECG signal during BMA is presumed to be an additive mix of signals due to cardiac activities, motion artifacts and sensor noise. A particular class of BMA is characterized by applying eigen decomposition on the corresponding ECG data. The classification accuracies range from 70% to 98% for various class combinations of BMAs depending on their uniqueness based on this technique. The above classification is also useful for analysis of P and T waves in the presence of BMA.     

Combination of canonical correlation analysis and empirical mode decomposition applied to denoising the labor electrohysterogram

The electrohysterogram (EHG) is often corrupted by electronic and electromagnetic noise as well as movement artifacts, skeletal electromyogram, and ECGs from both mother and fetus. The interfering signals are sporadic and/or have spectra overlapping the spectra of the signals of interest rendering classical filtering ineffective. In the absence of efficient methods for denoising the monopolar EHG signal, bipolar methods are usually used. In this paper, we propose a novel combination of blind source separation using canonical correlation analysis (BSS_CCA) and empirical mode decomposition (EMD) methods to denoise monopolar EHG. We first extract the uterine bursts by using BSS_CCA then the biggest part of any residual noise is removed from the bursts by EMD. Our algorithm, called CCA_EMD, was compared with wavelet filtering and independent component analysis. We also compared CCA_EMD with the corresponding bipolar signals to demonstrate that the new method gives signals that have not been degraded by the new method. The proposed method successfully removed artifacts from the signal without altering the underlying uterine activity as observed by bipolar methods. The CCA_EMD algorithm performed considerably better than the comparison methods.     

Multi-component filtering of magnetic signals

This paper presents two treatments of real signals recorded by an immersed passive array during an underwater magnetic detection experiment. In this situation, the signal-to-noise ratio is very low, and the noise characteristics are close to the characteristics of the signal that is to be detected. Methods that estimate the signal by minimizing the mean-square distance between the signal and its estimate (using a Wiener filter) are optimal, but some a priori knowledge must first be available before these methods can be used. In this specific application, this a priori knowledge can only be acquired by analyzing the data, and a straightforward application of the optimal method may be misleading. Two filtering methods are therefore implemented to estimate a vector valued signal in a noisy observation. These filtering methods are simple and robust, and they give good results on real data. One method is an adaptation of Wiener filtering. The other consists of global filtering using the eigen-structure of the spectral density matrix of the received signals. Results for each method are presented and compared.     

Main flattening directions and Quadtree decomposition for multi-way Wiener filtering

Previous studies have shown that multi-way Wiener filtering improves the restoration of tensors impaired by an additive white Gaussian noise. Multi-way Wiener filtering is based on the distinction between noise and signal subspaces. In this paper, we show that the lower is the signal subspace dimension, the better is the restored tensor. To reduce the signal subspace dimension, we propose a method based on array processing technique to estimate main orientations in a flattened tensor. The rotation of a tensor of its main orientation values permits to concentrate the information along either rows or columns of the flattened tensor. We show that multi-way Wiener filtering performed on the rotated noisy tensor enables an improved recovery of signal tensor. Moreover, we propose in this paper a quadtree decomposition to avoid a blurry effect in the recovered tensor by multi-way Wiener filtering. We show that this proposed block processing reduces the whole blur and restores local characteristics of the signal tensor. Thus, we show that multi-way Wiener filtering is significantly improved thanks to rotations of the estimated main orientations of tensors and a block processing approach.     

Epilepsy analysis using open source EDF tools for information science and data analytics

Electroencephalogram (EEG) is the signals that measure the electrical variances of brain using metal electrodes. We observe the EEG signals by using European Data Format (EDF) BROWSER and EEG STUDIO. By using EDF BROWSER, we can get the mean and frequency from the filtered output signal using band‐pass filter. Using EDF BROWSER, we can also perform Root Mean Square (RMS) and signal samples. Using EEG STUDIO, we can analyze the average frequency and standard deviation. Epileptic seizure prediction and detection are done by spike detection, frequency domain analysis, and nonlinear methods. EEG signal contains different artifacts like electrooculography (EOG), EKG, and electrocardiogram (ECG). ECG signals are produced by heart. EOG signals are produced by eyes. EMG signals are produced by muscle coordination.     

Noise Removal From Hyperspectral Images by Multidimensional Filtering

A generalized multidimensional Wiener filter for denoising is adapted to hyperspectral images (HSIs). Commonly, multidimensional data filtering is based on data vectorization or matricization. Few new approaches have been proposed to deal with multidimensional data. Multidimensional Wiener filtering (MWF) is one of these techniques. It considers a multidimensional data set as a third-order tensor. It also relies on the separability between a signal subspace and a noise subspace. Using multilinear algebra, MWF needs to flatten the tensor. However, flattening is always orthogonally performed, which may not be adapted to data. In fact, as a Tucker-based filtering, MWF only considers the useful signal subspace. When the signal subspace and the noise subspace are very close, it is difficult to extract all the useful information. This may lead to artifacts and loss of spatial resolution in the restored HSI. Our proposed method estimates the relevant directions of tensor flattening that may not be parallel either to rows or columns. When rearranging data so that flattening can be performed in the estimated directions, the signal subspace dimension is reduced, and the signal-to-noise ratio is improved. We adapt the bidimensional straight-line detection algorithm that estimates the HSI main directions, which are used to flatten the HSI tensor. We also generalize the quadtree partitioning to tensors in order to adapt the filtering to the image discontinuities. Comparative studies with MWF, wavelet thresholding, and channel-by-channel Wiener filtering show that our algorithm provides better performance while restoring impaired HYDICE HSIs.     

CFE条件下基于循环维纳滤波的单通道分离算法

在航天测控数传（C＆T）信号中,频谱混叠的现象非常普遍,因此对频谱严重混叠的信号进行单通道盲分离成为信号处理领域中研究的热点和难点。基于线性-共轭-线性频移（LCL-FRESH）滤波的基本概念,考虑到在实际非合作通信应用中基于循环平稳的LCL-FRESH滤波依赖于较高的循环频率精度,然而循环频率误差（CFE）总是无可避免,提出了一种改进的CFE校正算法,简单分析了误差条件下滤波器性能下降的原因。仿真表明,所研究的算法可以有效分离存在CFE下时频重叠的数传通信信号。     

A spatially adaptive nonparametric regression image deblurring.

We propose a novel nonparametric regression metthod for deblurring noisy images. The method is based on the local polynomial approximation (LPA) of the image and the paradigm of intersecting confidence intervals (ICI) that is applied to define the adaptive varying scales (window sizes) of the LPA estimators. The LPA-ICI algorithm is nonlinear and spatially adaptive with respect to smoothness and irregularities of the image corrupted by additive noise. Multiresolution wavelet algorithms produce estimates which are combined from different scale projections. In contrast to them, the proposed ICI algorithm gives a varying scale adaptive estimate defining a single best scale for each pixel. In the new algorithm, the actual filtering is performed in signal domain while frequency domain Fourier transform operations are applied only for calculation of convolutions. The regularized inverse and Wiener inverse filters serve as deblurring operators used jointly with the LPA-design directional kernel filters. Experiments demonstrate the state-of-art performance of the new estimators which visually and quantitatively outperform some of the best existing methods.     

A kind of denoising method for airborne waveform sampling LiDAR

This paper describes a full waveform sampling LiDAR system applying stripe principle. A kind of denoising method based on edge detection of original stripe signal is proposed. This method is compared with other denoising methods, such as Wiener filtering, mean filtering and median filtering. It is found that the proposed denoising method is much more effective for dealing with the waveform signals.     

Music/Voice Separation Based on Kernel Back‐Fitting Using Weighted β‐Order MMSE Estimation

Recent developments in the field of separation of mixed signals into music/voice components have attracted the attention of many researchers. Recently, iterative kernel back‐fitting, also known as kernel additive modeling, was proposed to achieve good results for music/voice separation. To obtain minimum mean square error (MMSE) estimates of short‐time Fourier transforms of sources, generalized spatial Wiener filtering (GW) is typically used. In this paper, we propose an advanced music/voice separation method that utilizes a generalized weighted β‐order MMSE estimation (WbE) based on iterative kernel back‐fitting (KBF). In the proposed method, WbE is used for the step of mixed music signal separation, while KBF permits kernel spectrogram model fitting at each iteration. Experimental results show that the proposed method achieves better separation performance than GW and existing Bayesian estimators.     

红外点目标图像的两种滤波技术

通过对点目标图像预处理方法的研究,采用基于形态学滤波与维纳差分滤波背景抑制算法,在matlab中利用这2种算法进行红外点目标图像的处理仿真。仿真结果表明形态学滤波有更好的滤波效果,实时性好,易于硬件实现,更适合作为红外图像弱小目标检测的预处理手段。     

Wiener predistorter for nonlinear satellite downlinks

The nonlinear power amplifier and the analogue output channel filter with nonconstant group delay add nonlinear and linear distortions to the satellite transmitted signal, respectively. Recently, there has been growing interest in using Wiener predistorters, consisting of separate nonlinear and linear components, to compensate for these distortions in the satellite downlinks. The Wiener predistorter has been shown to effectively deal with signal distortions and has lower complexity compared to other state-of-the-art methods. In this paper, we argue that fully compensating the nonconstant group-delay distortion degrades the overall performance of the Wiener predistortion systems. This is primarily due to the increased peak-to-average power ratio of the signal at the output of the linear predistorter. We show that the overall performance of the Wiener predistorter can be improved by undercompensating the group-delay distortion. We propose two optimizations to address the PAPR growth problem and show using simulations that our approach leads to significant improvement in predistortion performance. Using our Wiener predistorter, the total degradation gap to the ideal limiter channel can be reduced to merely 0.34 dB for 64-APSK modulated signals.