Alin Achim

Novel Bayesian multiscale method for speckle removal in medical ultrasound images

A novel speckle suppression method for medical ultrasound images is presented. First, the logarithmic transform of the original image is analyzed into the multiscale wavelet domain. The authors show that the subband decompositions of ultrasound images have significantly non-Gaussian statistics that are best described by families of heavy-tailed distributions such as the alpha-stable. Then, the authors design a Bayesian estimator that exploits these statistics. They use the alpha-stable model to develop a blind noise-removal processor that performs a nonlinear operation on the data. Finally, the authors compare their technique with current state-of-the-art soft and hard thresholding methods applied on actual ultrasound medical images and they quantify the achieved performance improvement.

SAR image denoising via Bayesian wavelet shrinkage based on heavy-tailed modeling

Synthetic aperture radar (SAR) images are inherently affected by multiplicative speckle noise, which is due to the coherent nature of the scattering phenomenon. This paper proposes a novel Bayesian-based algorithm within the framework of wavelet analysis, which reduces speckle in SAR images while preserving the structural features and textural information of the scene. First, we show that the subband decompositions of logarithmically transformed SAR images are accurately modeled by alpha-stable distributions, a family of heavy-tailed densities. Consequently, we exploit this a priori information by designing a maximum a posteriori (MAP) estimator. We use the alpha-stable model to develop a blind speckle-suppression processor that performs a nonlinear operation on the data and we relate this nonlinearity to the degree of non-Gaussianity of the data. Finally, we compare our proposed method to current state-of-the-art soft thresholding techniques applied on real SAR imagery and we quantify the achieved performance improvement.

Image denoising using bivariate α-stable distributions in the complex wavelet domain

Recently, the dual-tree complex wavelet transform has been proposed as an analysis tool featuring near shift-invariance and improved directional selectivity compared to the standard wavelet transform. Within this framework, we describe a novel technique for removing noise from digital images. We design a bivariate maximum a posteriori estimator, which relies on the family of isotropic α-stable distributions. Using this relatively new statistical model we are able to better capture the heavy-tailed nature of the data as well as the interscale dependencies of wavelet coefficients. We test our algorithm for the Cauchy case, in comparison with several recently published methods. The simulation results show that our proposed technique achieves state-of-the-art performance in terms of root mean squared (RMS) error.

SAR image filtering based on the heavy-tailed Rayleigh model

Synthetic aperture radar (SAR) images are inherently affected by a signal dependent noise known as speckle, which is due to the radar wave coherence. In this paper, we propose a novel adaptive despeckling filter and derive a maximum a posteriori (MAP) estimator for the radar cross section (RCS). We first employ a logarithmic transformation to change the multiplicative speckle into additive noise. We model the RCS using the recently introduced heavy-tailed Rayleigh density function, which was derived based on the assumption that the real and imaginary parts of the received complex signal are best described using the alpha-stable family of distribution. We estimate model parameters from noisy observations by means of second-kind statistics theory, which relies on the Mellin transform. Finally, we compare the proposed algorithm with several classical speckle filters applied on actual SAR images. Experimental results show that the homomorphic MAP filter based on the heavy-tailed Rayleigh prior for the RCS is among the best for speckle removal

SVM-based texture classification in Optical Coherence Tomography

This paper describes a new method for automated texture classification for glaucoma detection using high resolution retinal Optical Coherence Tomography (OCT). OCT is a non-invasive technique that produces cross-sectional imagery of ocular tissue. Here, we exploit information from OCT images, specifically the inner retinal layer thickness and speckle patterns, to detect glaucoma. The proposed method relies on support vector machines (SVM), while principal component analysis (PCA) is also employed to improve classification performance. Results show that texture features can improve classification accuracy over what is achieved using only layer thickness as existing methods currently do.

Segmentation-Driven Image Fusion Based on Alpha-Stable Modeling of Wavelet Coefficients

A novel region-based image fusion framework based on multiscale image segmentation and statistical feature extraction is proposed. A dual-tree complex wavelet transform (DT-CWT) and a statistical region merging algorithm are used to produce a region map of the source images. The input images are partitioned into meaningful regions containing salient information via symmetric alpha-stable (S alphaS) distributions. The region features are then modeled using bivariate alpha-stable (B alphaS) distributions, and the statistical measure of similarity between corresponding regions of the source images is calculated as the Kullback-Leibler distance (KLD) between the estimated B alphaS models. Finally, a segmentation-driven approach is used to fuse the images, region by region, in the complex wavelet domain. A novel decision method is introduced by considering the local statistical properties within the regions, which significantly improves the reliability of the feature selection and fusion processes. Simulation results demonstrate that the bivariate alpha-stable model outperforms the univariate alpha-stable and generalized Gaussian densities by not only capturing the heavy-tailed behavior of the subband marginal distribution, but also the strong statistical dependencies between wavelet coefficients at different scales. The experiments show that our algorithm achieves better performance in comparison with previously proposed pixel and region-level fusion approaches in both subjective and objective evaluation tests.

A novel system for robust lane detection and tracking

This paper presents a lane detection and tracking system based on a novel lane feature extraction approach and the Gaussian Sum Particle filter (GSPF). The proposed feature extraction approach is based on the fact that by zooming into the vanishing point of the lanes, the lane markings/boundaries will only move on the same straight lines they are on. Objects other than the lanes in the frame do not share this property and can be ignored during the model parameter estimation. This algorithm is able to iteratively refine various traditional feature maps and to operate with curved roads. The tracking part of the system is initialised by a deformable template matching algorithm. Three types of tracking algorithms are compared in our study: the original Sequential Importance Resampling (SIR) particle filter, the Gaussian Particle Filter (GPF) and the Gaussian Sum Particles Filter (GSPF). The GSPF achieves the best performance by integrating a novel likelihood function and an intuitive parameter selection process. Both the GSPF and GPF provide improved tracking performance and require less computational power than the SIR. It has also been found that the detection and tracking performance is enhanced significantly by incorporating the refined feature map.

Non-Gaussian model-based fusion of noisy images in the wavelet domain

This paper describes a new methodology for multimodal image fusion based on non-Gaussian statistical modelling of wavelet coefficients. Special emphasis is placed on the fusion of noisy images. The use of families of generalised Gaussian and alpha-stable distributions for modelling image wavelet coefficients is investigated and methods for estimating distribution parameters are proposed. Improved techniques for image fusion are developed, by incorporating these models into a weighted average image fusion algorithm. The proposed method has been shown to perform very well with both noisy and noise-free images from multimodal datasets, outperforming conventional methods in terms of fusion quality and noise reduction in the fused output.

Compressive image fusion

Compressive sensing (CS) has received a lot of interest due to its compression capability and lack of complexity on the sensor side. In this paper, we present a study of three sampling patterns and investigate their performance on CS reconstruction. We then propose a new image fusion algorithm in the compressive domain by using an improved sampling pattern. There are few studies regarding the applicability of CS to image fusion. The main purpose of this work is to explore the properties of compressive measurements through different sampling patterns and their potential use in image fusion. The study demonstrates that CS-based image fusion has a number of perceived advantages in comparison with image fusion in the multiresolution (MR) domain. The simulations show that the proposed CS-based image fusion algorithm provides promising results.

Compressive sensing for ultrasound RF echoes using a-Stable Distributions

This paper introduces a novel framework for compressive sensing of biomedical ultrasonic signals based on modelling data with stable distributions. We propose an approach to ℓp norm minimisation that employs the iteratively reweighted least squares (IRLS) algorithm but in which the parameter p is judiciously chosen by relating it to the characteristic exponent of the underlying alpha-stable distributed data. Our results show that the proposed algorithm, which we prefer to call S±S-IRLS, outperforms previously proposed ℓ1 minimisation algorithms, such as basis pursuit or orthogonal matching pursuit, both visually and in terms of PSNR.