# Shearlets

In the following we will give a short introduction into the theory of shearlets. The shearlet transform is unlike the traditional wavelet transform which does not posses the ability to detect directionality, since it is merely associated with two parameters, the scaling parameter *a* and the the translation parameter *t*. The idea now is to define a transform, which overcomes this vice, while retaining most aspects of the mathematical framework of wavelets, e.g., the fact that

- the associated system forms an affine system,
- the transform can be regarded as matrix coefficients of a unitary representation of a special group,
- there is an MRA-structure associated with the systems.

The shearlets satisfy all these properties in addition to showing optimal behavior with respect to the detection of directional information.

The basic idea for the definition of *continuous shearlets* is the usage of a 2-parameter dilation group, which consists of products of parabolic scaling matrices and shear matrices. Hence the continuous shearlets depend on three parameters, the scaling parameter a > 0, the shear parameter s ∈ *R* and the translation parameter t ∈ *R*^{2}, and they are defined by

_{a,s,t}(x)=a

^{-3/4}ψ ((D

_{a,s}

^{-1 }(x-t)), where D

_{a,s}= [a,-a

^{1/2}s;0,a

^{1/2}].

The mother shearlet function ψ is defined almost like a tensor product by

_{1},ξ

_{2}) =ψ

_{1}(ξ

_{1})ψ

_{2}(ξ

_{2}/ξ

_{1}),

where ψ_{1} is a wavelet and ψ_{2} is a bump function. The figure on the right hand side illustrates the behavior of the continuous shearlets in frequency domain assuming that ψ_{1} and ψ_{2} are chosen to be compactly supported in frequency domain.

The associated *continuous shearlet transform* again depends on the scaling parameter a, the shear parameter s and the translation parameter t, and is defined by

*SH*(a,s,t) = <f,ψ

_{f}_{a,s,t}>

This transform can also be regarded as matrix coefficients of the unitary representation

_{a,s,t}(x)=a

^{-3/4}ψ ((D

_{a,s}

^{-1 }(x-t))

of the *shearlet group* *S*=*R*^{+} x *R* x *R*^{2} with multiplication given by

^{1/2}, t + D

_{a,s}t').

The continuous shearlet transform precisely resolves the wavefront set, since the location parameter t does detect the location of singularities, whereas the shear parameter shows the direction perpendicular to the direction of the singularity.

To illustrate the ability of shearlets to detect directions, the following images show ψ_{a,s,t} in the time and frequency domain for a=0.3, s=0 and t=0.

The next images illustrate the behavior of shearlets with respect to the shearing operator. The image on the left-hand-side shows ψ_{a,s,t} for a=0.3, s=-0.5 and t=0, whereas the image on the right-hand-side shows its Fourier transform.

The directionality of the shearlets in time domain is evident.

By sampling the continuous shearlet transform on an appropriate discrete set of the scaling, shear, and translation parameters, it is possible to obtain a frame or even a Parseval frame for L^{2}(*R*). We will briefly discuss two possibilities to discretize the continuous shearlet transform.

To obtain the *discrete shearlets*, we sample the three parameters as a_{j}=2^{j} (j ∈ Z) , s_{j,k}=k a_{j}^{1/2} =k2^{j/2} (k ∈ Z), and t_{j,k,m}=D_{aj,sj,k} (m ∈ Z^{2}). We choose the mother shearlet function ψ in a similar fashion as in the continuous case, i.e., we now choose ψ_{1} to be a discrete wavelet and ψ_{2} to be bump function with certain weak additional properties. The tiling of the frequency plane is illustrated in the figure on the left hand side.

This system forms a Parseval frame for L^{2}(*R*), and they are optimally sparse. Furthermore, they are associated with a generalized MRA-structure, where the scaling space is not only translation invariant but also invariant under the shear operator.