Integro-Differential Equations Based on (BV, L 1 ) Image Decomposition
|
|
- Buddy Moses Watkins
- 5 years ago
- Views:
Transcription
1 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 A novel approach for multiscale image processing based on integro-differential equations (IDEs was proposed in [E Tadmor and P Athavale, Inverse Probl Imaging, 3 (29, pp 693 7] These IDEs, which stem naturally from multiscale (BV,L 2 hierarchical decompositions, yield inverse scale representations of images in the sense that the BV -dual norms of their residuals are inversely proportional to the scaling parameters Motivated by the fact that (BV,L decomposition is more suitable for extracting local scale-space features than (BV,L 2, we introduce here the IDEs which arise from multiscale (BV,L hierarchical decompositions We study several variants of this (BV,L -based IDE model, depending on modifications to the curvature term Key words hierarchical image decompositions, total variation, denoising, deblurring, inverse scale, variational problem, integro-differential equation, energy decomposition AMS subject classifications 26B3, 65C2, 68U, 92C5 DOI 37/79554 Introduction In [32] we introduced a class of integro-differential equations (IDEs to model images, based on (BV,L 2 hierarchical decompositions [33, 34] In this paper we study integro-differential equations arising from hierarchical (BV,L image decomposition As in [32], we work here with grayscale images which can be realized as a graph of a discrete function f :Ω R 2 R The values of this function, f(x, denote the intensity of the image at the discrete points x Ω: the function f attains its maximum value at the brightest spots in the image and its minimum value of zero at the darkest spots The graph of an image consists of discrete pixels, which, for mathematical analysis, is postulated as an L 2 (Ω function The image f, which is obtained by any image capturing device, may be blurred and may contain some additive noise η Thus, the observed image, f, could be written as f = KU+η, whereu is the clean image sought without blurring and noise and K : L 2 (Ω L 2 (Ω is some blurring operator The problem of obtaining the clean image U from the observed image f is an ill-posed problem of image restoration, which can be addressed by several inverse problem solvers We mention in this context variational techniques using Tikhonov-like regularization, PDE-based methods, filtering, stochastic modeling, and wavelet-based techniques that were developed for solving these image processing problems [5, 9, 4, 5, 7, 8, 2, 2, 24, 25, 29, 3, 35] Received by the editors May 7, 2; accepted for publication (in revised form January 3, 2; published electronically March 7, 2 This research was supported in part by NSF grants DMS and DMS Institute of Pure and Applied Mathematics, University of California, Los Angeles, CA 995 (prashant@math uclaedu Department of Mathematics, Institute for Physical Science & Technology and Center of Scientific Computation and Mathematical Modeling (CSCAMM, University of Maryland, College Park, MD 2742 (tadmor@cscammumd edu 3
2 (BV, L INTEGRO-DIFFERENTIAL EQUATIONS FOR IMAGES 3 Image restoration leads to image decomposition For example, any denoising of an observed image f results in the decomposing of type f = U λ + η λ,whereη λ is interpreted as a noise by the denoising method Here, λ is an algorithm-specific scaling parameter: inthecaseof Gaussian smoothing, for example, the variance of the Gaussian kernel may serve as such a scaling parameter Small-scale features, categorized as noise, are then forced into η λ, resulting in a larger-scale version, U λ, of the original image f PDE-based, variational, and hierarchical decompositions Two main approaches for multiscale representations of images are PDE-based methods and variational methods The basic PDE-based method is the heat equation, u t =Δu, where,startingwiththegiven image as initial condition, u( := f, it produces a multiscale representation, {u(t} t>,forthe image f The heat equation removes noise through blurring, but being an isotropic diffusion it also blurs edges which define main features in the image This drawback can be removed by using instead the Perona Malik [29] nonlinear equation based on nonisotropic diffusion, u t =div(g( u u The equation is ill-posed, however, and as already noted by these authors, the Perona Malik diffusion can lead to false detection of edges in the presence of noise To resolve this issue, Catté etal[] proposed using a modified nonisotropic and nonlocal diffusive amplitude of the form g( G σ u ; this class of methods was studied in [9, 36] All these methods give rise to a forward multiscale representation of the image f, in the sense that u(t begins with the finest scale, u( := f, and then progressively blurs into coarser representations of f as t The class of variational methods is a widely used alternative to PDE-based methods in image processing Methods such as Mumford Shah segmentation [25, 26] and Rudin Osher Fatemi (ROF decomposition [3] fall under a general category of Tikhonov regularization [35, 3] Here one attempts to find a close approximation to a function f Y,inaspace X Y, which is an appropriate space adapted to measure edges and textures sought in u This leads to the following minimization problem: { u X + λ v Y } f = u λ + v λ, [u λ,v λ ]:=arginf f=u+v The term u X is a regularizing term, and u λ + v λ is a multiscale decomposition of f which varies with the positive scaling parameter, λ For references on this class of u + v methods we refer the reader to [] In the case of the ROF model [3], for example, edges are sought in the space of bounded variations, X = BV (Ω and f Y = L 2 ; see, eg, [6] This yields the (BV,L 2 -decomposition: { u BV + λ v 2 L 2 }, ( f = u λ + v λ, [u λ,v λ ]:=arginf f=u+v where u BV := Ω u is the BV -seminorm For small values of λ, the minimizer u λ is a large-scale image which consists of only main features and prominent edges in f As λ becomes larger, u λ captures the small-scale details of f The point of view advocated in [33] is to replace the fixed scaling parameter λ by a varying sequence of inverse-scale parameters, λ <λ 2 <, which dictate the scaling of the iterative refinement step, (2a [u λj+,v λj+ ]:= arginf v λj =u+v { u BV + λ j+ v 2 L 2 },
3 32 PRASHANT ATHAVALE AND EITAN TADMOR which in turn generates a multiscale (BV,L 2 hierarchical representation, (2b f j u λj 2 A novel (BV, L -based integro-differential equation (IDE In [32] we introduced the integro-differential equation t (2 u(x, s ds = f(x+ ( u(x, t 2λ(t div u, = u(x, t n Ω This novel IDE was based on a hierarchical (BV,L 2 decomposition (2 which was introduced by Tadmor, Nezzar, and Vese [33, 34], from which it inherits many important properties Once we were dealing with the framework of IDEs, several extensions of (2 wereproposedin[32] to deal with denoising and deblurring of images, which were not strictly associated with any variational problem The ROF decomposition ( usesthel 2 -fidelity term f u 2 L Analternativevariational problem based on an L -fidelity term, f u L, was proposed earlier by Alliney [4] 2 and was studied extensively by Chan, Esedoḡlu, and Nikolova (see [2, 27, 3]: { } (22 f = u λ + v λ, [u λ,v λ ]:=arginf u BV + λ v L f=u+v The resulting (BV,L minimization differs from the (BV,L 2 modelinseveralimportant aspects which have attracted considerable attention in recent years; see [, 2, 3, 6] We shall mention two such aspects The (BV,L minimization is contrast invariant, asopposedto the (BV,L 2 minimization Chan and Esedoḡlu [2] alsoshowedthemorelocal geometric aspects of the (BV,L minimization In particular, recall that if supp(f B R (, then both the (BV,L andthe(bv,l 2 minimizers vanish, namely, u λ andv λ = f for small λ s below the minimal threshold λ<λ L R On the other hand, if the f is a characteristic function χ Σ, then the (BV,L minimizer of (22 admits a maximal threshold, λ H, such that for all λ>λ H we have u λ = f and v λ, in contrast to the leaking phenomenon with the (BV,L 2 -based ROF model, where v λ R [22] 2 Quantization To proceed with the derivation of a (BV,L -based IDE, we first need to address a proper quantization of the image intensity If we let τ denote the small intensity quanta, then we can rescale the coarse image in (22, u λ,inτ-units With this in mind, we decompose the given image f using the (BV,L scheme, starting with initial scale λ : f = u λ τ + v λ, u λ := arginf u { u BV + λ τ f uτ L } The term div ( u(,t, which often arises in image processing, is the curvature of isolevel curves of function u(,t[8, 7] In the present context, it arises as the first variation of the BV -norm, u BV,whichin practice is replaced by its regularized form, eg, u BV = u 2 L + ε2, thus removing the ambiguity when u = and the possible lack of uniqueness in the corresponding (BV,L minimization in (22 [, 3]
4 (BV, L INTEGRO-DIFFERENTIAL EQUATIONS FOR IMAGES 33 The residual image v λ := f u λ τ can be further decomposed into a smaller scale with λ >λ : { v λ = u λ τ + v λ, u λ := arginf u BV + λ } u τ v λ uτ L We can continue this process for λ <λ <λ 2 < : v λk = u λk τ + v λk, u λk := arginf u { u BV + λ } k τ v λ k uτ L A telescoping sum of the refinement step v λk v λk = u λk τ yields the following hierarchical (BV,L decomposition of f: f = u λ τ + v λ = u λ τ + u λ τ + v λ = = u λ τ + u λ τ + + u λn τ + v λn Remark 2 We note that in order to begin with only the main features of f, the initial scale needs to be small but not too small: to capture a nontrivial minimizer, λ needs to be larger than a minimal threshold, which is quantified in terms of the dual -norm specified in (27 below[34, eq (25], [22]; ie, (23 λ >λ L := / sgn(f This yields a hierarchical (BV,L multiscale image decomposition, (24 f = N k=k u λk τ + v λn, where k signals the initial scale λ k to yield the first nontrivial minimizer u λ Consider the Nth step in the (BV,L scheme, ( u λn =arginf u + λ N v λn uτ u Ω τ Ω The Euler Lagrange equation associated with this minimizer reads (with the usual understanding of a regularized curvature term as (25 sgn (u λn τ v λn = ( uλn div λ N u λn From (24 we have N v λn = f u λk τ, k=k
5 34 PRASHANT ATHAVALE AND EITAN TADMOR and using this above expression in (25 we find N sgn u λk τ f k=k = λ N div ( uλn u λn Letting τ, we arrive at the following (BV,L IDE: (26 sgn u(x, s ds f(x = ( u(x, t s=t λ(t div u, =, t t, u(x, t n Ω subject to the initial condition u(,t = The scaling function λ(t is any monotone increasing function at our disposal We discuss the role of this function in the next section 22 On the scaling function λ(t It is argued in [22] that the dual norm, w, ϕ ϕ (27 w := sup L 2, ϕ = ϕ(xdx, ϕ BV ϕ BV Ω Ω is a proper norm to measure texture (equivalently, one can characterize this dual norm as w := inf g { g L :divg = w w} To understand the critical role of the scaling function λ(t intheidemodel(26 and its relationship with the star-norm, we first prove the following useful lemma Lemma 2 For u BV,letκ u be the curvature term, κ u := div ( u, whenever u Then κ u = Proof For ϕ BV we have the following: ( (28 (κ u,ϕ L 2 = div ( u u,ϕ L 2 ϕ BV Thus, we have κ u Letting ϕ = u in (28, we obtain κ u,u L 2 = u BV and the result follows The following theorem is a direct consequence of this lemma Theorem 22 Consider the IDE model (26 sgn (U(x, t f(x = λ(t div ( u(x, t u(x, t t, U(x, t := u(x, s ds t Then, the size of the residual, V (,t:=f t t u(,s ds, is dictated by the scaling function λ(t: sgn(v (,t = λ(t The star-norm measures oscillations [22] Thus, the above equality could be interpreted as saying that the oscillations in the residual image V (,t are inversely proportional to the scaling function λ(t In particular, the oscillations in the residual function diminish in a precise manner controlled by increasing λ(t
6 (BV, L INTEGRO-DIFFERENTIAL EQUATIONS FOR IMAGES 35 Remark 22 Recall that t in Theorem 22 is determined as the minimal effective scale in the sense that (23 holds; ie, λ(t >λ L =/ sgn(f One does not need, however, to have a priori knowledge on the amount of star-norm oscillation in f, sinceu(,t as long as λ(t is underestimated, ie, λ<λ L Eventually, the increasing λ(t will pick up the minimal effective scale 23 L -energy decomposition The increase of λ(t is directly linked to the energy of the underlying image For the (BV,L 2 -based IDE (2, we proved in [32] thel 2 -energy decomposition t s=t λ(s u(,s BV ds + V (,t 2 L = f 2 2 L 2 For the (BV,L -based IDE we have the corresponding L -energy decomposition Proposition 23 Consider the IDE (26 and let V (,t denote the residual, V (,t := t t u(,s ds f Then we have the following L -energy decomposition: t t λ(s u(,s BV ds + V (,t L = f L Proof We denote by U(,t:= t t u(,s ds the cumulative image so that the IDE (26 takes the form ( u(x, t λ(t div = sgn(u(x, t f(x u(x, t Integrating the above against u d dt (U f, we get λ(t u(,t BV = sgn(u(,t f d U(,t f dx Ω dt( = d sgn(u(,t f(u(,t f dx = d U(,t f dx, dt dt Ω and the result follows by time integration 3 IDEs for images: Numerical results and extensions 3 Comparing the (BV, L -and(bv, L 2 -based IDEs The multiscale image representation generated by the IDE (26 is shown in Figure It is instructive to look at the image representation produced by the (BV,L -based IDE (26 withthe(bv,l 2 -based IDE (2 In the first row of Figure 2 we see the multiscale representation using the new (26 Note that the larger blocks appear before the small blocks, irrespective of their intensity levels We observe that the appearance of these blocks is abrupt, as opposed to the results obtained with (2 shown in the second row of Figure 2, where these blocks appear gradually and the multiscale representation produced by (2 depends on the intensity levels rather than the size of the blocks This feature makes (BV,L IDE(26 more suitable for multiscale representation than (BV,L 2 IDE(2 The difference between the (BV,L -based IDE and the (BV,L 2 -based IDE is demonstrated clearly in Figure 3 We note that, for the same time, the (BV,L IDE gives better Ω
7 36 PRASHANT ATHAVALE AND EITAN TADMOR t = t = t =3 t =25 Figure The images U(t = t u(,s ds of the IDE (26 at t =,, 3, 25 Here, λ(t = 2t t = t =3 t =5 t =7 t = t =4 t =6 t = Figure 2 First row: the images U(t = t u(,s ds of the (BV,L -based IDE (26 at t =, 3, 5, 7 Here, λ(t =5 2 t Second row: the images U(t = t u(,s ds of the (BV,L2 -based IDE (2 at t =, 4, 6, Here, λ(t = 2 t results at the edges than the (BV,L 2 IDE The(BV,L IDE is seen to be contrast invariant and depends only on the scale of the blocks, whereas the (BV,L 2 IDEisnotcontrast invariant, as expected 32 IDE with filtered diffusion Recall that one of the drawbacks in using the heat equation for denoising is that it results in an isotropic diffusion The Perona Malik model removes this drawback by introducing a diffusion controlling function that controls the diffusion near prominent edges in a given image We propose a similar modification to our IDE model, seeking u(x, t :Ω R + R such that (3a sgn u(x, s dx f(x = g( G ( σ u(x, t u(x, t div, t > t, s=t λ(t u(x, t
8 (BV, L INTEGRO-DIFFERENTIAL EQUATIONS FOR IMAGES 37 (a (b (c Figure 3 In (a, (b, and (c we see the mesh views of the original image from Figure 2, the result U(t = t u(,s ds for the (BV,L2 -based IDE (2 at t =, and the image U(t = t u(,s ds for the (BV,L -based IDE (26, alsoatt =, respectively subject to normal boundary conditions u n Ω = As before, the IDE (3a is activated after the minimal time scale, t t, such that u(,t is the first slice of f to capture nontrivial features at that scale Similar to the Perona Malik-type models, we can choose the prefactor function g so that it vanishes at infinity to control the diffusion at prominent edges in the image Thus, the function g acts here as a high-pass filter which retains prominent edges in the image t u(x, s ds without diffusing them As choices for such a g-filter, Figure 5 displays the results of the modified IDE (3a with (3b g(s = +(s/β 2 Here, the constant β determines the extent to which edges are preserved: for small β s, relevant edges are preserved, whereas for large β s, they are diffused Detailed discussion of the numerical scheme for the filtered diffusion model (3 is given in section 4 Comparing the results of the filtered IDE (3a showninfigure5, we observe that edges, which are diffused by the basic IDE (26 as depicted in Figure 4, are preserved in Figure 5 As inthecaseof(2 (see[32], this phenomenon is more apparent for smaller values of t due to the fact that as t increases, U(,t in both models approaches f and, consequently, suffers from less diffusion of the edges The usefulness of the filtered diffusion IDE model becomes apparent when certain edges are required in the scale space for smaller values of t For example, in Figure 4, the edges are blurred for smaller values of t with the standard IDE (26, but with the filtered diffusion IDE (3a we retain relevant edges, as shown in Figure 5 We remark that this edge preserving filtered diffusion could be used in landmark-based hierarchical image registration algorithms [23, 28] 33 IDE with tangential smoothing The approach of using the diffusion controlling function works well with natural images with moderate gradients With other images, however, such as those which often arise in computer vision and industrial applications, the boundaries of their internal objects are marked with large, sharp gradients, for example, characteristic function χ D,whereD Ω In such cases, we can choose to smooth only in the tangential
9 38 PRASHANT ATHAVALE AND EITAN TADMOR t =4 t =42 t =44 t =6 Figure 4 The images U(t = t u(,s ds of the IDE (26 at t =4, 42, 44, 6 Here, λ(t = 2t t =2 t =38 t =4 t =6 Figure 5 The images U(t = t u(,s ds of the IDE (3a at t =2, 38, 4, 6 Here, λ(t = 2t direction to the boundaries of the objects, eg, [5] To this end, write Δu := u TT + u NN, where u TT and u NN are the tangential and normal diffusion components, ie, ( u u u TT =Δu u NN = u div, u NN = u u, 2 u u u If we restrict the diffusion in our IDE model to tangential directions, we obtain modified IDEs with tangential smoothing, (32 sgn u(x, s ds f(x = ( u(x, t u u(x, t div, =, t λ(t u(x, t n Ω and with tangential smoothing and filtering, (33 sgn u(x, s ds f(x t = g( G ( σ u(x, t u(x, t u(x, t div λ(t u(x, t As before, u :Ω R + R evolves in inverse scale space using the conditions u(, and u n Ω = Numerical experiments are shown in Figures 6 and 7 Comparethe(BV,L IDE results (26 shown in Figure 6 with the tangential smoothing (33 shown in Figure 7: the point here is that the tangential diffusion model preserves the edges, while denoising the rest of the image at a much faster rate than in the standard IDE model
10 (BV, L INTEGRO-DIFFERENTIAL EQUATIONS FOR IMAGES 39 f at t = t = t =4 t =8 Figure 6 A given noisy image f and the IDE images t u(,s ds of (26 at t =, 4, 8 Here, the scaling function is λ(t = 2 t Most of the noise is present at scale t =8 f at t = t = t =4 t =8 Figure 7 A given noisy image f and the IDE images t u(,s ds of (33 at t =, 4, 8 Here, the scaling function is λ(t = 2 t We observe that at t =8the noise is suppressed and at the same time the boundary of the circle is kept intact 4 Appendix: Numerical discretizations In this appendix we describe the numerical implementation of (26 and its variants First let us concentrate on the basic (BV,L IDE model (26, rewritten here for convenience (for simplicity we set the initial scale at λ(t at t =, (4 sgn u(x, s ds f(x = ( u(x, t λ(t div u(x, t In practice we approximate the BV -seminorm u BV with Ω u 2 + ɛ 2, which leads to the following IDE: ( (42 sgn u(x, s ds f(x = λ(t div u(x, t ɛ 2 + u(x, t 2 In our numerical experiments we used the regularization parameter ɛ = 3 Then, as usual, U(t := t u(x, s ds is the exact solution Let Δt bethetimestep,andletun+ denote the corresponding computed solution at t n+ =(n +Δt: U n+ = U n + W n+, W n+ Wi,j n+ := u n+ i,j Δt, where u n+ i,j u n+ (ih, jh is the approximate solution of the IDE at grid point (ih, jh We introduce α n := U n f, and we compute the increment W n+ by fixed-point iterations {ω k }
11 3 PRASHANT ATHAVALE AND EITAN TADMOR With this, the IDE (42 is discretized at t = t n+ : (43 U n i,j + ω k+ i,j + + = f i,j ωi+,j k ωk+ i,j ωi,j k+ ωi,j k ε 2 +(D +x ω ki,j 2 +(D y ω ki,j 2 ε 2 +(D x ωi,j k 2 +(D y ωi,j k 2 ωi,j+ k ωk+ i,j ωi,j k+ ωi,j k ε 2 +(D x ω ki,j 2 +(D +y ω ki,j 2 ε 2 +(D x ωi,j k 2 +(D y ωi,j k 2 α n λ (n+ h 2 α n λ (n+ h 2 The nonlinear system (43 is solved using Jacobi iterations, which leads to the fixed-point iterations for computing ω k+ : (44a i,j = λ(n+ h 2 (f i,j Ui,j n +αn (c E ωi+,j k + c W ωi,j k + c Sωi,j+ k + c N ωi,j k λ (n+ h 2 + α n (c E + c W + c S + c N ω k+ Here, λ (n+ = λ(t n+ are the discrete scaling parameters, and c E,c W,c S,c N are the discrete coefficients c E := c S := ε 2 +(D +x ω ki,j 2 +(D y ω ki,j 2, c W := ε 2 +(D x ω ki,j 2 +(D +y ω ki,j 2, c N := ε 2 +(D x ω ki,j 2 +(D y ω ki,j 2, ε 2 +(D x ω ki,j 2 +(D y ω ki,j 2 In the computations above we set h = To minimize the grid effects, we alternate the directions in which the above iterations are carried out, starting at the top-left corner position (, Fixing i =,wevaryj =toj max (east-south direction, initiating the next iteration at the top-right corner, and so on The fixed-point iterations (44a yield ω k k W n+ u n+ Δt, and we can then update the computed image U, (44b U n+ = U n + W n+ Next, we consider the filtered IDE (3a, which is rewritten here for convenience as u(x, s ds f(x = g( G ( σ u(x, t u(x, t div λ(t u(x, t The only difference here is the additional diffusion controlling function g( G σ u(x, t, where G σ is the two-dimensional Gaussian smoothing with standard deviation σ The function g(s = with β = 5 is used in our numerical experiments We approximate +(s/β 2 ( Gσ g( G σ u(x, t g ωn i,j, Δt and the expression on the right enters into the right-hand side of (43 Weendupwiththe same discrete IDE scheme (44 withλ (n λ (n/ g( G σ ω n i,j /Δt
12 (BV, L INTEGRO-DIFFERENTIAL EQUATIONS FOR IMAGES 3 REFERENCES [] W K Allard, Total variation regularization for image denoising, I Geometric theory, SIAMJ Math Anal, 39 (27, pp 5 9 [2] W K Allard, Total variation regularization for image denoising, II Examples, SIAMJ ImagingSci, (28, pp 4 47 [3] W K Allard, Total variation regularization for image denoising, III Examples, SIAM J Imaging Sci, 2 (29, pp [4] S Alliney, Digital filters as absolute norm regularizers, IEEE Trans Signal Process, 4 (992, pp [5] L Alvarez, P-L Lions, and J-M Morel, Image selective smoothing and edge detection by nonlinear diffusion II, SIAM J Numer Anal, 29 (992, pp [6] L Ambrosio, N Fusco, and D Pallara, Functions of Bounded Variation and Free Discontinuity Problems, Oxford Math Monogr, The Clarendon Press, Oxford University Press, New York, 2 [7] F Andreu-Vaillo, V Casalles, and J Mazón, Parabolic Quasilinear Equations Minimizing Linear Growth Functionals, Progr Math 223, Birkhäuser, Basel, Switzerland, 24 [8] G Aubert and P Kornprobst, Mathematical Problems in Image Processing Partial Differential Equations and the Calculus of Variations, 2nd ed, Appl Math Sci 47, Springer-Verlag, New York, 26 [9] G Aubert and L Vese, A variational method in image recovery, SIAM J Numer Anal, 34 (997, pp [] A Buades, T Le, J-M Morel, and L Vese, Fast cartoon + texture image filters, IEEE Trans Image Process, 9 (2, pp [] F Catté, P-L Lions, J-M Morel, and T Coll, Image selective smoothing and edge detection by nonlinear diffusion, SIAM J Numer Anal, 29 (992, pp [2] T F Chan and S Esedoḡlu, Aspects of total variation regularized L function approximation, SIAM J Appl Math, 65 (25, pp [3] T F Chan, S Esedoḡlu, and M Nikolova, Algorithms for finding global minimizers of image segmentation and denoising models, SIAM J Appl Math, 66 (26, pp [4] T F Chan and J Shen, Image Processing and Analysis: Variational, PDE, Wavelet, and Stochastic Methods, SIAM, Philadelphia, 25 [5] T F Chan, J Shen, and L Vese, Variational PDE models in image processing, Notices Amer Math Soc, 5 (23, pp 4 26 [6] V Duval, J-F Aujol, and Y Gousseau, The TVL model: A geometric point of view, Multiscale Model Simul, 8 (29, pp [7] S Esedoḡlu and S J Osher, Decomposition of images by the anisotropic Rudin-Osher-Fatemi model, Comm Pure Appl Math, 57 (24, pp [8] S Geman and D Geman, Stochastic relaxation, Gibbs distributions, and the Bayesian restoration of images, IEEE Trans Pattern Anal Mach Intell, 6 (984, pp [9] P Guidotti, A new well-posed nonlinear nonlocal diffusion, Nonlinear Anal, 72 (2, pp [2] S Li, Markov Random Field Modeling in Computer Vision, Springer-Verlag, London, 995 [2] S Mallat, A Wavelet Tour of Signal Processing, Academic Press, San Diego, CA, 998 [22] Y Meyer, Oscillating Patterns in Image Processing and Nonlinear Evolution Equations The Fifteenth Dean Jacqueline B Lewis Memorial Lectures, Univ Lecture Ser 22, American Mathematical Society, Providence, RI, 2 [23] J Modersitzki, Numerical Methods for Image Registration, Oxford University Press, New York, 24 [24] V Morozov, Regularization Methods for Ill-Posed Problems, CRC Press, Boca Raton, FL, 993 Translated from the 987 Russian original [25] D Mumford and J Shah, Boundary detection by minimizing functionals, in Proceedings of the IEEE Computer Society Conference on Computer Vision and Pattern Recognition, San Francisco, CA, 985, pp [26] D Mumford and J Shah, Optimal approximations by piecewise smooth functions and associated variational problems, Comm Pure Appl Math, 42 (989, pp [27] M Nikolova, Minimizers of cost-functions involving nonsmooth data-fidelity terms Application to the processing of outliers, SIAM J Numer Anal, 4 (22, pp
13 32 PRASHANT ATHAVALE AND EITAN TADMOR [28] D Paquin, D Levy, E Schreibmann, and L Xing, Multiscale image registration, Math Biosci Eng, 3 (26, pp [29] P Perona and J Malik, Scale-space and edge detection using anisotropic diffusion, IEEE Trans Pattern Anal Mach Intell, 2 (99, pp [3] L Rudin, S Osher, and E Fatemi, Nonlinear total variation based noise removal algorithms, Phys D, 6 (992, pp [3] O Scherzer, M Grasmair, H Grossauer, M Haltmeier, and F Lenzen, Variational Methods in Imaging, Appl Math Sci 67, Springer-Verlag, New York, 29 [32] E Tadmor and P Athavale, Multiscale image representation using novel integro-differential equations, Inverse Probl Imaging, 3 (29, pp [33] E Tadmor, S Nezzar, and L Vese, A multiscale image representation using hierarchical (BV,L 2 decompositions, Multiscale Model Simul, 2 (24, pp [34] E Tadmor, S Nezzar, and L Vese, Multiscale hierarchical decomposition of images with applications to deblurring, denoising and segmentation, Commun Math Sci, 6 (28, pp 26 [35] A Tikhonov and V Arsenin, Solutions of Ill-Posed Problems, Scripta Series in Mathematics, V H Winston & Sons, Washington, DC, John Wiley & Sons, New York, 977 Translated from the Russian [36] J Weickert, Anisotropic Diffusion in Image Processing, Teubner-Verlag, Stuttgart, Germany, 998
Novel 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 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 informationOn a multiscale representation of images as hierarchy of edges. Eitan Tadmor. University of Maryland
On a multiscale representation of images as hierarchy of edges Eitan Tadmor Center for Scientific Computation and Mathematical Modeling (CSCAMM) Department of Mathematics and Institute for Physical Science
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 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 informationImage Cartoon-Texture Decomposition and Feature Selection using the Total Variation Regularized L 1 Functional
Image Cartoon-Texture Decomposition and Feature Selection using the Total Variation Regularized L 1 Functional Wotao Yin 1, Donald Goldfarb 1, and Stanley Osher 2 1 Department of Industrial Engineering
More informationAn image decomposition model using the total variation and the infinity Laplacian
An image decomposition model using the total variation and the inity Laplacian Christopher Elion a and Luminita A. Vese a a Department of Mathematics, University of California Los Angeles, 405 Hilgard
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 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 informationFundamentals of Non-local Total Variation Spectral Theory
Fundamentals of Non-local Total Variation Spectral Theory Jean-François Aujol 1,2, Guy Gilboa 3, Nicolas Papadakis 1,2 1 Univ. Bordeaux, IMB, UMR 5251, F-33400 Talence, France 2 CNRS, IMB, UMR 5251, F-33400
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 informationConvex Hodge Decomposition and Regularization of Image Flows
Convex Hodge Decomposition and Regularization of Image Flows Jing Yuan, Christoph Schnörr, Gabriele Steidl April 14, 2008 Abstract The total variation (TV) measure is a key concept in the field of variational
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 informationA Variational Approach to Reconstructing Images Corrupted by Poisson Noise
J Math Imaging Vis c 27 Springer Science + Business Media, LLC. Manufactured in The Netherlands. DOI: 1.7/s1851-7-652-y A Variational Approach to Reconstructing Images Corrupted by Poisson Noise TRIET
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 informationNonlinear Diffusion. Journal Club Presentation. Xiaowei Zhou
1 / 41 Journal Club Presentation Xiaowei Zhou Department of Electronic and Computer Engineering The Hong Kong University of Science and Technology 2009-12-11 2 / 41 Outline 1 Motivation Diffusion process
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 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 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 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 informationNonlinear Flows for Displacement Correction and Applications in Tomography
Nonlinear Flows for Displacement Correction and Applications in Tomography Guozhi Dong 1 and Otmar Scherzer 1,2 1 Computational Science Center, University of Vienna, Oskar-Morgenstern-Platz 1, 1090 Wien,
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 informationLecture 4 Colorization and Segmentation
Lecture 4 Colorization and Segmentation Summer School Mathematics in Imaging Science University of Bologna, Itay June 1st 2018 Friday 11:15-13:15 Sung Ha Kang School of Mathematics Georgia Institute of
More informationNonlinear regularizations of TV based PDEs for image processing
Nonlinear regularizations of TV based PDEs for image processing Andrea Bertozzi, John Greer, Stanley Osher, and Kevin Vixie ABSTRACT. We consider both second order and fourth order TV-based PDEs for image
More informationc 2006 Society for Industrial and Applied Mathematics
SIAM J. APPL. MATH. Vol. 66, No. 5, pp. 1632 1648 c 26 Society for Industrial and Applied Mathematics ALGORITHMS FOR FINING GLOBAL MINIMIZERS OF IMAGE SEGMENTATION AN ENOISING MOELS TONY F. CHAN, SELIM
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 informationADAPTIVE DIFFUSION CONSTRAINED TOTAL VARIATION SCHEME WITH APPLICATION TO CARTOON + TEXTURE + EDGE IMAGE DECOMPOSITION
Pré-Publicações do Departamento de Matemática Universidade de Coimbra Preprint Number 13 54 ADAPTIVE DIFFUSION CONSTRAINED TOTAL VARIATION SCHEME WITH APPLICATION TO CARTOON + TEXTURE + EDGE IMAGE DECOMPOSITION
More informationIMAGE RECOVERY USING FUNCTIONS OF BOUNDED VARIATION AND SOBOLEV SPACES OF NEGATIVE DIFFERENTIABILITY. Yunho Kim and Luminita A.
Inverse Problems and Imaging Volume 3, No. 1, 2009, 43 68 doi:10.3934/ipi.2009.3.43 IMAGE RECOVERY USING FUNCTIONS OF BOUNDED VARIATION AND SOBOLEV SPACES OF NEGATIVE DIFFERENTIABILITY Yunho Kim and Luminita
More informationAn Energy Regularization Method for the Backward Diffusion Problem and its Applications to Image Deblurring
COMMUNICATIONS IN COMPUTATIONAL PHYSICS Vol. 4, No. 1, pp. 177-194 Commun. Comput. Phys. July 8 An Energy Regularization Method for the Backward Diffusion Problem and its Applications to Image Deblurring
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. 402 Iterative TV-regularization of Grey-scale Images Martin Fuchs and Joachim Weickert Saarbrücken
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 informationKey words: denoising, higher-order regularization, stability, weak convergence, Brezis-Lieb condition. AMS subject classifications: 49N45, 94A08.
A FEW REMARKS ON VARIATIONAL MODELS FOR DENOISING RUSTUM CHOKSI, IRENE FONSECA, AND BARBARA ZWICKNAGL Abstract. Variational models for image and signal denoising are based on the minimization of energy
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 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 informationc 2008 International Press
COMMUN. MATH. SCI. Vol. 6, No. 2, pp. 281 307. c 2008 International Press MULTISCALE HIERARCHICAL DECOMPOSITION OF IMAGES WITH APPLICATIONS TO DEBLURRING, DENOISING AND SEGMENTATION EITAN TADMOR, SUZANNE
More informationIMA Preprint Series # 2098
LEVEL SET BASED BIMODAL SEGMENTATION WITH STATIONARY GLOBAL MINIMUM By Suk-Ho Lee and Jin Keun Seo IMA Preprint Series # 9 ( February ) INSTITUTE FOR MATHEMATICS AND ITS APPLICATIONS UNIVERSITY OF MINNESOTA
More informationA fast primal-dual method for generalized Total Variation denoising
Appl. Math. Inf. Sci. 6, No. 3, 401-409 (2012) 401 Applied Mathematics & Information Sciences An International Journal A fast primal-dual method for generalized Total Variation denoising Francesc Aràndiga,
More informationVariational Methods in Signal and Image Processing
Variational Methods in Signal and Image Processing XU WANG Texas A&M University Dept. of Electrical & Computer Eng. College Station, Texas United States xu.wang@tamu.edu ERCHIN SERPEDIN Texas A&M University
More informationTWO-PHASE APPROACH FOR DEBLURRING IMAGES CORRUPTED BY IMPULSE PLUS GAUSSIAN NOISE. Jian-Feng Cai. Raymond H. Chan. Mila Nikolova
Manuscript submitted to Website: http://aimsciences.org AIMS Journals Volume 00, Number 0, Xxxx XXXX pp. 000 000 TWO-PHASE APPROACH FOR DEBLURRING IMAGES CORRUPTED BY IMPULSE PLUS GAUSSIAN NOISE Jian-Feng
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 informationImage Deblurring in the Presence of Impulsive Noise
Image Deblurring in the Presence of Impulsive Noise Leah Bar Nir Sochen Nahum Kiryati School of Electrical Engineering Dept. of Applied Mathematics Tel Aviv University Tel Aviv, 69978, Israel April 18,
More informationError Analysis for H 1 Based Wavelet Interpolations
Error Analysis for H 1 Based Wavelet Interpolations Tony F. Chan Hao-Min Zhou Tie Zhou Abstract We rigorously study the error bound for the H 1 wavelet interpolation problem, which aims to recover missing
More informationVariational Image Restoration
Variational Image Restoration Yuling Jiao yljiaostatistics@znufe.edu.cn School of and Statistics and Mathematics ZNUFE Dec 30, 2014 Outline 1 1 Classical Variational Restoration Models and Algorithms 1.1
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 informationConvergence rates of convex variational regularization
INSTITUTE OF PHYSICS PUBLISHING Inverse Problems 20 (2004) 1411 1421 INVERSE PROBLEMS PII: S0266-5611(04)77358-X Convergence rates of convex variational regularization Martin Burger 1 and Stanley Osher
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 informationNON-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 informationANALYSIS OF THE TV REGULARIZATION AND H 1 FIDELITY MODEL FOR DECOMPOSING AN IMAGE INTO CARTOON PLUS TEXTURE. C.M. Elliott and S.A.
COMMUNICATIONS ON Website: http://aimsciences.org PURE AND APPLIED ANALYSIS Volume 6, Number 4, December 27 pp. 917 936 ANALYSIS OF THE TV REGULARIZATION AND H 1 FIDELITY MODEL FOR DECOMPOSING AN IMAGE
More informationITERATIVELY SOLVING LINEAR INVERSE PROBLEMS UNDER GENERAL CONVEX CONSTRAINTS. Ingrid Daubechies. Gerd Teschke. Luminita Vese
Inverse Problems and Imaging Volume 1, No. 1, 2007, 29 46 Web site: http://www.aimsciences.org ITERATIVELY SOLVING LINEAR INVERSE PROBLEMS UNDER GENERAL CONVEX CONSTRAINTS Ingrid Daubechies Princeton University,
More informationA Convex Relaxation of the Ambrosio Tortorelli Elliptic Functionals for the Mumford Shah Functional
A Convex Relaxation of the Ambrosio Tortorelli Elliptic Functionals for the Mumford Shah Functional Youngwook Kee and Junmo Kim KAIST, South Korea Abstract In this paper, we revisit the phase-field approximation
More informationImage Deblurring in the Presence of Impulsive Noise
Image Deblurring in the Presence of Impulsive Noise Leah Bar Nir Sochen Nahum Kiryati School of Electrical Engineering Dept. of Applied Mathematics Tel Aviv University Tel Aviv, 69978, Israel Abstract
More informationA FAMILY OF NONLINEAR DIFFUSIONS CONNECTING PERONA-MALIK TO STANDARD DIFFUSION. Patrick Guidotti
DISCRETE AND CONTINUOUS DYNAMICAL SYSTEMS SERIES S Volume 5, Number 3, June 202 doi:0.3934/dcdss.202.5.xx pp. X XX A FAMILY OF NONLINEAR DIFFUSIONS CONNECTING PERONA-MALIK TO STANDARD DIFFUSION Patrick
More informationImage enhancement by non-local reverse heat equation
Image enhancement by non-local reverse heat equation Antoni Buades 1, Bartomeu Coll 2, and Jean-Michel Morel 3 1 Dpt Matematiques i Informatica Univ. Illes Balears Ctra Valldemossa km 7.5, 07122 Palma
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 informationNonlinear diffusion filtering on extended neighborhood
Applied Numerical Mathematics 5 005) 1 11 www.elsevier.com/locate/apnum Nonlinear diffusion filtering on extended neighborhood Danny Barash Genome Diversity Center, Institute of Evolution, University of
More informationIteratively Solving Linear Inverse Problems under General Convex Constraints
Konrad-Zuse-Zentrum fu r Informationstechnik Berlin Takustraße 7 D-14195 Berlin-Dahlem Germany I NGRID DAUBECHIES, G ERD T ESCHKE, L UMINITA V ESE Iteratively Solving Linear Inverse Problems under General
More informationELLIPTIC RECONSTRUCTION AND A POSTERIORI ERROR ESTIMATES FOR PARABOLIC PROBLEMS
ELLIPTIC RECONSTRUCTION AND A POSTERIORI ERROR ESTIMATES FOR PARABOLIC PROBLEMS CHARALAMBOS MAKRIDAKIS AND RICARDO H. NOCHETTO Abstract. It is known that the energy technique for a posteriori error analysis
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 informationStability Properties of Perona-Malik Scheme
Stability Properties of Perona-Malik Scheme Selim Esedoglu Institute for Mathematics and its Applications December 1 Abstract The Perona-Malik scheme is a numerical technique for de-noising digital images
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 informationModeling oscillatory components with the homogeneous spaces BMO
Modeling oscillatory components with the homogeneous spaces BMO α and Ẇ α,p John B. Garnett, Peter W. Jones, Triet M. Le, and Luminita A. Vese UCLA CAM Report 07-2, 2007 Abstract This paper is devoted
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 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 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 informationregularization parameter choice
A duality-based splitting method for l 1 -TV image restoration with automatic regularization parameter choice Christian Clason Bangti Jin Karl Kunisch November 30, 2009 A novel splitting method is presented
More informationDue Giorni di Algebra Lineare Numerica (2GALN) Febbraio 2016, Como. Iterative regularization in variable exponent Lebesgue spaces
Due Giorni di Algebra Lineare Numerica (2GALN) 16 17 Febbraio 2016, Como Iterative regularization in variable exponent Lebesgue spaces Claudio Estatico 1 Joint work with: Brigida Bonino 1, Fabio Di Benedetto
More informationA posteriori error control for the binary Mumford Shah model
A posteriori error control for the binary Mumford Shah model Benjamin Berkels 1, Alexander Effland 2, Martin Rumpf 2 1 AICES Graduate School, RWTH Aachen University, Germany 2 Institute for Numerical Simulation,
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 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 informationThe heat equation in time dependent domains with Neumann boundary conditions
The heat equation in time dependent domains with Neumann boundary conditions Chris Burdzy Zhen-Qing Chen John Sylvester Abstract We study the heat equation in domains in R n with insulated fast moving
More informationOptimal stopping time formulation of adaptive image filtering
Optimal stopping time formulation of adaptive image filtering I. Capuzzo Dolcetta, R. Ferretti 19.04.2000 Abstract This paper presents an approach to image filtering based on an optimal stopping time problem
More informationA primal-dual approach for a total variation Wasserstein flow
A primal-dual approach for a total variation Wasserstein flow Martin Benning 1, Luca Calatroni 2, Bertram Düring 3, Carola-Bibiane Schönlieb 4 1 Magnetic Resonance Research Centre, University of Cambridge,
More informationControl of Interface Evolution in Multi-Phase Fluid Flows
Control of Interface Evolution in Multi-Phase Fluid Flows Markus Klein Department of Mathematics University of Tübingen Workshop on Numerical Methods for Optimal Control and Inverse Problems Garching,
More informationMathematical modeling of textures : application to color image decomposition with a projected gradient algorithm
Mathematical modeling of textures : application to color image decomposition with a projected gradient algorithm Vincent Duval, Jean-François Aujol, Luminita Vese To cite this version: Vincent Duval, Jean-François
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 informationAn unconstrained multiphase thresholding approach for image segmentation
An unconstrained multiphase thresholding approach for image segmentation Benjamin Berkels Institut für Numerische Simulation, Rheinische Friedrich-Wilhelms-Universität Bonn, Nussallee 15, 53115 Bonn, Germany
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 informationParameter Identification in Partial Differential Equations
Parameter Identification in Partial Differential Equations Differentiation of data Not strictly a parameter identification problem, but good motivation. Appears often as a subproblem. Given noisy observation
More informationA Theoretical Framework for the Regularization of Poisson Likelihood Estimation Problems
c de Gruyter 2007 J. Inv. Ill-Posed Problems 15 (2007), 12 8 DOI 10.1515 / JIP.2007.002 A Theoretical Framework for the Regularization of Poisson Likelihood Estimation Problems Johnathan M. Bardsley Communicated
More informationREFINED UPPER BOUNDS ON THE COARSENING RATE OF DISCRETE, ILL-POSED DIFFUSION EQUATIONS
REFINED UPPER BOUNDS ON THE COARSENING RATE OF DISCRETE, ILL-POSED DIFFUSION EQUATIONS SELIM ESEDOḠLU AND DEJAN SLEPČEV Abstract We study coarsening phenomena observed in discrete, ill-posed diffusion
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 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 informationIN this chapter, we introduce the Stabilized Inverse Diffusion Equations (SIDEs), as
Chapter Image Segmentation with Stabilized Inverse Diffusion Equations. Introduction. IN this chapter, we introduce the Stabilized Inverse Diffusion Equations (SIDEs), as well as illustrate their speed
More informationSemi-Blind Image Restoration via Mumford-Shah Regularization
Semi-Blind Image Restoration via Mumford-Shah Regularization L. Bar N. Sochen N. Kiryati School of Electrical Engineering Dept. of Applied Mathematics Tel Aviv University Tel Aviv, 69978, Israel Abstract
More informationMathematical modeling of textures : application to color image decomposition with a projected gradient algorithm
Mathematical modeling of textures : application to color image decomposition with a projected gradient algorithm Vincent Duval 1 & Jean-François Aujol 2 & Luminita A. Vese 3 1 Institut Telecom, Telecom
More informationResearch Article Topological Optimization with the p-laplacian Operator and an Application in Image Processing
Hindawi Publishing Corporation Boundary Value Problems Volume 9, Article ID 896813, 20 pages doi:10.1155/9/896813 Research Article Topological Optimization with the p-laplacian Operator and an Application
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 informationTwo-parameter regularization method for determining the heat source
Global Journal of Pure and Applied Mathematics. ISSN 0973-1768 Volume 13, Number 8 (017), pp. 3937-3950 Research India Publications http://www.ripublication.com Two-parameter regularization method for
More informationAsymmetric Cheeger cut and application to multi-class unsupervised clustering
Asymmetric Cheeger cut and application to multi-class unsupervised clustering Xavier Bresson Thomas Laurent April 8, 0 Abstract Cheeger cut has recently been shown to provide excellent classification results
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 informationLINEARIZED BREGMAN ITERATIONS FOR FRAME-BASED IMAGE DEBLURRING
LINEARIZED BREGMAN ITERATIONS FOR FRAME-BASED IMAGE DEBLURRING JIAN-FENG CAI, STANLEY OSHER, AND ZUOWEI SHEN Abstract. Real images usually have sparse approximations under some tight frame systems derived
More informationOn non negative solutions of some quasilinear elliptic inequalities
On non negative solutions of some quasilinear elliptic inequalities Lorenzo D Ambrosio and Enzo Mitidieri September 28 2006 Abstract Let f : R R be a continuous function. We prove that under some additional
More informationRegular article Nonlinear multigrid methods for total variation image denoising
Comput Visual Sci 7: 199 06 (004) Digital Object Identifier (DOI) 10.1007/s00791-004-0150-3 Computing and Visualization in Science Regular article Nonlinear multigrid methods for total variation image
More informationRegularization methods for large-scale, ill-posed, linear, discrete, inverse problems
Regularization methods for large-scale, ill-posed, linear, discrete, inverse problems Silvia Gazzola Dipartimento di Matematica - Università di Padova January 10, 2012 Seminario ex-studenti 2 Silvia Gazzola
More informationIMA Preprint Series # 2057
PIECEWISE H 1 + H 0 + H 1 IMAGES AND THE MUMFORD-SHAH-SOBOLEV MODEL FOR SEGMENTED IMAGE DECOMPOSITION By Jianhong (Jackie) Shen IMA Preprint Series # 2057 ( July 2005 ) INSTITUTE FOR MATHEMATICS AND ITS
More informationVariational image restoration by means of wavelets: simultaneous decomposition, deblurring and denoising
Variational image restoration by means of wavelets: simultaneous decomposition, deblurring and denoising I. Daubechies and G. Teschke December 2, 2004 Abstract Inspired by papers of Vese Osher [20] and
More informationE. YABLONOVITCH photonic crystals by using level set methods
Appl. Phys. B 81, 235 244 (2005) DOI: 10.1007/s00340-005-1877-3 Applied Physics B Lasers and Optics C.Y. KAO 1, Maximizing band gaps in two-dimensional S. OSHER 2 E. YABLONOVITCH photonic crystals by using
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 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 informationFast Multilevel Algorithm for a Minimization Problem in Impulse Noise Removal
Fast Multilevel Algorithm for a Minimization Problem in Impulse Noise Removal Raymond H. Chan and Ke Chen September 25, 7 Abstract An effective 2-phase method for removing impulse noise was recently proposed.
More information