We construct frames of wavepackets produced by parabolic dilation, rotation and translation of (a finite sum of) Gaussians and give asymptotics on the analogue of Daubechies frame criterion. We show that the coefficients in the corresponding approximate expansion decay fast away from the wavefront set of the original data.
While the differences between the various parabolic wavepackets are very relevant in practice, their asymptotic properties are similar. The Gaussian wavepackets in are an example of from @cite_22 and therefore they provide the same sparse approximation properties as curvelets for functions with discontinuities along smooth edges @cite_5 . The notion of curvelet molecules and the related notion of in @cite_4 are in turn particular cases of the notion of , as recently introduced in @cite_10 . In that article it is shown that all expansions with adequate systems of parabolic molecules share similar asymptotic properties. As a consequence, our results on the decay of frame coefficients away from the wavefront set can be transferred to other parabolic systems.
@cite_5 @cite_4 @cite_22 @cite_10
2069912449 2014172718 2168141504
This paper introduces new tight frames of curvelets to address the problem of finding optimally sparse representations of objects with discontinuities along piecewise C 2 edges. Conceptually, the curvelet transform is a multiscale pyramid with many directions and positions at each length scale, and needle-shaped elements at fine scales. These elements have many useful geometric multiscale features that set them apart from classical multiscale representations such as wavelets. For instance, curvelets obey a parabolic scaling relation which says that at scale 2 -j , each element has an envelope that is aligned along a ridge of length 2 -j 2 and width 2 -j . We prove that curvelets provide an essentially optimal representation of typical objects f that are C 2 except for discontinuities along piecewise C 2 curves. Such representations are nearly as sparse as if f were not singular and turn out to be far more sparse than the wavelet decomposition of the object. For instance, the n-term partial reconstruction f C n obtained by selecting the n largest terms in the curvelet series obeys ∥f - f C n ∥ 2 L2 ≤ C . n -2 . (log n) 3 , n → ∞. This rate of convergence holds uniformly over a class of functions that are C 2 except for discontinuities along piecewise C 2 curves and is essentially optimal. In comparison, the squared error of n-term wavelet approximations only converges as n -1 as n → ∞, which is considerably worse than the optimal behavior. Traditional methods of time-frequency and multiscale analysis, such as wavelets and Gabor frames, have been successfully employed for representing most classes of pseudodifferential operators. However, these methods are not equally effective in dealing with Fourier Integral Operators in general. In this article, we show that the shearlets, recently introduced by the authors and their collaborators, provide very efficient representations for a large class of Fourier Integral Operators. The shearlets are an affine-like system of well-localized waveforms at various scales, locations and orientations, which are particularly efficient in representing anisotropic functions. Using this approach, we prove that the matrix representation of a Fourier Integral Operator with respect to a Parseval frame of shearlets is sparse and well-organized. This fact recovers a similar result recently obtained by Candes and Demanet using curvelets, which illustrates the benefits of directional multiscale representations (such as curvelets and shearlets) in the study of those functions and operators where traditional multiscale methods are unable to provide the appropriate geometric analysis in the phase space. This paper argues that curvelets provide a powerful tool for representing very general linear symmetric systems of hyperbolic differential equations. Curvelets are a recently developed multiscale system [7, 9] in which the elements are highly anisotropic at fine scales, with effective support shaped according to the parabolic scaling principle width length2 at fine scales. We prove that for a wide class of linear hyperbolic differential equations, the curvelet representation of the solution operator is both optimally sparse and well organized. * It is sparse in the sense that the matrix entries decay nearly exponentially fast (i.e., faster than any negative polynomial) * and well organized in the sense that the very few nonnegligible entries occur near a few shifted diagonals. Indeed, we show that the wave group maps each curvelet onto a sum of curveletlike waveforms whose locations and orientations are obtained by following the different Hamiltonian flows - hence the diagonal shifts in the curvelet representation. A physical interpretation of this result is that curvelets may be viewed as coherent waveforms with enough frequency localization so that they behave like waves but at the same time, with enough spatial localization so that they simultaneously behave like particles.