Edgar R. Arce-Santana

Automatic sleep staging based on ballistocardiographic signals recorded through bed sensors

This study presents different methods for automatic sleep classification based on heart rate variability (HRV), respiration and movement signals recorded through bed sensors. Two methods for feature extraction have been implemented: time variant-autoregressive model (TVAM) and wavelet discrete transform (WDT); the obtained features are fed into two classifiers: Quadratic (QD) and Linear (LD) discriminant for staging sleep in REM, nonREM and WAKE periods. The performances of all the possible combinations of feature extractors and classifiers are compared in terms of accuracy and kappa index, using clinical polysomographyc evaluation as golden standard. 17 recordings from healthy subjects, including also polisomnography, were used to train and test the algorithms. When automatic classification is compared. QD-TVAM algorithm achieved a total accuracy of 76.81 ± 7.51 % and kappa index of 0.55 ± 0.10, while LD-WDT achieved a total accuracy of 79 ± 10% and kappa index of 0.51 ± 0.17. The results suggest that a good sleep evaluation can be achieved through non-conventional recording systems that could be used outside sleep centers.

Image registration using Markov random coefficient and geometric transformation fields

Image registration is central to different applications such as medical analysis, biomedical systems, and image guidance. In this paper we propose a new algorithm for multimodal image registration. A Bayesian formulation is presented in which a likelihood term is defined using an observation model based on coefficient and geometric fields. These coefficients, which represent the local intensity polynomial transformations, as the local geometric transformations, are modeled as prior information by means of Markov random fields. This probabilistic approach allows one to find optimal estimators by minimizing an energy function in terms of both fields, making the registration between the images possible.

Phase correlation based image alignment with subpixel accuracy

The phase correlation method is a well-known image alignment technique with broad applications in medical image processing, image stitching, and computer vision. This method relies on estimating the maximum of the phase-only correlation (POC) function, which is defined as the inverse Fourier transform of the normalized cross-spectrum between two images. The coordinates of the maximum correspond to the translation between the two images. One of the main drawbacks of this method, in its basic form, is that the location of the maximum can only be obtained with integer accuracy. In this paper, we propose a new technique to estimate the location with subpixel accuracy, by minimizing the magnitude of gradient of the POC function around a point near the maximum. We also present some experimental results where the proposed method shows an increased accuracy of at least one order of magnitude with respect to the base method. Finally, we illustrate the application of the proposed algorithm to the rigid registration of digital images.

Evaluation of multimodal medical image registration based on Particle Filter

This paper presents a performance evaluation of a new multimodal image registration algorithm which is based on Bayesian estimation theory, specifically on Particle Filters. The results point to an efficient, easy to implement and robust to noise algorithm. The registration method showed good performance when using partial data, and it was compared with an algorithm based on maximization of mutual information and a Hyperplanes optimization method. Finally, we showed that the algorithm may be parallelizable, so that it is possible to reduce the computation time for image registration.

Phase correlation with sub-pixel accuracy

Six methods for the accurate estimation of phase-correlation maxima are evaluated.Methods are tested under noise, extreme transformations, incomplete data, and for real cases with unknown transformations.Sinc function fitting provides the best average accuracy.Local Center of Mass, and Minimization of the POC gradient provide good balance between accuracy and efficiency. Six methods for the accurate estimation of the phase-correlation maxima are discussed and evaluated in this article for one- and two-dimensional signals. The evaluation was carried out under a rigid image registration framework, where artificially generated transformations were used in order to perform a quantitative assessment of the accuracy of each method and its robustness in the presence of noise, incomplete data, or extreme transformations. Another round of tests were performed with real cases where the true transformation is unknown, and not necessarily rigid; for these tests, quantitative evaluation was achieved by means of the root mean square error of the overlapping area between the two aligned images. While most methods behaved similarly under difficult conditions, three of the methods under study displayed clear advantages under mild levels of noise, low transformation complexity, and small percentages of missing data. These methods are the local center of mass, sinc function fitting, and minimization of the POC gradient magnitude. The other tested methods included quadratic fitting, linear fitting in the frequency domain, and up-sampling; however, these methods did not perform consistently well.

Blind End-Member and Abundance Extraction for Multispectral Fluorescence Lifetime Imaging Microscopy Data

This paper proposes a new blind end-member and abundance extraction (BEAE) method for multispectral fluorescence lifetime imaging microscopy (m-FLIM) data. The chemometrical analysis relies on an iterative estimation of the fluorescence decay end-members and their abundances. The proposed method is based on a linear mixture model with positivity and sum-to-one restrictions on the abundances and end-members to compensate for signature variability. The synthesis procedure depends on a quadratic optimization problem, which is solved by an alternating least-squares structure over convex sets. The BEAE strategy only assumes that the number of components in the analyzed sample is known a spriori. The proposed method is first validated by using synthetic m-FLIM datasets at 15, 20, and 25 dB signal-to-noise ratios. The samples simulate the mixed response of tissue containing multiple fluorescent intensity decays. Furthermore, the results were also validated with six m-FLIM datasets from fresh postmortem human coronary atherosclerotic plaques. A quantitative evaluation of the BEAE was made against two popular techniques: minimum volume constrained nonnegative matrix factorization (MVC-NMF) and multivariate curve resolution-alternating least-squares (MCR-ALS). Our proposed method (BEAE) was able to provide more accurate estimations of the end-members: 0.32% minimum relative error and 13.82% worst-case scenario, despite different initial conditions in the iterative optimization procedure and noise effect. Meanwhile, MVC-NMF and MCR-ALS presented more variability in estimating the end-members: 0.35% and 0.34% for minimum errors and 15.31% and 13.25% in the worst-case scenarios, respectively. This tendency was also maintained for the abundances, where BEAE obtained 0.05 as the minimum absolute error and 0.12 in the worst-case scenario; MCR-ALS and MVC-NMF achieved 0.04 and 0.06 for the minimum absolute errors, and 0.15 and 0.17 under the worst-case conditions, respectively. In addition, the average computation time was evaluated for the synthetic datasets, where MVC-NMF achieved the fastest time, followed by BEAE and finally MCR-ALS. Consequently, BEAE improved MVC-NMF in convergence to a local optimal solution and robustness against signal variability, and it is roughly 3.6 time faster than MCR-ALS.

Sleep-wake detection based on respiratory signal acquired through a Pressure Bed Sensor

This study proposes an automatic method for the sleep-wake staging in normal and pathologic sleep based only on respiratory effort acquired from a Pressure Bed Sensor (PBS). Motion and respiratory movements were obtained through a PBS and sleep-wake staging was evaluated from those time series. 20 all night polysomnographies, with annotations, used as gold standard and the time series coming from the PBS were used to develop and to evaluate the automatic wake-sleep staging. The database was built up by: 10 healthy subjects and 10 patients with severe sleep apnea. The agreement of the statistical measures between the automatic classification and the human scoring were: 83.59 ± 6.79 of sensitivity, 83.60 ± 15.13 of specificity and 81.91 ± 6.36 of accuracy. These results suggest that some important indexes, such as sleep efficiency, could be computed through a contactless technique.

A non-rigid multimodal image registration method based on particle filter and optical flow

Image Registration is a central task to many medical image analysis applications. In this paper, we present a novel iterative algorithm composed of two main steps: a global affine image registration based on particle filter, and a local refinement obtained from a linear optical flow approximation. The key idea is to iteratively apply these simple and robust steps to efficiently solve complex non-rigid multimodal or unimodal image registrations. Finally, we present a set of evaluation experiments demonstrating the accuracy and applicability of the method to medical images.

Evaluation of the sleep quality based on bed sensor signals: Time-variant analysis

Automatic detection of the sleep macrostructure (Wake, NREM -non Rapid Eye Movement- and REM -Rapid Eye Movement-) based on bed sensor signals is presented. This study assesses the feasibility of different methodologies to evaluate the sleep quality out of sleep centers. The study compares a) the features extracted from time-variant autoregressive modeling (TVAM) and Wavelet Decomposition (WD) and b) the performance of K-Nearest Neighbor (KNN) and Feed Forward Neural Networks (FFNN) classifiers. In the current analysis, 17 full polysomnography recordings from healthy subjects were used. The best agreement for Wake-NREM-REM with respect to the gold standard was 71.95 ± 7.47% of accuracy and 0.42 ± 0.10 of kappa index for TVAM-LD while WD-FFNN shows 67.17 ± 11.88% of accuracy and 0.39 ± 0.13 of kappa index. The results suggest that the sleep quality assessment out of sleep centers could be possible and as consequence more people could be beneficiated.

Phase-correlation guided area matching for realtime vision and video encoding

In computer vision and video encoding applications, one of the first and most important steps is to establish a pixel-to-pixel correspondence between two images of the same scene obtained at slightly different times or points of view. One of the most popular methods to find these correspondences, known as Area Matching, consists in performing a computationally intensive search for each pixel in the first image, around a neighborhood of the same pixel in the second image. In this work we propose a method which significantly reduces the search space to only a few candidates, and permits the implementation of real-time vision and video encoding algorithms which do not require specialized hardware such as GPU’s or FPGA’s. Theoretical and experimental support for this method is provided. Specifically, we present results from the application of the method to the realtime video compression and transmission, as well as the realtime estimation of dense optical flow and stereo disparity maps, where a basic implementation achieves up to 100 fps in a typical dual-core PC.