Sparse tomography. Samuli Siltanen. Department of Mathematics and Statistics University of Helsinki, Finland
|
|
- Rose Montgomery
- 6 years ago
- Views:
Transcription
1 Sparse tomography Samuli Siltanen Department of Mathematics and Statistics University of Helsinki, Finland Minisymposium: Fourier analytic methods in tomographic image reconstruction Applied Inverse Problems Conference Daejeon, Korea, July 1, 2013
2 Finland
3 This is a joint work with Keijo Hämäläinen, University of Helsinki Aki Kallonen, University of Helsinki Ville Kolehmainen, University of Eastern Finland Matti Lassas, University of Helsinki Esa Niemi, University of Helsinki Kati Niinimäki, University of Eastern Finland Eero Saksman, University of Helsinki
4 Outline Sparse sampling and tomography Total variation regularization Discretization-invariance Besov space regularization Parameter choice: the S-curve method
5 Outline Sparse sampling and tomography Total variation regularization Discretization-invariance Besov space regularization Parameter choice: the S-curve method
6 Let us study a simple two-dimensional example of tomographic imaging
7 Tomography is based on measuring densities of matter using X-ray attenuation data Detector X-ray source (=4+4+5)
8 A projection image is produced by parallel X-rays and several detector pixels (here three pixels) Detector (=4+4+5) 8 (=1+3+4) 3 (=1+0+2)
9 For tomographic imaging it is essential to record projection images from different directions
10 The length of X-rays traveling inside each pixel is important, thus here the square roots
11 The direct problem of tomography is to find the projection images from known tissue
12 The inverse problem of tomography is to reconstruct the interior from X-ray data ?????????
13 We write the reconstruction problem in matrix form and assume Gaussian noise m 3 m 2 m 1 f 1 f 4 f 7 f 2 f 5 f 8 f 3 f 6 f 9 f 2 f 3 m 1 f 4 m 2 f = f 5, m = m 3 f 6 m 4, f 7 m 5 f 8 m 6 f 9 Our measurement model is m = Af + ε with independently distributed Gaussian noise (white noise) with standard deviation σ > 0. m 4 m 5 m 6 f 1 m = Af
14 Outline Sparse sampling and tomography Total variation regularization Discretization-invariance Besov space regularization Parameter choice: the S-curve method
15 We reconstruct the unknown pixel values using total variation regularization and non-negativity We consider the minimizer of the anisotropic TV functional Af m α { L H f 1 + L V f 1 } where L H and L H are horizontal and vertical first-order difference matrices. [Rudin, Osher and Fatemi 1992] Primal-dual algorithms Chambolle, Chan, Chen, Esser, Golub, Mulet, Nesterov, Zhang Thresholding Candès, Chambolle, Chaux, Combettes, Daubechies, Defrise, DeMol, Donoho, Pesquet, Starck, Teschke, Vese, Wajs Bregman iteration Cai, Burger, Darbon, Dong, Goldfarb, Mao, Osher, Shen, Xu, Yin, Zhang Splitting approaches Chan, Esser, Fornasier, Goldstein, Langer, Osher, Schönlieb, Setzer, Wajs Nonlocal TV Bertozzi, Bresson, Burger, Chan, Lou, Osher, Zhang
16 We use quadratic programming (QP) to compute the non-negative minimizer of the TV functional The minimizer of the functional { } arg min Af m α L H f 1 + α L V f 1 f R n + can be transformed into the standard form { 1 2 zt Qz + c T z arg min z R 5n }, z 0, Ez = b, where Q is symmetric and E implements equality constraints. Large-scale primal-dual interior point QP method was developed in Kolehmainen, Lassas, Niinimäki & S (2012) and Hämäläinen, Kallonen, Kolehmainen, Lassas, Niinimäki & S (2013).
17 Reduction to arg min z R 5n { 1 2 zt Qz + c T z } Denote horizontal and vertical differences by where u ± H, u ± V arg min f R n + L H f = u + H u H and L V f = u + V u V, 0. TV minimization is now { f T A T Af 2f T A T m + α1 T (u + H + u H + u + V + u V ) }, where 1 R n is vector of all ones. Further, we denote 1 f A T A u + σ 2 H z = u H u +, Q = V , c = u V A T m α1 α1 α1 α1.
18 We collected X-ray projection data of a walnut from 1200 directions The data was collected by Keijo Hämäläinen and Aki Kallonen at University of Helsinki.
19 This is the reconstruction using all 1200 projections and filtered back-projection
20 When only few projection angles are available, TV regularization performs better than FBP 23 angles 15 angles 10 angles FBP TV These images were computed by Kati Niinimäki.
21 Outline Sparse sampling and tomography Total variation regularization Discretization-invariance Besov space regularization Parameter choice: the S-curve method
22 X-ray tomography: Continuum model
23 X-ray tomography: Practical measurement model m R k k = 4
24 X-ray tomography: Computational model m = Af m R k f R n k = 4 n = 48
25 The numbers k and n are independent m = Af m R k f R n k = 8 n = 48
26 The numbers k and n are independent m = Af m R k f R n k = 8 n = 156
27 The numbers k and n are independent m = Af m R k f R n k = 8 n = 440
28 The numbers k and n are independent m = Af m R k f R n k = 16 n = 440
29 The numbers k and n are independent m = Af m R k f R n k = 24 n = 440
30 Bayesian inversion using total variation prior does not have a well-defined continuous limit Theorem (Lassas and S 2004) Total variation prior is not discretization-invariant. Sketch of proof: Apply a variant of the central limit theorem to the independent, identically distributed random consecutive differences. New numerical experiments are reported in Kolehmainen, Lassas, Niinimäki and S (2012) and Lucka (2012).
31 Bayesian inversion using Besov space priors has a well-defined continuous limit Theorem (Lassas, Saksman and S 2009) Besov space priors are discretization-invariant. Sketch of proof: Construction of well-defined Bayesian inversion theory in infinite-dimensional Besov spaces that allow wavelet bases. Discretizations are achieved by truncating the wavelet expansion. Numerical experiments are reported in Kolehmainen, Lassas, Niinimäki and S (2012). Deterministic Besov space regularization was first introduced in Daubechies, Defrise and De Mol (2004); this applies to Bayesian MAP estimates.
32 Outline Sparse sampling and tomography Total variation regularization Discretization-invariance Besov space regularization Parameter choice: the S-curve method
33 The wavelet transform divides an image into three types of details at different scales
34 This is how the wavelet decomposition is defined for two-dimensional periodic images We can represent periodic functions by wavelet expansion f (x, y) = 2 J 0 1 k 1 =0 2 J 0 1 k 2 =0 J 1 c J0 k φ J0, k (x, y)+ 3 2 j 1 2 j 1 j=j 0 l=1 k 1 =0 k 2 =0 w j kl ψ l j, (x, y), k where k = (k 1, k 2 ) and the coefficients c J0 k and w j kl are defined by c J0 k = f, φ J0 k = f (x, y)φ j k (x, y)dxdy, T 2 w j kl = f, ψ l j = f (x, y)ψ k T l 2 j (x, y)dxdy. k
35 Besov space norm can be defined in terms of the wavelet coefficients A function f belongs to B s pq(t 2 ) if and only if the following expression is finite: 2 J 0 1 k 1 =0 2 J 0 1 k 2 =0 c J0 k p + 1 p j=j 0 2 jq(s+1 2 p ) We focus on the case p = 1 = q and s = 1: 3 2 j 1 2 j 1 l=1 k 1 =0 k 2 =0 w j kl p q p 1 q f B 1 11 (T 2 ) = k c J0 k + j,l, k w j kl.
36 We promote sparsity by minimizing a sum of l 2 and l 1 norms In the case of Besov space penalty we minimize f MAP = arg min f R n + { 1 2σ 2 Af m 2 l 2 + α f B 1 11 (T 2 ) where ψ j are the wavelet basis functions. The above kind of minimization task has been studied in Chambolle, DeVore, Lee & Lucier 1998, Daubechies, Defrise & De Mol 2004, Candès, Romberg & Tao 2006, Klann, Maaß & Ramlau 2006, Grasmair, Haltmeier & Scherzer 2008, Kolehmainen, Lassas, Niinimäki & S }, We use constrained quadratic programming for the minimization.
37 History of multiresolution tomography 1994 Olson and DeStefano 1994 Sahiner and Yagle 1995 Delaney and Bresler 1996 Berenstein and Walnut 1996 Bhatia, Karl and Willsky 1997 Rashid-Farrokhi, Liu, Berenstein and Walnut 1997 Zhao, Welland and Wang 1998 Louis, Maaß and Rieder 1999 Candès and Donoho 1999 Madych 2000 Smith and Adhami 2000 Bonnet, Peyrin, Turjman and Prost 2002 Das and Sastry 2002 Frese, Bouman and Sauer 2004 Soleski and Walter 2004 Zhong, Ning and Conover 2006 Sastry and Das 2006 Soleski and Walter 2006 Lee and Lucier 2006 Rantala, Vänskä, Järvenpää, Kalke, Lassas, Moberg and S 2007 Niinimäki, S and Kolehmainen 2008 Soussen and Idier 2009 Vänskä, Lassas and S 2011 Klann, Ramlau and Reichel 2011 Terzija and McCann 2011 Frikel
38 Besov regularized sparse-data reconstructions compared to filtered back-projection
39 Outline Sparse sampling and tomography Total variation regularization Discretization-invariance Besov space regularization Parameter choice: the S-curve method
40 We took photographs of walnuts cut in half These photos are used for estimating the expected number of nonzero wavelet coefficients in a two-dimensional tomographic reconstruction. Special thanks go to Esa Niemi for his careful job in sawing the walnuts.
41 The S-curve method determines a regularization parameter value giving the right sparsity level Number of nonzero wavelet coefficients Regularization parameter α Kolehmainen, Lassas, Niinimäki & S 2012
42 Something needs to be done to deal with the erratic behaviour of the raw S-curve
43 Reconstruction from 30 projections using interpolated S-curve α = 0.019
44 Interpolated S-curve method for TV These unpublished images were computed by Kati Niinimäki.
45 All Matlab codes freely available on a website! Part I: Linear Inverse Problems 1 Introduction 2 Naïve reconstructions and inverse crimes 3 Ill-Posedness in Inverse Problems 4 Truncated singular value decomposition 5 Tikhonov regularization 6 Total variation regularization 7 Besov space regularization using wavelets 8 Discretization-invariance 9 Practical X-ray tomography with limited data 10 Projects Part II: Nonlinear Inverse Problems 11 Nonlinear inversion 12 Electrical impedance tomography 13 Simulation of noisy EIT data 14 Complex geometrical optics solutions 15 A regularized D-bar method for direct EIT 16 Other direct solution methods for EIT 17 Projects
46 Inverse Days, Dec 11 13, 2013, Inari, Finland Organizers: Maarten de Hoop Matti Lassas Markku Lehtinen Lassi Roininen S. S. Gunther Uhlmann
47 Thank you for your attention! Preprints available at
Discretization-invariant Bayesian inversion and Besov space priors. Samuli Siltanen RICAM Tampere University of Technology
Discretization-invariant Bayesian inversion and Besov space priors Samuli Siltanen RICAM 28.10.2008 Tampere University of Technology http://math.tkk.fi/inverse-coe/ This is a joint work with Matti Lassas
More informationRecent Advances in Bayesian Inference for Inverse Problems
Recent Advances in Bayesian Inference for Inverse Problems Felix Lucka University College London, UK f.lucka@ucl.ac.uk Applied Inverse Problems Helsinki, May 25, 2015 Bayesian Inference for Inverse Problems
More informationarxiv: v2 [math.na] 11 Aug 2015
MULTI-RESOLUTION PARAMETER CHOICE METHOD FOR TOTAL VARIATION REGULARIZED TOMOGRAPHY KATI NIINIMÄKI, MATTI LASSAS, KEIJO HÄMÄLÄINEN, AKI KALLONEN, VILLE KOLEHMAINEN, ESA NIEMI, AND SAMULI SILTANEN arxiv:1407.2386v2
More informationStrengthened Sobolev inequalities for a random subspace of functions
Strengthened Sobolev inequalities for a random subspace of functions Rachel Ward University of Texas at Austin April 2013 2 Discrete Sobolev inequalities Proposition (Sobolev inequality for discrete images)
More informationMCMC Sampling for Bayesian Inference using L1-type Priors
MÜNSTER MCMC Sampling for Bayesian Inference using L1-type Priors (what I do whenever the ill-posedness of EEG/MEG is just not frustrating enough!) AG Imaging Seminar Felix Lucka 26.06.2012 , MÜNSTER Sampling
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 informationc 2011 International Press Vol. 18, No. 1, pp , March DENNIS TREDE
METHODS AND APPLICATIONS OF ANALYSIS. c 2011 International Press Vol. 18, No. 1, pp. 105 110, March 2011 007 EXACT SUPPORT RECOVERY FOR LINEAR INVERSE PROBLEMS WITH SPARSITY CONSTRAINTS DENNIS TREDE Abstract.
More informationEnhanced Compressive Sensing and More
Enhanced Compressive Sensing and More Yin Zhang Department of Computational and Applied Mathematics Rice University, Houston, Texas, U.S.A. Nonlinear Approximation Techniques Using L1 Texas A & M University
More informationSolving DC Programs that Promote Group 1-Sparsity
Solving DC Programs that Promote Group 1-Sparsity Ernie Esser Contains joint work with Xiaoqun Zhang, Yifei Lou and Jack Xin SIAM Conference on Imaging Science Hong Kong Baptist University May 14 2014
More informationSignal Restoration with Overcomplete Wavelet Transforms: Comparison of Analysis and Synthesis Priors
Signal Restoration with Overcomplete Wavelet Transforms: Comparison of Analysis and Synthesis Priors Ivan W. Selesnick a and Mário A. T. Figueiredo b a Polytechnic Institute of New York University, Brooklyn,
More informationCompressive Sensing Theory and L1-Related Optimization Algorithms
Compressive Sensing Theory and L1-Related Optimization Algorithms Yin Zhang Department of Computational and Applied Mathematics Rice University, Houston, Texas, USA CAAM Colloquium January 26, 2009 Outline:
More informationConvergence rates in l 1 -regularization when the basis is not smooth enough
Convergence rates in l 1 -regularization when the basis is not smooth enough Jens Flemming, Markus Hegland November 29, 2013 Abstract Sparsity promoting regularization is an important technique for signal
More informationResolving the White Noise Paradox in the Regularisation of Inverse Problems
1 / 32 Resolving the White Noise Paradox in the Regularisation of Inverse Problems Hanne Kekkonen joint work with Matti Lassas and Samuli Siltanen Department of Mathematics and Statistics University of
More informationSparsity Regularization
Sparsity Regularization Bangti Jin Course Inverse Problems & Imaging 1 / 41 Outline 1 Motivation: sparsity? 2 Mathematical preliminaries 3 l 1 solvers 2 / 41 problem setup finite-dimensional formulation
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 informationFeature Reconstruction in Tomography
Feature Reconstruction in Tomography Alfred K. Louis Institut für Angewandte Mathematik Universität des Saarlandes 66041 Saarbrücken http://www.num.uni-sb.de louis@num.uni-sb.de Wien, July 20, 2009 Images
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 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 informationwissen leben WWU Münster
MÜNSTER Sparsity Constraints in Bayesian Inversion Inverse Days conference in Jyväskylä, Finland. Felix Lucka 19.12.2012 MÜNSTER 2 /41 Sparsity Constraints in Inverse Problems Current trend in high dimensional
More informationProximal tools for image reconstruction in dynamic Positron Emission Tomography
Proximal tools for image reconstruction in dynamic Positron Emission Tomography Nelly Pustelnik 1 joint work with Caroline Chaux 2, Jean-Christophe Pesquet 3, and Claude Comtat 4 1 Laboratoire de Physique,
More informationCompressive Sensing (CS)
Compressive Sensing (CS) Luminita Vese & Ming Yan lvese@math.ucla.edu yanm@math.ucla.edu Department of Mathematics University of California, Los Angeles The UCLA Advanced Neuroimaging Summer Program (2014)
More informationIterative Methods for Inverse & Ill-posed Problems
Iterative Methods for Inverse & Ill-posed Problems Gerd Teschke Konrad Zuse Institute Research Group: Inverse Problems in Science and Technology http://www.zib.de/ag InverseProblems ESI, Wien, December
More information(4D) Variational Models Preserving Sharp Edges. Martin Burger. Institute for Computational and Applied Mathematics
(4D) Variational Models Preserving Sharp Edges Institute for Computational and Applied Mathematics Intensity (cnt) Mathematical Imaging Workgroup @WWU 2 0.65 0.60 DNA Akrosom Flagellum Glass 0.55 0.50
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 informationPrimal-dual fixed point algorithms for separable minimization problems and their applications in imaging
1/38 Primal-dual fixed point algorithms for separable minimization problems and their applications in imaging Xiaoqun Zhang Department of Mathematics and Institute of Natural Sciences Shanghai Jiao Tong
More informationLINEARIZED BREGMAN ITERATIONS FOR COMPRESSED SENSING
LINEARIZED BREGMAN ITERATIONS FOR COMPRESSED SENSING JIAN-FENG CAI, STANLEY OSHER, AND ZUOWEI SHEN Abstract. Finding a solution of a linear equation Au = f with various minimization properties arises from
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 informationCompressed Sensing with Quantized Measurements
Compressed Sensing with Quantized Measurements Argyrios Zymnis, Stephen Boyd, and Emmanuel Candès April 28, 2009 Abstract We consider the problem of estimating a sparse signal from a set of quantized,
More informationarxiv: v1 [math.st] 19 Mar 2016
arxiv:63.635v [math.st] 9 Mar 6 Cauchy difference priors for edge-preserving Bayesian inversion with an application to X-ray tomography Markku Markkanen, Lassi Roininen, Janne M J Huttunen 3 and Sari Lasanen
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 informationMathematisches Forschungsinstitut Oberwolfach. Computational Inverse Problems
Mathematisches Forschungsinstitut Oberwolfach Report No. 51/2012 DOI: 10.4171/OWR/2012/51 Computational Inverse Problems Organised by Habib Ammari, Paris Liliana Borcea, Houston Thorsten Hohage, Göttingen
More informationLarge-Scale L1-Related Minimization in Compressive Sensing and Beyond
Large-Scale L1-Related Minimization in Compressive Sensing and Beyond Yin Zhang Department of Computational and Applied Mathematics Rice University, Houston, Texas, U.S.A. Arizona State University March
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 informationInverse problems and sparse models (6/6) Rémi Gribonval INRIA Rennes - Bretagne Atlantique, France.
Inverse problems and sparse models (6/6) Rémi Gribonval INRIA Rennes - Bretagne Atlantique, France remi.gribonval@inria.fr Overview of the course Introduction sparsity & data compression inverse problems
More informationAdaptive Primal Dual Optimization for Image Processing and Learning
Adaptive Primal Dual Optimization for Image Processing and Learning Tom Goldstein Rice University tag7@rice.edu Ernie Esser University of British Columbia eesser@eos.ubc.ca Richard Baraniuk Rice University
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 informationOptimization Algorithms for Compressed Sensing
Optimization Algorithms for Compressed Sensing Stephen Wright University of Wisconsin-Madison SIAM Gator Student Conference, Gainesville, March 2009 Stephen Wright (UW-Madison) Optimization and Compressed
More informationContraction Methods for Convex Optimization and Monotone Variational Inequalities No.11
XI - 1 Contraction Methods for Convex Optimization and Monotone Variational Inequalities No.11 Alternating direction methods of multipliers for separable convex programming Bingsheng He Department of Mathematics
More informationThe convex algebraic geometry of rank minimization
The convex algebraic geometry of rank minimization Pablo A. Parrilo Laboratory for Information and Decision Systems Massachusetts Institute of Technology International Symposium on Mathematical Programming
More informationA Primal-dual Three-operator Splitting Scheme
Noname manuscript No. (will be inserted by the editor) A Primal-dual Three-operator Splitting Scheme Ming Yan Received: date / Accepted: date Abstract In this paper, we propose a new primal-dual algorithm
More informationSparse Optimization: Algorithms and Applications. Formulating Sparse Optimization. Motivation. Stephen Wright. Caltech, 21 April 2007
Sparse Optimization: Algorithms and Applications Stephen Wright 1 Motivation and Introduction 2 Compressed Sensing Algorithms University of Wisconsin-Madison Caltech, 21 April 2007 3 Image Processing +Mario
More informationHierarchical Matérn fields and Cauchy difference priors for Bayesian inversion
Hierarchical Matérn fields and Cauchy difference priors for Bayesian inversion Lassi Roininen Bayesian Computation for High-Dimensional Statistical Models National University of Singapore, 27 August 2018
More informationImage Noise: Detection, Measurement and Removal Techniques. Zhifei Zhang
Image Noise: Detection, Measurement and Removal Techniques Zhifei Zhang Outline Noise measurement Filter-based Block-based Wavelet-based Noise removal Spatial domain Transform domain Non-local methods
More informationUser s Guide to Compressive Sensing
A Walter B. Richardson, Jr. University of Texas at San Antonio walter.richardson@utsa.edu Engineering November 18, 2011 Abstract During the past decade, Donoho, Candes, and others have developed a framework
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 informationRecovering overcomplete sparse representations from structured sensing
Recovering overcomplete sparse representations from structured sensing Deanna Needell Claremont McKenna College Feb. 2015 Support: Alfred P. Sloan Foundation and NSF CAREER #1348721. Joint work with Felix
More informationStable Signal Recovery from Incomplete and Inaccurate Measurements
Stable Signal Recovery from Incomplete and Inaccurate Measurements EMMANUEL J. CANDÈS California Institute of Technology JUSTIN K. ROMBERG California Institute of Technology AND TERENCE TAO University
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 informationApproximate Message Passing for Bilinear Models
Approximate Message Passing for Bilinear Models Volkan Cevher Laboratory for Informa4on and Inference Systems LIONS / EPFL h,p://lions.epfl.ch & Idiap Research Ins=tute joint work with Mitra Fatemi @Idiap
More informationA FISTA-like scheme to accelerate GISTA?
A FISTA-like scheme to accelerate GISTA? C. Cloquet 1, I. Loris 2, C. Verhoeven 2 and M. Defrise 1 1 Dept. of Nuclear Medicine, Vrije Universiteit Brussel 2 Dept. of Mathematics, Université Libre de Bruxelles
More informationBias-free Sparse Regression with Guaranteed Consistency
Bias-free Sparse Regression with Guaranteed Consistency Wotao Yin (UCLA Math) joint with: Stanley Osher, Ming Yan (UCLA) Feng Ruan, Jiechao Xiong, Yuan Yao (Peking U) UC Riverside, STATS Department March
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 informationRegularizing inverse problems using sparsity-based signal models
Regularizing inverse problems using sparsity-based signal models Jeffrey A. Fessler William L. Root Professor of EECS EECS Dept., BME Dept., Dept. of Radiology University of Michigan http://web.eecs.umich.edu/
More informationRecent developments on sparse representation
Recent developments on sparse representation Zeng Tieyong Department of Mathematics, Hong Kong Baptist University Email: zeng@hkbu.edu.hk Hong Kong Baptist University Dec. 8, 2008 First Previous Next Last
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 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 informationVariational image restoration by means of wavelets: Simultaneous decomposition, deblurring, and denoising
Appl. Comput. Harmon. Anal. 19 (2005) 1 16 www.elsevier.com/locate/acha Variational image restoration by means of wavelets: Simultaneous decomposition, deblurring, and denoising I. Daubechies a,g.teschke
More informationCalderón s inverse problem in 2D Electrical Impedance Tomography
Calderón s inverse problem in 2D Electrical Impedance Tomography Kari Astala (University of Helsinki) Joint work with: Matti Lassas, Lassi Päivärinta, Samuli Siltanen, Jennifer Mueller and Alan Perämäki
More informationInterpolation procedure in filtered backprojection algorithm for the limited-angle tomography
Interpolation procedure in filtered backprojection algorithm for the limited-angle tomography Aleksander Denisiuk 1 Abstract A new data completion procedure for the limited-angle tomography is proposed.
More informationGauge optimization and duality
1 / 54 Gauge optimization and duality Junfeng Yang Department of Mathematics Nanjing University Joint with Shiqian Ma, CUHK September, 2015 2 / 54 Outline Introduction Duality Lagrange duality Fenchel
More informationAdaptive Corrected Procedure for TVL1 Image Deblurring under Impulsive Noise
Adaptive Corrected Procedure for TVL1 Image Deblurring under Impulsive Noise Minru Bai(x T) College of Mathematics and Econometrics Hunan University Joint work with Xiongjun Zhang, Qianqian Shao June 30,
More informationAN INTRODUCTION TO COMPRESSIVE SENSING
AN INTRODUCTION TO COMPRESSIVE SENSING Rodrigo B. Platte School of Mathematical and Statistical Sciences APM/EEE598 Reverse Engineering of Complex Dynamical Networks OUTLINE 1 INTRODUCTION 2 INCOHERENCE
More informationCompressed Sensing and Linear Codes over Real Numbers
Compressed Sensing and Linear Codes over Real Numbers Henry D. Pfister (joint with Fan Zhang) Texas A&M University College Station Information Theory and Applications Workshop UC San Diego January 31st,
More informationCOMPRESSIVE OPTICAL DEFLECTOMETRIC TOMOGRAPHY
COMPRESSIVE OPTICAL DEFLECTOMETRIC TOMOGRAPHY Adriana González ICTEAM/UCL March 26th, 214 1 ISPGroup - ICTEAM - UCL http://sites.uclouvain.be/ispgroup Université catholique de Louvain, Louvain-la-Neuve,
More informationMultiscale Geometric Analysis: Thoughts and Applications (a summary)
Multiscale Geometric Analysis: Thoughts and Applications (a summary) Anestis Antoniadis, University Joseph Fourier Assimage 2005,Chamrousse, February 2005 Classical Multiscale Analysis Wavelets: Enormous
More informationSparsity in Underdetermined Systems
Sparsity in Underdetermined Systems Department of Statistics Stanford University August 19, 2005 Classical Linear Regression Problem X n y p n 1 > Given predictors and response, y Xβ ε = + ε N( 0, σ 2
More informationConvergence rates for Morozov s Discrepancy Principle using Variational Inequalities
Convergence rates for Morozov s Discrepancy Principle using Variational Inequalities Stephan W Anzengruber Ronny Ramlau Abstract We derive convergence rates for Tikhonov-type regularization with conve
More informationMorozov s discrepancy principle for Tikhonov-type functionals with non-linear operators
Morozov s discrepancy principle for Tikhonov-type functionals with non-linear operators Stephan W Anzengruber 1 and Ronny Ramlau 1,2 1 Johann Radon Institute for Computational and Applied Mathematics,
More informationSparse linear models
Sparse linear models Optimization-Based Data Analysis http://www.cims.nyu.edu/~cfgranda/pages/obda_spring16 Carlos Fernandez-Granda 2/22/2016 Introduction Linear transforms Frequency representation Short-time
More informationStatistical and Computational Inverse Problems with Applications Part 2: Introduction to inverse problems and example applications
Statistical and Computational Inverse Problems with Applications Part 2: Introduction to inverse problems and example applications Aku Seppänen Inverse Problems Group Department of Applied Physics University
More informationA REWEIGHTED l 2 METHOD FOR IMAGE RESTORATION WITH POISSON AND MIXED POISSON-GAUSSIAN NOISE. Jia Li. Zuowei Shen. Rujie Yin.
Volume X, No. 0X, 0xx, X XX doi:10.3934/ipi.xx.xx.xx A REWEIGHTED l METHOD FOR IMAGE RESTORATION WITH POISSON AND MIXED POISSON-GAUSSIAN NOISE JIA LI AND ZUOWEI SHEN AND RUJIE YIN AND XIAOQUN ZHANG Jia
More informationACCELERATED LINEARIZED BREGMAN METHOD. June 21, Introduction. In this paper, we are interested in the following optimization problem.
ACCELERATED LINEARIZED BREGMAN METHOD BO HUANG, SHIQIAN MA, AND DONALD GOLDFARB June 21, 2011 Abstract. In this paper, we propose and analyze an accelerated linearized Bregman (A) method for solving the
More informationThe impact of a curious type of smoothness conditions on convergence rates in l 1 -regularization
The impact of a curious type of smoothness conditions on convergence rates in l 1 -regularization Radu Ioan Boț and Bernd Hofmann March 1, 2013 Abstract Tikhonov-type regularization of linear and nonlinear
More informationSparse & Redundant Signal Representation, and its Role in Image Processing
Sparse & Redundant Signal Representation, and its Role in Michael Elad The CS Department The Technion Israel Institute of technology Haifa 3000, Israel Wave 006 Wavelet and Applications Ecole Polytechnique
More informationLecture Notes 9: Constrained Optimization
Optimization-based data analysis Fall 017 Lecture Notes 9: Constrained Optimization 1 Compressed sensing 1.1 Underdetermined linear inverse problems Linear inverse problems model measurements of the form
More informationInverse problems in statistics
Inverse problems in statistics Laurent Cavalier (Université Aix-Marseille 1, France) Yale, May 2 2011 p. 1/35 Introduction There exist many fields where inverse problems appear Astronomy (Hubble satellite).
More informationNear Optimal Signal Recovery from Random Projections
1 Near Optimal Signal Recovery from Random Projections Emmanuel Candès, California Institute of Technology Multiscale Geometric Analysis in High Dimensions: Workshop # 2 IPAM, UCLA, October 2004 Collaborators:
More informationBindel, Fall 2016 Matrix Computations (CS 6210) Notes for
1 A cautionary tale Notes for 2016-10-05 You have been dropped on a desert island with a laptop with a magic battery of infinite life, a MATLAB license, and a complete lack of knowledge of basic geometry.
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 informationReconstructing conductivities with boundary corrected D-bar method
Reconstructing conductivities with boundary corrected D-bar method Janne Tamminen June 24, 2011 Short introduction to EIT The Boundary correction procedure The D-bar method Simulation of measurement data,
More informationCompressed Sensing and Neural Networks
and Jan Vybíral (Charles University & Czech Technical University Prague, Czech Republic) NOMAD Summer Berlin, September 25-29, 2017 1 / 31 Outline Lasso & Introduction Notation Training the network Applications
More informationSPARSITY-ASSISTED SIGNAL SMOOTHING (REVISITED) Ivan Selesnick
SPARSITY-ASSISTED SIGNAL SMOOTHING (REVISITED) Ivan Selesnick Electrical and Computer Engineering Tandon School of Engineering, New York University Brooklyn, New York, USA ABSTRACT This paper proposes
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 informationInverse problems and sparse models (1/2) Rémi Gribonval INRIA Rennes - Bretagne Atlantique, France
Inverse problems and sparse models (1/2) Rémi Gribonval INRIA Rennes - Bretagne Atlantique, France remi.gribonval@inria.fr Structure of the tutorial Session 1: Introduction to inverse problems & sparse
More informationarxiv: v2 [math.fa] 18 Jul 2008
Optimal Convergence Rates for Tikhonov Regularization in Besov Scales arxiv:0806095v [mathfa] 8 Jul 008 Introduction D A Lorenz and D Trede Zentrum für Technomathematik University of Bremen D 8334 Bremen
More informationA primal-dual fixed point algorithm for multi-block convex minimization *
Journal of Computational Mathematics Vol.xx, No.x, 201x, 1 16. http://www.global-sci.org/jcm doi:?? A primal-dual fixed point algorithm for multi-block convex minimization * Peijun Chen School of Mathematical
More informationA General Framework for a Class of Primal-Dual Algorithms for TV Minimization
A General Framework for a Class of Primal-Dual Algorithms for TV Minimization Ernie Esser UCLA 1 Outline A Model Convex Minimization Problem Main Idea Behind the Primal Dual Hybrid Gradient (PDHG) Method
More informationNoisy Signal Recovery via Iterative Reweighted L1-Minimization
Noisy Signal Recovery via Iterative Reweighted L1-Minimization Deanna Needell UC Davis / Stanford University Asilomar SSC, November 2009 Problem Background Setup 1 Suppose x is an unknown signal in R d.
More informationSparse Approximation: from Image Restoration to High Dimensional Classification
Sparse Approximation: from Image Restoration to High Dimensional Classification Bin Dong Beijing International Center for Mathematical Research Beijing Institute of Big Data Research Peking University
More informationCompressed Sensing in Astronomy
Compressed Sensing in Astronomy J.-L. Starck CEA, IRFU, Service d'astrophysique, France jstarck@cea.fr http://jstarck.free.fr Collaborators: J. Bobin, CEA. Introduction: Compressed Sensing (CS) Sparse
More informationA Primal-Dual Method for Total Variation-Based. Wavelet Domain Inpainting
A Primal-Dual Method for Total Variation-Based 1 Wavelet Domain Inpainting You-Wei Wen, Raymond H. Chan, Andy M. Yip Abstract Loss of information in a wavelet domain can occur during storage or transmission
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 informationMinimizing Isotropic Total Variation without Subiterations
MITSUBISHI ELECTRIC RESEARCH LABORATORIES http://www.merl.com Minimizing Isotropic Total Variation without Subiterations Kamilov, U. S. TR206-09 August 206 Abstract Total variation (TV) is one of the most
More informationMinimizing the Difference of L 1 and L 2 Norms with Applications
1/36 Minimizing the Difference of L 1 and L 2 Norms with Department of Mathematical Sciences University of Texas Dallas May 31, 2017 Partially supported by NSF DMS 1522786 2/36 Outline 1 A nonconvex approach:
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 informationLearning MMSE Optimal Thresholds for FISTA
MITSUBISHI ELECTRIC RESEARCH LABORATORIES http://www.merl.com Learning MMSE Optimal Thresholds for FISTA Kamilov, U.; Mansour, H. TR2016-111 August 2016 Abstract Fast iterative shrinkage/thresholding algorithm
More informationRegularized Multiplicative Algorithms for Nonnegative Matrix Factorization
Regularized Multiplicative Algorithms for Nonnegative Matrix Factorization Christine De Mol (joint work with Loïc Lecharlier) Université Libre de Bruxelles Dept Math. and ECARES MAHI 2013 Workshop Methodological
More informationA Linearly Convergent First-order Algorithm for Total Variation Minimization in Image Processing
A Linearly Convergent First-order Algorithm for Total Variation Minimization in Image Processing Cong D. Dang Kaiyu Dai Guanghui Lan October 9, 0 Abstract We introduce a new formulation for total variation
More informationA Unified Primal-Dual Algorithm Framework Based on Bregman Iteration
J Sci Comput (2011) 46: 20 46 DOI 10.1007/s10915-010-9408-8 A Unified Primal-Dual Algorithm Framework Based on Bregman Iteration Xiaoqun Zhang Martin Burger Stanley Osher Received: 23 November 2009 / Revised:
More informationBasic Formulation. minimize. subject to Ax + By = b. functions f, h can take the extended value, can be nonsmooth. x,y
Basic Formulation minimize x,y f(x) + h(y) subject to Ax + By = b functions f, h can take the extended value, can be nonsmooth ADMM Define the augmented Lagrangian L β (x, y; w) = f(x) + h(y) + w, Ax +
More information