Wilfried Philips

A versatile wavelet domain noise filtration technique for medical imaging

We propose a robust wavelet domain method for noise filtering in medical images. The proposed method adapts itself to various types of image noise as well as to the preference of the medical expert; a single parameter can be used to balance the preservation of (expert-dependent) relevant details against the degree of noise reduction. The algorithm exploits generally valid knowledge about the correlation of significant image features across the resolution scales to perform a preliminary coefficient classification. This preliminary coefficient classification is used to empirically estimate the statistical distributions of the coefficients that represent useful image features on the one hand and mainly noise on the other. The adaptation to the spatial context in the image is achieved by using a wavelet domain indicator of the local spatial activity. The proposed method is of low complexity, both in its implementation and execution time. The results demonstrate its usefulness for noise suppression in medical ultrasound and magnetic resonance imaging. In these applications, the proposed method clearly outperforms single-resolution spatially adaptive algorithms, in terms of quantitative performance measures as well as in terms of visual quality of the images.

Noise reduction by fuzzy image filtering

A new fuzzy filter is presented for the noise reduction of images corrupted with additive noise. The filter consists of two stages. The first stage computes a fuzzy derivative for eight different directions. The second stage uses these fuzzy derivatives to perform fuzzy smoothing by weighting the contributions of neighboring pixel values. Both stages are based on fuzzy rules which make use of membership functions. The filter can be applied iteratively to effectively reduce heavy noise. In particular, the shape of the membership functions is adapted according to the remaining noise level after each iteration, making use of the distribution of the homogeneity in the image. A statistical model for the noise distribution can be incorporated to relate the homogeneity to the adaptation scheme of the membership functions. Experimental results are obtained to show the feasibility of the proposed approach. These results are also compared to other filters by numerical measures and visual inspection.

An imaging system with calibrated color image acquisition for use in dermatology

The authors propose a novel imaging system useful in dermatology, more precisely, for the follow-up of patients with an increased risk of skin cancer. The system consists of a Pentium PC equipped with an RGB frame grabber, a three-chip charge coupled devices (CCD) camera controlled by the serial port and equipped with a zoom lens and a halogen annular light source. Calibration of the imaging system provides a way to transform the acquired images, which are defined in an unknown color space, to a standard, well-defined color space called SRGB, sRGB has a known relation to the CIE/sup 1/ XYZ and CIE L*a*b* colorimetric spaces. These CIE color spaces are based on the human vision, and they allow the computation of a color difference metric called CIE /spl Delta/E/sub ab/*, which is proportional to the color difference, as seen by a human observer. Several types of polynomial RGB to sRGB transforms will be tried, including some optimized in perceptually uniform color spaces. The use of a standard and well-defined color space also allows meaningful exchange of images, e.g., in teledermatology. The calibration procedure is based on 24 patches with known color properties, and it takes about 5 minutes to perform. It results in a number of settings called a profile that remains valid for tens of hours of operation. Such a profile is checked before acquiring images using just one color patch, and is adjusted on the fly to compensate for short-term drift in the response of the imaging system. Precision or reproducibility of subsequent color measurements is very good with =0.3 and /spl Delta/E/sub ab/*<1.2. Accuracy compared with spectrophotometric measurements is fair with =6.2 and /spl Delta/E/sub ab/*><13.3.

Hyperspectral and LiDAR Data Fusion: Outcome of the 2013 GRSS Data Fusion Contest

The 2013 Data Fusion Contest organized by the Data Fusion Technical Committee (DFTC) of the IEEE Geoscience and Remote Sensing Society aimed at investigating the synergistic use of hyperspectral and Light Detection And Ranging (LiDAR) data. The data sets distributed to the participants during the Contest, a hyperspectral imagery and the corresponding LiDAR-derived digital surface model (DSM), were acquired by the NSF-funded Center for Airborne Laser Mapping over the University of Houston campus and its neighboring area in the summer of 2012. This paper highlights the two awarded research contributions, which investigated different approaches for the fusion of hyperspectral and LiDAR data, including a combined unsupervised and supervised classification scheme, and a graph-based method for the fusion of spectral, spatial, and elevation information.

Improved Classification of VHR Images of Urban Areas Using Directional Morphological Profiles

Meter to submeter resolution satellite images have generated new interests in extracting man-made structures in the urban area. However, classification accuracies for such purposes are far from satisfactory. Spectral characteristics of urban land cover classes are so similar that they cannot be separated using only spectral information. As a result, there is an increased interest in incorporating geometrical information. One possible approach is the use of morphological profiles (MPs). In this paper, we introduce two improvements on the use of MPs. Current approaches use disk-shaped structuring elements (SEs) to derive an MP. This profile contains information about the minimum dimension of objects. In this paper, we extend this approach by using linear SEs. This results in a profile containing information about the maximum object dimension. We show that the addition of the line-based MP gives a substantial improvement of the classification result. A second improvement is achieved by using ldquopartial morphological reconstructionrdquo instead of the normal morphological reconstruction. Morphological reconstruction is commonly used to better preserve the shape of objects. However, we show that this leads to ldquoover-reconstructionrdquo in typical remote sensing images and a decreased classification performance. With ldquopartial reconstruction,rdquo we are able to overcome this problem and still preserve the shape of objects.

Wavelet-Domain Video Denoising Based on Reliability Measures

This paper proposes a novel video denoising method based on nondecimated wavelet band filtering. In the proposed method, motion estimation and adaptive recursive temporal filtering are performed in a closed loop, followed by an intra-frame spatially adaptive filter. All processing occurs in the wavelet domain. The paper introduces new wavelet-based motion reliability measures. We make a difference between motion reliability per orientation and reliability per wavelet band. These two reliability measures are employed in different stages of the proposed denoising scheme. The reliability per orientation (horizontal and vertical) measure is used in the proposed motion estimation scheme while the reliability of the estimated motion vectors (MVs) per wavelet band is utilized for subsequent adaptive temporal and spatial filtering. We propose a novel cost function for motion estimation which takes into account the spatial orientation of image structures and their motion matching values. Our motion estimation approach is a novel wavelet-domain three-step scheme, where the refinement of MVs in each step is determined based on the proposed motion reliabilities per orientation. The temporal filtering is performed separately in each wavelet band along the estimated motion trajectory and the parameters of the temporal filter depend on the motion reliabilities per wavelet band. The final spatial filtering step employs an adaptive smoothing of wavelet coefficients that yields a stronger filtering at the positions where the temporal filter was less effective. The results on various grayscale sequences demonstrate that the proposed filter outperforms several state-of-the-art filters visually (as judged by a small test panel) as well as in terms of peak signal-to-noise ratio

Semisupervised Local Discriminant Analysis for Feature Extraction in Hyperspectral Images

We propose a novel semisupervised local discriminant analysis method for feature extraction in hyperspectral remote sensing imagery, with improved performance in both ill-posed and poor-posed conditions. The proposed method combines unsupervised methods (local linear feature extraction methods and supervised method (linear discriminant analysis) in a novel framework without any free parameters. The underlying idea is to design an optimal projection matrix, which preserves the local neighborhood information inferred from unlabeled samples, while simultaneously maximizing the class discrimination of the data inferred from the labeled samples. Experimental results on four real hyperspectral images demonstrate that the proposed method compares favorably with conventional feature extraction methods.

MRI Segmentation of the Human Brain: Challenges, Methods, and Applications

Image segmentation is one of the most important tasks in medical image analysis and is often the first and the most critical step in many clinical applications. In brain MRI analysis, image segmentation is commonly used for measuring and visualizing the brains anatomical structures, for analyzing brain changes, for delineating pathological regions, and for surgical planning and image-guided interventions. In the last few decades, various segmentation techniques of different accuracy and degree of complexity have been developed and reported in the literature. In this paper we review the most popular methods commonly used for brain MRI segmentation. We highlight differences between them and discuss their capabilities, advantages, and limitations. To address the complexity and challenges of the brain MRI segmentation problem, we first introduce the basic concepts of image segmentation. Then, we explain different MRI preprocessing steps including image registration, bias field correction, and removal of nonbrain tissue. Finally, after reviewing different brain MRI segmentation methods, we discuss the validation problem in brain MRI segmentation.

Hex-splines: a novel spline family for hexagonal lattices

This paper proposes a new family of bivariate, nonseparable splines, called hex-splines, especially designed for hexagonal lattices. The starting point of the construction is the indicator function of the Voronoi cell, which is used to define in a natural way the first-order hex-spline. Higher order hex-splines are obtained by successive convolutions. A mathematical analysis of this new bivariate spline family is presented. In particular, we derive a closed form for a hex-spline of arbitrary order. We also discuss important properties, such as their Fourier transform and the fact they form a Riesz basis. We also highlight the approximation order. For conventional rectangular lattices, hex-splines revert to classical separable tensor-product B-splines. Finally, some prototypical applications and experimental results demonstrate the usefulness of hex-splines for handling hexagonally sampled data.

Rosette Tracker: an open source image analysis tool for automatic quantification of genotype effects

Image analysis of Arabidopsis (Arabidopsis thaliana) rosettes is an important nondestructive method for studying plant growth. Some work on automatic rosette measurement using image analysis has been proposed in the past but is generally restricted to be used only in combination with specific high-throughput monitoring systems. We introduce Rosette Tracker, a new open source image analysis tool for evaluation of plant-shoot phenotypes. This tool is not constrained by one specific monitoring system, can be adapted to different low-budget imaging setups, and requires minimal user input. In contrast with previously described monitoring tools, Rosette Tracker allows us to simultaneously quantify plant growth, photosynthesis, and leaf temperature-related parameters through the analysis of visual, chlorophyll fluorescence, and/or thermal infrared time-lapse sequences. Freely available, Rosette Tracker facilitates the rapid understanding of Arabidopsis genotype effects.