Nonlinear Diffusion. Journal Club Presentation. Xiaowei Zhou
|
|
- Leonard Lucas
- 6 years ago
- Views:
Transcription
1 1 / 41 Journal Club Presentation Xiaowei Zhou Department of Electronic and Computer Engineering The Hong Kong University of Science and Technology
2 2 / 41 Outline 1 Motivation Diffusion process Linear diffusion 2 3 Other applications of PDE based methods
3 3 / 41 Outline Motivation Diffusion process Linear diffusion 1 Motivation Diffusion process Linear diffusion 2 3 Other applications of PDE based methods
4 4 / 41 What are the problems? Motivation Diffusion process Linear diffusion Figure: Left: a noisy picture; Right: ultrasound image to examine the heart of zebrafish
5 5 / 41 Why diffusion Motivation Diffusion process Linear diffusion How to reduce noises and restore the original image contents. Traditional linear filtering is unable to preserve "semantically meaningful" features while reducing the noise. Key idea of diffusion: evolve the image as a dynamic system to reduce system randomness. Some applications of diffusion Image restoration Scale space theory feature enhancement...
6 6 / 41 Outline Motivation Diffusion process Linear diffusion 1 Motivation Diffusion process Linear diffusion 2 3 Other applications of PDE based methods
7 Motivation Diffusion process Linear diffusion Definition of diffusion Diffusion In physics, diffusion is a process that equilibrates concentration differences without creating or destroying mass. Examples: heat diffusion. molecular motion. An animation of molecular diffusion In image processing, the idea of diffusion is borrowed to reduce the random variation and restore the image. Space: 2D image plane Concentration: Image intensity 7 / 41
8 8 / 41 Laws for diffusion Motivation Diffusion process Linear diffusion Fick s Law Intensity gradient u generates flux j. j = D u (1) u: Intensity (mass density, temperature, image intensity...) D: diffusion tensor, symmetric positive definite if D = ki: flux j is always parallel to u else: flux j is not always parallel to u Continuity Equation no mass generated or destroyed. tu = div(j) (2)
9 9 / 41 Diffusion Equation Motivation Diffusion process Linear diffusion Fick s law and continuity equation yeild: Diffusion Equation j tu = div(d u) u(x, 0) = f (x) The diffusion process is determined by intensity gradient and the tensor field. In image processing: boundary condition is needed: < D u, n >= 0 (no normal flow on boundary) D is critical to design a diffusion method. For linear diffusion, D is constant over image. For nonlinear diffusion, D is adapted to local image structure. (3)
10 10 / 41 Outline Motivation Diffusion process Linear diffusion 1 Motivation Diffusion process Linear diffusion 2 3 Other applications of PDE based methods
11 11 / 41 Solution for linear diffusion Motivation Diffusion process Linear diffusion Solve the simplest case: D = ki { t u = k div( u) = k u u(x, 0) = f (x) (4) D is independent with the image structure and isotropic. Solution: u(x, t) = φ(x, t) f (x) φ(x, t). = 1 4πkt n exp( x T x 4kt ) (5) φ is a Gaussian function with σ = 2kt! That means the solution of linear diffusion process is Gaussian filtering with σ increasing with time.
12 12 / 41 Drawbacks of linear diffusion Motivation Diffusion process Linear diffusion the Gaussian process will distort the image structure Figure: Left: linear diffusion result at different time; Right: corresponding zero-crossing points. Dislocation of edges occurs at coarse scales.
13 13 / 41 Outline 1 Motivation Diffusion process Linear diffusion 2 3 Other applications of PDE based methods
14 14 / 41 Perona and Malik modified the diffusion coefficient D to be adapted { to the image: large, when u is small inside a homogeneous region D = small, when u is large near boundaries between regions Perona and Malik s Anisotropic Diffusion j tu = div(g( u 2 ) u) u(x, 0) = f (x) where, g( ) is a nonnegative monotonically decreasing function with g(0) = 1 P. Perona, J. Malik Scale space and edge detection using anisotropic diffusion IEEE Trans. Pattern Anal. Mach. Intell., Vol. 12, 629 C639, (6)
15 15 / 41 Diffusion coefficient Requirements on the diffusion coefficient function: nonnegative: satisfy the maximum principle to avoid generating spurious features in coarse levels of the scale space. monotonically decreasing: discourage the diffusion over edges where the image gradient is large. A such function popularly used: g(s 2 ) = s 2 /λ 2 (7)
16 16 / 41 Diffusivity and Flux function g(s 2 ) is usually called diffusivity, describing the diffusion or conduction property. s g(s 2 ) is named flux function, describing magnitude of flux. Figure: Left: diffusivity in Eq.(7); Right: the corresponding flux function
17 17 / 41 Example result of anisotropic diffusion Figure: Left: original noisy image; Middle: after linear diffusion; Right: after nonlinear diffusion
18 18 / 41 Example result of anisotropic diffusion Figure: Left: original noisy image; Middle: after 100 iterations; Right: after 1000 iterations
19 19 / 41 Edge enhancement property We can analyze the property of the nonlinear diffusion in Eq.(6) first from 1-D case. Define the flux function: φ(s) = s g(s 2 ). For certain type of g(s 2 ) (Page 16): In 1-D space, Eq.(6) reads: φ (s) 0, s λ φ (s) < 0, s > λ t u = φ (u x )u xx (8) If we check the derivative of the magnitude of gradient with respective to time: t (u 2 x ) = 2u x x (u t ) = 2φ (u x )u x u 2 xx + 2φ (u x )u x u xxx (9)
20 20 / 41 Edge enhancement property cont. At a edge in a bounded image, we have following assumputions: u xx = 0 u x u xxx < 0 From Eq.(9), it can be found that: if u x > λ,φ (u x ) < 0, t (u 2 x ) > 0 if u x < λ,φ (u x ) > 0, t (u 2 x ) < 0 It means the gradient at edges where u x > λ will get larger. The edge is enhanced.
21 21 / 41 Edge enhancement property cont. In 2D case: t u = div(g( u 2 ) u) = g( u 2 ) u + g( u 2 ) T u = g( u 2 ) u + 2g ( u 2 )u ηη = g( u 2 )u ξξ + [g( u 2 ) + 2 u 2 g ( u 2 )]u ηη = g( u 2 )u ξξ + φ ( u )u ηη (10) where u ξξ and u ηη are the second derivative along directions perpendicular and parallel to u x respectively: H u is the Hessian matrix. u ξξ = ξ T H u ξ u ηη = η T H u η
22 22 / 41 Edge enhancement property cont. Eq.(10) gives us more insights into 2D diffusion: The 2D diffusion can be decomposed to process along two orthogonal directions along the local edge and normal to the edge. if u x << λ, the diffusion is the ordinary isotropic diffusion driven by Laplacian. if u x λ, φ ( u ) 0, the diffusion u ηη parallel to the image gradient (normal to the edge) is suppressed. if u x >> λ, φ ( u ) < 0, backward diffusion happens normal to the edge!
23 23 / 41 Spatial regularization The original Perona-Malik model has the following problem: in the beginning of the diffusion, the diffusion coefficient can t be estimated robustly since the high variation in a noisy image. This issue can be numerically solved by adding a spatial filtering when estimating the gradient, and the diffusion coefficient is modified as: where G σ is a Gaussian mask. D = g( G σ u 2 ) (11)
24 24 / 41 Outline 1 Motivation Diffusion process Linear diffusion 2 3 Other applications of PDE based methods
25 25 / 41 An example: edge-enhancing Anisotropic Diffusion Using a anisotropic tensor D, such that D reflect the edge structure: Eigenvectors: v 1 u σ, v 2 u σ Eigenvalues: damped diffusion across edge: λ 1 = g( u σ 2 ) full diffusion parallel to edge: λ 2 = 1 J. Weickert Theoretical foundations of anisotropic diffusion in image processing Computing Supplement; Vol. 11 Proceedings of the 7th TFCV on Theoretical Foundations of Computer Vision, 1996.
26 26 / 41 result of edge-enhancing anisotropic diffusion Figure: Left: original noisy image; Right: after edge-enhancing anisotropic diffusion. Adopted from Weickert 1996
27 27 / 41 Another example J. Weickert Coherence-Enhancing Diffusion Filtering International Journal of Computer Vision, Make use of the image structure tensor K ρ ( u σ u T σ )
28 28 / 41 Outline 1 Motivation Diffusion process Linear diffusion 2 3 Other applications of PDE based methods
29 29 / 41 SRAD Yu, Y. and Acton, ST Speckle reducing anisotropic diffusion IEEE Transactions on Image Processing, 2002, vol. 11, issue 11, pp Basic formulation: the same with Eq.(6) modification:the diffusion coefficient is modified to be: g(s 2 ) = [s 2 s 2 0 ]/[s2 0 (1 + s2 0 )] (12) s 0 and s are local statistics of the homogeneous region and the region of interest respectively. s 2 = var(u) ū 2 (13)
30 30 / 41 Comments on SRAD SRAD is intrinsically the same with the original model of Perona and Malik. the gradient magnitude estimator G σ u 2 in the Perona-Malik model is replaced by local statistics. This idea is borrowed from the Lee and Frost filters, which are popular adaptive filters in the radar community to reduce speckle noises in radar images. In the Lee and Frost filters, this statistics is used to determine the coefficient of smoothing. To determine s 0 in Eq.(12), a homogeneous region needs to be chosen. For fully automatic implementation, in this paper s 0 is set to be decreasing exponentially with the diffusion process going on. s 0 (t) = s 0 (0) exp( ρt).
31 31 / 41 Example results of SRAD For the phantom data: Figure: Left: original speckle image; Middle: result of SRAD; Right: result of the Perona-Malik model
32 32 / 41 Example results of SRAD For the zebrafish data: Figure: Left: original speckle image; Middle: result of SRAD; Right: result of the Perona-Malik model
33 33 / 41 Outline 1 Motivation Diffusion process Linear diffusion 2 3 Other applications of PDE based methods
34 34 / 41 Energy minimization The diffusion process can be viewed as minimizing a potential function evaluating the total variation of the image. If we can find a potential function ψ( ) such that: ( u) (ψ( u )) = g( u 2 ) u (14) ( u) means gradient respect to u. Then minimizing the energy functional: E(u) =. ψ( u )dxdy (15) Ω with the Euler-Lagrange equation leads to the gradient descent: t u = div( u (ψ( u ))) = div(g( u 2 ) u) (16)
35 35 / 41 diffusivity vs. potential Linear diffusion: g(s 2 ) = 1, ψ(s) = s2 2 the penalty function is exactly 2 norm. Perona-Malik model: g(s 2 ) = 1 1+s 2 /λ 2, ψ(s) = λ2 2 log(1 + ( s λ) 2) the penalty function is a robust estimator. Generally, only need to find ψ (s) = sg(s 2 ) Thus, diffusion is minimizing the total variation by gradient descent. The diffusion coefficient is determined by the penalty function (estimator) for the variation.
36 36 / 41 Some other comments The time of diffusion controls the fidelity to the original image. It is corresponding to the coefficient balancing the data term and regularization in ordinary regularized energy functions. Instead aiming at the objective function, designing a diffusion function (tensor) is more direct and easy to understand in the view of physics. (e.g. the diffusion with anisotropic tensor is difficult to be expressed as an energy minimization form) Why in robust estimation theory ψ (s) is usually called the flux function? Because it is exactly the flux function of diffusion defined in Page 16.
37 37 / 41 Outline Other applications of PDE based methods 1 Motivation Diffusion process Linear diffusion 2 3 Other applications of PDE based methods
38 38 / 41 Segmentation Other applications of PDE based methods Chan and Vese s active contour without edges: E CV (c 1, c 2, φ) = λ 1 c 1 I 0 (x, y) 2 H(φ(x, y))dxdy Ω +λ 2 Ω c 2 I 0 (x, y) 2 H( φ(x, y))dxdy +µ H(φ(x, y)) dxdy + ν H(φ(x, y))dxdy (17) Ω The gradient descent flow: [ ( ) ] φ φ = δ ε (φ) µ div ν λ 1 (u 0 c 1 ) 2 + λ 2 (u 0 c 2 ) 2 t φ φ(0, x, y) = φ 0 (x, y) (18) Ω
39 39 / 41 Motion tracking Other applications of PDE based methods Brox, T. et al. High accuracy optical flow estimation based on a theory for warping Proc. 8th European Conference on Computer Vision The energy functional with respect to the displacement field: E(d = (u, v)) = ρ( I 2 (x + d) I 1 (x) 2 + γ I 2 (x + d) I 1 (x) 2 )dx Ω +α ρ( u 2 + v 2 )dx (19) Ω The gradient descent flow (two complicated, please refer to the paper): u t = x E ux + y E uy E u v t = x E vx + y E vy E v
40 40 / 41 Outline Other applications of PDE based methods 1 Motivation Diffusion process Linear diffusion 2 3 Other applications of PDE based methods
41 41 / 41 Other applications of PDE based methods Why PDE based methods model and explore the image structure mathematical foundation, concise representation physically meaningful and understandable numerically implementable...
NON-LINEAR DIFFUSION FILTERING
NON-LINEAR DIFFUSION FILTERING Chalmers University of Technology Page 1 Summary Introduction Linear vs Nonlinear Diffusion Non-Linear Diffusion Theory Applications Implementation References Page 2 Introduction
More informationNONLINEAR DIFFUSION PDES
NONLINEAR DIFFUSION PDES Erkut Erdem Hacettepe University March 5 th, 0 CONTENTS Perona-Malik Type Nonlinear Diffusion Edge Enhancing Diffusion 5 References 7 PERONA-MALIK TYPE NONLINEAR DIFFUSION The
More informationNonlinear Diffusion. 1 Introduction: Motivation for non-standard diffusion
Nonlinear Diffusion These notes summarize the way I present this material, for my benefit. But everything in here is said in more detail, and better, in Weickert s paper. 1 Introduction: Motivation for
More informationImage enhancement. Why image enhancement? Why image enhancement? Why image enhancement? Example of artifacts caused by image encoding
13 Why image enhancement? Image enhancement Example of artifacts caused by image encoding Computer Vision, Lecture 14 Michael Felsberg Computer Vision Laboratory Department of Electrical Engineering 12
More informationLinear Diffusion and Image Processing. Outline
Outline Linear Diffusion and Image Processing Fourier Transform Convolution Image Restoration: Linear Filtering Diffusion Processes for Noise Filtering linear scale space theory Gauss-Laplace pyramid for
More informationLinear Diffusion. E9 242 STIP- R. Venkatesh Babu IISc
Linear Diffusion Derivation of Heat equation Consider a 2D hot plate with Initial temperature profile I 0 (x, y) Uniform (isotropic) conduction coefficient c Unit thickness (along z) Problem: What is temperature
More informationTensor-based Image Diffusions Derived from Generalizations of the Total Variation and Beltrami Functionals
Generalizations of the Total Variation and Beltrami Functionals 29 September 2010 International Conference on Image Processing 2010, Hong Kong Anastasios Roussos and Petros Maragos Computer Vision, Speech
More informationITK Filters. Thresholding Edge Detection Gradients Second Order Derivatives Neighborhood Filters Smoothing Filters Distance Map Image Transforms
ITK Filters Thresholding Edge Detection Gradients Second Order Derivatives Neighborhood Filters Smoothing Filters Distance Map Image Transforms ITCS 6010:Biomedical Imaging and Visualization 1 ITK Filters:
More informationMedical Image Analysis
Medical Image Analysis CS 593 / 791 Computer Science and Electrical Engineering Dept. West Virginia University 23rd January 2006 Outline 1 Recap 2 Edge Enhancement 3 Experimental Results 4 The rest of
More informationRIGOROUS MATHEMATICAL INVESTIGATION OF A NONLINEAR ANISOTROPIC DIFFUSION-BASED IMAGE RESTORATION MODEL
Electronic Journal of Differential Equations, Vol. 2014 (2014), No. 129, pp. 1 9. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu RIGOROUS MATHEMATICAL
More informationA Nonlocal p-laplacian Equation With Variable Exponent For Image Restoration
A Nonlocal p-laplacian Equation With Variable Exponent For Image Restoration EST Essaouira-Cadi Ayyad University SADIK Khadija Work with : Lamia ZIAD Supervised by : Fahd KARAMI Driss MESKINE Premier congrès
More informationErkut Erdem. Hacettepe University February 24 th, Linear Diffusion 1. 2 Appendix - The Calculus of Variations 5.
LINEAR DIFFUSION Erkut Erdem Hacettepe University February 24 th, 2012 CONTENTS 1 Linear Diffusion 1 2 Appendix - The Calculus of Variations 5 References 6 1 LINEAR DIFFUSION The linear diffusion (heat)
More informationIntroduction to Nonlinear Image Processing
Introduction to Nonlinear Image Processing 1 IPAM Summer School on Computer Vision July 22, 2013 Iasonas Kokkinos Center for Visual Computing Ecole Centrale Paris / INRIA Saclay Mean and median 2 Observations
More informationMedical Image Analysis
Medical Image Analysis CS 593 / 791 Computer Science and Electrical Engineering Dept. West Virginia University 20th January 2006 Outline 1 Discretizing the heat equation 2 Outline 1 Discretizing the heat
More informationEnergy-Based Image Simplification with Nonlocal Data and Smoothness Terms
Energy-Based Image Simplification with Nonlocal Data and Smoothness Terms Stephan Didas 1, Pavel Mrázek 2, and Joachim Weickert 1 1 Mathematical Image Analysis Group Faculty of Mathematics and Computer
More informationConsistent Positive Directional Splitting of Anisotropic Diffusion
In: Boštjan Likar (ed.): Proc. of Computer Vision Winter Workshop, Bled, Slovenia, Feb. 7-9, 2001, pp. 37-48. Consistent Positive Directional Splitting of Anisotropic Diffusion Pavel Mrázek and Mirko Navara
More informationTotal Variation Image Edge Detection
Total Variation Image Edge Detection PETER NDAJAH Graduate School of Science and Technology, Niigata University, 8050, Ikarashi 2-no-cho, Nishi-ku, Niigata, 950-28, JAPAN ndajah@telecom0.eng.niigata-u.ac.jp
More informationA Partial Differential Equation Approach to Image Zoom
A Partial Differential Equation Approach to Image Zoom Abdelmounim Belahmidi and Frédéric Guichard January 2004 Abstract We propose a new model for zooming digital image. This model, driven by a partial
More informationA Riemannian Framework for Denoising Diffusion Tensor Images
A Riemannian Framework for Denoising Diffusion Tensor Images Manasi Datar No Institute Given Abstract. Diffusion Tensor Imaging (DTI) is a relatively new imaging modality that has been extensively used
More informationModulation-Feature based Textured Image Segmentation Using Curve Evolution
Modulation-Feature based Textured Image Segmentation Using Curve Evolution Iasonas Kokkinos, Giorgos Evangelopoulos and Petros Maragos Computer Vision, Speech Communication and Signal Processing Group
More informationENERGY METHODS IN IMAGE PROCESSING WITH EDGE ENHANCEMENT
ENERGY METHODS IN IMAGE PROCESSING WITH EDGE ENHANCEMENT PRASHANT ATHAVALE Abstract. Digital images are can be realized as L 2 (R 2 objects. Noise is introduced in a digital image due to various reasons.
More informationMotion Estimation (I) Ce Liu Microsoft Research New England
Motion Estimation (I) Ce Liu celiu@microsoft.com Microsoft Research New England We live in a moving world Perceiving, understanding and predicting motion is an important part of our daily lives Motion
More informationPartial Differential Equations and Image Processing. Eshed Ohn-Bar
Partial Differential Equations and Image Processing Eshed Ohn-Bar OBJECTIVES In this presentation you will 1) Learn what partial differential equations are and where do they arise 2) Learn how to discretize
More informationAnalysis of Numerical Methods for Level Set Based Image Segmentation
Analysis of Numerical Methods for Level Set Based Image Segmentation Björn Scheuermann and Bodo Rosenhahn Institut für Informationsverarbeitung Leibnitz Universität Hannover {scheuerm,rosenhahn}@tnt.uni-hannover.de
More informationScale Space Analysis by Stabilized Inverse Diffusion Equations
Scale Space Analysis by Stabilized Inverse Diffusion Equations Ilya Pollak, Alan S. Willsky and Hamid Krim Laboratory for Information and Decision Systems, Massachusetts Institute of Technology, 77 Massachusetts
More informationA BAYESIAN APPROACH TO NONLINEAR DIFFUSION BASED ON A LAPLACIAN PRIOR FOR IDEAL IMAGE GRADIENT
A BAYESIAN APPROACH TO NONLINEAR DIFFUSION BASED ON A LAPLACIAN PRIOR FOR IDEAL IMAGE GRADIENT Aleksandra Pižurica, Iris Vanhamel, Hichem Sahli, Wilfried Philips and Antonis Katartzis Ghent University,
More informationT H E S I S. Computer Engineering May Smoothing of Matrix-Valued Data. Thomas Brox
D I P L O M A Computer Engineering May 2002 T H E S I S Smoothing of Matrix-Valued Data Thomas Brox Computer Vision, Graphics, and Pattern Recognition Group Department of Mathematics and Computer Science
More informationEdge detection and noise removal by use of a partial differential equation with automatic selection of parameters
Volume 24, N. 1, pp. 131 150, 2005 Copyright 2005 SBMAC ISSN 0101-8205 www.scielo.br/cam Edge detection and noise removal by use of a partial differential equation with automatic selection of parameters
More informationNovel integro-differential equations in image processing and its applications
Novel integro-differential equations in image processing and its applications Prashant Athavale a and Eitan Tadmor b a Institute of Pure and Applied Mathematics, University of California, Los Angeles,
More informationRestoration of images with rotated shapes
Numer Algor (008) 48:49 66 DOI 0.007/s075-008-98-y ORIGINAL PAPER Restoration of images with rotated shapes S. Setzer G. Steidl T. Teuber Received: 6 September 007 / Accepted: 3 January 008 / Published
More informationA Tutorial on Active Contours
Elham Sakhaee March 4, 14 1 Introduction I prepared this technical report as part of my preparation for Computer Vision PhD qualifying exam. Here we discuss derivation of curve evolution function for some
More informationGradient Descent and Implementation Solving the Euler-Lagrange Equations in Practice
1 Lecture Notes, HCI, 4.1.211 Chapter 2 Gradient Descent and Implementation Solving the Euler-Lagrange Equations in Practice Bastian Goldlücke Computer Vision Group Technical University of Munich 2 Bastian
More informationImage Segmentation Using A New Scalar Texture Descriptor based on Extended Structure Tensor
Journal of Information Hiding and Multimedia Signal Processing c 2016 ISSN 2073-4212 Ubiquitous International Volume 7, Number 5, September 2016 Image Segmentation Using A New Scalar Texture Descriptor
More informationKillingFusion: Non-rigid 3D Reconstruction without Correspondences. Supplementary Material
KillingFusion: Non-rigid 3D Reconstruction without Correspondences Supplementary Material Miroslava Slavcheva 1, Maximilian Baust 1 Daniel Cremers 1 Slobodan Ilic 1, {mira.slavcheva,maximilian.baust,cremers}@tum.de,
More informationLEAST-SQUARES FINITE ELEMENT MODELS
LEAST-SQUARES FINITE ELEMENT MODELS General idea of the least-squares formulation applied to an abstract boundary-value problem Works of our group Application to Poisson s equation Application to flows
More informationEnergy Minimization via Graph Cuts
Energy Minimization via Graph Cuts Xiaowei Zhou, June 11, 2010, Journal Club Presentation 1 outline Introduction MAP formulation for vision problems Min-cut and Max-flow Problem Energy Minimization via
More informationLecture No 1 Introduction to Diffusion equations The heat equat
Lecture No 1 Introduction to Diffusion equations The heat equation Columbia University IAS summer program June, 2009 Outline of the lectures We will discuss some basic models of diffusion equations and
More informationMotion Estimation (I)
Motion Estimation (I) Ce Liu celiu@microsoft.com Microsoft Research New England We live in a moving world Perceiving, understanding and predicting motion is an important part of our daily lives Motion
More informationOn Mean Curvature Diusion in Nonlinear Image Filtering. Adel I. El-Fallah and Gary E. Ford. University of California, Davis. Davis, CA
On Mean Curvature Diusion in Nonlinear Image Filtering Adel I. El-Fallah and Gary E. Ford CIPIC, Center for Image Processing and Integrated Computing University of California, Davis Davis, CA 95616 Abstract
More informationFraunhofer Institute for Computer Graphics Research Interactive Graphics Systems Group, TU Darmstadt Fraunhoferstrasse 5, Darmstadt, Germany
Scale Space and PDE methods in image analysis and processing Arjan Kuijper Fraunhofer Institute for Computer Graphics Research Interactive Graphics Systems Group, TU Darmstadt Fraunhoferstrasse 5, 64283
More informationUniversität des Saarlandes. Fachrichtung 6.1 Mathematik
Universität des Saarlandes U N I V E R S I T A S S A R A V I E N I S S Fachrichtung 6.1 Mathematik Preprint Nr. 143 PDEs for Tensor Image Processing Joachim Weickert, Christian Feddern, Martin Welk, Bernhard
More informationTotal Variation Theory and Its Applications
Total Variation Theory and Its Applications 2nd UCC Annual Research Conference, Kingston, Jamaica Peter Ndajah University of the Commonwealth Caribbean, Kingston, Jamaica September 27, 2018 Peter Ndajah
More informationEfficient Beltrami Filtering of Color Images via Vector Extrapolation
Efficient Beltrami Filtering of Color Images via Vector Extrapolation Lorina Dascal, Guy Rosman, and Ron Kimmel Computer Science Department, Technion, Institute of Technology, Haifa 32000, Israel Abstract.
More informationCarotid Lumen Segmentation Based on Tubular Anisotropy and Contours Without Edges Release 0.00
Carotid Lumen Segmentation Based on Tubular Anisotropy and Contours Without Edges Release 0.00 Julien Mille, Fethallah Benmansour and Laurent D. Cohen July 20, 2009 CEREMADE, UMR CNRS 7534, Université
More informationEdge Detection PSY 5018H: Math Models Hum Behavior, Prof. Paul Schrater, Spring 2005
Edge Detection PSY 5018H: Math Models Hum Behavior, Prof. Paul Schrater, Spring 2005 Gradients and edges Points of sharp change in an image are interesting: change in reflectance change in object change
More informationNOISE ENHANCED ANISOTROPIC DIFFUSION FOR SCALAR IMAGE RESTORATION. 6, avenue du Ponceau Cergy-Pontoise, France.
NOISE ENHANCED ANISOTROPIC DIFFUSION FOR SCALAR IMAGE RESTORATION Aymeric HISTACE 1, David ROUSSEAU 2 1 Equipe en Traitement d'image et du Signal UMR CNRS 8051 6, avenue du Ponceau 95014 Cergy-Pontoise,
More informationEdge Detection. Computer Vision P. Schrater Spring 2003
Edge Detection Computer Vision P. Schrater Spring 2003 Simplest Model: (Canny) Edge(x) = a U(x) + n(x) U(x)? x=0 Convolve image with U and find points with high magnitude. Choose value by comparing with
More informationCN780 Final Lecture. Low-Level Vision, Scale-Space, and Polyakov Action. Neil I. Weisenfeld
CN780 Final Lecture Low-Level Vision, Scale-Space, and Polyakov Action Neil I. Weisenfeld Department of Cognitive and Neural Systems Boston University chapter 14.2-14.3 May 9, 2005 p.1/25 Introduction
More informationAn improved active contour model based on level set method
1 015 1 ( ) Journal of East China Normal University (Natural Science) No. 1 Jan. 015 : 1000-5641(015)01-0161-11, (, 0006) :,., (SPF),.,,.,,,,. : ; ; ; : O948 : A DOI: 10.3969/j.issn.1000-5641.015.01.00
More informationPDEs in Image Processing, Tutorials
PDEs in Image Processing, Tutorials Markus Grasmair Vienna, Winter Term 2010 2011 Direct Methods Let X be a topological space and R: X R {+ } some functional. following definitions: The mapping R is lower
More informationOutline Introduction Edge Detection A t c i ti ve C Contours
Edge Detection and Active Contours Elsa Angelini Department TSI, Telecom ParisTech elsa.angelini@telecom-paristech.fr 2008 Outline Introduction Edge Detection Active Contours Introduction The segmentation
More informationA Tensor Variational Formulation of Gradient Energy Total Variation
A Tensor Variational Formulation of Gradient Energy Total Variation Freddie Åström, George Baravdish and Michael Felsberg Linköping University Post Print N.B.: When citing this work, cite the original
More informationA Four-Pixel Scheme for Singular Differential Equations
A Four-Pixel Scheme for Singular Differential Equations Martin Welk 1, Joachim Weickert 1, and Gabriele Steidl 1 Mathematical Image Analysis Group Faculty of Mathematics and Computer Science, Bldg. 7 Saarland
More informationConvex Hodge Decomposition of Image Flows
Convex Hodge Decomposition of Image Flows Jing Yuan 1, Gabriele Steidl 2, Christoph Schnörr 1 1 Image and Pattern Analysis Group, Heidelberg Collaboratory for Image Processing, University of Heidelberg,
More informationUniversität des Saarlandes. Fachrichtung 6.1 Mathematik
Universität des Saarlandes U N I V E R S I T A S S A R A V I E N I S S Fachrichtung 6.1 Mathematik Preprint Nr. 141 Adaptive Structure Tensors and their Applications Thomas Brox, Rein van den Boomgaard,
More informationPDE-based image restoration, I: Anti-staircasing and anti-diffusion
PDE-based image restoration, I: Anti-staircasing and anti-diffusion Kisee Joo and Seongjai Kim May 16, 2003 Abstract This article is concerned with simulation issues arising in the PDE-based image restoration
More informationIMAGE RESTORATION: TOTAL VARIATION, WAVELET FRAMES, AND BEYOND
IMAGE RESTORATION: TOTAL VARIATION, WAVELET FRAMES, AND BEYOND JIAN-FENG CAI, BIN DONG, STANLEY OSHER, AND ZUOWEI SHEN Abstract. The variational techniques (e.g., the total variation based method []) are
More informationNONLOCAL DIFFUSION EQUATIONS
NONLOCAL DIFFUSION EQUATIONS JULIO D. ROSSI (ALICANTE, SPAIN AND BUENOS AIRES, ARGENTINA) jrossi@dm.uba.ar http://mate.dm.uba.ar/ jrossi 2011 Non-local diffusion. The function J. Let J : R N R, nonnegative,
More informationNeighborhood filters and PDE s
Neighborhood filters and PDE s Antoni Buades, Bartomeu Coll, Jean-Michel Morel To cite this version: Antoni Buades, Bartomeu Coll, Jean-Michel Morel. Neighborhood filters and PDE s. Numerische Mathematik,
More informationFinal: Solutions Math 118A, Fall 2013
Final: Solutions Math 118A, Fall 2013 1. [20 pts] For each of the following PDEs for u(x, y), give their order and say if they are nonlinear or linear. If they are linear, say if they are homogeneous or
More informationUniversitat Autònoma de Barcelona. Geometric Differential Operators for Shape Modelling
Universitat Autònoma de Barcelona Geometric Differential Operators for Shape Modelling A dissertation submitted by Debora Gil Resina at Universitat Autònoma de Barcelona to fulfil the degree of Doctor
More informationA New Scheme for Anisotropic Diffusion: Performance Evaluation of Rotation Invariance, Dissipativity and Efficiency no author given no address given A
A New Scheme for Anisotropic Diffusion: Performance Evaluation of Rotation Invariance, Dissipativity and Efficiency Hanno Scharr 1 and Joachim Weickert 1 Interdisciplinary Center for Scientific Computing,
More informationDiffusion of a density in a static fluid
Diffusion of a density in a static fluid u(x, y, z, t), density (M/L 3 ) of a substance (dye). Diffusion: motion of particles from places where the density is higher to places where it is lower, due to
More informationWavelet Inpainting Based on Tensor Diffusion
Vol. 39, No. 7 ACTA AUTOMATICA SINICA July, 03 Wavelet Inpainting Based on Tensor Diffusion YANG Xiu-Hong GUO Bao-Long WU Xian-Xiang Abstract Due to the lossy transmission in the JPEG000 image compression
More informationAlexei F. Cheviakov. University of Saskatchewan, Saskatoon, Canada. INPL seminar June 09, 2011
Direct Method of Construction of Conservation Laws for Nonlinear Differential Equations, its Relation with Noether s Theorem, Applications, and Symbolic Software Alexei F. Cheviakov University of Saskatchewan,
More informationSOLUTIONS OF FOURTH-ORDER PARTIAL DIFFERENTIAL EQUATIONS IN A NOISE REMOVAL MODEL
Electronic Journal of Differential Equations, Vol. 7(7, No., pp.. ISSN: 7-669. URL: ttp://ejde.mat.txstate.edu or ttp://ejde.mat.unt.edu ftp ejde.mat.txstate.edu (login: ftp SOLUTIONS OF FOURTH-ORDER PARTIAL
More informationMathematical Problems in Image Processing
Gilles Aubert Pierre Kornprobst Mathematical Problems in Image Processing Partial Differential Equations and the Calculus of Variations Second Edition Springer Foreword Preface to the Second Edition Preface
More informationScientific Data Computing: Lecture 3
Scientific Data Computing: Lecture 3 Benson Muite benson.muite@ut.ee 23 April 2018 Outline Monday 10-12, Liivi 2-207 Monday 12-14, Liivi 2-205 Topics Introduction, statistical methods and their applications
More informationSpectral Processing. Misha Kazhdan
Spectral Processing Misha Kazhdan [Taubin, 1995] A Signal Processing Approach to Fair Surface Design [Desbrun, et al., 1999] Implicit Fairing of Arbitrary Meshes [Vallet and Levy, 2008] Spectral Geometry
More informationSelection of Optimal Stopping Time for Nonlinear Diffusion Filtering
International Journal of Computer Vision 52(2/3), 189 23, 23 c 23 Kluwer Academic Publishers. Manufactured in The Netherlands. Selection of Optimal Stopping Time for Nonlinear Diffusion Filtering PAVEL
More informationA High Order Conservative Semi-Lagrangian Discontinuous Galerkin Method for Two-Dimensional Transport Simulations
Motivation Numerical methods Numerical tests Conclusions A High Order Conservative Semi-Lagrangian Discontinuous Galerkin Method for Two-Dimensional Transport Simulations Xiaofeng Cai Department of Mathematics
More informationPROBABILISTIC CONTINUOUS EDGE DETECTION USING LOCAL SYMMETRY
PROBABILISTIC CONTINUOUS EDGE DETECTION USING LOCAL SYMMETRY Gerald Mwangi, Paul Fieguth, Christoph S. Garbe : University of Heidelberg, Germany, : University of Waterloo, Canada ABSTRACT We describe a
More informationCONVERGENCE THEORY. G. ALLAIRE CMAP, Ecole Polytechnique. 1. Maximum principle. 2. Oscillating test function. 3. Two-scale convergence
1 CONVERGENCE THEOR G. ALLAIRE CMAP, Ecole Polytechnique 1. Maximum principle 2. Oscillating test function 3. Two-scale convergence 4. Application to homogenization 5. General theory H-convergence) 6.
More informationSINCE THE elegant formulation of anisotropic diffusion
IEEE TRANSACTIONS ON IMAGE PROCESSING, VOL. 7, NO. 3, MARCH 1998 421 Robust Anisotropic Diffusion Michael J. Black, Member, IEEE, Guillermo Sapiro, Member, IEEE, David H. Marimont, Member, IEEE, and David
More informationA Localized Linearized ROF Model for Surface Denoising
1 2 3 4 A Localized Linearized ROF Model for Surface Denoising Shingyu Leung August 7, 2008 5 Abstract 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 1 Introduction CT/MRI scan becomes a very
More informationEdge-preserving wavelet thresholding for image denoising
Journal of Computational and Applied Mathematics ( ) www.elsevier.com/locate/cam Edge-preserving wavelet thresholding for image denoising D. Lazzaro, L.B. Montefusco Departement of Mathematics, University
More informationA Generative Model Based Approach to Motion Segmentation
A Generative Model Based Approach to Motion Segmentation Daniel Cremers 1 and Alan Yuille 2 1 Department of Computer Science University of California at Los Angeles 2 Department of Statistics and Psychology
More informationDigital Signal Processing
JID:YDSPR AID:53 /FLA [m5g; v1.9; Prn:17/07/2017; 11:02] P.1 (1-13) Digital Signal Processing ( ) 1 Contents lists available at ScienceDirect 2 68 3 69 Digital Signal Processing 6 72 7 www.elsevier.com/locate/dsp
More informationEE 367 / CS 448I Computational Imaging and Display Notes: Image Deconvolution (lecture 6)
EE 367 / CS 448I Computational Imaging and Display Notes: Image Deconvolution (lecture 6) Gordon Wetzstein gordon.wetzstein@stanford.edu This document serves as a supplement to the material discussed in
More informationHigh Order Semi-Lagrangian WENO scheme for Vlasov Equations
High Order WENO scheme for Equations Department of Mathematical and Computer Science Colorado School of Mines joint work w/ Andrew Christlieb Supported by AFOSR. Computational Mathematics Seminar, UC Boulder
More informationChapter 3 Second Order Linear Equations
Partial Differential Equations (Math 3303) A Ë@ Õæ Aë áöß @. X. @ 2015-2014 ú GA JË@ É Ë@ Chapter 3 Second Order Linear Equations Second-order partial differential equations for an known function u(x,
More informationNumerical Approximation of Phase Field Models
Numerical Approximation of Phase Field Models Lecture 2: Allen Cahn and Cahn Hilliard Equations with Smooth Potentials Robert Nürnberg Department of Mathematics Imperial College London TUM Summer School
More informationIterative Methods for Ill-Posed Problems
Iterative Methods for Ill-Posed Problems Based on joint work with: Serena Morigi Fiorella Sgallari Andriy Shyshkov Salt Lake City, May, 2007 Outline: Inverse and ill-posed problems Tikhonov regularization
More informationNonlinear stabilization via a linear observability
via a linear observability Kaïs Ammari Department of Mathematics University of Monastir Joint work with Fathia Alabau-Boussouira Collocated feedback stabilization Outline 1 Introduction and main result
More informationInverse problems Total Variation Regularization Mark van Kraaij Casa seminar 23 May 2007 Technische Universiteit Eindh ove n University of Technology
Inverse problems Total Variation Regularization Mark van Kraaij Casa seminar 23 May 27 Introduction Fredholm first kind integral equation of convolution type in one space dimension: g(x) = 1 k(x x )f(x
More informationVariational Methods in Image Denoising
Variational Methods in Image Denoising Jamylle Carter Postdoctoral Fellow Mathematical Sciences Research Institute (MSRI) MSRI Workshop for Women in Mathematics: Introduction to Image Analysis 22 January
More informationMATH 425, HOMEWORK 3 SOLUTIONS
MATH 425, HOMEWORK 3 SOLUTIONS Exercise. (The differentiation property of the heat equation In this exercise, we will use the fact that the derivative of a solution to the heat equation again solves the
More informationMIT (Spring 2014)
18.311 MIT (Spring 014) Rodolfo R. Rosales February 13, 014. Problem Set # 01. Due: Mon. February 4. IMPORTANT: Turn in the regular and the special problems stapled in two SEPARATE packages. Print your
More informationGeneralized Newton-Type Method for Energy Formulations in Image Processing
Generalized Newton-Type Method for Energy Formulations in Image Processing Leah Bar and Guillermo Sapiro Department of Electrical and Computer Engineering University of Minnesota Outline Optimization in
More informationRECENT years have seen a great number of exciting developments
484 IEEE TRANSACTIONS ON SIGNAL PROCESSING, VOL. 53, NO. 2, FEBRUARY 2005 Nonlinear Evolution Equations as Fast and Exact Solvers of Estimation Problems Ilya Pollak, Alan S. Willsky, Fellow, IEEE, and
More informationDense Optical Flow Estimation from the Monogenic Curvature Tensor
Dense Optical Flow Estimation from the Monogenic Curvature Tensor Di Zang 1, Lennart Wietzke 1, Christian Schmaltz 2, and Gerald Sommer 1 1 Department of Computer Science, Christian-Albrechts-University
More informationLucio Demeio Dipartimento di Ingegneria Industriale e delle Scienze Matematiche
Scuola di Dottorato THE WAVE EQUATION Lucio Demeio Dipartimento di Ingegneria Industriale e delle Scienze Matematiche Lucio Demeio - DIISM wave equation 1 / 44 1 The Vibrating String Equation 2 Second
More informationSparse Regularization via Convex Analysis
Sparse Regularization via Convex Analysis Ivan Selesnick Electrical and Computer Engineering Tandon School of Engineering New York University Brooklyn, New York, USA 29 / 66 Convex or non-convex: Which
More informationIntroduction to motion correspondence
Introduction to motion correspondence 1 IPAM - UCLA July 24, 2013 Iasonas Kokkinos Center for Visual Computing Ecole Centrale Paris / INRIA Saclay Why estimate visual motion? 2 Tracking Segmentation Structure
More informationLEVEL SET BASED MULTISPECTRAL SEGMENTATION WITH CORNERS
LEVEL SET BASED MULTISPECTRAL SEGMENTATION WITH CORNERS WENHUA GAO AND ANDREA BERTOZZI Abstract. In this paper we propose an active contour model for segmentation based on the Chan-Vese model. The new
More informationENO and WENO schemes. Further topics and time Integration
ENO and WENO schemes. Further topics and time Integration Tefa Kaisara CASA Seminar 29 November, 2006 Outline 1 Short review ENO/WENO 2 Further topics Subcell resolution Other building blocks 3 Time Integration
More informationIntegro-Differential Equations Based on (BV, L 1 ) Image Decomposition
SIAM J IMAGING SCIENCES Vol 4, No, pp 3 32 c 2 Society for Industrial and Applied Mathematics Integro-Differential Equations Based on (BV, L Image Decomposition Prashant Athavale and Eitan Tadmor Abstract
More informationAllen Cahn Equation in Two Spatial Dimension
Allen Cahn Equation in Two Spatial Dimension Yoichiro Mori April 25, 216 Consider the Allen Cahn equation in two spatial dimension: ɛ u = ɛ2 u + fu) 1) where ɛ > is a small parameter and fu) is of cubic
More informationSeparation of Variables in Linear PDE: One-Dimensional Problems
Separation of Variables in Linear PDE: One-Dimensional Problems Now we apply the theory of Hilbert spaces to linear differential equations with partial derivatives (PDE). We start with a particular example,
More informationOn Multigrid for Phase Field
On Multigrid for Phase Field Carsten Gräser (FU Berlin), Ralf Kornhuber (FU Berlin), Rolf Krause (Uni Bonn), and Vanessa Styles (University of Sussex) Interphase 04 Rome, September, 13-16, 2004 Synopsis
More information