code: def sheardec2D(X, shearletsystem):
      
      coeffs = np.zeros(shearletsystem.shearlets.shape, dtype=complex)
      Xfreq = fftshift(fft2(ifftshift(X)))
      for i in range(shearletsystem.nShearlets):
          coeffs[:, :, i] = fftshift(ifft2(ifftshift(Xfreq * np.conj(
                                     shearletsystem.shearlets[:, :, i]))))
      return coeffs.real
docs: Shearlet Decomposition function.