Maximum-likelihood estimation of Rician distribution parameters

The problem of parameter estimation from Rician distributed data (e.g., magnitude magnetic resonance images) is addressed. The properties of conventional estimation methods are discussed and compared to maximum-likelihood (ML) estimation which is known to yield optimal results asymptotically. In contrast to previously proposed methods, ML estimation is demonstrated to be unbiased for high signal-to-noise ratio (SNR) and to yield physical relevant results for low SNR     

Watershed-based segmentation of 3D MR data for volume quantization

The aim of this work is the development of a semiautomatic segmentation technique for efficient and accurate volume quantization of Magnetic Resonance (MR) data. The proposed technique uses a 3D variant of Vincent and Soilles immersion-based watershed algorithm that is applied to the gradient magnitude of the MR data and that produces small volume primitives. The known drawback of the watershed algorithm, oversegmentation, is strongly reduced by a priori application of a 3D adaptive anisotropic diffusion filter to the MR data. Furthermore, oversegmentation is a posteriori reduced by properly merging small volume primitives that have similar gray level distributions. The outcome of the proceeding image processing steps is presented to the user for manual segmentation. Through selection of volume primitives, the user quickly segments of first slice, which contains the object of interest. Afterwards, the subsequent slices are automatically segmented by extrapolation. Segmentation results are contingently manually corrected. The proposed segmentation technique is tested on phantom objects, where segmentation errors less than 2% are observed. In addition, the technique is demonstrated on 3D MR data of the mouse head from which the cerebellum is extracted. Volumes of the mouse cerebellum and the mouse brains in toto are calculated.     

Wavelet denoising of multicomponent images using gaussian scale mixture models and a noise-free image as priors.

In this paper, a Bayesian wavelet-based denoising procedure for multicomponent images is proposed. A denoising procedure is constructed that (1) fully accounts for the multicomponent image covariances, (2) makes use of Gaussian scale mixtures as prior models that approximate the marginal distributions of the wavelet coefficients well, and (3) makes use of a noise-free image as extra prior information. It is shown that such prior information is available with specific multicomponent image data of, e.g., remote sensing and biomedical imaging. Experiments are conducted in these two domains, in both simulated and real noisy conditions.     

On the Geometry of Multivariate Generalized Gaussian Models

This paper concerns the geometry of the zero-mean multivariate generalized Gaussian distribution (MGGD) and the calculation of geodesic distances on the MGGD manifold. The MGGD is a suitable distribution for the modeling of multivariate (color, multispectral, vector and tensor images, etc.) image wavelet statistics. Expressions are derived for the Fisher-Rao metric for the zero-mean MGGD model. A closed-form expression is obtained for the geodesic distance on the submanifolds characterized by a fixed MGGD shape parameter. Suitable approximate solutions to the geodesic equations are presented in the case of MGGDs with varying shape parameters. An application to image texture similarity measurement in the wavelet domain is briefly discussed, comparing the performance of the geodesic distance and the Kullback-Leibler divergence.     

Fusion and merging of multispectral images with use of multiscale fundamental forms

A new multispectral image wavelet representation is introduced, based on multiscale fundamental forms. This representation describes gradient information of multispectral images in a multiresolution framework. The representation is, in particular, extremely suited for fusion and merging of multispectral images. For fusion as well as for merging, a strategy is described. Experiments are performed on multispectral images, where Landsat Thematic Mapper images are fused and merged with SPOT Panchromatic images. The proposed techniques are compared with wavelet-based techniques described in the literature.     

An orthogonal wavelet representation of multivalued images

A new orthogonal wavelet representation of multivalued images is presented. The idea for this representation is based on the concept of maximal gradient of multivalued images. This concept is generalized from gradients toward linear vector operators in the image plane with equal components along rows and columns. Using this generalization, the pyramidal dyadic wavelet transform algorithm using quadrature mirror filters is modified to be applied to multivalued images. This results in a representation of a single image, containing multiscale detail information from all component images involved. This representation leads to multiple applications ranging from multispectral image fusion to color and multivalued image enhancement, denoising and segmentation. In this paper, the representation is applied for fusion of images. More in particular, we introduce a scheme to merge high spatial resolution greylevel images with low spatial resolution multivalued images to improve spatial resolution of the latter while preserving spectral resolution. Two applications are studied: demosaicing of color images and merging of multispectral remote sensing images.     

Classification of heathland vegetation in a hierarchical contextual framework

Heathlands in Western Europe have shown dramatic declines over the last century and therefore have been given a high conservation priority in the Habitats Directive of the European Union (EU). Accurate surveying and monitoring of heathland habitats is essential for appropriate conservation management, but the large heterogeneity of vegetation types within habitats as well as the occurrence of similar vegetation across habitat types hinders a straightforward, automated mapping based on aerial images. In such a case, a context-dependent classification algorithm is expected to be superior to traditional classification techniques. This article presents a novel approach to map the conservation status of heathland vegetation by using a hierarchical classification scheme that describes the structural dependencies in the field between the basic vegetation and the land-cover types that habitats are composed of. These dependency relationships are included as contextual information in the classification process, using a tree-structured Markov random field (TS-MRF) technique with a tree that reflects the hierarchy of the classification scheme. Results of this approach for a heathland area in Belgium were compared with results from more conventional classification approaches. Validation of the results showed that the structure of the scheme contained important spatial relationships, which were further reinforced by using the contextual classification strategy, especially for the most detailed level of the classification scheme. Accuracy increased and the classification results were more suitable for visual interpretation.     

Geodesics on the Manifold of Multivariate Generalized Gaussian Distributions with an Application to Multicomponent Texture Discrimination

We consider the Rao geodesic distance (GD) based on the Fisher information as a similarity measure on the manifold of zero-mean multivariate generalized Gaussian distributions (MGGD). The MGGD is shown to be an adequate model for the heavy-tailed wavelet statistics in multicomponent images, such as color or multispectral images. We discuss the estimation of MGGD parameters using various methods. We apply the GD between MGGDs to color texture discrimination in several classification experiments, taking into account the correlation structure between the spectral bands in the wavelet domain. We compare the performance, both in terms of texture discrimination capability and computational load, of the GD and the Kullback-Leibler divergence (KLD). Likewise, both uni- and multivariate generalized Gaussian models are evaluated, characterized by a fixed or a variable shape parameter. The modeling of the interband correlation significantly improves classification efficiency, while the GD is shown to consistently outperform the KLD as a similarity measure.     

Statistical texture characterization from discrete waveletrepresentations

We conjecture that texture can be characterized by the statistics of the wavelet detail coefficients and therefore introduce two feature sets: (1) the wavelet histogram signatures which capture all first order statistics using a model based approach and (2) the wavelet co-occurrence signatures, which reflect the coefficients' second-order statistics. The introduced feature sets outperform the traditionally used energy. Best performance is achieved by combining histogram and co-occurrence signatures