@article {vandenberg2009ACMstf,
	title = {Sparco: a testing framework for sparse reconstruction},
	journal = {ACM Transactions on Mathematical Software},
	volume = {35},
	number = {4},
	year = {2009},
	month = {02},
	pages = {1-16},
	abstract = {Sparco is a framework for testing and benchmarking algorithms for sparse reconstruction. It includes a large collection of sparse reconstruction problems drawn from the imaging, compressed sensing, and geophysics literature. Sparco is also a framework for implementing new test problems and can be used as a tool for reproducible research. Sparco is implemented entirely in Matlab, and is released as open-source software under the GNU Public License.},
	keywords = {compressed sensing, linear operators, Sparse recovery},
	url = {http://doi.acm.org/10.1145/1462173.1462178},
	author = {Ewout van den Berg and Michael P. Friedlander and Gilles Hennenfent and Felix J. Herrmann and Rayan Saab and Ozgur Yilmaz}
}
@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}
}
@conference {johnson2008SINBADsdi,
	title = {Seismic data interpolation with symmetry},
	booktitle = {SINBAD 2008},
	year = {2008},
	abstract = {Due to the physics of reciprocity seismic data sets are symmetric in the source and receiver coordinates. Often seismic data sets are incomplete and the missing data must be interpolated. Typically, missing traces do not occur symmetrically. The purpose of this project is to extend the current formulation for solving the seismic interpolation problems in such a way that they enforce reciprocity. The method decomposes the seismic data volume into symmetric and antisymmetric parts. This decomposition leads to an augmented system of equations for the L1-solver that promotes sparsity in the curvelet domain. Interpolation is carried out on the entire system during which the asymmetric component of the volume is forced to zero, while the symmetric part of the data volume is matched to the measured data.},
	keywords = {Presentation, SINBAD, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2008/johnson2008SINBADsdi/johnson2008SINBADsdi.pdf},
	author = {James Johnson and Gilles Hennenfent}
}
@article {hennenfent2008GEOPsdw,
	title = {Simply denoise: wavefield reconstruction via jittered undersampling},
	journal = {Geophysics},
	volume = {73},
	number = {3},
	year = {2008},
	month = {05},
	pages = {V19-V28},
	publisher = {SEG},
	abstract = {In this paper, we present a new discrete undersampling scheme designed to favor wavefield reconstruction by sparsity-promoting inversion with transform elements that are localized in the Fourier domain. Our work is motivated by empirical observations in the seismic community, corroborated by recent results from compressive sampling, which indicate favorable (wavefield) reconstructions from random as opposed to regular undersampling. As predicted by theory, random undersampling renders coherent aliases into harmless incoherent random noise, effectively turning the interpolation problem into a much simpler denoising problem. A practical requirement of wavefield reconstruction with localized sparsifying transforms is the control on the maximum gap size. Unfortunately, random undersampling does not provide such a control and the main purpose of this paper is to introduce a sampling scheme, coined jittered undersampling, that shares the benefits of random sampling, while offering control on the maximum gap size. Our contribution of jittered sub-Nyquist sampling proves to be key in the formulation of a versatile wavefield sparsity-promoting recovery scheme that follows the principles of compressive sampling. After studying the behavior of the jittered undersampling scheme in the Fourier domain, its performance is studied for curvelet recovery by sparsity-promoting inversion (CRSI). Our findings on synthetic and real seismic data indicate an improvement of several decibels over recovery from regularly-undersampled data for the same amount of data collected.},
	keywords = {Acquisition, Compressive Sensing, Geophysics, Optimization, Processing, sampling, SLIM},
	doi = {10.1190/1.2841038},
	url = {https://slim.gatech.edu/Publications/Public/Journals/Geophysics/2008/hennenfent08GEOsdw/hennenfent08GEOsdw.pdf},
	html_version = {https://slim.gatech.edu/Publications/Public/Journals/Geophysics/2008/hennenfent08GEOsdw/paper_html/paper.html},
	author = {Gilles Hennenfent and Felix J. Herrmann}
}
@conference {hennenfent2008SINBADsdw2,
	title = {Simply denoise: wavefield reconstruction via jittered undersampling},
	booktitle = {SINBAD 2008},
	year = {2008},
	abstract = {We present a new discrete undersampling scheme designed to favor wavefield reconstruction by sparsity-promoting inversion with transform elements that are localized in the Fourier domain. Our work is motivated by empirical observations in the seismic community, corroborated by recent results from compressive sampling, which indicate favorable (wavefield) reconstructions from random as opposed to regular undersampling. As predicted by theory, random undersampling renders coherent aliases into harmless incoherent random noise, effectively turning the interpolation problem into a much simpler denoising problem. A practical requirement of wavefield reconstruction with localized sparsifying transforms is the control on the maximum gap size. Unfortunately, random undersampling does not provide such a control and the main purpose of this paper is to introduce a sampling scheme, coined jittered undersampling, that shares the benefits of random sampling, while offering control on the maximum gap size. Our contribution of jittered sub-Nyquist sampling proofs to be key in the formulation of a versatile wavefield sparsity-promoting recovery scheme that follows the principles of compressive sampling. After studying the behavior of the jittered-undersampling scheme in the Fourier domain, its performance is studied for curvelet recovery by sparsity-promoting inversion (CRSI). Our findings on synthetic and real seismic data indicate an improvement of several decibels over recovery from regularly-undersampled data for the same amount of data collected.},
	keywords = {Presentation, SINBAD, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2008/hennenfent2008SINBADsdw2/hennenfent2008SINBADsdw2.pdf},
	author = {Gilles Hennenfent}
}
@conference {herrmann2007SEGsdp,
	title = {Seismic data processing with curvelets: a multiscale and nonlinear approach},
	booktitle = {SEG Technical Program Expanded Abstracts},
	volume = {26},
	year = {2007},
	pages = {2220-2224},
	publisher = {SEG},
	organization = {SEG},
	abstract = {In this abstract, we present a nonlinear curvelet-based sparsity-promoting formulation of a seismic processing flow, consisting of the following steps: seismic data regularization and the restoration of migration amplitudes. We show that the curvelet{\textquoteright}s wavefront detection capability and invariance under the migration-demigration operator lead to a formulation that is stable under noise and missing data. \copyright2007 Society of Exploration Geophysicists},
	keywords = {Presentation, SEG, SLIM},
	doi = {10.1190/1.2792927},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2007/herrmann07SEGsdp/herrmann07SEGsdp.pdf},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2007/herrmann07SEGsdp/herrmann07SEGsdp_pres.pdf},
	author = {Felix J. Herrmann and Deli Wang and Gilles Hennenfent and Peyman P. Moghaddam}
}
@conference {herrmann2007EAGEsia,
	title = {Seismic imaging and processing with curvelets},
	booktitle = {EAGE Annual Conference Proceedings},
	year = {2007},
	month = {06},
	abstract = {In this paper, we present a nonlinear curvelet-based sparsity-promoting formulation for three problems in seismic processing and imaging namely, seismic data regularization from data with large percentages of traces missing; seismic amplitude recovery for sub-salt images obtained by reverse-time migration and primary-multiple separation, given an inaccurate multiple prediction. We argue why these nonlinear formulations are beneficial.},
	keywords = {EAGE, Presentation, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/EAGE/2007/herrmann07EAGEsia/herrmann07EAGEsia.pdf},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/EAGE/2007/herrmann07EAGEsia/herrmann07EAGEsia_pres.pdf},
	url2 = {http://earthdoc.eage.org/publication/publicationdetails/?publication=7075},
	author = {Felix J. Herrmann and Gilles Hennenfent and Peyman P. Moghaddam}
}
@conference {challa2007EAGEsrf,
	title = {Signal reconstruction from incomplete and misplaced measurements},
	booktitle = {EAGE Annual Conference Proceedings},
	year = {2007},
	month = {06},
	abstract = {Constrained by practical and economical considerations, one often uses seismic data with missing traces. The use of such data results in image artifacts and poor spatial resolution. Sometimes due to practical limitations, measurements may be available on a perturbed grid, instead of on the designated grid. Due to algorithmic requirements, when such measurements are viewed as those on the designated grid, the recovery procedures may result in additional artifacts. This paper interpolates incomplete data onto regular grid via the Fourier domain, using a recently developed greedy algorithm. The basic objective is to study experimentally as to what could be the size of the perturbation in measurement coordinates that allows for the measurements on the perturbed grid to be considered as on the designated grid for faithful recovery. Our experimental work shows that for compressible signals, a uniformly distributed perturbation can be offset with slightly more number of measurements.},
	keywords = {EAGE, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/EAGE/2007/challa07EAGEsrf/challa07EAGEsrf.pdf},
	url2 = {http://earthdoc.eage.org/publication/publicationdetails/?publication=6917},
	author = {Challa S. Sastry and Gilles Hennenfent and Felix J. Herrmann}
}
@article {hennenfent2006CiSEsdn,
	title = {Seismic denoising with nonuniformly sampled curvelets},
	journal = {Computing in Science \& Engineering},
	volume = {8},
	number = {3},
	year = {2006},
	month = {05},
	pages = {16-25},
	abstract = {The authors present an extension of the fast discrete curvelet transform (FDCT) to nonuniformly sampled data. This extension not only restores curvelet compression rates for nonuniformly sampled data but also removes noise and maps the data to a regular grid.},
	keywords = {CiSE, Processing},
	doi = {10.1109/MCSE.2006.49},
	url = {https://slim.gatech.edu/Publications/Public/Journals/CiSE/2006/hennenfent06CiSEsdn/hennenfent06CiSEsdn.pdf},
	author = {Gilles Hennenfent and Felix J. Herrmann}
}
@conference {hennenfent2005EAGEsdr,
	title = {Seismic deconvolution revisited with curvelet frames},
	booktitle = {EAGE Annual Conference Proceedings},
	year = {2005},
	month = {06},
	abstract = {We propose an efficient iterative curvelet-regularized deconvolution algorithm that exploits continuity along reflectors in seismic images. Curvelets are a new multiscale transform that provides sparse representations for images (such as seismic images) that comprise smooth objects separated by piece-wise smooth discontinuities. Our technique combines conjugate gradient-based convolution operator inversion with noise regularization that is performed using non-linear curvelet coefficient shrinkage (thresholding). The shrinkage operation leverages the sparsity of curvelets representations. Simulations demonstrate that our algorithm provides improved resolution compared to the traditional Wiener-based deconvolution approach.},
	keywords = {EAGE, Presentation, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/EAGE/2005/Hennenfent05EAGEsdr/Hennenfent05EAGEsdr.pdf},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/EAGE/2005/Hennenfent05EAGEsdr/Hennenfent05EAGEsdr_poster.pdf},
	url2 = {http://earthdoc.eage.org/publication/publicationdetails/?publication=1383},
	author = {Gilles Hennenfent and R. Neelamani and Felix J. Herrmann}
}
@conference {hennenfent2005SEGscd,
	title = {Sparseness-constrained data continuation with frames: applications to missing traces and aliased signals in 2/3-D},
	booktitle = {SEG Technical Program Expanded Abstracts},
	year = {2005},
	pages = {2162-2165},
	publisher = {SEG},
	organization = {SEG},
	abstract = {We present a robust iterative sparseness-constrained interpolation algorithm using 2-/3-D curvelet frames and Fourier-like transforms that exploits continuity along reflectors in seismic data. By choosing generic transforms, we circumvent the necessity to make parametric assumptions (e.g. through linear/parabolic Radon or demigration) regarding the shape of events in seismic data. Simulation and real data examples for data with moderately sized gaps demonstrate that our algorithm provides interpolated traces that accurately reproduce the wavelet shape as well as the AVO behavior. Our method also shows good results for de-aliasing judged by the behavior of the ($f-k$)-spectrum before and after regularization. \copyright2005 Society of Exploration Geophysicists},
	keywords = {Presentation, SEG, SLIM},
	doi = {10.1190/1.2148142},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2005/Hennenfent05SEGscd/Hennenfent05SEGscd.pdf},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2005/Hennenfent05SEGscd/Hennenfent05SEGscd_pres.pdf},
	author = {Gilles Hennenfent and Felix J. Herrmann}
}
@conference {hennenfent2005CSEGscs,
	title = {Sparseness-constrained seismic deconvolution with curvelets},
	booktitle = {CSEG Annual Conference Proceedings},
	year = {2005},
	month = {05},
	publisher = {CSEG},
	organization = {CSEG},
	abstract = {Continuity along reflectors in seismic images is used via Curvelet representation to stabilize the convolution operator inversion. The Curvelet transform is a new multiscale transform that provides sparse representations for images that comprise smooth objects separated by piece-wise smooth discontinuities (e.g. seismic images). Our iterative Curvelet-regularized deconvolution algorithm combines conjugate gradient-based inversion with noise regularization performed using non-linear Curvelet coefficient thresholding. The thresholding operation enhances the sparsity of Curvelet representations. We show on a synthetic example that our algorithm provides improved resolution and continuity along reflectors as well as reduced ringing effect compared to the iterative Wiener-based deconvolution approach.},
	keywords = {Presentation, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/CSEG/2005/Hennenfent05CSEGscs/Hennenfent05CSEGscs.pdf},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/CSEG/2005/Hennenfent05CSEGscs/Hennenfent05CSEGscs_pres.pdf},
	author = {Gilles Hennenfent and Felix J. Herrmann and R. Neelamani}
}
