Wavelet-Galerkin solution of boundary value problems

International audienceIn this paper we review the application of wavelets to the solution of partial differential equations. We consider in detail both the single scale and the multiscale Wavelet Galerkin method. The theory of wavelets is described here using the language and mathematics of signal processing. We show a method of adapting wavelets to an interval using an extrapolation technique called Wavelet Extrapolation. Wavelets on an interval allow boundary conditions to be enforced in partial differential equations and image boundary problems to be overcome in image processing. Finally, we discuss the fast inversion of matrices arising from differential operators by preconditioning the multiscale wavelet matrix. Wavelet preconditioning is shown to limit the growth of the matrix condition number, such that Krylov subspace iteration methods can accomplish fast inversion of large matrices

Multiscale Astronomical Image Processing Based on Nonlinear Partial Differential Equations

Astronomical applications of recent advances in the field of nonastronomical image processing are presented. These innovative methods, applied to multiscale astronomical images, increase signal-to-noise ratio, do not smear point sources or extended diffuse structures, and are thus a highly useful preliminary step for detection of different features including point sources, smoothing of clumpy data, and removal of contaminants from background maps. We show how the new methods, combined with other algorithms of image processing, unveil fine diffuse structures while at the same time enhance detection of localized objects, thus facilitating interactive morphology studies and paving the way for the automated recognition and classification of different features. We have also developed a new application framework for astronomical image processing that implements some recent advances made in computer vision and modern image processing, along with original algorithms based on nonlinear partial differential equations. The framework enables the user to easily set up and customize an image-processing pipeline interactively; it has various common and new visualization features and provides access to many astronomy data archives. Altogether, the results presented here demonstrate the first implementation of a novel synergistic approach based on integration of image processing, image visualization, and image quality assessment

Cross-diffusion based filtering as pre-processing step for remote sensing procedures

A new methodology combining 2 × 2 cross-diffusion systems of nonlinear partial differential equations (CDS) with classical image classification procedures is proposed in the present paper. Such a kind of mathematical models (CDS) have been theoretically studied in previous works in the context of image processing, however here they are tested and stressed in very practical instances. In particular, the main contribution of this paper is the improvement of the classification of satellite images when they are previously filtered by means of a CDS model. This conclusion is based on a wide and costly experimentation with satellite images of areas damaged by forest fires and surface coal mining, all of them located in Mediterranean areas. The efficiency of our metho- dology is not only in terms of the classification improvement but also in terms of the runtime saving since CDS based filtering is much less costly than other classical partial differential equations based filtering mathematical models as for example anisotropic models or higher order ones, always within the framework of nonlinear partial differential equations

Image Processing Application for Cognition: IPAC Architecture and Implementation in Java

An application framework for advanced image processing and visualization is presented. It provides common two-dimensional operators and implements recent developments in the field of image processing as well as original algorithms based on nonlinear partial differential equations (PDEs). It is platform independent and has the capability of extensibility. This objective is achieved by exploiting the object-oriented paradigm. A graphical user interface (GUI) provides processing, analysis and visualization in a highly integrated, easy to use environment. Applications of the developed system to images obtained by the Spitzer Space Telescope are demonstrated

ADI splitting schemes for a fourth-order nonlinear partial differential equation from image processing

We present directional operator splitting schemes for the numerical solution of a fourth-order, nonlinear partial differential evolution equation which arises in image processing. This equation constitutes the H−1-gradient flow of the total variation and represents a prototype of higher-order equations of similar type which are popular in imaging for denoising, deblurring and inpainting problems. The efficient numerical solution of this equation is very challenging due to the stiffness of most numerical schemes. We show that the combination of directional splitting schemes with implicit time-stepping provides a stable and computationally cheap numerical realisation of the equation

Introduction to the Special Issue on Partial Differential Equations and Geometry-Driven Diffusion in Image Processing and Analysis

©1998 IEEE. Personal use of this material is permitted. However, permission to reprint/republish this material for advertising or promotional purposes or for creating new collective works for resale or distribution to servers or lists, or to reuse any copyrighted component of this work in other works must be obtained from the IEEE. This material is presented to ensure timely dissemination of scholarly and technical work. Copyright and all rights therein are retained by authors or by other copyright holders. All persons copying this information are expected to adhere to the terms and constraints invoked by each author's copyright. In most cases, these works may not be reposted without the explicit permission of the copyright holder.DOI: 10.1109/TIP.1998.66117

Exact asymptotics of the uniform error of interpolation by multilinear splines

The question of adaptive mesh generation for approximation by splines has been studied for a number of years by various authors. The results have numerous applications in computational and discrete geometry, computer aided geometric design, finite element methods for numerical solutions of partial differential equations, image processing, and mesh generation for computer graphics, among others. In this paper we will investigate the questions regarding adaptive approximation of C2 functions with arbitrary but fixed throughout the domain signature by multilinear splines. In particular, we will study the asymptotic behavior of the optimal error of the weighted uniform approximation by interpolating and quasi-interpolating multilinear splines

Anisotropic Diffusion Partial Differential Equations in Multi-Channel Image Processing : Framework and Applications

We review recent methods based on diffusion PDE's (Partial Differential Equations) for the purpose of multi-channel image regularization. Such methods have the ability to smooth multi-channel images anisotropically and can preserve then image contours while removing noise or other undesired local artifacts. We point out the pros and cons of the existing equations, providing at each time a local geometric interpretation of the corresponding processes. We focus then on an alternate and generic tensor-driven formulation, able to regularize images while specifically taking the curvatures of local image structures into account. This particular diffusion PDE variant is actually well suited for the preservation of thin structures and gives regularization results where important image features can be particularly well preserved compared to its competitors. A direct link between this curvature-preserving equation and a continuous formulation of the Line Integral Convolution technique (Cabral and Leedom, 1993) is demonstrated. It allows the design of a very fast and stable numerical scheme which implements the multi-valued regularization method by successive integrations of the pixel values along curved integral lines. Besides, the proposed implementation, based on a fourth-order Runge Kutta numerical integration, can be applied with a subpixel accuracy and preserves then thin image structures much better than classical finite-differences discretizations, usually chosen to implement PDE-based diffusions. We finally illustrate the efficiency of this diffusion PDE's for multi-channel image regularization - in terms of speed and visual quality - with various applications and results on color images, including image denoising, inpainting and edge-preserving interpolation