Akram Aldroubi

In Vivo Fiber Tractography Using DT-MRI Data

Fiber tract trajectories in coherently organized brain white matter pathways were computed from in vivo diffusion tensor magnetic resonance imaging (DT‐MRI) data. First, a continuous diffusion tensor field is constructed from this discrete, noisy, measured DT‐MRI data. Then a Frenet equation, describing the evolution of a fiber tract, was solved. This approach was validated using synthesized, noisy DT‐MRI data. Corpus callosum and pyramidal tract trajectories were constructed and found to be consistent with known anatomy. The methods reliability, however, degrades where the distribution of fiber tract directions is nonuniform. Moreover, background noise in diffusion‐weighted MRIs can cause a computed trajectory to hop from tract to tract. Still, this method can provide quantitative information with which to visualize and study connectivity and continuity of neural pathways in the central and peripheral nervous systems in vivo, and holds promise for elucidating architectural features in other fibrous tissues and ordered media. Magn Reson Med 44:625–632, 2000. Published 2000 Wiley‐Liss, Inc.

B-spline signal processing. I. Theory

The use of continuous B-spline representations for signal processing applications such as interpolation, differentiation, filtering, noise reduction, and data compressions is considered. The B-spline coefficients are obtained through a linear transformation, which unlike other commonly used transforms is space invariant and can be implemented efficiently by linear filtering. The same property also applies for the indirect B-spline transform as well as for the evaluation of approximating representations using smoothing or least squares splines. The filters associated with these operations are fully characterized by explicitly evaluating their transfer functions for B-splines of any order. Applications to differentiation, filtering, smoothing, and least-squares approximation are examined. The extension of such operators for higher-dimensional signals such as digital images is considered. >

B-spline signal processing. II. Efficiency design and applications

For pt.I see ibid., vol.41, no.2, p.821-33 (1993). A class of recursive filtering algorithms for the efficient implementation of B-spline interpolation and approximation techniques is described. In terms of simplicity of realization and reduction of computational complexity, these algorithms compare favorably with conventional matrix approaches. A filtering interpretation (low-pass filter followed by an exact polynomial spline interpolator) of smoothing spline and least-squares approximation methods is proposed. These techniques are applied to the design of digital filters for cubic spline signal processing. An efficient implementation of a smoothing spline edge detector is proposed. It is also shown how to construct a cubic spline image pyramid that minimizes the loss of information in passage from one resolution level to the next. In terms of common measures of fidelity, this data structure appears to be superior to the Gaussian/Laplacian pyramid. >

A review of wavelets in biomedical applications

We present an overview of the various uses of the wavelet transform (WT) in medicine and biology. We start by describing the wavelet properties that are the most important for biomedical applications. In particular we provide an interpretation of the the continuous wavelet transform (CWT) as a prewhitening multiscale matched filter. We also briefly indicate the analogy between the WT and some of the the biological processing that occurs in the early components of the auditory and visual system. We then review the uses of the WT for the analysis of 1-D physiological signals obtained by phonocardiography, electrocardiography (ECG), mid electroencephalography (EEG), including evoked response potentials. Next, we provide a survey of wavelet developments in medical imaging. These include biomedical image processing algorithms (e.g., noise reduction, image enhancement, and detection of microcalcifications in mammograms), image reconstruction and acquisition schemes (tomography, and magnetic resonance imaging (MRI)), and multiresolution methods for the registration and statistical analysis of functional images of the brain (positron emission tomography (PET) and functional MRI (fMRI)). In each case, we provide the reader with same general background information and a brief explanation of how the methods work.

Nonuniform Sampling and Reconstruction in Shift-Invariant Spaces

This article discusses modern techniques for nonuniform sampling and reconstruction of functions in shift-invariant spaces. It is a survey as well as a research paper and provides a unified framework for uniform and nonuniform sampling and reconstruction in shift-invariant subspaces by bringing together wavelet theory, frame theory, reproducing kernel Hilbert spaces, approximation theory, amalgam spaces, and sampling. Inspired by applications taken from communication, astronomy, and medicine, the following aspects will be emphasized: (a) The sampling problem is well defined within the setting of shift-invariant spaces. (b) The general theory works in arbitrary dimension and for a broad class of generators. (c) The reconstruction of a function from any sufficiently dense nonuniform sampling set is obtained by efficient iterative algorithms. These algorithms converge geometrically and are robust in the presence of noise. (d) To model the natural decay conditions of real signals and images, the sampling theory is developed in weighted L p-spaces.

Fast B-spline transforms for continuous image representation and interpolation

Efficient algorithms for the continuous representation of a discrete signal in terms of B-splines (direct B-spline transform) and for interpolative signal reconstruction (indirect B-spline transform) with an expansion factor m are described. Expressions for the z-transforms of the sampled B-spline functions are determined and a convolution property of these kernels is established. It is shown that both the direct and indirect spline transforms involve linear operators that are space invariant and are implemented efficiently by linear filtering. Fast computational algorithms based on the recursive implementations of these filters are proposed. A B-spline interpolator can also be characterized in terms of its transfer function and its global impulse response (cardinal spline of order n). The case of the cubic spline is treated in greater detail. The present approach is compared with previous methods that are reexamined from a critical point of view. It is concluded that B-spline interpolation correctly applied does not result in a loss of image resolution and that this type of interpolation can be performed in a very efficient manner. >

The adaptive bases algorithm for intensity-based nonrigid image registration

Nonrigid registration of medical images is important for a number of applications such as the creation of population averages, atlas-based segmentation, or geometric correction of functional magnetic resonance imaging (IMRI) images to name a few. In recent years, a number of methods have been proposed to solve this problem, one class of which involves maximizing a mutual information (Ml)-based objective function over a regular grid of splines. This approach has produced good results but its computational complexity is proportional to the compliance of the transformation required to register the smallest structures in the image. Here, we propose a method that permits the spatial adaptation of the transformations compliance. This spatial adaptation allows us to reduce the number of degrees of freedom in the overall transformation, thus speeding up the process and improving its convergence properties. To develop this method, we introduce several novelties: 1) we rely on radially symmetric basis functions rather than B-splines traditionally used to model the deformation field; 2) we propose a metric to identify regions that are poorly registered and over which the transformation needs to be improved; 3) we partition the global registration problem into several smaller ones; and 4) we introduce a new constraint scheme that allows us to produce transformations that are topologically correct. We compare the approach we propose to more traditional ones and show that our new algorithm compares favorably to those in current use.

On the asymptotic convergence of B-spline wavelets to Gabor functions

A family of nonorthogonal polynomial spline wavelet transforms is considered. These transforms are fully reversible and can be implemented efficiently. The corresponding wavelet functions have a compact support. It is proven that these B-spline wavelets converge to Gabor functions (modulated Gaussian) pointwise and in all L/sub p/-norms with 1 >

A family of polynomial spline wavelet transforms

Abstract This paper presents an extension of the family of orthogonal Battle/Lemarie spline wavelet transforms with emphasis on filter bank implementation. Spline wavelets that are not necessarily orthogonal within the same resoluton level, are constructed by linear combination of polynomial spline wavelets of compact support, the natural counterpart of classical B-spline functions. Mallats fast wavelet transform algorithm is extended to deal with these non-orthogonal basis functions. The impulse and frequency responses of the corresponding analysis and synthesis filters are derived explicitly for polynomial splines of any order n (n odd). The link with the general framework of biorthogonal wavelet transforms is also made explicit. The special cases of orthogonal, B-spline, cardinal and dual wavelets are considered in greater detail. The B-spline (respectively dual) representation is associated with simple FIR binomial synthesis (respectively analysis) filters and recursive analysis (respectively synthesis) filters. The cardinal representation provides a sampled representation of the underlying continuous functions (interpolation property). The distinction between cardinal and orthogonal representation vanishes as the order of the spline is increased; both wavelets tend asymptotically to the bandlimited sinc-wavelet. The distinctive features of these various representations are discussed and illustrated with a texture analysis example.

Sampling procedures in function spaces and asymptotic equivalence with shannon's sampling theory

We view Shannons sampling procedure as a problem of approximation in the space S = {s: s (x) = (c * sinc)(x)c e l 2}. We show that under suitable conditions on a generating function λ e L 2, the approximation problem onto the space V = {v:v(x) = (c * λ)(x)c e l 2} produces a sampling procedure similar to the classical one. It consists of an optimal prefiltering, a pure jitter-stable sampling, and a postfiltering for the reconstruction. We describe equivalent signal representations using generic, dual, cardinal, and orthogonal basis functions and give the expression of the corresponding filters. We then consider sequences λn, where λn denotes the n-fold convolution of λ. They provide a sequence of increasingly regular sampling schemes as the value of n increases. We show that the cardinal and orthogonal pre- and postfilters associated with these sequences asymptotically converge to the ideal lowpass filter of Shannon. The theory is illustrated using several eamples.