THESE/codes/wave/allcode/cplxdual2D_plots.m

   1 % cplxdual2D_plots
   2 % DISPLAY 2D WAVELETS OF cplxdual2D.M
   3
   4 J = 4;
   5 L = 3*2^(J+1);
   6 N = L/2^J;
   7 [Faf, Fsf] = FSfarras;
   8 [af, sf] = dualfilt1;
   9 x = zeros(2*L,6*L);
  10 w = cplxdual2D(x, J, Faf, af);
  11 w{J}{1}{2}{2}(N/2,N/2+0*N) = 1;
  12 w{J}{1}{1}{3}(N/2,N/2+1*N) = 1;
  13 w{J}{1}{2}{1}(N/2,N/2+2*N) = 1;
  14 w{J}{1}{1}{1}(N/2,N/2+3*N) = 1;
  15 w{J}{1}{2}{3}(N/2,N/2+4*N) = 1;
  16 w{J}{1}{1}{2}(N/2,N/2+5*N) = 1;
  17 w{J}{2}{2}{2}(N/2+N,N/2+0*N) = 1;
  18 w{J}{2}{1}{3}(N/2+N,N/2+1*N) = 1;
  19 w{J}{2}{2}{1}(N/2+N,N/2+2*N) = 1;
  20 w{J}{2}{1}{1}(N/2+N,N/2+3*N) = 1;
  21 w{J}{2}{2}{3}(N/2+N,N/2+4*N) = 1;
  22 w{J}{2}{1}{2}(N/2+N,N/2+5*N) = 1;
  23 y = icplxdual2D(w, J, Fsf, sf);
  24 y = [y; sqrt(y(1:L,:).^2+y(L+[1:L],:).^2)];
  25 figure(1)
  26 clf
  27 imagesc(y);
  28 title('2D Dual-Tree Complex Wavelets')
  29 axis image
  30 axis off
  31 colormap(gray(128))
  32 print -djpeg95 cplxdual2D_plots