# zbMATH — the first resource for mathematics

Nonlinear regularization techniques for seismic tomography. (English) Zbl 1182.86003
Summary: The effects of several nonlinear regularization techniques are discussed in the framework of 3D seismic tomography. Traditional, linear, $$\ell _{2}$$ penalties are compared to so-called sparsity promoting $$\ell _{1}$$ and $$\ell _{0}$$ penalties, and a total variation penalty. Which of these algorithms is judged optimal depends on the specific requirements of the scientific experiment. If the correct reproduction of model amplitudes is important, classical damping towards a smooth model using an $$\ell _{2}$$ norm works almost as well as minimizing the total variation but is much more efficient. If gradients (edges of anomalies) should be resolved with a minimum of distortion, we prefer $$\ell _{1}$$ damping of Daubechies-4 wavelet coefficients. It has the additional advantage of yielding a noiseless reconstruction, contrary to simple $$\ell _{2}$$ minimization (‘Tikhonov regularization’) which should be avoided. In some of our examples, the $$\ell _{0}$$ method produced notable artifacts. In addition we show how nonlinear $$\ell _{1}$$ methods for finding sparse models can be competitive in speed with the widely used $$\ell _{2}$$ methods, certainly under noisy conditions, so that there is no need to shun $$\ell _{1}$$ penalizations.

##### MSC:
 86A15 Seismology (including tsunami modeling), earthquakes 86A22 Inverse problems in geophysics 42C40 Nontrigonometric harmonic analysis involving wavelets and other special systems
##### Keywords:
inverse problem; one-norm; sparsity; tomography; wavelets; regularization
##### Software:
DT-CWT; SPGL1; TwIST
Full Text:
##### References:
 [1] Backus, G.; Gilbert, J., Numerical applications of a formalism for geophysical inverse problems, Geophysical journal of the royal astronomical society, 13, 247-276, (1967) [2] Backus, G.; Gilbert, J., Uniqueness on the inversion of inaccurate Gross Earth data, Philosophical transactions of the royal society of London, A266, 123-192, (1970) [3] Loris, I.; Nolet, G.; Daubechies, I.; Dahlen, F.A., Tomographic inversion using $$\ell_1$$-norm regularization of wavelet coefficients, Geophysical journal international, 170, 1, 359-370, (2007) [4] Blumensath, T.; Davies, M., Iterative thresholding for sparse approximations, Journal of Fourier analysis and applications, 14, 5, 629-654, (2008) · Zbl 1175.94060 [5] Blumensath, T.; Davies, M.E., Iterative hard thresholding for compressed sensing, Applied and computational harmonic analysis, 27, 3, 265-274, (2009) · Zbl 1174.94008 [6] Chevrot, S.; Zhao, L., Multiscale finite frequency Rayleigh wave tomography of the kaapvaal craton, Geophysical journal international, 169, 201-215, (2007) [7] Scales, J.; Gersztenkorn, A.; Treitel, S., Fast $$\ell_p$$ solutions of large sparse linear systems: application to seismic travel time tomography, Journal of computational physics, 75, 314-333, (1988) · Zbl 0637.65035 [8] Nolet, G., Seismic wave propagation and seismic tomography, (), 1-23 [9] Taylor, H.L.; Banks, S.C.; McCoy, J.F., Deconvolution with the $$\ell_1$$ norm, Geophysics, 44, 39-52, (1979) [10] Oldenburg, D.; Scheuer, T.; Levy, S., Recovery of the acoustic impedance from reflection seismograms, Geophysics, 48, 1318-1337, (1983) [11] Santosa, F.; Symes, W.W., Linear inversion of band-limited reflection seismograms, SIAM journal on scientific and statistical computing, 7, 1307-1330, (1986) · Zbl 0602.73113 [12] Sacchi, M.D.; Ulrych, T.J.; Walker, C., Interpolation and extrapolation using a high-resolution discrete Fourier transform, IEEE transactions on signal processing, 46, 31-38, (1998) [13] Candès, E.J.; Romberg, J.K.; Tao, T., Stable signal recovery from incomplete and inaccurate measurements, Communications on pure and applied mathematics, 59, 8, 1207-1223, (2006) · Zbl 1098.94009 [14] Donoho, D., Compressed sensing, IEEE transactions on information theory, 52, 4, 1289-1306, (2006) · Zbl 1288.94016 [15] Donoho, D.L., For most large underdetermined systems of linear equations the minimal $$\ell_1$$-norm solution is also the sparsest solution, Communications on pure and applied mathematics, 59, 797-829, (2006) · Zbl 1113.15004 [16] Herrmann, F.J.; Hennenfent, G., Non-parametric seismic data recovery with curvelet frames, Geophysical journal international, 173, 1, 233-248, (2008) [17] Hennenfent, G.; van den Berg, E.; Friedlander, M.P.; Herrmann, F.J., New insights into one-norm solvers from the Pareto curve, Geophysics, 73, 4, A23-A26, (2008) [18] van den Berg, E.; Friedlander, M.P., Probing the Pareto frontier for basis pursuit solutions, SIAM journal on scientific computing, 31, 2, 890-912, (2008) · Zbl 1193.49033 [19] Favier, N.; Chevrot, S., Sensitivity kernels for shear wave splitting in transverse isotropic media, Geophysical journal international, 153, 213-228, (2003) [20] Nolet, G., A breviary of seismic tomography, (2008), Cambridge University Press · Zbl 1166.86001 [21] Dahlen, F.A., Resolution limit of traveltime tomography, Geophysical journal international, 157, 1, 315-331, (2004) [22] Tikhonov, A., Solution of incorrectly formulated problems and the regularization method, Soviet mathematics – doklady, 4, 1035-1038, (1963) · Zbl 0141.11001 [23] Daubechies, I.; Defrise, M.; De Mol, C., An iterative thresholding algorithm for linear inverse problems with a sparsity constraint, Communications on pure and applied mathematics, 57, 11, 1413-1457, (2004) · Zbl 1077.65055 [24] I. Daubechies, Ten Lectures on Wavelets, Society for Industrial and Applied Mathematics, 1992. · Zbl 0776.42018 [25] Beck, A.; Teboulle, M., A fast iterative shrinkage-threshold algorithm for linear inverse problems, SIAM journal on imaging sciences, 2, 183-202, (2009) · Zbl 1175.94009 [26] Nesterov, Y.E., A method for solving the convex programming problem with convergence rate $$\mathcal{O}(1 / k^2)$$, Doklady academii nauk SSSR, 269, 543-547, (1983), (in Russian) [27] Nesterov, Y.E., A method for solving a convex programming problem with convergence rate $$\mathcal{O}(1 / k^2)$$, Soviet mathematics – doklady, 27, 372-376, (1983) · Zbl 0535.90071 [28] Nesterov, Y., Introductory lectures on convex optimization, (2004), Kluwer Academic Publishers · Zbl 1086.90045 [29] Oliveira, J.; Bioucas-Dias, J.; Figueiredo, M., Adaptive total variation image deblurring: a majorization – minimization approach, Signal processing, 89, 9, 1683-1693, (2009) · Zbl 1178.94029 [30] Bect, J.; Blanc-Féraud, L.; Aubert, G.; Chambolle, A., A $$\ell^1$$-unified variational framework for image restoration, (), 1-13 · Zbl 1098.68728 [31] Bioucas-Dias, J.; Figueiredo, M., A new twist: two-step iterative shrinkage/thresholding algorithms for image restoration, IEEE transactions on image processing, 16, 2980-2991, (2007) [32] Haar, A., Zur theorie der orthogonalen funktionensysteme, Mathematische annalen, 69, 331-371, (1910) · JFM 41.0469.03 [33] Cohen, A.; Daubechies, I.; Jawerth, B.; Vial, P., Multiresolution analysis, wavelets and fast algorithms on the interval, Comptes rendus de l’académie des sciences Paris, 316, 417-421, (1992) · Zbl 0768.42015 [34] Cohen, A.; Daubechies, I.; Vial, P., Wavelets on the interval and fast wavelet transforms, Journal for applied and computational harmonic analysis, 1, 54-81, (1993) · Zbl 0795.42018 [35] Kingsbury, N., Image processing with complex wavelets, Philosophical transactions of the royal society A: mathematical, physical and engineering sciences, 357, 1760, 2543-2560, (1999) · Zbl 0976.68527 [36] Kingsbury, N., Complex wavelets for shift invariant analysis and filtering of signals, Applied and computational harmonic analysis, 10, 234-253, (2001) · Zbl 0990.94005 [37] D. Labate, W. Lim, G. Kutyniok, G. Weiss, Sparse multidimensional representation using shearlets, in: SPIE Proc. 5914, SPIE, Bellingham, WA, 2005, pp. 254-262, Wavelets XI (San Diego, CA, 2005). [38] Kutyniok, G.; Labate, D., Construction of regular and irregular shearlets, Journal of wavelet theory and applications, 1, 1-10, (2007) [39] Candes, E.J.; Demanet, L.; Ying, D.L.D.L., Fast discrete curvelet transforms, SIAM multiscale modeling and simulation, 5, 3, 861-899, (2006) · Zbl 1122.65134 [40] M.A.T. Figueiredo, R.D. Nowak, S.J. Wright, Gradient projection for sparse reconstruction: application to compressed sensing and other inverse problems, IEEE Journal of Selected Topics in Signal Processing (Special Issue on Convex Optimization Methods for Signal Processing) 1 (2007) 586-598. doi:10.1109/JSTSP.2007.910281. [41] Elad, M.; Matalon, B.; Zibulevsky, M., Coordinate and subspace optimization methods for linear least squares with non-quadratic regularization, Applied and computational harmonic analysis, 23, 3, 346-367, (2007) · Zbl 1133.65022 [42] Kim, S.-J.; Koh, K.; Lustig, M.; Boyd, S.; Gorinevsky, D., An interior-point method for large-scale $$\ell_1$$-regularized least squares, IEEE journal of selected topics in signal processing, 1, 4, 606-617, (2007) [43] Hale, E.; Yin, W.; Zhang, Y., Fixed-point continuation for $$\ell_1$$-minimization: methodology and convergence, SIAM journal on optimization, 19, 3, 1107-1130, (2008) · Zbl 1180.65076 [44] Yin, W.; Osher, S.; Goldfarb, D.; Darbon, J., Bregman iterative algorithms for $$\ell_1$$-minimization with applications to compressed sensing, SIAM journal on imaging sciences, 1, 1, 143-168, (2008) · Zbl 1203.90153 [45] Loris, I., On the performance of algorithms for the minimization of $$\ell_1$$-penalized functionals, Inverse problems, 25, 035008, (2009) · Zbl 1162.65333 [46] Defrise, M.; De Mol, C., A note on stopping rules for iterative regularization methods and filtered SVD, (), 261-268
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. It attempts to reflect the references listed in the original paper as accurately as possible without claiming the completeness or perfect precision of the matching.