On nonstationary preconditioned iterative regularization methods for image deblurring
|
|
- Belinda McGee
- 5 years ago
- Views:
Transcription
1 On nonstationary preconditioned iterative regularization methods for image deblurring Alessandro Buccini joint work with Prof. Marco Donatelli University of Insubria Department of Science and High Technology Cetraro Summer School Cetraro, Italy June 22-26, 2015
2 Outline 1 Image Deblur 2 Tikhonov Regularization 3 4 Conclusions and Future Work Alessandro Buccini (University of Insubria) Methods for image deblurring 1 / 27
3 Image Deblur Outline 1 Image Deblur 2 Tikhonov Regularization 3 4 Conclusions and Future Work Alessandro Buccini (University of Insubria) Methods for image deblurring 2 / 27
4 Image Deblur Image Deblur The image deblur problem is a deconvolution problem. g(x,y) = K(s x,t y)f(s,t)dsdt = K f, (1) K is called PSF (Point Spread Function). Alessandro Buccini (University of Insubria) Methods for image deblurring 3 / 27
5 Image Deblur Image Deblur The image deblur problem is a deconvolution problem. g(x,y) = K(s x,t y)f(s,t)dsdt = K f, (1) K is called PSF (Point Spread Function). Discretizing (1) we obtain a linear system Ax = b. (2) Alessandro Buccini (University of Insubria) Methods for image deblurring 3 / 27
6 Image Deblur Image Deblur The image deblur problem is a deconvolution problem. g(x,y) = K(s x,t y)f(s,t)dsdt = K f, (1) K is called PSF (Point Spread Function). Discretizing (1) we obtain a linear system It is impossible to avoid noise so Ax = b. (2) b = b true +η. Alessandro Buccini (University of Insubria) Methods for image deblurring 3 / 27
7 Image Deblur Image Deblur The image deblur problem is a deconvolution problem. g(x,y) = K(s x,t y)f(s,t)dsdt = K f, (1) K is called PSF (Point Spread Function). Discretizing (1) we obtain a linear system It is impossible to avoid noise so Ax = b. (2) b = b true +η. We can not directly solve this linear system for x because A is severely ill-conditioned, we have to use regularization methods. Alessandro Buccini (University of Insubria) Methods for image deblurring 3 / 27
8 Image Deblur In practical application we only have finite data, so a problem arises near the boundary. In particular we need to make assumptions on what is outside the field of view. Alessandro Buccini (University of Insubria) Methods for image deblurring 4 / 27
9 Image Deblur (a) Zero (b) Periodic (c) Reflexive (d) Antireflexive Reflexive: [Ng et al., SISC1999]; Antireflexive: [S. Serra-Capizzano, SISC2003]. Alessandro Buccini (University of Insubria) Methods for image deblurring 5 / 27
10 Tikhonov Regularization Outline 1 Image Deblur 2 Tikhonov Regularization 3 4 Conclusions and Future Work Alessandro Buccini (University of Insubria) Methods for image deblurring 6 / 27
11 Tikhonov Regularization Tikhonov Regularization The general Tikhonov regularization has the form with N(L) N(A) = {0}. x α = argmin x Ax b 2 2 +α Lx 2 2, (3) Alessandro Buccini (University of Insubria) Methods for image deblurring 7 / 27
12 Tikhonov Regularization Tikhonov Regularization The general Tikhonov regularization has the form with N(L) N(A) = {0}. The solution of (3) is x α = argmin x Ax b 2 2 +α Lx 2 2, (3) x α = (A A+αL L) 1 A b. Alessandro Buccini (University of Insubria) Methods for image deblurring 7 / 27
13 Tikhonov Regularization Tikhonov Regularization The general Tikhonov regularization has the form with N(L) N(A) = {0}. The solution of (3) is x α = argmin x Ax b 2 2 +α Lx 2 2, (3) x α = (A A+αL L) 1 A b. By applying a refinement technique we obtain the Iterated Tikhonov method x k+1 = x k +(A A+αL L) 1 A (b Ax k ), which can be seen as a preconditioned Landweber. Alessandro Buccini (University of Insubria) Methods for image deblurring 7 / 27
14 Outline 1 Image Deblur 2 Tikhonov Regularization 3 4 Conclusions and Future Work Alessandro Buccini (University of Insubria) Methods for image deblurring 8 / 27
15 General Idea Iterated Tikhonov works on the error equation Ae k r k, (4) where e k = x x k and r k = b Ax k. Alessandro Buccini (University of Insubria) Methods for image deblurring 9 / 27
16 General Idea Iterated Tikhonov works on the error equation where e k = x x k and r k = b Ax k. Let C be such that for some 0 < ρ < 1 2. Ae k r k, (4) (C A)z 2 ρ Az 2, z, (5) Alessandro Buccini (University of Insubria) Methods for image deblurring 9 / 27
17 General Idea Iterated Tikhonov works on the error equation where e k = x x k and r k = b Ax k. Let C be such that for some 0 < ρ < 1 2. We approximate equation (4) with Ae k r k, (4) (C A)z 2 ρ Az 2, z, (5) Ce k r k. Alessandro Buccini (University of Insubria) Methods for image deblurring 9 / 27
18 AIT Algorithm (AIT, [M. Donatelli and M. Hanke, IP2013]) Let x 0 be fixed and set r 0 = b Ax 0. Choose τ = 1+2ρ 1 2ρ with ρ as in (5), and fix q [2ρ,1]. { While r k 2 > τ η 2, let τ k = r k 2 / η 2 and q k = max q,2ρ+ 1+ρ τ k }, compute h k = C (CC +α k I) 1 r k, where α k is such that and update r k Ch k 2 = q k r k 2, x k+1 = x k +h k, r k+1 = b Ax k+1. Alessandro Buccini (University of Insubria) Methods for image deblurring 10 / 27
19 AIT Remark In AIT, L = I. Alessandro Buccini (University of Insubria) Methods for image deblurring 11 / 27
20 AIT Remark In AIT, L = I. Proposition While r k 2 > τ η 2 the norm of the reconstruction error e k = x x k decreases monotonically. Alessandro Buccini (University of Insubria) Methods for image deblurring 11 / 27
21 AIT Remark In AIT, L = I. Proposition While r k 2 > τ η 2 the norm of the reconstruction error e k = x x k decreases monotonically. Theorem Let δ b δ a function such that δ b δ b 2 true < δ. Let C,τ, and q as in AIT and denote with x δ the approximation obtained with AIT with right hand side b δ. For δ 0 we have that x δ x 0, which is the nearest solution to x 0 of the system. Alessandro Buccini (University of Insubria) Methods for image deblurring 11 / 27
22 ARIT Algorithm (ARIT, [A. B. and M. Donatelli, submitted 2015]) Let x 0 be fixed and set r 0 = b Ax 0. Choose τ = 1+2ρ 1 2ρ with ρ as in (5), and fix q [2ρ,1]. { While r k 2 > τ η 2, let τ k = r k 2 / η 2 and q k = max q,2ρ+ 1+ρ τ k }, compute h k = C (CC +α k LL ) 1 r k, where α k is such that and update r k Ch k 2 = q k r k 2, x k+1 = x k +h k, r k+1 = b Ax k+1. Alessandro Buccini (University of Insubria) Methods for image deblurring 12 / 27
23 ARIT Remark For ARIT to work we need to make the same assumptions of AIT on C and A and the following N(A) N(L) = {0}; C N(L) = A N(L) ; L and C are diagonalized by the same unitary transformation. Alessandro Buccini (University of Insubria) Methods for image deblurring 13 / 27
24 ARIT Remark For ARIT to work we need to make the same assumptions of AIT on C and A and the following N(A) N(L) = {0}; C N(L) = A N(L) ; L and C are diagonalized by the same unitary transformation. Proposition While r k 2 > τ η 2 we have L 2 Le k L 2 Le k+1 0. Alessandro Buccini (University of Insubria) Methods for image deblurring 13 / 27
25 ARIT Remark For ARIT to work we need to make the same assumptions of AIT on C and A and the following N(A) N(L) = {0}; C N(L) = A N(L) ; L and C are diagonalized by the same unitary transformation. Proposition While r k 2 > τ η 2 we have L 2 Le k L 2 Le k+1 0. Theorem As δ 0 we have that x δ x 0, with the same notation of AIT. Alessandro Buccini (University of Insubria) Methods for image deblurring 13 / 27
26 APIT Algorithm (APIT, [A. B. and M. Donatelli, submitted 2015]) Let x 0 be fixed and set r 0 = b Ax 0. Choose τ = 1+2ρ 1 2ρ with ρ as in (5), and fix q [2ρ,1]. Let Ω be a closed and convex set such that x true Ω and P Ω the metric projection over it. { While r k 2 > τ η 2, let τ k = r k 2 / η 2 and q k = max q,2ρ+ 1+ρ τ k }, compute h k = C (CC +α k I) 1 r k, where α k is such that r k Ch k 2 = q k r k 2, and update x k+1 = P Ω (x k +h k ), r k+1 = b Ax k+1. Alessandro Buccini (University of Insubria) Methods for image deblurring 14 / 27
27 APIT Remark In APIT, L = I. Alessandro Buccini (University of Insubria) Methods for image deblurring 15 / 27
28 APIT Remark In APIT, L = I. Proposition While r k 2 > τ η 2 the norm of the reconstruction error e k = x x k decreases monotonically. Alessandro Buccini (University of Insubria) Methods for image deblurring 15 / 27
29 APIT Remark In APIT, L = I. Proposition While r k 2 > τ η 2 the norm of the reconstruction error e k = x x k decreases monotonically. Theorem As δ 0 we have that x δ x 0, with the same notation of AIT. Alessandro Buccini (University of Insubria) Methods for image deblurring 15 / 27
30 APRIT Algorithm (APRIT, [A. B. and M. Donatelli, submitted 2015]) Let x 0 be fixed and set r 0 = b Ax 0. Choose τ = 1+2ρ 1 2ρ with ρ as in (5), and fix q [2ρ,1]. Let Ω be a closed and convex set such that x true Ω and P Ω the metric projection over it. { While r k 2 > τ η 2, let τ k = r k 2 / η 2 and q k = max q,2ρ+ 1+ρ τ k }, compute h k = C (CC +α k LL ) 1 r k, where α k is such that r k Ch k 2 = q k r k 2, and update x k+1 = P Ω (x k +h k ), r k+1 = b Ax k+1. Alessandro Buccini (University of Insubria) Methods for image deblurring 16 / 27
31 Numerical Experiments Numerical Experiments A is the blurring matrix with the desired boundary conditions; Alessandro Buccini (University of Insubria) Methods for image deblurring 17 / 27
32 Numerical Experiments Numerical Experiments A is the blurring matrix with the desired boundary conditions; C is the blurring matrix with the same PSF of A but with periodic boundary conditions; Alessandro Buccini (University of Insubria) Methods for image deblurring 17 / 27
33 Numerical Experiments Numerical Experiments A is the blurring matrix with the desired boundary conditions; C is the blurring matrix with the same PSF of A but with periodic boundary conditions; L = L 1 I +I L 1, where L 1 is defined as L 1 = ; Alessandro Buccini (University of Insubria) Methods for image deblurring 17 / 27
34 Numerical Experiments Numerical Experiments A is the blurring matrix with the desired boundary conditions; C is the blurring matrix with the same PSF of A but with periodic boundary conditions; L = L 1 I +I L 1, where L 1 is defined as L 1 = ; Ω = {x R M : i = 1,...,M, x i 0}; Alessandro Buccini (University of Insubria) Methods for image deblurring 17 / 27
35 Numerical Experiments Numerical Experiments A is the blurring matrix with the desired boundary conditions; C is the blurring matrix with the same PSF of A but with periodic boundary conditions; L = L 1 I +I L 1, where L 1 is defined as L 1 = ; Ω = {x R M : i = 1,...,M, x i 0}; ρ = 10 3, q = 0.7; Alessandro Buccini (University of Insubria) Methods for image deblurring 17 / 27
36 Numerical Experiments Numerical Experiments A is the blurring matrix with the desired boundary conditions; C is the blurring matrix with the same PSF of A but with periodic boundary conditions; L = L 1 I +I L 1, where L 1 is defined as L 1 = ; Ω = {x R M : i = 1,...,M, x i 0}; ρ = 10 3, q = 0.7; Note that N(L) N(A) = {0}. Alessandro Buccini (University of Insubria) Methods for image deblurring 17 / 27
37 Numerical Experiments (a) (b) (c) Figure: Barbara Test case: (a) Barbara Test image, (b) Diagonal motion PSF, (c) Blurred image, η 2 = 0.03 b true 2. Since the image is generic we use the antireflexive boundary conditions. Alessandro Buccini (University of Insubria) Methods for image deblurring 18 / 27
38 Numerical Experiments AIT APIT ARIT APRIT Figure: Barbara Test Case: RRE History Alessandro Buccini (University of Insubria) Methods for image deblurring 19 / 27
39 Numerical Experiments Method RRE Computational Time AIT sec. ARIT sec. APIT sec. APRIT sec. Hybrid (Optimal: ) sec. TwIST sec. FlexiAT sec. RRAT sec. NN-ReStart-GAT sec. Table: Barbara: Relative errors comparisons Hybrid: [J. Chung, J. Nagy and, D. O Leary, ETNA2008]; TwIST: [J. Bioucas-Dias and M. Figueiredo, IEEE2007]; FlexiAT: [S. Gazzola and J. Nagy, SISC2014]; RRAT: [A. Neuman, L. Reichel and, H. Sadok, NA2012]; NN-ReStart-GAT: [S. Gazzola and J. Nagy, SISC2014]. Alessandro Buccini (University of Insubria) Methods for image deblurring 20 / 27
40 Numerical Experiments (a) (b) (c) (d) Figure: Barbara Reconstructions: (a) ARIT, (b) APIT, (c) Hybrid, (d) TwIST. Alessandro Buccini (University of Insubria) Methods for image deblurring 21 / 27
41 Numerical Experiments (a) (b) (c) Figure: Satellite Test case: (a) Test image, (b) Astronomic PSF, (c) Blurred image, η 2 = 0.01 b true 2. Since the image is all black around the boundaries we use the zero boundary conditions. Alessandro Buccini (University of Insubria) Methods for image deblurring 22 / 27
42 Numerical Experiments AIT APIT ARIT APRIT Figure: Satellite Test Case: RRE History Alessandro Buccini (University of Insubria) Methods for image deblurring 23 / 27
43 Numerical Experiments Method RRE Computational Time AIT sec. ARIT sec. APIT sec. APRIT sec. Hybrid sec. TwIST sec. FlexiAT sec. RRAT sec. NN-ReStart-GAT sec. Table: Satellite: Relative errors comparisons Alessandro Buccini (University of Insubria) Methods for image deblurring 24 / 27
44 Numerical Experiments (a) (b) (c) (d) Figure: Satellite Reconstructions: (a) APIT, (b) APRIT, (c) RRAT, (d) NN-Restart-GAT. Alessandro Buccini (University of Insubria) Methods for image deblurring 25 / 27
45 Conclusions and Future Work Outline 1 Image Deblur 2 Tikhonov Regularization 3 4 Conclusions and Future Work Alessandro Buccini (University of Insubria) Methods for image deblurring 26 / 27
46 Conclusions and Future Work Conclusions and Future Work We proposed four methods with good properties Alessandro Buccini (University of Insubria) Methods for image deblurring 27 / 27
47 Conclusions and Future Work Conclusions and Future Work We proposed four methods with good properties Almost parameter free; Alessandro Buccini (University of Insubria) Methods for image deblurring 27 / 27
48 Conclusions and Future Work Conclusions and Future Work We proposed four methods with good properties Almost parameter free; Capable of dealing with high level of noise; Alessandro Buccini (University of Insubria) Methods for image deblurring 27 / 27
49 Conclusions and Future Work Conclusions and Future Work We proposed four methods with good properties Almost parameter free; Capable of dealing with high level of noise; Any boundary condition can be easily implemented. Alessandro Buccini (University of Insubria) Methods for image deblurring 27 / 27
50 Conclusions and Future Work Conclusions and Future Work We proposed four methods with good properties Almost parameter free; Capable of dealing with high level of noise; Any boundary condition can be easily implemented. Alessandro Buccini (University of Insubria) Methods for image deblurring 27 / 27
51 Conclusions and Future Work Conclusions and Future Work We proposed four methods with good properties Almost parameter free; Capable of dealing with high level of noise; Any boundary condition can be easily implemented. Possible future improvements Alessandro Buccini (University of Insubria) Methods for image deblurring 27 / 27
52 Conclusions and Future Work Conclusions and Future Work We proposed four methods with good properties Almost parameter free; Capable of dealing with high level of noise; Any boundary condition can be easily implemented. Possible future improvements Non spatially invariant blur; Alessandro Buccini (University of Insubria) Methods for image deblurring 27 / 27
53 Conclusions and Future Work Conclusions and Future Work We proposed four methods with good properties Almost parameter free; Capable of dealing with high level of noise; Any boundary condition can be easily implemented. Possible future improvements Non spatially invariant blur; More advanced regularization term (p norms, TV,...); Alessandro Buccini (University of Insubria) Methods for image deblurring 27 / 27
54 Conclusions and Future Work Conclusions and Future Work We proposed four methods with good properties Almost parameter free; Capable of dealing with high level of noise; Any boundary condition can be easily implemented. Possible future improvements Non spatially invariant blur; More advanced regularization term (p norms, TV,...); Blind deconvolution. Alessandro Buccini (University of Insubria) Methods for image deblurring 27 / 27
55 Thank you for your attention.
A fast nonstationary preconditioning strategy for ill-posed problems, with application to image deblurring
A fast nonstationary preconditioning strategy for ill-posed problems, with application to image deblurring Marco Donatelli Dept. of Science and High Tecnology U. Insubria (Italy) Joint work with M. Hanke
More informationAn l 2 -l q regularization method for large discrete ill-posed problems
Noname manuscript No. (will be inserted by the editor) An l -l q regularization method for large discrete ill-posed problems Alessandro Buccini Lothar Reichel the date of receipt and acceptance should
More informationIterative Krylov Subspace Methods for Sparse Reconstruction
Iterative Krylov Subspace Methods for Sparse Reconstruction James Nagy Mathematics and Computer Science Emory University Atlanta, GA USA Joint work with Silvia Gazzola University of Padova, Italy Outline
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 informationPreconditioning. Noisy, Ill-Conditioned Linear Systems
Preconditioning Noisy, Ill-Conditioned Linear Systems James G. Nagy Emory University Atlanta, GA Outline 1. The Basic Problem 2. Regularization / Iterative Methods 3. Preconditioning 4. Example: Image
More informationIterated Tikhonov regularization with a general penalty term
NUMERICAL LINEAR ALGEBRA WITH APPLICATIONS Numer. Linear Algebra Appl. ; 00:1 19 Published online in Wiley InterScience (www.interscience.wiley.com). Iterated Tihonov regularization with a general penalty
More informationTikhonov-type iterative regularization methods for ill-posed inverse problems: theoretical aspects and applications
UNIVERSITY OF INSUBRIA Department of Science and High Technology Ph.D. course in Computer Science and Computational Mathematics (XXIX cycle Tikhonov-type iterative regularization methods for ill-posed
More informationPreconditioning. Noisy, Ill-Conditioned Linear Systems
Preconditioning Noisy, Ill-Conditioned Linear Systems James G. Nagy Emory University Atlanta, GA Outline 1. The Basic Problem 2. Regularization / Iterative Methods 3. Preconditioning 4. Example: Image
More informationAn iterative multigrid regularization method for Toeplitz discrete ill-posed problems
NUMERICAL MATHEMATICS: Theory, Methods and Applications Numer. Math. Theor. Meth. Appl., Vol. xx, No. x, pp. 1-18 (200x) An iterative multigrid regularization method for Toeplitz discrete ill-posed problems
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 informationInverse Ill Posed Problems in Image Processing
Inverse Ill Posed Problems in Image Processing Image Deblurring I. Hnětynková 1,M.Plešinger 2,Z.Strakoš 3 hnetynko@karlin.mff.cuni.cz, martin.plesinger@tul.cz, strakos@cs.cas.cz 1,3 Faculty of Mathematics
More informationStatistically-Based Regularization Parameter Estimation for Large Scale Problems
Statistically-Based Regularization Parameter Estimation for Large Scale Problems Rosemary Renaut Joint work with Jodi Mead and Iveta Hnetynkova March 1, 2010 National Science Foundation: Division of Computational
More informationarxiv: v1 [math.na] 15 Jun 2009
Noname manuscript No. (will be inserted by the editor) Fast transforms for high order boundary conditions Marco Donatelli arxiv:0906.2704v1 [math.na] 15 Jun 2009 the date of receipt and acceptance should
More informationDealing with edge effects in least-squares image deconvolution problems
Astronomy & Astrophysics manuscript no. bc May 11, 05 (DOI: will be inserted by hand later) Dealing with edge effects in least-squares image deconvolution problems R. Vio 1 J. Bardsley 2, M. Donatelli
More informationBlind image restoration as a convex optimization problem
Int. J. Simul. Multidisci.Des. Optim. 4, 33-38 (2010) c ASMDO 2010 DOI: 10.1051/ijsmdo/ 2010005 Available online at: http://www.ijsmdo.org Blind image restoration as a convex optimization problem A. Bouhamidi
More informationKrylov subspace iterative methods for nonsymmetric discrete ill-posed problems in image restoration
Krylov subspace iterative methods for nonsymmetric discrete ill-posed problems in image restoration D. Calvetti a, B. Lewis b and L. Reichel c a Department of Mathematics, Case Western Reserve University,
More informationLecture 1: Numerical Issues from Inverse Problems (Parameter Estimation, Regularization Theory, and Parallel Algorithms)
Lecture 1: Numerical Issues from Inverse Problems (Parameter Estimation, Regularization Theory, and Parallel Algorithms) Youzuo Lin 1 Joint work with: Rosemary A. Renaut 2 Brendt Wohlberg 1 Hongbin Guo
More informationThe Chi-squared Distribution of the Regularized Least Squares Functional for Regularization Parameter Estimation
The Chi-squared Distribution of the Regularized Least Squares Functional for Regularization Parameter Estimation Rosemary Renaut Collaborators: Jodi Mead and Iveta Hnetynkova DEPARTMENT OF MATHEMATICS
More informationA Hybrid LSQR Regularization Parameter Estimation Algorithm for Large Scale Problems
A Hybrid LSQR Regularization Parameter Estimation Algorithm for Large Scale Problems Rosemary Renaut Joint work with Jodi Mead and Iveta Hnetynkova SIAM Annual Meeting July 10, 2009 National Science Foundation:
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 informationThe Chi-squared Distribution of the Regularized Least Squares Functional for Regularization Parameter Estimation
The Chi-squared Distribution of the Regularized Least Squares Functional for Regularization Parameter Estimation Rosemary Renaut DEPARTMENT OF MATHEMATICS AND STATISTICS Prague 2008 MATHEMATICS AND STATISTICS
More informationIterative regularization of nonlinear ill-posed problems in Banach space
Iterative regularization of nonlinear ill-posed problems in Banach space Barbara Kaltenbacher, University of Klagenfurt joint work with Bernd Hofmann, Technical University of Chemnitz, Frank Schöpfer and
More informationNumerical Methods for Separable Nonlinear Inverse Problems with Constraint and Low Rank
Numerical Methods for Separable Nonlinear Inverse Problems with Constraint and Low Rank Taewon Cho Thesis submitted to the Faculty of the Virginia Polytechnic Institute and State University in partial
More information1. Introduction. We are concerned with the approximate solution of linear least-squares problems
FRACTIONAL TIKHONOV REGULARIZATION WITH A NONLINEAR PENALTY TERM SERENA MORIGI, LOTHAR REICHEL, AND FIORELLA SGALLARI Abstract. Tikhonov regularization is one of the most popular methods for solving linear
More informationAMS classification scheme numbers: 65F10, 65F15, 65Y20
Improved image deblurring with anti-reflective boundary conditions and re-blurring (This is a preprint of an article published in Inverse Problems, 22 (06) pp. 35-53.) M. Donatelli, C. Estatico, A. Martinelli,
More informationResearch Article Vector Extrapolation Based Landweber Method for Discrete Ill-Posed Problems
Hindawi Mathematical Problems in Engineering Volume 2017, Article ID 1375716, 8 pages https://doi.org/10.1155/2017/1375716 Research Article Vector Extrapolation Based Landweber Method for Discrete Ill-Posed
More informationA Double Regularization Approach for Inverse Problems with Noisy Data and Inexact Operator
A Double Regularization Approach for Inverse Problems with Noisy Data and Inexact Operator Ismael Rodrigo Bleyer Prof. Dr. Ronny Ramlau Johannes Kepler Universität - Linz Florianópolis - September, 2011.
More informationIterative Methods for Smooth Objective Functions
Optimization Iterative Methods for Smooth Objective Functions Quadratic Objective Functions Stationary Iterative Methods (first/second order) Steepest Descent Method Landweber/Projected Landweber Methods
More informationA Fast Alternating Minimization Algorithm for Total Variation Deblurring Without Boundary Artifacts
A Fast Alternating Minimization Algorithm for Total Variation Deblurring Without Boundary Artifacts Zheng-Jian Bai Daniele Cassani Marco Donatelli Stefano Serra-Capizzano November 7, 13 Abstract Recently,
More informationSquare Smoothing Regularization Matrices with Accurate Boundary Conditions
Square Smoothing Regularization Matrices with Accurate Boundary Conditions M. Donatelli a,, L. Reichel b a Dipartimento di Scienza e Alta Tecnologia, Università dell Insubria, Via Valleggio 11, 22100 Como,
More informationSPARSE SIGNAL RESTORATION. 1. Introduction
SPARSE SIGNAL RESTORATION IVAN W. SELESNICK 1. Introduction These notes describe an approach for the restoration of degraded signals using sparsity. This approach, which has become quite popular, is useful
More informationAn Interior-Point Trust-Region-Based Method for Large-Scale Nonnegative Regularization
An Interior-Point Trust-Region-Based Method for Large-Scale Nonnegative Regularization Marielba Rojas Trond Steihaug July 6, 2001 (Revised December 19, 2001) CERFACS Technical Report TR/PA/01/11 Abstract
More informationAn Alternating Projections Algorithm for Sparse Blind Deconvolution
An Alternating Projections Algorithm for Sparse Blind Deconvolution Aniruddha Adiga Department of Electrical Engineering Indian Institute of Science, Bangalore Email: aniruddha@ee.iisc.ernet.in Advisor:
More informationA Double Regularization Approach for Inverse Problems with Noisy Data and Inexact Operator
A Double Regularization Approach for Inverse Problems with Noisy Data and Inexact Operator Ismael Rodrigo Bleyer Prof. Dr. Ronny Ramlau Johannes Kepler Universität - Linz Cambridge - July 28, 211. Doctoral
More informationThe Global Krylov subspace methods and Tikhonov regularization for image restoration
The Global Krylov subspace methods and Tikhonov regularization for image restoration Abderrahman BOUHAMIDI (joint work with Khalide Jbilou) Université du Littoral Côte d Opale LMPA, CALAIS-FRANCE bouhamidi@lmpa.univ-littoral.fr
More informationETNA Kent State University
Electronic Transactions on Numerical Analysis. Volume 31, pp. 204-220, 2008. Copyright 2008,. ISSN 1068-9613. ETNA NOISE PROPAGATION IN REGULARIZING ITERATIONS FOR IMAGE DEBLURRING PER CHRISTIAN HANSEN
More informationNewton s method for obtaining the Regularization Parameter for Regularized Least Squares
Newton s method for obtaining the Regularization Parameter for Regularized Least Squares Rosemary Renaut, Jodi Mead Arizona State and Boise State International Linear Algebra Society, Cancun Renaut and
More informationLevenberg-Marquardt method in Banach spaces with general convex regularization terms
Levenberg-Marquardt method in Banach spaces with general convex regularization terms Qinian Jin Hongqi Yang Abstract We propose a Levenberg-Marquardt method with general uniformly convex regularization
More informationScaled gradient projection methods in image deblurring and denoising
Scaled gradient projection methods in image deblurring and denoising Mario Bertero 1 Patrizia Boccacci 1 Silvia Bonettini 2 Riccardo Zanella 3 Luca Zanni 3 1 Dipartmento di Matematica, Università di Genova
More informationLinear Inverse Problems
Linear Inverse Problems Ajinkya Kadu Utrecht University, The Netherlands February 26, 2018 Outline Introduction Least-squares Reconstruction Methods Examples Summary Introduction 2 What are inverse problems?
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 informationRapid, Robust, and Reliable Blind Deconvolution via Nonconvex Optimization
Rapid, Robust, and Reliable Blind Deconvolution via Nonconvex Optimization Shuyang Ling Department of Mathematics, UC Davis Oct.18th, 2016 Shuyang Ling (UC Davis) 16w5136, Oaxaca, Mexico Oct.18th, 2016
More informationInverse problem and optimization
Inverse problem and optimization Laurent Condat, Nelly Pustelnik CNRS, Gipsa-lab CNRS, Laboratoire de Physique de l ENS de Lyon Decembre, 15th 2016 Inverse problem and optimization 2/36 Plan 1. Examples
More informationOn the regularization properties of some spectral gradient methods
On the regularization properties of some spectral gradient methods Daniela di Serafino Department of Mathematics and Physics, Second University of Naples daniela.diserafino@unina2.it contributions from
More informationNewton s Method for Estimating the Regularization Parameter for Least Squares: Using the Chi-curve
Newton s Method for Estimating the Regularization Parameter for Least Squares: Using the Chi-curve Rosemary Renaut, Jodi Mead Arizona State and Boise State Copper Mountain Conference on Iterative Methods
More informationTowards Improved Sensitivity in Feature Extraction from Signals: one and two dimensional
Towards Improved Sensitivity in Feature Extraction from Signals: one and two dimensional Rosemary Renaut, Hongbin Guo, Jodi Mead and Wolfgang Stefan Supported by Arizona Alzheimer s Research Center and
More informationarxiv: v1 [math.na] 3 Jan 2019
STRUCTURED FISTA FOR IMAGE RESTORATION ZIXUAN CHEN, JAMES G. NAGY, YUANZHE XI, AND BO YU arxiv:9.93v [math.na] 3 Jan 29 Abstract. In this paper, we propose an efficient numerical scheme for solving some
More informationKey words. Boundary conditions, fast transforms, matrix algebras and Toeplitz matrices, Tikhonov regularization, regularizing iterative methods.
REGULARIZATION OF IMAGE RESTORATION PROBLEMS WITH ANTI-REFLECTIVE BOUNDARY CONDITIONS MARCO DONATELLI, CLAUDIO ESTATICO, AND STEFANO SERRA-CAPIZZANO Abstract. Anti-reflective boundary conditions have been
More informationINVERSE SUBSPACE PROBLEMS WITH APPLICATIONS
INVERSE SUBSPACE PROBLEMS WITH APPLICATIONS SILVIA NOSCHESE AND LOTHAR REICHEL Abstract. Given a square matrix A, the inverse subspace problem is concerned with determining a closest matrix to A with a
More informationETNA Kent State University
Electronic Transactions on Numerical Analysis. Volume 28, pp. 49-67, 28. Copyright 28,. ISSN 68-963. A WEIGHTED-GCV METHOD FOR LANCZOS-HYBRID REGULARIZATION JULIANNE CHUNG, JAMES G. NAGY, AND DIANNE P.
More informationNUMERICAL OPTIMIZATION METHODS FOR BLIND DECONVOLUTION
NUMERICAL OPTIMIZATION METHODS FOR BLIND DECONVOLUTION ANASTASIA CORNELIO, ELENA LOLI PICCOLOMINI, AND JAMES G. NAGY Abstract. This paper describes a nonlinear least squares framework to solve a separable
More informationMathematical Beer Goggles or The Mathematics of Image Processing
How Mathematical Beer Goggles or The Mathematics of Image Processing Department of Mathematical Sciences University of Bath Postgraduate Seminar Series University of Bath 12th February 2008 1 How 2 How
More informationSelf-Calibration and Biconvex Compressive Sensing
Self-Calibration and Biconvex Compressive Sensing Shuyang Ling Department of Mathematics, UC Davis July 12, 2017 Shuyang Ling (UC Davis) SIAM Annual Meeting, 2017, Pittsburgh July 12, 2017 1 / 22 Acknowledgements
More informationDownloaded 06/11/15 to Redistribution subject to SIAM license or copyright; see
SIAM J. SCI. COMPUT. Vol. 37, No. 2, pp. B332 B359 c 2015 Society for Industrial and Applied Mathematics A FRAMEWORK FOR REGULARIZATION VIA OPERATOR APPROXIMATION JULIANNE M. CHUNG, MISHA E. KILMER, AND
More informationProgetto di Ricerca GNCS 2016 PING Problemi Inversi in Geofisica Firenze, 6 aprile Regularized nonconvex minimization for image restoration
Progetto di Ricerca GNCS 2016 PING Problemi Inversi in Geofisica Firenze, 6 aprile 2016 Regularized nonconvex minimization for image restoration Claudio Estatico Joint work with: Fabio Di Benedetto, Flavia
More informationAMS526: Numerical Analysis I (Numerical Linear Algebra) Lecture 23: GMRES and Other Krylov Subspace Methods; Preconditioning
AMS526: Numerical Analysis I (Numerical Linear Algebra) Lecture 23: GMRES and Other Krylov Subspace Methods; Preconditioning Xiangmin Jiao SUNY Stony Brook Xiangmin Jiao Numerical Analysis I 1 / 18 Outline
More informationSIGNAL AND IMAGE RESTORATION: SOLVING
1 / 55 SIGNAL AND IMAGE RESTORATION: SOLVING ILL-POSED INVERSE PROBLEMS - ESTIMATING PARAMETERS Rosemary Renaut http://math.asu.edu/ rosie CORNELL MAY 10, 2013 2 / 55 Outline Background Parameter Estimation
More informationArnoldi-Tikhonov regularization methods
Arnoldi-Tikhonov regularization methods Bryan Lewis a, Lothar Reichel b,,1 a Rocketcalc LLC, 100 W. Crain Ave., Kent, OH 44240, USA. b Department of Mathematical Sciences, Kent State University, Kent,
More informationAlternating Krylov Subspace Image Restoration Methods
Alternating Krylov Subspace Image Restoration Methods J. O. Abad a, S. Morigi b, L. Reichel c, F. Sgallari d, a Instituto de Ciencia y Tecnologia de Materiales (IMRE), Universidad de La Habana, Habana,
More informationStochastic dynamical modeling:
Stochastic dynamical modeling: Structured matrix completion of partially available statistics Armin Zare www-bcf.usc.edu/ arminzar Joint work with: Yongxin Chen Mihailo R. Jovanovic Tryphon T. Georgiou
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 informationAn example of Bayesian reasoning Consider the one-dimensional deconvolution problem with various degrees of prior information.
An example of Bayesian reasoning Consider the one-dimensional deconvolution problem with various degrees of prior information. Model: where g(t) = a(t s)f(s)ds + e(t), a(t) t = (rapidly). The problem,
More informationMathematics and Computer Science
Technical Report TR-2004-012 Kronecker Product Approximation for Three-Dimensional Imaging Applications by MIsha Kilmer, James Nagy Mathematics and Computer Science EMORY UNIVERSITY Kronecker Product Approximation
More informationDEBLURRING AND SPARSE UNMIXING OF HYPERSPECTRAL IMAGES USING MULTIPLE POINT SPREAD FUNCTIONS G = XM + N,
DEBLURRING AND SPARSE UNMIXING OF HYPERSPECTRAL IMAGES USING MULTIPLE POINT SPREAD FUNCTIONS SEBASTIAN BERISHA, JAMES G NAGY, AND ROBERT J PLEMMONS Abstract This paper is concerned with deblurring and
More informationNonstationary Iterated Thresholding Algorithms for Image Deblurring
Nonstationary Iterated Thresholding Algorithms for Image Deblurring J. Huang, M. Donatelli, and R. Chan Abstract We propose iterative thresholding algorithms based on the iterated Tikhonov method for image
More informationRegularization Parameter Estimation for Least Squares: A Newton method using the χ 2 -distribution
Regularization Parameter Estimation for Least Squares: A Newton method using the χ 2 -distribution Rosemary Renaut, Jodi Mead Arizona State and Boise State September 2007 Renaut and Mead (ASU/Boise) Scalar
More informationAccelerated Gradient Methods for Constrained Image Deblurring
Accelerated Gradient Methods for Constrained Image Deblurring S Bonettini 1, R Zanella 2, L Zanni 2, M Bertero 3 1 Dipartimento di Matematica, Università di Ferrara, Via Saragat 1, Building B, I-44100
More informationGeneralized Local Regularization for Ill-Posed Problems
Generalized Local Regularization for Ill-Posed Problems Patricia K. Lamm Department of Mathematics Michigan State University AIP29 July 22, 29 Based on joint work with Cara Brooks, Zhewei Dai, and Xiaoyue
More informationLanczos tridigonalization and Golub - Kahan bidiagonalization: Ideas, connections and impact
Lanczos tridigonalization and Golub - Kahan bidiagonalization: Ideas, connections and impact Zdeněk Strakoš Academy of Sciences and Charles University, Prague http://www.cs.cas.cz/ strakos Hong Kong, February
More informationRegularization in Banach Space
Regularization in Banach Space Barbara Kaltenbacher, Alpen-Adria-Universität Klagenfurt joint work with Uno Hämarik, University of Tartu Bernd Hofmann, Technical University of Chemnitz Urve Kangro, University
More informationBlind Deconvolution Using Convex Programming. Jiaming Cao
Blind Deconvolution Using Convex Programming Jiaming Cao Problem Statement The basic problem Consider that the received signal is the circular convolution of two vectors w and x, both of length L. How
More informationConjugate Gradient Tutorial
Conjugate Gradient Tutorial Prof. Chung-Kuan Cheng Computer Science and Engineering Department University of California, San Diego ckcheng@ucsd.edu December 1, 2015 Prof. Chung-Kuan Cheng (UC San Diego)
More informationAdaptive discretization and first-order methods for nonsmooth inverse problems for PDEs
Adaptive discretization and first-order methods for nonsmooth inverse problems for PDEs Christian Clason Faculty of Mathematics, Universität Duisburg-Essen joint work with Barbara Kaltenbacher, Tuomo Valkonen,
More informationDeblurring Jupiter (sampling in GLIP faster than regularized inversion) Colin Fox Richard A. Norton, J.
Deblurring Jupiter (sampling in GLIP faster than regularized inversion) Colin Fox fox@physics.otago.ac.nz Richard A. Norton, J. Andrés Christen Topics... Backstory (?) Sampling in linear-gaussian hierarchical
More informationCOMPUTATION OF FOURIER TRANSFORMS FOR NOISY B FOR NOISY BANDLIMITED SIGNALS
COMPUTATION OF FOURIER TRANSFORMS FOR NOISY BANDLIMITED SIGNALS October 22, 2011 I. Introduction Definition of Fourier transform: F [f ](ω) := ˆf (ω) := + f (t)e iωt dt, ω R (1) I. Introduction Definition
More informationCOMPUTATIONAL ISSUES RELATING TO INVERSION OF PRACTICAL DATA: WHERE IS THE UNCERTAINTY? CAN WE SOLVE Ax = b?
COMPUTATIONAL ISSUES RELATING TO INVERSION OF PRACTICAL DATA: WHERE IS THE UNCERTAINTY? CAN WE SOLVE Ax = b? Rosemary Renaut http://math.asu.edu/ rosie BRIDGING THE GAP? OCT 2, 2012 Discussion Yuen: Solve
More informationSIMPLIFIED GSVD COMPUTATIONS FOR THE SOLUTION OF LINEAR DISCRETE ILL-POSED PROBLEMS (1.1) x R Rm n, x R n, b R m, m n,
SIMPLIFIED GSVD COMPUTATIONS FOR THE SOLUTION OF LINEAR DISCRETE ILL-POSED PROBLEMS L. DYKES AND L. REICHEL Abstract. The generalized singular value decomposition (GSVD) often is used to solve Tikhonov
More informationParallel Cimmino-type methods for ill-posed problems
Parallel Cimmino-type methods for ill-posed problems Cao Van Chung Seminar of Centro de Modelización Matemática Escuela Politécnica Naciónal Quito ModeMat, EPN, Quito Ecuador cao.vanchung@epn.edu.ec, cvanchung@gmail.com
More informationNumerical Methods for Parameter Estimation in Poisson Data Inversion
DOI 10.1007/s10851-014-0553-9 Numerical Methods for Parameter Estimation in Poisson Data Inversion Luca Zanni Alessandro Benfenati Mario Bertero Valeria Ruggiero Received: 28 June 2014 / Accepted: 11 December
More informationREGULARIZATION MATRICES FOR DISCRETE ILL-POSED PROBLEMS IN SEVERAL SPACE-DIMENSIONS
REGULARIZATION MATRICES FOR DISCRETE ILL-POSED PROBLEMS IN SEVERAL SPACE-DIMENSIONS LAURA DYKES, GUANGXIN HUANG, SILVIA NOSCHESE, AND LOTHAR REICHEL Abstract. Many applications in science and engineering
More informationTikhonov Regularization of Large Symmetric Problems
NUMERICAL LINEAR ALGEBRA WITH APPLICATIONS Numer. Linear Algebra Appl. 2000; 00:1 11 [Version: 2000/03/22 v1.0] Tihonov Regularization of Large Symmetric Problems D. Calvetti 1, L. Reichel 2 and A. Shuibi
More informationLPA-ICI Applications in Image Processing
LPA-ICI Applications in Image Processing Denoising Deblurring Derivative estimation Edge detection Inverse halftoning Denoising Consider z (x) =y (x)+η (x), wherey is noise-free image and η is noise. assume
More informationRegolarizzazione adattiva in spazi di Lebesgue a esponente variabile. Claudio Estatico Dipartimento di Matematica, Università di Genova, Italy
GNCS 2018 Convegno Biennale del Gruppo Nazionale per il Calcolo Scientifico Montecatini Terme, 14-16 febbraio 2018 Regolarizzazione adattiva in spazi di Lebesgue a esponente variabile Claudio Estatico
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 informationChapter 7 Iterative Techniques in Matrix Algebra
Chapter 7 Iterative Techniques in Matrix Algebra Per-Olof Persson persson@berkeley.edu Department of Mathematics University of California, Berkeley Math 128B Numerical Analysis Vector Norms Definition
More informationETNA Kent State University
Electronic Transactions on Numerical Analysis. Volume 39, pp. 476-57,. Copyright,. ISSN 68-963. ETNA ON THE MINIMIZATION OF A TIKHONOV FUNCTIONAL WITH A NON-CONVEX SPARSITY CONSTRAINT RONNY RAMLAU AND
More informationGlobal Golub Kahan bidiagonalization applied to large discrete ill-posed problems
Global Golub Kahan bidiagonalization applied to large discrete ill-posed problems A. H. Bentbib Laboratory LAMAI, FSTG, University of Cadi Ayyad, Marrakesh, Morocco. M. El Guide Laboratory LAMAI, FSTG,
More informationTRACKING SOLUTIONS OF TIME VARYING LINEAR INVERSE PROBLEMS
TRACKING SOLUTIONS OF TIME VARYING LINEAR INVERSE PROBLEMS Martin Kleinsteuber and Simon Hawe Department of Electrical Engineering and Information Technology, Technische Universität München, München, Arcistraße
More informationEE 367 / CS 448I Computational Imaging and Display Notes: Noise, Denoising, and Image Reconstruction with Noise (lecture 10)
EE 367 / CS 448I Computational Imaging and Display Notes: Noise, Denoising, and Image Reconstruction with Noise (lecture 0) Gordon Wetzstein gordon.wetzstein@stanford.edu This document serves as a supplement
More informationLaplace-distributed increments, the Laplace prior, and edge-preserving regularization
J. Inverse Ill-Posed Probl.? (????), 1 15 DOI 1.1515/JIIP.????.??? de Gruyter???? Laplace-distributed increments, the Laplace prior, and edge-preserving regularization Johnathan M. Bardsley Abstract. For
More informationShiqian Ma, MAT-258A: Numerical Optimization 1. Chapter 9. Alternating Direction Method of Multipliers
Shiqian Ma, MAT-258A: Numerical Optimization 1 Chapter 9 Alternating Direction Method of Multipliers Shiqian Ma, MAT-258A: Numerical Optimization 2 Separable convex optimization a special case is min f(x)
More informationComputer Vision & Digital Image Processing
Computer Vision & Digital Image Processing Image Restoration and Reconstruction I Dr. D. J. Jackson Lecture 11-1 Image restoration Restoration is an objective process that attempts to recover an image
More informationInexact Alternating Direction Method of Multipliers for Separable Convex Optimization
Inexact Alternating Direction Method of Multipliers for Separable Convex Optimization Hongchao Zhang hozhang@math.lsu.edu Department of Mathematics Center for Computation and Technology Louisiana State
More informationA MODIFIED TSVD METHOD FOR DISCRETE ILL-POSED PROBLEMS
A MODIFIED TSVD METHOD FOR DISCRETE ILL-POSED PROBLEMS SILVIA NOSCHESE AND LOTHAR REICHEL Abstract. Truncated singular value decomposition (TSVD) is a popular method for solving linear discrete ill-posed
More informationExtension of GKB-FP algorithm to large-scale general-form Tikhonov regularization
Extension of GKB-FP algorithm to large-scale general-form Tikhonov regularization Fermín S. Viloche Bazán, Maria C. C. Cunha and Leonardo S. Borges Department of Mathematics, Federal University of Santa
More informationarxiv: v1 [math.na] 21 Aug 2014 Barbara Kaltenbacher
ENHANCED CHOICE OF THE PARAMETERS IN AN ITERATIVELY REGULARIZED NEWTON- LANDWEBER ITERATION IN BANACH SPACE arxiv:48.526v [math.na] 2 Aug 24 Barbara Kaltenbacher Alpen-Adria-Universität Klagenfurt Universitätstrasse
More informationNoise-Blind Image Deblurring Supplementary Material
Noise-Blind Image Deblurring Supplementary Material Meiguang Jin University of Bern Switzerland Stefan Roth TU Darmstadt Germany Paolo Favaro University of Bern Switzerland A. Upper and Lower Bounds Our
More informationComparison of A-Posteriori Parameter Choice Rules for Linear Discrete Ill-Posed Problems
Comparison of A-Posteriori Parameter Choice Rules for Linear Discrete Ill-Posed Problems Alessandro Buccini a, Yonggi Park a, Lothar Reichel a a Department of Mathematical Sciences, Kent State University,
More informationProjected Nonstationary Iterated Tikhonov Regularization
BIT manuscript No. (will be inserted by the editor) Projected Nonstationary Iterated Tihonov Regularization Guangxin Huang Lothar Reichel Feng Yin Received: date / Accepted: date Dedicated to Heinrich
More information10. Multi-objective least squares
L Vandenberghe ECE133A (Winter 2018) 10 Multi-objective least squares multi-objective least squares regularized data fitting control estimation and inversion 10-1 Multi-objective least squares we have
More information