Sparse wavelet expansions for seismic tomography: Methods and algorithms
|
|
- Sharyl McCormick
- 5 years ago
- Views:
Transcription
1 Sparse wavelet expansions for seismic tomography: Methods and algorithms Ignace Loris Université Libre de Bruxelles International symposium on geophysical imaging with localized waves July 2011 (Joint work with: G. Nolet, I.Daubechies, F.J.Simons, C. Verhoeven et al.)
2 Theme Discussion of: fewer data than unknowns -problem recent advances in regularization methods in context of global seismic tomography How to exploit such insight in seismic inversions? Philosophy: Ockham s razor the simplest explanation is usually the best In applied mathematics: principle of sparsity (here sparsity refers to a property of the model, NOT of a matrix)
3 Overview Mathematical developments: How can model sparsity help in solving an inverse problem? Illustration on toy example Algorithms Application to seismic tomography: Good model parametrization: so-called wavelets Mathematical microscope : analyze seismic heterogeneity at different length scales Computationally efficient (fast algorithm) Can express seismic heterogeneity efficiently Formulate solution method
4 Mathematical framework Solve linear relations between unknown x and measurement data y: Kx = y Here: K =matrix with sensitivity kernels (known) y =data vector (known) x =model vector (coefficients representing the wave speed anomaly) Problems: insufficient data, inconsistent data, ill-conditioning of K No solution or no unique solution Way out: minimize a penalized least-squares functional: x rec = arg min Kx y 2 + penalty x
5 Penalization strategy Traditional: quadratic penalties penalty x 2 or penalty x 2 = smoothing operator Advantage: leads to linear equations K T Kx + λx = K T y or K T Kx + λ T x = K T y However: model can be too smooth... Can better reconstructions be found by using other information? Increased interest in using sparsity as a priori information (sparsity solution described by few nonzero coefficients)
6 Penalization strategy How to enforce sparsity on solutions of Kx = y? Minimize a penalized least-squares functional: x rec 1 = arg min Kx x 2 y 2 + penalty with a judiciously chosen sparsity promoting penalty. A trade-off between sparsity promotion and tractability (convexity) D. Donoho (2006) DOI: /TIT A. Bruckstein et al. (2009) DOI: /
7 Penalization strategy How to enforce sparsity on solutions of Kx = y? Minimize a penalized least-squares functional: x rec 1 = arg min Kx x 2 y 2 + penalty with a judiciously chosen sparsity promoting penalty. A trade-off between sparsity promotion and tractability (convexity) D. Donoho (2006) DOI: /TIT A. Bruckstein et al. (2009) DOI: /
8 Penalization strategy How to enforce sparsity on solutions of Kx = y? Minimize a penalized least-squares functional: x rec 1 = arg min Kx x 2 y 2 + penalty with a judiciously chosen sparsity promoting penalty. A trade-off between sparsity promotion and tractability (convexity) x 2 Kx = y x 1 # nonzero 1 combinatorial D. Donoho (2006) DOI: /TIT A. Bruckstein et al. (2009) DOI: /
9 Penalization strategy How to enforce sparsity on solutions of Kx = y? Minimize a penalized least-squares functional: x rec 1 = arg min Kx x 2 y 2 + penalty with a judiciously chosen sparsity promoting penalty. A trade-off between sparsity promotion and tractability (convexity) x 2 Kx = y x 2 Kx = y x 1 x 1 # nonzero 1 x1 + x 2 1 combinatorial algebraic D. Donoho (2006) DOI: /TIT A. Bruckstein et al. (2009) DOI: /
10 Penalization strategy How to enforce sparsity on solutions of Kx = y? Minimize a penalized least-squares functional: x rec 1 = arg min Kx x 2 y 2 + penalty with a judiciously chosen sparsity promoting penalty. A trade-off between sparsity promotion and tractability (convexity) x 2 Kx = y x 2 Kx = y x 2 Kx = y x 1 x 1 x 1 # nonzero 1 x1 + x 2 1 x 1 + x 2 1 combinatorial algebraic convex D. Donoho (2006) DOI: /TIT A. Bruckstein et al. (2009) DOI: /
11 Penalization strategy How to enforce sparsity on solutions of Kx = y? Minimize a penalized least-squares functional: x rec 1 = arg min Kx x 2 y 2 + penalty with a judiciously chosen sparsity promoting penalty. A trade-off between sparsity promotion and tractability (convexity) x 2 Kx = y x 2 Kx = y x 2 Kx = y x 2 Kx = y x 1 x 1 x 1 x 1 # nonzero 1 x1 + x 2 1 x 1 + x 2 1 x1 2 + x2 2 r 2 combinatorial algebraic convex linear D. Donoho (2006) DOI: /TIT A. Bruckstein et al. (2009) DOI: /
12 Penalization strategy How to enforce sparsity on solutions of Kx = y? Minimize a penalized least-squares functional: x rec 1 = arg min Kx x 2 y 2 + penalty with a judiciously chosen sparsity promoting penalty. A trade-off between sparsity promotion and tractability (convexity) x 2 Kx = y x 2 Kx = y x 2 Kx = y x 2 Kx = y x 1 x 1 x 1 x 1 # nonzero 1 combinatorial algebraic x1 + x 2 1 convex x 1 + x 2 1 linear x x2 2 r 2 l 1 norm penalty promotes sparsity and is tractable ( x 1 = i x i ) D. Donoho (2006) DOI: /TIT A. Bruckstein et al. (2009) DOI: /
13 Sparsity and l 1 norm Unit l 1 ball in 2-D 1 2
14 Sparsity and l 1 norm Unit l 1 ball in N-D 1 N
15 Sparsity and l 1 norm Unit l 1 ball in N-D 1 N looks like (when N is large)
16 Illustration of the efficiency of sparse recovery x input K
17 Illustration of the efficiency of sparse recovery x input K Synthetic noisy data: y = Kx input + n
18 Illustration of the efficiency of sparse recovery 30 K Synthetic noisy data: y = Kx input + n
19 Illustration of the efficiency of sparse recovery 30 K Synthetic noisy data: y = Kx input + n min x Kx y 2 + λ x 2 min x Kx y 2 + 2λ x 1
20 Illustration of the efficiency of sparse recovery 30 K Synthetic noisy data: y = Kx input + n min Kx y 2 + λ x 2 min x x Kx y 2 + 2λ x
21 Illustration of the efficiency of sparse recovery x input K Synthetic noisy data: y = Kx input + n min Kx y 2 + λ x 2 min x x Kx y 2 + 2λ x
22 Analysis sparsity vs. synthesis sparsity Analysis-style sparsity: Require that certain linear combinations Ax of unknown x are sparse: x with many (Ax) i = 0 Here A is explicitly known. Synthesis-style sparsity: Express unknown x as a sparse linear combination of a set of known basis/frame vectors: x = Su with many u i = 0 If AS = SA = 1 then: synthesis sparsity x rec = analysis sparsity x rec (algorithms may still be different)
23 Analysis sparsity vs. synthesis sparsity Analysis-style sparsity: Require that certain linear combinations Ax of unknown x are sparse: x with many (Ax) i = 0 Here A is explicitly known. Synthesis-style sparsity: Express unknown x as a sparse linear combination of a set of known basis/frame vectors: x = Su with many u i = 0 If AS = SA = 1 then: synthesis sparsity x rec = analysis sparsity x rec (algorithms may still be different)
24 Analysis sparsity vs. synthesis sparsity Analysis-style sparsity: Require that certain linear combinations Ax of unknown x are sparse: x with many (Ax) i = 0 Here A is explicitly known. Synthesis-style sparsity: Express unknown x as a sparse linear combination of a set of known basis/frame vectors: x = Su with many u i = 0 If AS = SA = 1 then: synthesis sparsity x rec = analysis sparsity x rec (algorithms may still be different)
25 Analysis sparsity vs. synthesis sparsity Analysis-style sparsity: Require that certain linear combinations Ax of unknown x are sparse: x with many (Ax) i = 0 Here A is explicitly known. Synthesis-style sparsity: Express unknown x as a sparse linear combination of a set of known basis/frame vectors: x = Su with many u i = 0 If AS = SA = 1 then: synthesis sparsity x rec = analysis sparsity x rec (algorithms may still be different)
26 Cost functions for analysis and synthesis sparsity Analysis-style sparsity (require that Ax is sparse): x rec 1 = arg min Kx x 2 y 2 + λ Ax 1 (1) Primary example of (1) is total variation (TV) penalty: A =grad Synthesis-style sparsity (express x = Su with u sparse): u rec 1 = arg min KSu u 2 y 2 + λ u 1 and x rec = Su rec (2) Example of (2): sparse combination of wavelets (2) is a special case of (1): A 1 and K KS Variational equations are nonlinear but problem is convex: convex minimization techniques can be applied
27 Cost functions for analysis and synthesis sparsity Analysis-style sparsity (require that Ax is sparse): x rec 1 = arg min Kx x 2 y 2 + λ Ax 1 (1) Primary example of (1) is total variation (TV) penalty: A =grad Synthesis-style sparsity (express x = Su with u sparse): u rec 1 = arg min KSu u 2 y 2 + λ u 1 and x rec = Su rec (2) Example of (2): sparse combination of wavelets (2) is a special case of (1): A 1 and K KS Variational equations are nonlinear but problem is convex: convex minimization techniques can be applied
28 Cost functions for analysis and synthesis sparsity Analysis-style sparsity (require that Ax is sparse): x rec 1 = arg min Kx x 2 y 2 + λ Ax 1 (1) Primary example of (1) is total variation (TV) penalty: A =grad Synthesis-style sparsity (express x = Su with u sparse): u rec 1 = arg min KSu u 2 y 2 + λ u 1 and x rec = Su rec (2) Example of (2): sparse combination of wavelets (2) is a special case of (1): A 1 and K KS Variational equations are nonlinear but problem is convex: convex minimization techniques can be applied
29 Cost functions for analysis and synthesis sparsity Analysis-style sparsity (require that Ax is sparse): x rec 1 = arg min Kx x 2 y 2 + λ Ax 1 (1) Primary example of (1) is total variation (TV) penalty: A =grad Synthesis-style sparsity (express x = Su with u sparse): u rec 1 = arg min KSu u 2 y 2 + λ u 1 and x rec = Su rec (2) Example of (2): sparse combination of wavelets (2) is a special case of (1): A 1 and K KS Variational equations are nonlinear but problem is convex: convex minimization techniques can be applied
30 Cost functions for analysis and synthesis sparsity Analysis-style sparsity (require that Ax is sparse): x rec 1 = arg min Kx x 2 y 2 + λ Ax 1 (1) Primary example of (1) is total variation (TV) penalty: A =grad Synthesis-style sparsity (express x = Su with u sparse): u rec 1 = arg min KSu u 2 y 2 + λ u 1 and x rec = Su rec (2) Example of (2): sparse combination of wavelets (2) is a special case of (1): A 1 and K KS Variational equations are nonlinear but problem is convex: convex minimization techniques can be applied
31 Algorithm for analysis style sparsity To solve: 1 arg min Kx x 2 y 2 + λ Ax 1 Explicit iterative algorithm: x n+1 = x n + K T (y Kx n ) A T w n w n+1 ( = P λ w n + A x n+1) x n+1 = x n + K T (y Kx n ) A T w n+1, λ u u > λ with P λ (u) = u (component-wise) u u λ Explicit: uses only K, K T, A and A T at every step, No equation solving at every step No non-trivial sub-problem at every step No smoothing of u 1 with small ǫ: u i ui 2 + ǫ 2 I. Loris and C. Verhoeven (2011).
32 Algorithm for synthesis style sparsity When A = 1 then the previous algorithm reduces to iterative soft-thresholding : ) x n+1 = S λ (x n + K T (y Kx n 1 ) for arg min Kx x 2 y 2 + λ x 1 S λ (x) with S λ = soft-thresholding: Advantages: 1 very simple 2 F(x n ) F(x ) x 0 x 2 3 Can be accelerated 2n λ I. Daubechies et al. (2004). DOI: /cpa λ x n > 0
33 Example: Total variation penalty Total Variation: A =grad 1 arg min Kx x 2 y 2 + λ grad(x) 1 grad(x) is sparse grad(x) is piece-wise zero x is piece-wise constant
34 Toy TV regularization problem Seismic surface wave tomography (2D): T = receiver source (ray approximation) ds c receiver perturbation on c 1 δc δt = ds source c b c b y = Kx This example: 8490 data for unknowns. Sum of rows of K Min=0 Max= Sum of rows of K
35 (10 5 iterations of Generalized IST, 10% noise on data) Synthetic model and reference TV reconstruction Input noisy data Reconstruction (A =grad) Min= 1 Max= grad(input) Min=0 Max=
36 (10 5 iterations of Generalized IST, 10% noise on data) Synthetic model and reference TV reconstruction Input noisy data Reconstruction (A =grad) Min= 1 Max= grad(input) Min= Max= grad(reconstruction) Min=0 Max= Min=2.7207e 014 Max=
37 Alternatives for Total Variation (1) TV will impose piece-wise constant reconstruction Good: sharp edge will be preserved Bad: Smooth transition replaced by succession of sharp edges ( staircasing effect ) Improvements: generalized TV, Huber prior,... Other possibility: Try synthesis sparsity with an appropriate basis: express x = Su with u sparse S = basis or frame that can represent an overall smooth model with some sharp edges in a sparse manner S may be overcomplete
38 Alternatives for Total Variation (1) For synthesis style sparsity, faster algorithms exist Fast Iterative Soft-thresholding algorithm: with u n+1 = T(u n + β n (u n u n 1 )) 1 T(u) S λ [ u + S T K T (y KSu) ], 2 β n = n 1 n+2, 3 S λ = soft-thresholding 4 F(x n ) F(x ) 4 x 0 x 2 Sλ(u) λ (n + 1) 2 n > 0 A. Beck and M. Teboule (2009). DOI: / Yu E. Nesterov, Soviet Math. Dokl. 27 (1983) λ u
39 Application: Seismic tomography of Mantle Model parametrization is key to success of sparse recovery: Spherical harmonics well-adapted to spherical geometry good spectral localization poor spatial localization sparse spherical harmonic parametrization very smooth model Wavelet basis Constructed by local averaging/differencing and by subsampling (iterated filter bank) In other words: based on translating and resizing (scaling) of a single function Good compromise between spectral and spatial localization Sparse in wavelet basis = some discontinuities (edges) allowed Fast algorithms exist (O(N) for length N signal) Not well-adapted to spherical geometry Wavelet basis is a good candidate for setting-up a sparse regularization method for seismic tomography
40 Application: Seismic tomography of Mantle Model parametrization is key to success of sparse recovery: Spherical harmonics well-adapted to spherical geometry good spectral localization poor spatial localization sparse spherical harmonic parametrization very smooth model Wavelet basis Constructed by local averaging/differencing and by subsampling (iterated filter bank) In other words: based on translating and resizing (scaling) of a single function Good compromise between spectral and spatial localization Sparse in wavelet basis = some discontinuities (edges) allowed Fast algorithms exist (O(N) for length N signal) Not well-adapted to spherical geometry Wavelet basis is a good candidate for setting-up a sparse regularization method for seismic tomography
41 Application: Seismic tomography of Mantle Model parametrization is key to success of sparse recovery: Spherical harmonics well-adapted to spherical geometry good spectral localization poor spatial localization sparse spherical harmonic parametrization very smooth model Wavelet basis Constructed by local averaging/differencing and by subsampling (iterated filter bank) In other words: based on translating and resizing (scaling) of a single function Good compromise between spectral and spatial localization Sparse in wavelet basis = some discontinuities (edges) allowed Fast algorithms exist (O(N) for length N signal) Not well-adapted to spherical geometry Wavelet basis is a good candidate for setting-up a sparse regularization method for seismic tomography
42 Application: Seismic tomography of Mantle Model parametrization is key to success of sparse recovery: Spherical harmonics well-adapted to spherical geometry good spectral localization poor spatial localization sparse spherical harmonic parametrization very smooth model Wavelet basis Constructed by local averaging/differencing and by subsampling (iterated filter bank) In other words: based on translating and resizing (scaling) of a single function Good compromise between spectral and spatial localization Sparse in wavelet basis = some discontinuities (edges) allowed Fast algorithms exist (O(N) for length N signal) Not well-adapted to spherical geometry Wavelet basis is a good candidate for setting-up a sparse regularization method for seismic tomography
43 Wavelets for seismic tomography of Mantle Many wavelet families exist These will yield different reconstructions Important: the ability to easily switch between different wavelet bases Rather not use a custom-designed wavelet basis for the sphere Preference goes to using cartesian geometry (allows for off-the-shelve wavelets but also curvelets, shearlets etc) Needs a special parametrization of the sphere: cubed sphere (Ronchi et al, 1996)
44 Wavelets for seismic tomography of Mantle Many wavelet families exist These will yield different reconstructions Important: the ability to easily switch between different wavelet bases Rather not use a custom-designed wavelet basis for the sphere Preference goes to using cartesian geometry (allows for off-the-shelve wavelets but also curvelets, shearlets etc) Needs a special parametrization of the sphere: cubed sphere (Ronchi et al, 1996)
45 Wavelets for seismic tomography of Mantle Many wavelet families exist These will yield different reconstructions Important: the ability to easily switch between different wavelet bases Rather not use a custom-designed wavelet basis for the sphere Preference goes to using cartesian geometry (allows for off-the-shelve wavelets but also curvelets, shearlets etc) Needs a special parametrization of the sphere: cubed sphere (Ronchi et al, 1996)
46 The cubed sphere Sphere is divided into six identical regions (by projecting the sides of a circumscribed cube) All coordinate lines: arcs of great circles 4 subdivisions 8 subdivisions 16 subdivisions In radial dimension: simple subdivisions are used Unwraps to 6 Cartesian chunks with tolerable area distortion
47 The cubed sphere unfolded: coastline Min=0 Max=
48 The cubed sphere unfolded: topography
49 Wavelets on the cubed sphere (CDF 4-2) Min= Max= NB: Basis functions can cross edges smoothly
50 Wavelets on the cubed sphere (CDF 4-2) Min= Max=
51 Spherical harmonics Min= 1 Max=
52 Why wavelets for sparse seismic tomography? A model that is sparse in a wavelet basis is generally smooth but can have some sharp discontinuities Because wavelets are localized in space and scale, sharp discontinuities do not affect all wavelet coefficients Illustration: Developing an existing model and throwing away most of the coefficients still allows for a faithful model
53 An existing Earth model model: 100% nonzero coefficients, err=0%. Min= 1 Max= (Montelli, 2004)
54 An existing Earth model: sparsified Compressed model: 80% zero coefficients, error=3%. Min= Max= (Montelli, 2004)
55 An existing Earth model: sparsified Compressed model: 90% zero coefficients, error=5%. Min= Max= (Montelli, 2004)
56 An existing Earth model: sparsified Compressed model: 95% zero coefficients, error=10%. Min= Max= (Montelli, 2004)
57 An existing Earth model Los Angeles on top. montelli from129.mat Depth: 6371km n=(34 o N, 118 o E) x =20km y =20km (Montelli, 2004)
58 An existing Earth model: sparsified Los Angeles on top. Spars.: 20%, Err=1.967% Depth: 6371km n=(34 o N, 118 o E) x =20km y =20km (Montelli, 2004)
59 An existing Earth model: sparsified Los Angeles on top. Spars.: 10%, Err=3.8699% Depth: 6371km n=(34 o N, 118 o E) x =20km y =20km (Montelli, 2004)
60 An existing Earth model: sparsified Los Angeles on top. Spars.: 5%, Err=6.5745% Depth: 6371km n=(34 o N, 118 o E) x =20km y =20km (Montelli, 2004)
61 An existing Earth model: sparsified Los Angeles on top. Spars.: 2%, Err= % Depth: 6371km n=(34 o N, 118 o E) x =20km y =20km (Montelli, 2004)
62 Opportunity: Exploit model sparsity in reconstruction An accurate model can be obtained using few nonzero wavelet coefficients The location of these few nonzero coefficients is not, and does not have to be specified in advance As wavelets are localized in space, this is an automatic procedure of allowing more degrees of freedom where the data calls for it In contrast: With a standard spherical harmonics based reconstruction, an a priori choice is made of nonzero spherical harmonic coefficients (namely does belonging to the lowest orders)
63 Handling chunk boundaries (surface wave toy problem): Importance of using wavelet functions that cross chunk edges smoothly: Synthetic input model Ray density Min= 1 Min=0 Max=1595 Max= Synthetic output model (incorrect handling of boundaries) Synthetic output model (correct handling of boundaries) Min= Min= Max= Max=
64 Two pillars 1 Expected sparsity of the model: Seismic velocity anomalies can be sparsely described by wavelets Fast algorithms exist for wavelet transforms 2 Sparse recovery algorithms exist: l 1 penalty for recovering sparse solutions to linear equations Functional to be minimized: where Kx y 2 + 2λ u 1 1 x = Earth model 2 K = matrix containing sensitivity kernels 3 y = data 4 u = wavelet coefficients of x; so that x = Su 5 S = suitable basis or frame Together: better reconstructions than usual smoothing regularization
65 Conclusion Sparse reconstruction useful for seismic tomography (model is sparse in wavelet basis) Good parametrization essential Use l 1 penalty on wavelet coefficients to promote sparsity in wavelet basis Now: parallel implementation and inverting real data More: igloris I. Loris and C. Verhoeven (2011). F. J. Simons, I. Loris et al. (2011).
arxiv: v2 [physics.geo-ph] 6 Apr 2012
Iterative algorithms for total variation-like reconstructions in seismic tomography Ignace Loris and Caroline Verhoeven Université Libre de Bruxelles, Brussels, Belgium. August 27, 27 arxiv:23.445v2 [physics.geo-ph]
More informationDeconvolution with curvelet-domain sparsity Vishal Kumar, EOS-UBC and Felix J. Herrmann, EOS-UBC
Deconvolution with curvelet-domain sparsity Vishal Kumar, EOS-UBC and Felix J. Herrmann, EOS-UBC SUMMARY We use the recently introduced multiscale and multidirectional curvelet transform to exploit the
More informationELEG Compressive Sensing and Sparse Signal Representations
ELEG 867 - Compressive Sensing and Sparse Signal Representations Gonzalo R. Arce Depart. of Electrical and Computer Engineering University of Delaware Fall 211 Compressive Sensing G. Arce Fall, 211 1 /
More informationTotal variation tomographic inversion via the Alternating Direction Method of Multipliers
Total variation tomographic inversion via the Alternating Direction Method of Multipliers Landon Safron and Mauricio D. Sacchi Department of Physics, University of Alberta Summary Geophysical inverse problems
More informationIncoherent noise suppression with curvelet-domain sparsity Vishal Kumar, EOS-UBC and Felix J. Herrmann, EOS-UBC
Incoherent noise suppression with curvelet-domain sparsity Vishal Kumar, EOS-UBC and Felix J. Herrmann, EOS-UBC SUMMARY The separation of signal and noise is a key issue in seismic data processing. By
More informationG Practical Magnetic Resonance Imaging II Sackler Institute of Biomedical Sciences New York University School of Medicine. Compressed Sensing
G16.4428 Practical Magnetic Resonance Imaging II Sackler Institute of Biomedical Sciences New York University School of Medicine Compressed Sensing Ricardo Otazo, PhD ricardo.otazo@nyumc.org Compressed
More informationApproximate Constant Density Acoustic Inverse Scattering Using Dip-Dependent Scaling Rami Nammour and William Symes, Rice University
Approximate Constant Density Acoustic Inverse Scattering Using Dip-Dependent Scaling Rami Nammour and William Symes, Rice University SUMMARY This abstract presents a computationally efficient method to
More informationIntroduction to ANSYS DesignXplorer
Lecture 4 14. 5 Release Introduction to ANSYS DesignXplorer 1 2013 ANSYS, Inc. September 27, 2013 s are functions of different nature where the output parameters are described in terms of the input parameters
More informationIterative CT Reconstruction Using Curvelet-Based Regularization
Iterative CT Reconstruction Using Curvelet-Based Regularization Haibo Wu 1,2, Andreas Maier 1, Joachim Hornegger 1,2 1 Pattern Recognition Lab (LME), Department of Computer Science, 2 Graduate School in
More informationRobust Seismic Image Amplitude Recovery Using Curvelets
Robust Seismic Image Amplitude Recovery Using Curvelets Peyman P. Moghaddam Joint work with Felix J. Herrmann and Christiaan C. Stolk UNIVERSITY OF BRITISH COLUMBIA University of Twente References Curvelets
More informationLarge-scale workflows for wave-equation based inversion in Julia
Large-scale workflows for wave-equation based inversion in Julia Philipp A. Witte, Mathias Louboutin and Felix J. Herrmann SLIM University of British Columbia Motivation Use Geophysics to understand the
More informationOutline Introduction Problem Formulation Proposed Solution Applications Conclusion. Compressed Sensing. David L Donoho Presented by: Nitesh Shroff
Compressed Sensing David L Donoho Presented by: Nitesh Shroff University of Maryland Outline 1 Introduction Compressed Sensing 2 Problem Formulation Sparse Signal Problem Statement 3 Proposed Solution
More information3D Computer Vision. Dense 3D Reconstruction II. Prof. Didier Stricker. Christiano Gava
3D Computer Vision Dense 3D Reconstruction II Prof. Didier Stricker Christiano Gava Kaiserlautern University http://ags.cs.uni-kl.de/ DFKI Deutsches Forschungszentrum für Künstliche Intelligenz http://av.dfki.de
More informationCompressed Sensing for Rapid MR Imaging
Compressed Sensing for Rapid Imaging Michael Lustig1, Juan Santos1, David Donoho2 and John Pauly1 1 Electrical Engineering Department, Stanford University 2 Statistics Department, Stanford University rapid
More informationRandomized sampling strategies
Randomized sampling strategies Felix J. Herrmann SLIM Seismic Laboratory for Imaging and Modeling the University of British Columbia SLIM Drivers & Acquisition costs impediments Full-waveform inversion
More informationCompressed Sensing Reconstructions for Dynamic Contrast Enhanced MRI
1 Compressed Sensing Reconstructions for Dynamic Contrast Enhanced MRI Kevin T. Looby klooby@stanford.edu ABSTRACT The temporal resolution necessary for dynamic contrast enhanced (DCE) magnetic resonance
More informationCollaborative Sparsity and Compressive MRI
Modeling and Computation Seminar February 14, 2013 Table of Contents 1 T2 Estimation 2 Undersampling in MRI 3 Compressed Sensing 4 Model-Based Approach 5 From L1 to L0 6 Spatially Adaptive Sparsity MRI
More informationFull-waveform inversion using seislet regularization a
1 Full-waveform inversion using seislet regularization a a Published in Geophysics, 82, no. 5, A43-A49, (2017) Zhiguang Xue 1, Hejun Zhu 2 and Sergey Fomel 1 ABSTRACT Because of inaccurate, incomplete
More informationSUMMARY. min. m f (m) s.t. m C 1
Regularizing waveform inversion by projections onto convex sets application to the D Chevron synthetic blind-test dataset Bas Peters *, Zhilong Fang *, Brendan Smithyman #, Felix J. Herrmann * * Seismic
More informationReversible Wavelets for Embedded Image Compression. Sri Rama Prasanna Pavani Electrical and Computer Engineering, CU Boulder
Reversible Wavelets for Embedded Image Compression Sri Rama Prasanna Pavani Electrical and Computer Engineering, CU Boulder pavani@colorado.edu APPM 7400 - Wavelets and Imaging Prof. Gregory Beylkin -
More informationLearning Splines for Sparse Tomographic Reconstruction. Elham Sakhaee and Alireza Entezari University of Florida
Learning Splines for Sparse Tomographic Reconstruction Elham Sakhaee and Alireza Entezari University of Florida esakhaee@cise.ufl.edu 2 Tomographic Reconstruction Recover the image given X-ray measurements
More informationP257 Transform-domain Sparsity Regularization in Reconstruction of Channelized Facies
P257 Transform-domain Sparsity Regularization in Reconstruction of Channelized Facies. azemi* (University of Alberta) & H.R. Siahkoohi (University of Tehran) SUMMARY Petrophysical reservoir properties,
More informationDetecting Burnscar from Hyperspectral Imagery via Sparse Representation with Low-Rank Interference
Detecting Burnscar from Hyperspectral Imagery via Sparse Representation with Low-Rank Interference Minh Dao 1, Xiang Xiang 1, Bulent Ayhan 2, Chiman Kwan 2, Trac D. Tran 1 Johns Hopkins Univeristy, 3400
More informationTomographic reconstruction: the challenge of dark information. S. Roux
Tomographic reconstruction: the challenge of dark information S. Roux Meeting on Tomography and Applications, Politecnico di Milano, 20-22 April, 2015 Tomography A mature technique, providing an outstanding
More informationIntroduction to seismic body waves Tomography
Introduction to seismic body waves Tomography Presented by Han Zhang Prepared by group work of Justin Wilgus and Han Zhang Outlines 1. Velocity derived from wave equation 2. Two commonly used reference
More informationA Patch Prior for Dense 3D Reconstruction in Man-Made Environments
A Patch Prior for Dense 3D Reconstruction in Man-Made Environments Christian Häne 1, Christopher Zach 2, Bernhard Zeisl 1, Marc Pollefeys 1 1 ETH Zürich 2 MSR Cambridge October 14, 2012 A Patch Prior for
More informationG009 Scale and Direction-guided Interpolation of Aliased Seismic Data in the Curvelet Domain
G009 Scale and Direction-guided Interpolation of Aliased Seismic Data in the Curvelet Domain M. Naghizadeh* (University of Alberta) & M. Sacchi (University of Alberta) SUMMARY We propose a robust interpolation
More informationLecture 17 Sparse Convex Optimization
Lecture 17 Sparse Convex Optimization Compressed sensing A short introduction to Compressed Sensing An imaging perspective 10 Mega Pixels Scene Image compression Picture Why do we compress images? Introduction
More information(x, y, z) m 2. (x, y, z) ...] T. m 2. m = [m 1. m 3. Φ = r T V 1 r + λ 1. m T Wm. m T L T Lm + λ 2. m T Hm + λ 3. t(x, y, z) = m 1
Class 1: Joint Geophysical Inversions Wed, December 1, 29 Invert multiple types of data residuals simultaneously Apply soft mutual constraints: empirical, physical, statistical Deal with data in the same
More informationA NUMERICAL EXPERIMENT ON NONLINEAR IMAGE RECONSTRUCTION FROM FIRST-ARRIVAL TIMES FOR TWO-DIMENSIONAL ISLAND ARC STRUCTURE
Letter J. Phys. Earth, 34, 195-201, 1986 A NUMERICAL EXPERIMENT ON NONLINEAR IMAGE RECONSTRUCTION FROM FIRST-ARRIVAL TIMES FOR TWO-DIMENSIONAL ISLAND ARC STRUCTURE Ichiro NAKANISHI and Kaoru YAMAGUCHI
More informationEffects of multi-scale velocity heterogeneities on wave-equation migration Yong Ma and Paul Sava, Center for Wave Phenomena, Colorado School of Mines
Effects of multi-scale velocity heterogeneities on wave-equation migration Yong Ma and Paul Sava, Center for Wave Phenomena, Colorado School of Mines SUMMARY Velocity models used for wavefield-based seismic
More informationSparse Reconstruction / Compressive Sensing
Sparse Reconstruction / Compressive Sensing Namrata Vaswani Department of Electrical and Computer Engineering Iowa State University Namrata Vaswani Sparse Reconstruction / Compressive Sensing 1/ 20 The
More informationCurves and Surfaces Computer Graphics I Lecture 9
15-462 Computer Graphics I Lecture 9 Curves and Surfaces Parametric Representations Cubic Polynomial Forms Hermite Curves Bezier Curves and Surfaces [Angel 10.1-10.6] February 19, 2002 Frank Pfenning Carnegie
More informationInterval velocity estimation through convex optimization
Stanford Exploration Project, Report 125, January 16, 2007, pages 1?? Interval velocity estimation through convex optimization Ben Witten and Michael Grant ABSTRACT Convex optimization is an optimization
More informationEfficient 3D Gravity and Magnetic Modeling
Efficient 3D Gravity and Magnetic Modeling X. Li Fugro Gravity & Magnetic Services Inc., Houston, Texas, USA Summary There are many different spatial-domain algorithms for 3D gravity and magnetic forward
More informationGEOPHYS 242: Near Surface Geophysical Imaging. Class 8: Joint Geophysical Inversions Wed, April 20, 2011
GEOPHYS 4: Near Surface Geophysical Imaging Class 8: Joint Geophysical Inversions Wed, April, 11 Invert multiple types of data residuals simultaneously Apply soft mutual constraints: empirical, physical,
More informationThe Benefit of Tree Sparsity in Accelerated MRI
The Benefit of Tree Sparsity in Accelerated MRI Chen Chen and Junzhou Huang Department of Computer Science and Engineering, The University of Texas at Arlington, TX, USA 76019 Abstract. The wavelet coefficients
More informationMain Menu. Summary. sampled) f has a sparse representation transform domain S with. in certain. f S x, the relation becomes
Preliminary study on Dreamlet based compressive sensing data recovery Ru-Shan Wu*, Yu Geng 1 and Lingling Ye, Modeling and Imaging Lab, Earth & Planetary Sciences/IGPP, University of California, Santa
More informationOne-norm regularized inversion: learning from the Pareto curve
One-norm regularized inversion: learning from the Pareto curve Gilles Hennenfent and Felix J. Herrmann, Earth & Ocean Sciences Dept., the University of British Columbia ABSTRACT Geophysical inverse problems
More informationAdvanced phase retrieval: maximum likelihood technique with sparse regularization of phase and amplitude
Advanced phase retrieval: maximum likelihood technique with sparse regularization of phase and amplitude A. Migukin *, V. atkovnik and J. Astola Department of Signal Processing, Tampere University of Technology,
More informationCompressed Sensing for Electron Tomography
University of Maryland, College Park Department of Mathematics February 10, 2015 1/33 Outline I Introduction 1 Introduction 2 3 4 2/33 1 Introduction 2 3 4 3/33 Tomography Introduction Tomography - Producing
More informationTRAVELTIME TOMOGRAPHY (CONT)
30 March 005 MODE UNIQUENESS The forward model TRAVETIME TOMOGRAPHY (CONT di = A ik m k d = Am (1 Data vecto r Sensitivit y matrix Model vector states the linearized relationship between data and model
More informationSparse sampling in MRI: From basic theory to clinical application. R. Marc Lebel, PhD Department of Electrical Engineering Department of Radiology
Sparse sampling in MRI: From basic theory to clinical application R. Marc Lebel, PhD Department of Electrical Engineering Department of Radiology Objective Provide an intuitive overview of compressed sensing
More informationECE 8201: Low-dimensional Signal Models for High-dimensional Data Analysis
ECE 8201: Low-dimensional Signal Models for High-dimensional Data Analysis Yuejie Chi Departments of ECE and BMI The Ohio State University September 24, 2015 Time, location, and office hours Time: Tue/Thu
More informationInversion after depth imaging
Robin P. Fletcher *, Stewart Archer, Dave Nichols, and Weijian Mao, WesternGeco Summary In many areas, depth imaging of seismic data is required to construct an accurate view of the reservoir structure.
More informationPerception Viewed as an Inverse Problem
Perception Viewed as an Inverse Problem Fechnerian Causal Chain of Events - an evaluation Fechner s study of outer psychophysics assumes that the percept is a result of a causal chain of events: Distal
More informationMultiple-View Object Recognition in Band-Limited Distributed Camera Networks
in Band-Limited Distributed Camera Networks Allen Y. Yang, Subhransu Maji, Mario Christoudas, Kirak Hong, Posu Yan Trevor Darrell, Jitendra Malik, and Shankar Sastry Fusion, 2009 Classical Object Recognition
More informationTotal Variation Denoising with Overlapping Group Sparsity
1 Total Variation Denoising with Overlapping Group Sparsity Ivan W. Selesnick and Po-Yu Chen Polytechnic Institute of New York University Brooklyn, New York selesi@poly.edu 2 Abstract This paper describes
More informationWriting Kirchhoff migration/modelling in a matrix form
Writing Kirchhoff migration/modelling in a matrix form Abdolnaser Yousefzadeh and John C. Bancroft Kirchhoff migration matrix ABSTRACT Kirchhoff prestack migration and modelling are linear operators. Therefore,
More informationCurves and Surfaces Computer Graphics I Lecture 10
15-462 Computer Graphics I Lecture 10 Curves and Surfaces Parametric Representations Cubic Polynomial Forms Hermite Curves Bezier Curves and Surfaces [Angel 10.1-10.6] September 30, 2003 Doug James Carnegie
More informationSeismic velocity tomography with co-located soft data
Seismic velocity tomography with co-located soft data Bob Clapp Stanford University SCRF Affiliate Meeting May 2009 Motivation soft data constrains flow simulation model 2 Outline Introduction Seismic
More informationNoise-thresholding sparse-spike inversion with global convergence: calibration and applications
Noise-thresholding sparse-spike inversion with global convergence: calibration and applications Douglas S. Sassen and Michael Glinsky, Ion Geophysical 2013 SEG Douglas Sassen 9/25/13 1 Illustrative Abstract!
More informationIntroduction to Topics in Machine Learning
Introduction to Topics in Machine Learning Namrata Vaswani Department of Electrical and Computer Engineering Iowa State University Namrata Vaswani 1/ 27 Compressed Sensing / Sparse Recovery: Given y :=
More informationCombinatorial Selection and Least Absolute Shrinkage via The CLASH Operator
Combinatorial Selection and Least Absolute Shrinkage via The CLASH Operator Volkan Cevher Laboratory for Information and Inference Systems LIONS / EPFL http://lions.epfl.ch & Idiap Research Institute joint
More informationRandomized sampling without repetition in timelapse seismic surveys
Randomized sampling without repetition in timelapse seismic surveys Felix Oghenekohwo SLIM University of British Columbia Randomized sampling without repetition in timelapse seismic surveys Felix Oghenekohwo
More information1. Introduction. performance of numerical methods. complexity bounds. structural convex optimization. course goals and topics
1. Introduction EE 546, Univ of Washington, Spring 2016 performance of numerical methods complexity bounds structural convex optimization course goals and topics 1 1 Some course info Welcome to EE 546!
More informationLecture VIII. Global Approximation Methods: I
Lecture VIII Global Approximation Methods: I Gianluca Violante New York University Quantitative Macroeconomics G. Violante, Global Methods p. 1 /29 Global function approximation Global methods: function
More informationBilevel Sparse Coding
Adobe Research 345 Park Ave, San Jose, CA Mar 15, 2013 Outline 1 2 The learning model The learning algorithm 3 4 Sparse Modeling Many types of sensory data, e.g., images and audio, are in high-dimensional
More informationTime-jittered ocean bottom seismic acquisition Haneet Wason and Felix J. Herrmann
Time-jittered ocean bottom seismic acquisition Haneet Wason and Felix J. Herrmann SLIM University of British Columbia Challenges Need for full sampling - wave-equation based inversion (RTM & FWI) - SRME/EPSI
More informationIteratively Re-weighted Least Squares for Sums of Convex Functions
Iteratively Re-weighted Least Squares for Sums of Convex Functions James Burke University of Washington Jiashan Wang LinkedIn Frank Curtis Lehigh University Hao Wang Shanghai Tech University Daiwei He
More informationA Scalable Parallel LSQR Algorithm for Solving Large-Scale Linear System for Seismic Tomography
1 A Scalable Parallel LSQR Algorithm for Solving Large-Scale Linear System for Seismic Tomography He Huang, Liqiang Wang, Po Chen(University of Wyoming) John Dennis (NCAR) 2 LSQR in Seismic Tomography
More informationLeast squares Kirchhoff depth migration: important details
Least squares Kirchhoff depth migration: important details Daniel Trad CREWES-University of Calgary Summary Least squares migration has been an important research topic in the academia for about two decades,
More informationVariational Geometric Modeling with Wavelets
Variational Geometric Modeling with Wavelets Steven J. Gortler and Michael F. Cohen Microsoft Research Redmond, WA (excerpted from Hierarchical and Variational Geometric Modeling with Wavelets, by Steven
More informationHigher Degree Total Variation for 3-D Image Recovery
Higher Degree Total Variation for 3-D Image Recovery Greg Ongie*, Yue Hu, Mathews Jacob Computational Biomedical Imaging Group (CBIG) University of Iowa ISBI 2014 Beijing, China Motivation: Compressed
More informationA fast iterative thresholding algorithm for wavelet-regularized deconvolution
A fast iterative thresholding algorithm for wavelet-regularized deconvolution Cédric Vonesch and Michael Unser Biomedical Imaging Group, EPFL, Lausanne, Switzerland ABSTRACT We present an iterative deconvolution
More informationComputerLab: compressive sensing and application to MRI
Compressive Sensing, 207-8 ComputerLab: compressive sensing and application to MRI Aline Roumy This computer lab addresses the implementation and analysis of reconstruction algorithms for compressive sensing.
More informationImage Transformation Techniques Dr. Rajeev Srivastava Dept. of Computer Engineering, ITBHU, Varanasi
Image Transformation Techniques Dr. Rajeev Srivastava Dept. of Computer Engineering, ITBHU, Varanasi 1. Introduction The choice of a particular transform in a given application depends on the amount of
More informationReconstruction of CT Images from Sparse-View Polyenergetic Data Using Total Variation Minimization
1 Reconstruction of CT Images from Sparse-View Polyenergetic Data Using Total Variation Minimization T. Humphries and A. Faridani Abstract Recent work in CT image reconstruction has seen increasing interest
More informationAdaptive Waveform Inversion: Theory Mike Warner*, Imperial College London, and Lluís Guasch, Sub Salt Solutions Limited
Adaptive Waveform Inversion: Theory Mike Warner*, Imperial College London, and Lluís Guasch, Sub Salt Solutions Limited Summary We present a new method for performing full-waveform inversion that appears
More informationGeophysical Journal International
Geophysical Journal International Geophys. J. Int. (2012) 191, 1391 1402 doi: 10.1111/j.1365-246X.2012.05686.x Optimized discrete wavelet transforms in the cubed sphere with the lifting scheme implications
More informationCompressive Sensing Algorithms for Fast and Accurate Imaging
Compressive Sensing Algorithms for Fast and Accurate Imaging Wotao Yin Department of Computational and Applied Mathematics, Rice University SCIMM 10 ASU, Tempe, AZ Acknowledgements: results come in part
More informationOutline. Level Set Methods. For Inverse Obstacle Problems 4. Introduction. Introduction. Martin Burger
For Inverse Obstacle Problems Martin Burger Outline Introduction Optimal Geometries Inverse Obstacle Problems & Shape Optimization Sensitivity Analysis based on Gradient Flows Numerical Methods University
More informationSparsity and image processing
Sparsity and image processing Aurélie Boisbunon INRIA-SAM, AYIN March 6, Why sparsity? Main advantages Dimensionality reduction Fast computation Better interpretability Image processing pattern recognition
More informationTarget-oriented wave-equation inversion
Stanford Exploration Project, Report 120, May 3, 2005, pages 23 40 Target-oriented wave-equation inversion Alejandro A. Valenciano, Biondo Biondi, and Antoine Guitton 1 ABSTRACT A target-oriented strategy
More informationAim. Structure and matrix sparsity: Part 1 The simplex method: Exploiting sparsity. Structure and matrix sparsity: Overview
Aim Structure and matrix sparsity: Part 1 The simplex method: Exploiting sparsity Julian Hall School of Mathematics University of Edinburgh jajhall@ed.ac.uk What should a 2-hour PhD lecture on structure
More informationOptimal Sampling Geometries for TV-Norm Reconstruction of fmri Data
Optimal Sampling Geometries for TV-Norm Reconstruction of fmri Data Oliver M. Jeromin, Student Member, IEEE, Vince D. Calhoun, Senior Member, IEEE, and Marios S. Pattichis, Senior Member, IEEE Abstract
More informationEfficient Beam Velocity Model Building with Tomography Designed to Accept 3d Residuals Aligning Depth Offset Gathers
Efficient Beam Velocity Model Building with Tomography Designed to Accept 3d Residuals Aligning Depth Offset Gathers J.W.C. Sherwood* (PGS), K. Sherwood (PGS), H. Tieman (PGS), R. Mager (PGS) & C. Zhou
More informationAdaptive Reconstruction Methods for Low-Dose Computed Tomography
Adaptive Reconstruction Methods for Low-Dose Computed Tomography Joseph Shtok Ph.D. supervisors: Prof. Michael Elad, Dr. Michael Zibulevsky. Technion IIT, Israel, 011 Ph.D. Talk, Apr. 01 Contents of this
More informationLasso. November 14, 2017
Lasso November 14, 2017 Contents 1 Case Study: Least Absolute Shrinkage and Selection Operator (LASSO) 1 1.1 The Lasso Estimator.................................... 1 1.2 Computation of the Lasso Solution............................
More informationSome Blind Deconvolution Techniques in Image Processing
Some Blind Deconvolution Techniques in Image Processing Tony Chan Math Dept., UCLA Joint work with Frederick Park and Andy M. Yip IPAM Workshop on Mathematical Challenges in Astronomical Imaging July 26-30,
More informationMain Menu. Summary. Introduction
Noise-thresholding sparse-spike inversion with global convergence: calibration and applications Douglas S. Sassen* and Michael Glinsky, Ion Geophysical Corp. Summary We demonstrate an innovative inversion
More informationAdvanced Operations Research Techniques IE316. Quiz 1 Review. Dr. Ted Ralphs
Advanced Operations Research Techniques IE316 Quiz 1 Review Dr. Ted Ralphs IE316 Quiz 1 Review 1 Reading for The Quiz Material covered in detail in lecture. 1.1, 1.4, 2.1-2.6, 3.1-3.3, 3.5 Background material
More informationConvex Optimization. Stephen Boyd
Convex Optimization Stephen Boyd Electrical Engineering Computer Science Management Science and Engineering Institute for Computational Mathematics & Engineering Stanford University Institute for Advanced
More informationDiffuse Optical Tomography, Inverse Problems, and Optimization. Mary Katherine Huffman. Undergraduate Research Fall 2011 Spring 2012
Diffuse Optical Tomography, Inverse Problems, and Optimization Mary Katherine Huffman Undergraduate Research Fall 11 Spring 12 1. Introduction. This paper discusses research conducted in order to investigate
More informationIntroduction. Wavelets, Curvelets [4], Surfacelets [5].
Introduction Signal reconstruction from the smallest possible Fourier measurements has been a key motivation in the compressed sensing (CS) research [1]. Accurate reconstruction from partial Fourier data
More informationSparse coding for image classification
Sparse coding for image classification Columbia University Electrical Engineering: Kun Rong(kr2496@columbia.edu) Yongzhou Xiang(yx2211@columbia.edu) Yin Cui(yc2776@columbia.edu) Outline Background Introduction
More informationFace Recognition via Sparse Representation
Face Recognition via Sparse Representation John Wright, Allen Y. Yang, Arvind, S. Shankar Sastry and Yi Ma IEEE Trans. PAMI, March 2008 Research About Face Face Detection Face Alignment Face Recognition
More informationIMAGE DE-NOISING IN WAVELET DOMAIN
IMAGE DE-NOISING IN WAVELET DOMAIN Aaditya Verma a, Shrey Agarwal a a Department of Civil Engineering, Indian Institute of Technology, Kanpur, India - (aaditya, ashrey)@iitk.ac.in KEY WORDS: Wavelets,
More informationSource Estimation for Wavefield Reconstruction Inversion
Source Estimation for Wavefield Reconstruction Inversion Zhilong Fang * and Felix J. Herrmann * * Seismic Laboratory for Imaging and Modeling SLIM), University of British Columbia Abstract Wavefield reconstruction
More informationNonparametric Methods Recap
Nonparametric Methods Recap Aarti Singh Machine Learning 10-701/15-781 Oct 4, 2010 Nonparametric Methods Kernel Density estimate (also Histogram) Weighted frequency Classification - K-NN Classifier Majority
More informationMostafa Naghizadeh and Mauricio D. Sacchi
Ground-roll elimination by scale and direction guided curvelet transform Mostafa Naghizadeh and Mauricio D. Sacchi Summary We propose a curvelet domain strategy to subtract ground-roll from seismic records.
More informationFull-waveform inversion for reservoir characterization: A synthetic study
CWP-889 Full-waveform inversion for reservoir characterization: A synthetic study Nishant Kamath, Ilya Tsvankin & Ehsan Zabihi Naeini ABSTRACT Most current reservoir-characterization workflows are based
More informationMotion Estimation. There are three main types (or applications) of motion estimation:
Members: D91922016 朱威達 R93922010 林聖凱 R93922044 謝俊瑋 Motion Estimation There are three main types (or applications) of motion estimation: Parametric motion (image alignment) The main idea of parametric motion
More informationE044 Ray-based Tomography for Q Estimation and Q Compensation in Complex Media
E044 Ray-based Tomography for Q Estimation and Q Compensation in Complex Media M. Cavalca* (WesternGeco), I. Moore (WesternGeco), L. Zhang (WesternGeco), S.L. Ng (WesternGeco), R.P. Fletcher (WesternGeco)
More informationApplications of the generalized norm solver
Applications of the generalized norm solver Mandy Wong, Nader Moussa, and Mohammad Maysami ABSTRACT The application of a L1/L2 regression solver, termed the generalized norm solver, to two test cases,
More informationRevision of the SolidWorks Variable Pressure Simulation Tutorial J.E. Akin, Rice University, Mechanical Engineering. Introduction
Revision of the SolidWorks Variable Pressure Simulation Tutorial J.E. Akin, Rice University, Mechanical Engineering Introduction A SolidWorks simulation tutorial is just intended to illustrate where to
More informationInverse Continuous Wavelet Transform Deconvolution Summary ICWT deconvolution Introduction
Marcilio Castro de Matos*, Sismo Research&Consulting and AASPI/OU, and Kurt J. Marfurt, The University of Oklahoma (OU) Summary Most deconvolution algorithms try to transform the seismic wavelet into spikes
More information03 - Reconstruction. Acknowledgements: Olga Sorkine-Hornung. CSCI-GA Geometric Modeling - Spring 17 - Daniele Panozzo
3 - Reconstruction Acknowledgements: Olga Sorkine-Hornung Geometry Acquisition Pipeline Scanning: results in range images Registration: bring all range images to one coordinate system Stitching/ reconstruction:
More informationAlgorithms for convex optimization
Algorithms for convex optimization Michal Kočvara Institute of Information Theory and Automation Academy of Sciences of the Czech Republic and Czech Technical University kocvara@utia.cas.cz http://www.utia.cas.cz/kocvara
More informationLecture 2: Seismic Inversion and Imaging
Lecture 2: Seismic Inversion and Imaging Youzuo Lin 1 Joint work with: Lianjie Huang 1 Monica Maceira 2 Ellen M. Syracuse 1 Carene Larmat 1 1: Los Alamos National Laboratory 2: Oak Ridge National Laboratory
More information