software  /  papers  /  gallery  /  links  /  mailing list

Gallery
- click thumbnails for detail -







In this thought experiment, we consider the wave equation modeling acoustic, electromagnetic or elastic waves. If the initial condition is chosen localized and directional like a curvelet, then it will remain so for larger times as well---only it will travel along the light (or sound) rays. This is a result we proved in this paper, and remains true even if the index of refraction of the medium is smoothly varying. As a result the curvelet matrix of the Green's function exhibits an optimally sparse structure. This result does not hold for other transforms such as wavelets, Gabor, ridgelets, or contourlets.


This coherence of curvelets under the wave flow has far-reaching consequences. In image (2), we depict a wave front---where the underlying function is nonsmooth---initially in the shape of an ellipse. If we let it evolve under the wave equation, the ellipse will shrink to the fishtail pattern (the outgoing ellipse, if it exists, is not shown). A caustic has formed, and the wavefront seems to have lost some of its smoothness near the tips of the tail, named cusps. This is not a concern in phase-space, or in the curvelet domain, where the cusp pattern is in fact explained as a superposition of curvelets. Curvelets do not see caustics; keep in mind that the wave equation is linear and that the superposition principle is valid.

We have developed a numerical wave solver based on the same sparsity ideas.



Last modified 9 June 2006 - Maintained by Laurent Demanet - -