Abstract
Large-scale inverse problems arise in a variety of significant applications in image processing, and efficient regularization methods are needed to compute meaningful solutions. This chapter surveys three common mathematical models including a linear, a separable nonlinear, and a general nonlinear model. Techniques for regularization and large-scale implementations are considered, with particular focus on algorithms and computations that can exploit structure in the problem. Examples from image deconvolution, multi-frame blind deconvolution, and tomosynthesis illustrate the potential of these algorithms. Much progress has been made in the field of large-scale inverse problems, but many challenges still remain for future research.
Access this chapter
Tax calculation will be finalised at checkout
Purchases are for personal use only
Similar content being viewed by others
References and Further Reading
Andrews HC, Hunt BR (1977) Digital image restoration. Prentice-Hall, Englewood Cliffs
Bachmayr M, Burger M (2009) Iterative total variation schemes for nonlinear inverse problems. Inverse Prob 25:105004
Bardsley JM (2008) An efficient computational method for total variation-penalized Poisson likelihood estimation. Inverse Prob Imaging 2(2):167–185
Bardsley JM (2008) Stopping rules for a nonnegatively constrained iterative method for illposed Poisson imaging problems. BIT 48(4):651–664
Bardsley JM, Vogel CR (2003) A nonnegatively constrained convex programming method for image reconstruction. SIAM J Sci Comput 25(4):1326–1343
Barzilai J, Borwein JM (1988) Two-point step size gradient methods. IMA J Numer Anal 8(1):141–148
Björck Å (1988) A bidiagonalization algorithm for solving large and sparse ill-posed systems of linear equations. BIT, 28(3):659–670
Björck Å (1996) Numerical methods for least squares problems. SIAM, Philadelphia
Björck Å, Grimme E, van Dooren P (1994) An implicit shift bidiagonalization algorithm for ill-posed systems. BIT 34(4):510–534
Brakhage H (1987) On ill-posed problems and the method of conjugate gradients. In: Engl HW, Groetsch CW (eds) Inverse and ill-posed problems. Academic, Boston, pp 165–175
Calvetti D, Reichel L (2003) Tikhonov regularization of large linear problems. BIT 43(2):263–283
Calvetti D, Somersalo E (2007) Introduction to Bayesian scientific computing. Springer, New York
Candès EJ, Romberg JK, Tao T (2006) Stable signal recovery from incomplete and inaccurate measurements. Commun Pure Appl Math 59(8):1207–1223
Carasso AS (2001) Direct blind deconvolution. SIAM J Appl Math 61(6):1980–2007
Chadan K, Colton D, Päivärinta L, Rundell W (1997) An introduction to inverse scattering and inverse spectral problems. SIAM, Philadelphia
Chan TF, Shen J (2005) Image processing and analysis: variational, PDE, wavelet, and stochastic methods. SIAM, Philadelphia
Cheney M, Borden B (2009) Fundamentals of radar imaging. SIAM, Philadelphia
Chung J, Haber E, Nagy J (2006) Numerical methods for coupled super-resolution. Inverse Prob 22(4):1261–1272
Chung J, Nagy J (2010) An efficient iterative approach for large-scale separable nonlinear inverse problems. SIAM J Sci Comput 31(6):4654–4674
Chung J, Nagy J, Sechopoulos I (2010) Numerical algorithms for polyenergetic digital breast tomosynthesis reconstruction. SIAM J Imaging Sci 3(1):133–152
Chung J, Nagy JG, O’Leary DP (2008) A weighted GCV method for Lanczos hybrid regularization. Elec Trans Numer Anal 28: 149–167
Chung J, Sternberg P, Yang C (2010) High performance 3-d image reconstruction for molecular structure determination. Int J High Perform Comput Appl 24(2):117–135
De Man B, Nuyts J, Dupont P, Marchal G, Suetens P (2001) An iterative maximumlikelihood polychromatic algorithm for CT. IEEE Trans Med Imaging 20(10):999–1008
Diaspro A, Corosu M, Ramoino P, Robello M (1999) Two-photon excitation imaging based on a compact scanning head. IEEE Eng Med Biol 18(5):18–30
Dobbins JT III, Godfrey DJ (2003) Digital X-ray tomosynthesis: current state of the art and clinical potential. Phys Med Biol 48(19):R65–R106
Easley GR, Healy DM, Berenstein CA (2009) Image deconvolution using a general ridgelet and curvelet domain. SIAM J Imaging Sci 2(1):253–283
Elad M, Feuer A (1997) Restoration of a single superresolution image from several blurred, noisy, and undersampled measured images. IEEE Trans Image Process 6(12):1646–1658
Engl HW, Hanke M, Neubauer A (2000) Regularization of inverse problems. Kluwer, Dordrecht
Engl HW, Kügler P (2005) Nonlinear inverse problems: theoretical aspects and some industrial applications. In: Capasso V, Périaux J (eds) Multidisciplinary methods for analysis optimization and control of complex systems. Springer, Berlin, pp 3–48
Engl HW, Kunisch K, Neubauer A (1989) Convergence rates for Tikhonov regularisation of nonlinear ill-posed problems. Inverse Prob 5(4):523–540
Engl HW, Louis AK, Rundell W (eds) (1996) Inverse problems in geophysical applications. SIAM, Philadelphia
Eriksson J, Wedin P (2004) Truncated Gauss-Newton algorithms for ill-conditioned nonlinear least squares problems. Optim Meth Softw 19(6):721–737
Faber TL, Raghunath N, Tudorascu D, Votaw JR (2009) Motion correction of PET brain images through deconvolution: I. Theoretical development and analysis in software simulations. Phys Med Biol 54(3):797–811
Figueiredo MAT, Nowak RD, Wright SJ (2007) Gradient projection for sparse reconstruction: Application to compressed sensing and other inverse problems. IEEE J Sel Top Signal Process 1(4):586–597
Frank J (2006) Three-dimensional electron microscopy of macromolecular assemblies. Oxford University Press, New York
Golub GH, Heath M, Wahba G (1979) Generalized cross-validation as a method for choosing a good ridge parameter. Technometrics 21(2):215–223
Golub GH, Luk FT, Overton ML (1981) A block Lanczos method for computing the singular values and corresponding singular vectors of a matrix. ACM Trans Math Softw 7(2): 149–169
Golub GH, Pereyra V (1973) The differentiation of pseudo-inverses and nonlinear least squares problems whose variables separate. SIAM J Numer Anal 10(2):413–432
Golub GH, Pereyra V (2003) Separable nonlinear least squares: the variable projection method and its applications. Inverse Prob 19: R1–R26
Haber E, Ascher UM, Oldenburg D (2000) On optimization techniques for solving nonlinear inverse problems. Inverse Prob 16(5):1263–1280
Haber E, Oldenburg D (2000) A GCV based method for nonlinear ill-posed problems. Comput Geosci 4(1):41–63
Hammerstein GR, Miller DW, White DR, Masterson ME, Woodard HQ, Laughlin JS (1979) Absorbed radiation dose in mammography. Radiology 130(2):485–491
Hanke M (1995) Conjugate gradient type methods for ill-posed problems. Pitman research notes in mathematics, Longman Scientific & Technical, Harlow
Hanke M (1996) Limitations of the L-curve method in ill-posed problems. BIT 36(2):287–301
Hanke M (2001) On Lanczos based methods for the regularization of discrete ill-posed problems. BIT 41(5):1008–1018
Hansen PC (1992) Analysis of discrete ill-posed problems by means of the L-curve. SIAM Rev 34(4):561–580
Hansen PC (1992) Numerical tools for analysis and solution of Fredholm integral equations of the first kind. Inverse Prob 8(6):849–872
Hansen PC (1994) Regularization tools: a MATLAB package for analysis and solution of discrete ill-posed problems. Numer Algorithms 6(1):1–35
Hansen PC (1998) Rank-deficient and discrete ill-posed problems. SIAM, Philadelphia
Hansen PC (2010) Discrete inverse problems: insight and algorithms. SIAM, Philadelphia
Hansen PC, Nagy JG, O’Leary DP (2006) Deblurring images: matrices, spectra and filtering. SIAM, Philadelphia
Hansen PC, O’Leary DP (1993) The use of the L-curve in the regularization of discrete ill-posed problems. SIAM J Sci Comput 14(6):1487–1503
Hardy JW (1994) Adapt Opt Sci Am 270(6): 60–65
Hofmann B (1993) Regularization of nonlinear problems and the degree of ill-posedness. In: Anger G, Gorenflo R, Jochmann H, Moritz H, Webers W (eds) Inverse problems: principles and applications in geophysics, technology, and medicine. Akademie Verlag, Berlin
Hohn M, Tang G, Goodyear G, Baldwin PR, Huang Z, Penczek PA, Yang C, Glaeser RM, Adams PD, Ludtke SJ (2007) SPARX, a new environment for Cryo-EM image processing. J Struct Biol 157(1):47–55
Jain AK (1989) Fundamentals of digital image processing. Prentice-Hall, Englewood Cliffs
Kang MG, Chaudhuri S (2003) Super-resolution image reconstruction. IEEE Signal Process Mag 20(3):19–20
Kaufman L (1975) A variable projection method for solving separable nonlinear least squares problems. BIT 15(1):49–57
Kilmer ME, Hansen PC, Español MI (2007) A projection-based approach to general-form Tikhonov regularization. SIAM J Sci Comput 29(1):315–330
Kilmer ME, O’Leary DP (2001) Choosing regularization parameters in iterative methods for ill-posed problems. SIAM J Matrix Anal Appl 22(4):1204–1221
Landweber L (1951) An iteration formula for Fredholm integral equations of the first kind. Am J Math 73(3):615–624
Larsen RM (1998) Lanczos bidiagonalization with partial reorthogonalization. PhD thesis, Department of Computer Science, University of Aarhus, Denmark
Lawson CL, Hanson RJ (1995) Solving least squares problems. SIAM, Philadelphia
Löfdahl MG (2002) Multi-frame blind deconvolution with linear equality constraints. In: Bones PJ, Fiddy MA, Millane RP (eds) Image reconstruction from incomplete data II, vol 4792-21. SPIE, pp 146–155
Lohmann AW, Paris DP (1965) Space-variant image formation. J Opt Soc Am 55(8):1007–1013
Marabini R, Herman GT, Carazo JM (1998) 3D reconstruction in electron microscopy using ART with smooth spherically symmetric volume elements (blobs). Ultramicroscopy 72(1–2):53–65
Matson CL, Borelli K, Jefferies S, Beckner CC Jr, Hege EK, Lloyd-Hart M (2009) Fast and optimal multiframe blind deconvolution algorithm for high-resolution groundbased imaging of space objects. Appl Opt 48(1):A75–A92
McNown SR, Hunt BR (1994) Approximate shift-invariance by warping shift-variant systems. In: Hanisch RJ, White RL (eds) The restoration of HST images and spectra II. Space Telescope Science Institute, Baltimore, MD, pp 181–187
Miller K (1970) Least squares methods for ill-posed problems with a prescribed bound. SIAM J Math Anal 1(1):52–74
Modersitzki J (2004) Numerical methods for image registration. Oxford University Press, Oxford
Morozov VA (1966) On the solution of functional equations by the method of regularization. Sov Math Dokl 7:414–417
Nagy JG, O’Leary DP (1997) Fast iterative image restoration with a spatially varying PSF. In: Luk FT (ed) Advanced signal processing: algorithms, architectures, and implementations VII, vol 3162. SPIE, pp 388–399
Nagy JG, O’Leary DP (1998) Restoring images degraded by spatially-variant blur. SIAM J Sci Comput 19(4):1063–1082
Natterer F (2001) The mathematics of computerized tomography. SIAM, Philadelphia
Natterer F, Wübbeling F (2001) Mathematical methods in image reconstruction. SIAM, Philadelphia
Nguyen N, Milanfar P, Golub G (2001) Efficient generalized cross-validation with applications to parametric image restoration and resolution enhancement. IEEE Trans Image Process 10(9):1299–1308
Nocedal J, Wright S (1999) Numerical optimization. Springer, New York
O’Leary DP, Simmons JA (1981) A bidiagonali- zation-regularization procedure for large scale discretizations of ill-posed problems. SIAM J Sci Stat Comput 2(4):474–489
Osborne MR (2007) Separable least squares, variable projection, and the Gauss-Newton algorithm. Elec Trans Numer Anal 28:1–15
Paige CC, Saunders MA (1982) Algorithm 583 LSQR: Sparse linear equations and least squares problems. ACM Trans Math Softw 8(2): 195–209
Paige CC, Saunders MA (1982) LSQR: An algorithm for sparse linear equations and sparse least squares. ACM Trans Math Softw 8(1): 43–71
Penczek PA, Radermacher M, Frank J (1992) Three-dimensional reconstruction of single particles embedded in ice. Ultramicroscopy 40(1):33–53
Phillips DL (1962) A technique for the numerical solution of certain integral equations of the first kind. J Assoc Comput Mach 9(1): 84–97
Raghunath N, Faber TL, Suryanarayanan S, Votaw JR (2009) Motion correction of PET brain images through deconvolution: II. Practical implementation and algorithm optimization. Phys Med Biol 54(3):813–829
Rudin LI, Osher S, Fatemi E (1992) Nonlinear total variation based noise removal algorithms. Physica D 60:259–268
Ruhe A, Wedin P (1980) Algorithms for separable nonlinear least squares problems. SIAM Rev 22(3):318–337
Saad Y (1980) On the rates of convergence of the Lanczos and the block-Lanczos methods. SIAM J Numer Anal 17(5):687–706
Saban SD, Silvestry M, Nemerow GR, Stewart PL (2006) Visualization of α-helices in a 6-Ångstrom resolution cryoelectron microscopy structure of adenovirus allows refinement of capsid protein assignments. J Virol 80(24): 49–59
Tikhonov AN (1963) Regularization of incorrectly posed problems. Sov Math Dokl 4:1624–1627
Tikhonov AN (1963) Solution of incorrectly formulated problems and the regularization method. Sov Math Dokl 4:1035–1038
Tikhonov AN, Arsenin VY (1977) Solutions of ill-posed problems. Winston, Washington
Tikhonov AN, Leonov AS, Yagola AG (1998) Nonlinear ill-posed problems, vol 1–2. Chapman and Hall, London
Trussell HJ, Fogel S (1992) Identification and restoration of spatially variant motion blurs in sequential images. IEEE Trans Image Process 1(1):123–126
Tsaig Y, Donoho DL (2006) Extensions of compressed sensing. Signal Process 86(3): 549–571
Varah JM (1983) Pitfalls in the numerical solution of linear ill-posed problems. SIAM J Sci Stat Comput 4(2):164–176
Vogel CR (1986) Optimal choice of a truncation level for the truncated SVD solution of linear first kind integral equations when data are noisy. SIAM J Numer Anal 23(1):109–117
Vogel CR (1987) An overview of numerical methods for nonlinear ill-posed problems. In: Engl HW, Groetsch CW (eds) Inverse and ill-posed problems. Academic Press, Boston, pp 231–245
Vogel CR (1996) Non-convergence of the L-curve regularization parameter selection method. Inverse Prob 12(4):535–547
Vogel CR (2002) Computational methods for inverse problems. SIAM, Philadelphia
Wagner FC, Macovski A, Nishimura DG (1988) A characterization of the scatter pointspread-function in terms of air gaps. IEEE Trans Med Imaging 7(4):337–344
Editor information
Editors and Affiliations
Rights and permissions
Copyright information
© 2011 Springer Science+Business Media, LLC
About this entry
Cite this entry
Chung, J., Knepper, S., Nagy, J.G. (2011). Large-Scale Inverse Problems in Imaging. In: Scherzer, O. (eds) Handbook of Mathematical Methods in Imaging. Springer, New York, NY. https://doi.org/10.1007/978-0-387-92920-0_2
Download citation
DOI: https://doi.org/10.1007/978-0-387-92920-0_2
Publisher Name: Springer, New York, NY
Print ISBN: 978-0-387-92919-4
Online ISBN: 978-0-387-92920-0
eBook Packages: Mathematics and StatisticsReference Module Computer Science and Engineering