THESE/codes/wave/allcode/dwt2D_plots.m

   1 % dwt2D_plots
   2 % DISPLAY 2D WAVELETS OF DWT2D.M
   3
   4 [af, sf] = farras;
   5 J = 5;
   6 L = 3*2^(J+1);
   7 N = L/2^J;
   8 x = zeros(L,3*L);
   9 w = dwt2D(x,J,af);
  10 w{J}{1}(N/2,N/2+0*N) = 1;
  11 w{J}{2}(N/2,N/2+1*N) = 1;
  12 w{J}{3}(N/2,N/2+2*N) = 1;
  13 y = idwt2D(w,J,sf);
  14
  15 figure(1)
  16 clf
  17 imagesc(y);
  18 axis equal
  19 axis off
  20 colormap(gray(128))
  21