A trust-region approach to the regularization of large-scale discrete forms of ill-posed problems

SIAM Journal on Scientific Computing

Volumn 23, Issue 6, 2002, Pages 1842-1860

  Rojas, Marielba ^a   Sorensen, Danny C ^b  

^a CERFACS   (France)

^b Rice University   (United States)

Author keywords

Constrained quadratic optimization; Ill posed problems; Inverse problems; Lanczos method; Regularization; Seismic inversion; Trust region

Indexed keywords

INTERPOLATION; INVERSE PROBLEMS; LEAST SQUARES APPROXIMATIONS; MATRIX ALGEBRA; OPTIMIZATION; VECTORS;

LANCZOS METHODS;

PROBLEM SOLVING;

EID: 0036445949     PISSN: 10648275     EISSN: None     Source Type: Journal    
DOI: 10.1137/S1064827500378167     Document Type: Article

References (40)

1. Z. Ben-Avraham, G. Amit, A. Golan, And Z.B. Begin, The bathymetry of Lake Kinneret and its structural significance, Israel J. Earth Sci., 39 (1992), pp. 77-84.
2. M. Bertero And P. Bocacci, Introduction to Inverse Problems in Imaging, Institute of Physics, Bristol, UK, 1998.
3. Å. Björck, Numerical Methods for Least Squares Problems, SIAM, Philadelphia, 1996.
4. Å. Björck, E. Grimme, And P. Van Dooren, An implicit shift bidiagonalization algorithm for ill-posed systems, BIT, 34 (1994), pp. 510-534.
5. J.O. Blanch, A Study of Viscous Effects in Seismic Modeling, Imaging, and Inversion: Methodology, Computational Aspects; and Sensitivity, Ph.D. thesis, Department of Geology and Geophysics, Rice University, Houston, 1995.
6. J.O. Blanch, W.W. Sykes, And R.J. Versteeg, A numerical study of linear viscoacoustic inversion, in Comparison of Seismic Inversion Methods on a Single Real Data Set, R.G. Keys and D.J. Foster, eds., Society of Exploration Geophysicists, Tulsa, OK, 1998, pp. 13-44.
7. D. Calvetti, G.H. Golub, And L. Reichel, Estimation of the L-curve via Lanczos bidiagonalization, BIT, 39 (1999), pp. 603-619.
8. D. Calvetti, L. Reichel, And Q. Zhang, Iterative exponential filtering for large discrete ill-posed problems, Numar. Math., 83 (1999), pp. 535-556.
9. J. Clærbout, Geophysical Estimation by Example: Environmental Soundings Image Enhancement: Multidimensional Autoregression, available online from http://sepwww.stanford. edu/sep/jon/index.html.
10. L. Eldën, Algorithms for the regularization of ill-conditioned least squares problems, BIT, 17 (1977), pp. 134-145.
11. D.M. Gay, Computing optimal locally constrained steps, SIAM J. Sci. Statist. Comput., 2 (1981), pp. 186-197.
12. M.S. Gockenbach And W.W. Symes, Duality for inverse problems in wave propagation, in Large Scale Optimization, IMA Vol. Math, Appl. 92, L. Biegler, T. Coleman, F. Santosa, and A. Conn, eds., Springer-Verlag, New York, 1997, pp. 37-61.
13. G.H. Golub And C.F. Van Loan Matrix Computations, 3rd ed., Johns Hopkins University Press, Baltimore, 1996.
14. G.H. Golub And U. Von Matt, Quadratically constrained least squares and quadratic problems, Numer. Math., 59 (1991), pp. 561-580.
15. M. Hanke And P.C. Hansen, Regularization methods for large-scale problems, Surveys Math. Indust., 3 (1993), pp. 253-315.
16. P.C. Hansen, The discrete Picard condition for discrete ill-posed problems, BIT, 30 (1990), pp. 658-672.
17. P.G. Hansen, Analysis of discrete ill-posed problems by means of the L-curve, SIAM Rev., 34 (1992), pp. 561-580.
18. P.C. Hansen And D.P. O'Leary, The use of the L-curve in the regularization of discrete ill-posed problems, SIAM J. Sci. Comput., 14 (1993), pp. 1487-1503.
19. P.C. Hansen, Regularization Tools: A MATLAB package for analysis and solution of discrete ill-posed problems, Numar. Algorithms, 6 (1994), pp. 1-35.
20. P.O. Hansen, Rank-Deficient and Discrete Ill-Posed Problems, Doctoral Dissertation, UNI•C, Technical University of Denmark, Lyngby, Denmark, 1996.
21. P.O. Hansen, M. Jacobsen, J.M. Rasmussen, And H. Sørensen, The PP-TSVD algorithm for image restoration problems, in Methods and Applications of Inversion, Lecture Notes in Earth Sci. 92, P.C. Hansen, B.H. Jacobsen, and K. Mosegaard, eds., Springer-Verlag, Berlin, 2000, pp. 171-186.
22. L. Kaufman And A. Neumaier, Image Reconstruction through Regularization by Envelope Guided Conjugate Gradients, AT&T Bell Laboratories, 4-14, 1994; available online from http:cm.bells-labs.com/cm/cs/doe/94/4-14.ps.gz.
23. R.G. Keys And D.J. Foster, A data set for evaluating and comparing seismic inversion methods, in Comparison of Seismic Inversion Methods on a Single Real Data Set, R.G. Keys and D.J. Foster, eds., Society of Exploration Geophysicists, Tulsa, OK, 1998, pp. 1-12.
24. M.E. Kilmer And D.P. O'Leary, Pivoted Cauchy-like preconditioners for regularized solution of ill-posed problems, SIAM J. Sci. Comput, 21 (1999), pp. 88-110.
25. M.E. Kilmer And D.P. O'Leary, Choosing regularization parameters in iterative methods for ill-posed problems, SIAM J. Matrix Anal. Appl., 22 (2001), pp. 1204-1221.
26. C.L. Lawson And R.J. Hanson, Solving Least Squares Problems, Classics Appl. Math. 15, SIAM, Philadelphia, 1995.
27. R.B. Lehoucq, D.C. Sorensen, And C. Yang, ARPACK Users' Guide: Solution of Large-Scale Eigenvalue Problems by Implicitly Restarted Arnoldi Methods, SIAM, Philadelphia, 1998.
28. J.J. Moré And D.C. Sorensen, Computing a trust region step, SIAM J. Sci. Statist. Comput., 4 (1983), pp. 553-572.
29. J. Nagy, V. Pauga, R. Plemmons, And T. Torgersen, Space-varying restoration of optical images, Opt. Soc. Amer. A, 14 (1997), pp. 3162-3174.
30. F. Natterer, The Mathematics of Computerized Tomography, John Wiley, New York, 1986.
31. A. Neumaier, Solving ill-conditioned and singular linear systems: A tutorial on regularization, SIAM Rev., 40 (1998), pp. 636-666.
32. G. Nolet, Ed., Seismic Tomography, with Applications in Global Seismology and Exploration Geophysics, D. Reidel, Dordrecht, 1987.
33. B.N. Parlett, The Symmetric Eigenvalue Problem, Prentice-Hall, Englewood Cliffs, NJ, 1980.
34. M. Rojas, A Large-Scale Trust-Region Approach to the Regularization of Discrete Ill-Posed Problems, Ph.D. thesis, Technical Report TR98-19, Department of Computational and Applied Mathematics, Rice University, Houston, TX, 1998.
35. M. Rojas, S.A. Santos, And D.C. Sorensen, A new matrix-free algorithm for the large-scale trust-region subproblem, SIAM J. Optim., 11 (2000), pp. 611-646.
36. D.C. Sorensen, Newton's method with a model trust region modification, SIAM J. Numer. Anal., 19 (1982), pp. 409-425.
37. D.C. Sorensen, Implicit application of polynomial filters in a k-step Arnoldi method, SIAM J. Matrix Anal. Appl, 13 (1992), pp. 357-385.
38. D.C. Sorensen, Minimization of a large-scale quadratic function subject to a spherical constraint, SIAM J. Optim., 7 (1997), pp. 141-161.
39. W.W. Symes, A differential semblance criterion for inversion of multioffset seismic reflection data, J. Geophysical Research, 98 (1993), pp. 2061-2073.
40. A.N. Tikhonov, Regularization of incorrectly posed problems, Soviet Math., 4 (1963), pp. 1624-1627.