THESE/codes/wave/allcode/cplxdual3D_plots.m

   1 % cplxdual3D_plots
   2 % DISPLAY 3D WAVELETS OF CPLXDUAL3D.M
   3
   4 [Faf, Fsf] = FSfarras;
   5 [af, sf] = dualfilt1;
   6 J = 4;
   7 L = 2*2^(J+1);
   8 N = L/2^J;
   9 x = zeros(L,L,L);
  10 w = cplxdual3D(x, J, Faf, af);
  11 w{J}{1}{1}{1}{3}(N/2,N/2,N/2) = 1;
  12 y = icplxdual3D(w, J, Fsf, sf);
  13 figure(1)
  14 clf
  15 v = 1:L;
  16 S = 0.0010;
  17 p1 = patch(isosurface(v,v,v,y,S));
  18 isonormals(v,v,v,y,p1);
  19 set(p1,'FaceColor','red','EdgeColor','none');
  20 hold on
  21 p2 = patch(isosurface(v,v,v,y,-S));
  22 isonormals(v,v,v,y,p2);
  23 set(p2,'FaceColor','blue','EdgeColor','none');
  24 hold off
  25 daspect([1 1 1]);
  26 view(-30,30);
  27 camlight;
  28 lighting phong
  29 grid
  30 axis([12 38 12 38 12 38])
  31 set(gca,'fontsize',7)
  32 title('3-D WAVELET ISOSURFACE (REAL PART)')
  33 set(gcf,'paperposition',[0.5 0.5 0 0]+[0 0 3 3])
  34 print -djpeg95 cplxdual3D_plots_1
  35 print -depsc cplxdual3D_plots_1
  36 yr = y;
  37
  38 w{J}{1}{1}{1}{3}(N/2,N/2,N/2) = 0;
  39 w{J}{2}{2}{2}{3}(N/2,N/2,N/2) = 1;
  40 y = icplxdual3D(w, J, Fsf, sf);
  41 figure(2)
  42 clf
  43 v = 1:L;
  44 S = 0.0010;
  45 p1 = patch(isosurface(v,v,v,y,S));
  46 isonormals(v,v,v,y,p1);
  47 set(p1,'FaceColor','red','EdgeColor','none');
  48 hold on
  49 p2 = patch(isosurface(v,v,v,y,-S));
  50 isonormals(v,v,v,y,p2);
  51 set(p2,'FaceColor','blue','EdgeColor','none');
  52 hold off
  53 daspect([1 1 1]);
  54 view(-30,30);
  55 camlight;
  56 lighting phong
  57 grid
  58 axis([12 38 12 38 12 38])
  59 set(gca,'fontsize',7)
  60 title('3-D WAVELET ISOSURFACE (IMAGINARY PART)')
  61 set(gcf,'paperposition',[0.5 0.5 0 0]+[0 0 3 3])
  62 print -djpeg95 cplxdual3D_plots_2
  63 print -depsc cplxdual3D_plots_2
  64 yi = y;
  65
  66 y = sqrt(yr.^2 + yi.^2);
  67 figure(3)
  68 clf
  69 v = 1:L;
  70 S = 0.0010;
  71 p1 = patch(isosurface(v,v,v,y,S));
  72 isonormals(v,v,v,y,p1);
  73 set(p1,'FaceColor','red','EdgeColor','none');
  74 hold on
  75 p2 = patch(isosurface(v,v,v,y,-S));
  76 isonormals(v,v,v,y,p2);
  77 set(p2,'FaceColor','blue','EdgeColor','none');
  78 hold off
  79 daspect([1 1 1]);
  80 view(-30,30);
  81 camlight;
  82 lighting phong
  83 grid
  84 axis([12 38 12 38 12 38])
  85 set(gca,'fontsize',7)
  86 title('3-D WAVELET ISOSURFACE (MAGNITUDE)')
  87 set(gcf,'paperposition',[0.5 0.5 0 0]+[0 0 3 3])
  88 print -djpeg95 cplxdual3D_plots_3
  89 print -depsc cplxdual3D_plots_3
  90