THESE/codes/wave/allcode/dwt3D_plots.m

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