THESE/codes/wave/allcode/dualtree2D_plots.m

   1
   2 % dualtree2D_plots
   3 % DISPLAY 2D WAVELETS OF dualtree2D.M
   4
   5 J = 4;
   6 L = 3*2^(J+1);
   7 N = L/2^J;
   8 [Faf, Fsf] = FSfarras;
   9 [af, sf] = dualfilt1;
  10 x = zeros(2*L,3*L);
  11 w = dualtree2D(x, J, Faf, af);
  12 w{J}{1}{1}(N/2,N/2+0*N) = 1;
  13 w{J}{1}{2}(N/2,N/2+1*N) = 1;
  14 w{J}{1}{3}(N/2,N/2+2*N) = 1;
  15 w{J}{2}{1}(N/2+N,N/2+0*N) = 1;
  16 w{J}{2}{2}(N/2+N,N/2+1*N) = 1;
  17 w{J}{2}{3}(N/2+N,N/2+2*N) = 1;
  18 y = idualtree2D(w, J, Fsf, sf);
  19
  20 figure(1)
  21 clf
  22 imagesc(y);
  23 axis equal
  24 axis off
  25 colormap(gray(128))
  26