400 likes | 414 Views
This talk provides an overview of numerical methods for large-scale inverse problems, including regularization algorithms, Krylov subspace methods, preconditioning, and signal subspaces. Examples from various fields are discussed.
E N D
Large-Scale Methods in Inverse Problems • Per Christian Hansen • Informatics and Mathematical Modelling • Technical University of Denmark • With contributions from: • Michael Jacobsen, Toke Koldborg Jensen - PhD students • Line H. Clemmensen, Iben Kraglund, Kristine Horn,Jesper Pedersen, Marie-Louise H. Rasmussen - Master students Large-Scale Methods in Inverse Problems
Overview of Talk • A survey of numerical methods for large-scale inverse problems • Some examples. • The need for regularization algorithms. • Krylov subspace methods for large-scale problems. • Preconditioning for regularization problems. • Signal subspaces and (semi)norms. • GMRES as a regularization method. • Alternatives to spectral filtering. • Many details are skipped, to get the big picture!!! Large-Scale Methods in Inverse Problems
Related Work • Many people work on similar problems and algorithms: • Åke Björck, Lars Eldén, Tommy Elfving • Martin Hanke, James G. Nagy, Robert Plemmons • Misha E. Kilmer, Dianne P. Oleary • Daniela Calvetti, Lothar Reichel, Brian Lewis • Gene H. Golub, Urs von Matt • Uri Asher, Eldad Haber, Douglas Oldenburg • Jerry Eriksson, Mårten Gullikson, Per-Åke Wedin • Marielba Rojas, Trond Steihaug • Tony Chan, Stanley Osher, Curtis R. Vogel • Jesse Barlow, Raymond Chan, Michael Ng • Recent Matlab software packages: • Restore Tools (Nagy, Palmer, Perrone, 2004) • MOORe Tools (Jacobsen, 2004) • GeoTools (Pedersen, 2005) Large-Scale Methods in Inverse Problems
Inverse Geomagnetic Problems Large-Scale Methods in Inverse Problems
Inverse Acoustic Problems Oticon/ Rhinometrics Large-Scale Methods in Inverse Problems
Image Restoration Problems blurring deblurring Io (moon of Saturn) You cannot depend on your eyes when your imagination is out of focus – Mark Twain Large-Scale Methods in Inverse Problems
Model Problem and Discretization Vertical component of magnetic field from a dipole Large-Scale Methods in Inverse Problems
The Need for Regularization Regularization: keep the “good” SVD components and discard the noisy ones! Large-Scale Methods in Inverse Problems
Regularization – TSVD & Tikhonov Large-Scale Methods in Inverse Problems
Singular Vectors (Always) Oscillate Large-Scale Methods in Inverse Problems
Large-Scale Aspects (the easy case) Large-Scale Methods in Inverse Problems
Large-Scale Aspects (the real problems) Toeplitz matrix-vector multiplication flop count. Large-Scale Methods in Inverse Problems
Large-Scale Tikhonov Regularization Large-Scale Methods in Inverse Problems
Difficulties and Remedies I Large-Scale Methods in Inverse Problems
Difficulties and Remedies II Large-Scale Methods in Inverse Problems
The Art of Preconditioning Large-Scale Methods in Inverse Problems
Explicit Subspace Preconditiong Large-Scale Methods in Inverse Problems
Krylov Signal Subspaces Smiley Crater, Mars Large-Scale Methods in Inverse Problems
Pros and Cons of Regularizing Iterations Large-Scale Methods in Inverse Problems
Projection, then Regularization Large-Scale Methods in Inverse Problems
Bounds on “Everything” Large-Scale Methods in Inverse Problems
A Dilemma With Projection + Regular. Large-Scale Methods in Inverse Problems
Better Basis Vectors! Large-Scale Methods in Inverse Problems
Considerations in 2D … … Large-Scale Methods in Inverse Problems
Good Seminorms for 2D Problems Large-Scale Methods in Inverse Problems
Seminorms and Regularizing Iterations Large-Scale Methods in Inverse Problems
Krylov Implementation Large-Scale Methods in Inverse Problems
Avoiding the Transpose: GMRES Large-Scale Methods in Inverse Problems
GMRES and CGLS Basis Vectors Large-Scale Methods in Inverse Problems
CGLS and GMRES Solutions Large-Scale Methods in Inverse Problems
The “Freckles’’ DCT spectrum spatial domain Large-Scale Methods in Inverse Problems
Preconditioning for GMRES Large-Scale Methods in Inverse Problems
A New and Better Approach Large-Scale Methods in Inverse Problems
(P)CGLS and (P)GMRES Large-Scale Methods in Inverse Problems
Away From 2-Norms Io (moon of Saturn) q = 1.1 q = 2 Large-Scale Methods in Inverse Problems
Functionals Defined on Sols. to DIP Large-Scale Methods in Inverse Problems
Large-Scale Algorithm MLFIP Large-Scale Methods in Inverse Problems
Confidence Invervals with MLFIP Large-Scale Methods in Inverse Problems
Many Topics Not Covered … • Algorithms for other norms (p and q≠ 2). • In particular, total variation (TV). • Nonnegativity constraints. • General linear inequality constraints. • Compression of dense coefficient matrix A. • Color images (and color TV). • Implementation aspects and software. • The choice the of regularization parameter. Large-Scale Methods in Inverse Problems
“Conclusions and Further Work” • I hesitate to give any conclusion – • the work is ongoing; • there are many open problems, • lots of challenges (mathematical and numerical), • and a multitude of practical problems waiting to be solved. Large-Scale Methods in Inverse Problems