This paper introduces a Parseval frame of shearlets for the representation of threedimensional (3D) data, which is especially designed to handle geometric features such as discontinuous boundaries with very high efficiency. This system of 3D shearlets forms a multiscale pyramid of well-localized waveforms at various locations and orientations, which become increasingly thin and elongated at fine scales. We prove that this 3D shearlet construction provides essentially optimal sparse representations for functions on R3 which are C 2-regular away from discontinuities along C 2 surfaces. As a consequence, we show that, within this class of functions, the N-term approximation f SN obtained by selecting the N largest coefficients of the shearlet expansion of f satisfies the asymptotic estimate ||f-f SN ||22 N-1(logN)2 as N 8. This asymptotic behavior significantly outperforms wavelet and Fourier series approximations, which yield an approximation rate of only O(N-1/2) and O(N-1/3), respectively. This result extends to the 3D setting the (essentially) optimally sparse approximation results obtained by the authors using 2D shearlets and by Candès and Donoho using curvelets and is the first nonadaptive construction to provide provably (nearly) optimal representations for a large class of 3D data.



Document Type


Additional Information

The authors acknowledge support from NSF grant DMS 1008900/1008907; Demetrio Labate also acknowledges support from NSF grant DMS 1005799.



Rights Information

© Society for Industrial and Applied Mathematics


Affine systems, Curvelets, Nonlinear approximations, Shearlets, Sparsity, Wavelets

Publication Date


Journal Title

SIAM Journal on Mathematical Analysis