THESE/codes/wave/allcode/dualtree3D_plots.m

   1 % dualtree3D_plots
   2 % DISPLAY 3D WAVELETS OF DUALTREE3D.M
   3
   4
   5 [Faf, Fsf] = FSfarras;
   6 [af, sf] = dualfilt1;
   7 J = 4;
   8 L = 3*2^(J+1);
   9 N = L/2^J;
  10 x = zeros(L,L,L);
  11 w = dualtree3D(x, J, Faf, af);
  12 w{J}{4}{7}(N/2,N/2,N/2) = 1;
  13 y = idualtree3D(w, J, Fsf, sf);
  14 figure(1)
  15 clf
  16 v = 1:L;
  17 S = 0.0025;
  18 p1 = patch(isosurface(v,v,v,y,S));
  19 isonormals(v,v,v,y,p1);
  20 set(p1,'FaceColor','red','EdgeColor','none');
  21 hold on
  22 p2 = patch(isosurface(v,v,v,y,-S));
  23 isonormals(v,v,v,y,p2);
  24 set(p2,'FaceColor','blue','EdgeColor','none');
  25 hold off
  26 daspect([1 1 1]);
  27 view(-30,30);
  28 camlight;
  29 lighting phong
  30 grid
  31 axis([32 48 32 48 32 48])
  32 set(gca,'fontsize',7)
  33 title('3-D WAVELET ISOSURFACE (COMPLEX DT TRANSFORM)')
  34 set(gcf,'paperposition',[0.5 0.5 0 0]+[0 0 4 3])
  35 print -djpeg95 dualtree3D_plots