SHEARLETTRANSFORM

Digital Shearlet Transform (W sqrt(w) P)X W is the windowing, w is the weight operator, and P is the ppft operator

Contents

Desscription

SHX = SHEARLETTRANSFORM_XZ(X,R,BETA,WEIGHT,LEVEL) Input:

   X     - square N*N matrix  N dyadic
   R     - oversamping rate, 2,4,8,16.
   beta  - 2 or 4;
   w     - Wegithing Matrix, size: 2*(R*N+1)*(N+1) must be nonnegative
   level - decomposition level, level =0 means decompose to the lowest possible level

Output:

   ShX shearlet_coeff : shearlet coefficients (cell structure)

Examples

  tic
  img = imread('barbara.gif');
  img = double(img);
  N   = size(img,1); beta = 4;
  R   = 4; basisChoice  = 5;
  w   = generateW(N,R,basisChoice);
  shX = ShearletTransform(img, R, beta, w, 0);
  toc
Forward Shearlet Transform...
ppFT ...
Weighting ...
Windowing ...
Decomposition of high-frequency part...
scale = 5, Number of Tiles = 65
scale = 4, Number of Tiles = 33
scale = 3, Number of Tiles = 17
scale = 2, Number of Tiles = 9
scale = 1, Number of Tiles = 5
scale = 0, Number of Tiles = 3
scale = -1, Number of Tiles = 1
Decomposition of low frequency part...
Done!
Done with Forward Shearlet Transform!
Elapsed time is 3.434457 seconds.

See also

ADJSHEARLETTRANSFORM, INVSHEARLETTRANSFORM, PPFT, ADJPPFT, INVPPFTCG, WINDOWONPPGRID ADJWINDOWONPPGRID GENERATEW

Copyright

Copyright (C) 2011. Xiaosheng Zhuang, University of Osnabrueck