@mastersthesis {hennenfent08phd,
	title = {Sampling and reconstruction of seismic wavefields in the curvelet domain},
	year = {2008},
	note = {(PhD)},
	month = {05},
	school = {The University of British Columbia},
	type = {phd},
	address = {Vancouver},
	abstract = {Wavefield reconstruction is a crucial step in the seismic processing flow. For instance, unsuccessful interpolation leads to erroneous multiple predictions that adversely affect the performance of multiple elimination, and to imaging artifacts. We present a new non-parametric transform-based reconstruction method that exploits the compression of seismic data by the recently developed curvelet transform. The elements of this transform, called curvelets, are multi-dimensional, multi-scale, and multi-directional. They locally resemble wavefronts present in the data, which leads to a compressible representation for seismic data. This compression enables us to formulate a new curvelet-based seismic data recovery algorithm through sparsity-promoting inversion (CRSI). The concept of sparsity-promoting inversion is in itself not new to geophysics. However, the recent insights from the field of {\textquoteleft}{\textquoteleft}compressed sensing{\textquoteright}{\textquoteright} are new since they clearly identify the three main ingredients that go into a successful formulation of a reconstruction problem, namely a sparsifying transform, a sub-Nyquist sampling strategy that subdues coherent aliases in the sparsifying domain, and a data-consistent sparsity-promoting program. After a brief overview of the curvelet transform and our seismic-oriented extension to the fast discrete curvelet transform, we detail the CRSI formulation and illustrate its performance on synthetic and real datasets. Then, we introduce a sub-Nyquist sampling scheme, termed jittered undersampling, and show that, for the same amount of data acquired, jittered data are best interpolated using CRSI compared to regular or random undersampled data. We also discuss the large-scale one-norm solver involved in CRSI. Finally, we extend CRSI formulation to other geophysical applications and present results on multiple removal and migration-amplitude recovery.},
	keywords = {curvelet transform, reconstruction, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Thesis/2008/hennenfent08phd.pdf},
	author = {Gilles Hennenfent}
}
@article {herrmann2008ACHAsac,
	title = {Sparsity- and continuity-promoting seismic image recovery with curvelet frames},
	journal = {Applied and Computational Harmonic Analysis},
	volume = {24},
	number = {2},
	year = {2008},
	month = {03},
	pages = {150-173},
	abstract = {A nonlinear singularity-preserving solution to seismic image recovery with sparseness and continuity constraints is proposed. We observe that curvelets, as a directional frame expansion, lead to sparsity of seismic images and exhibit invariance under the normal operator of the linearized imaging problem. Based on this observation we derive a method for stable recovery of the migration amplitudes from noisy data. The method corrects the amplitudes during a post-processing step after migration, such that the main additional cost is one ap- plication of the normal operator, i.e. a modeling followed by a migration. Asymptotically this normal operator corresponds to a pseudodifferential operator, for which a convenient diagonal approximation in the curvelet domain is derived, including a bound for its error and a method for the estimation of the diagonal from a compound operator consisting of discrete implementations for the scattering operator and its adjoint the migration operator. The solution is formulated as a nonlinear optimization problem where sparsity in the curvelet domain as well as continuity along the imaged reflectors are jointly promoted. To enhance sparsity, the $ell_1$-norm on the curvelet coefficients is minimized, while continuity is promoted by minimizing an anisotropic diffusion norm on the image. The performance of the recovery scheme is evaluated with a time-reversed {\textquoteright}wave-equation{\textquoteright} migration code on synthetic datasets, including the complex SEG/EAGE AA salt model.},
	keywords = {curvelet transform, Imaging, Processing, SLIM},
	doi = {10.1016/j.acha.2007.06.007},
	url = {https://slim.gatech.edu/Publications/Public/Journals/ACHA/2008/herrmann2008ACHAsac/herrmann2008ACHAsac.pdf},
	author = {Felix J. Herrmann and Peyman P. Moghaddam and Christiaan C. Stolk}
}
