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 paper, we show that the shearlets, recently introduced by the authors and their collaborators...