wpplotcf.m
上传用户:haiyisale
上传日期:2013-01-09
资源大小:3246k
文件大小:8k
- function wpplotcf(Ts,Ds,colorMode,axe,nb_colors)
- %WPPLOTCF Plot wavelet packets colored coefficients.
- %
- % WPPLOTCF(TREE,DATA) OR
- % WPPLOTCF(TREE,DATA,COLORMODE) OR
- % WPPLOTCF(TREE,DATA,COLORMODE,AXE) OR
- % WPPLOTCF(TREE,DATA,COLORMODE,AXE,NB_COLORS)
- % M. Misiti, Y. Misiti, G. Oppenheim, J.M. Poggi 12-Mar-96.
- % Last Revision: 07-May-2003.
- % Copyright 1995-2004 The MathWorks, Inc.
- % $Revision: 1.15.4.2 $
- % Tag for image of Coefficients.
- %--------------------------------
- tag_img_cfs = 'Img_WPCfs';
- DefNBColor = 128;
- if nargin<3
- colorMode = 1; axe = gca; nb_colors = DefNBColor;
- elseif nargin<4
- axe = gca; nb_colors = DefNBColor;
- elseif nargin<5
- nb_colors = DefNBColor;
- end
- tab = wtreemgr('table',Ts,0)';
- nodes = tab(1,:)';
- sizes = wdatamgr('rsizes',Ds);
- nbtn = size(tab,2);
- order = wtreemgr('order',Ts);
- [depths,posis] = ind2depo(order,nodes);
- dmax = size(tab,1)-1;
- cfs = wdatamgr('rallcfs',Ds);
- switch order
- case 2 , flg_line = 6;
- case 4 , flg_line = 5;
- end
- if dmax==flg_line , lwidth = 0.5;
- elseif dmax==flg_line-1 , lwidth = 1;
- else , lwidth = 2;
- end
- switch order
- case 2
- switch colorMode
- case {1,2,3,4} , ord = wpfrqord(nodes);
- case {5,6,7,8} , ord = [1:length(nodes)];
- end
- switch colorMode
- case {1,2,5,6} , abs_val = 1;
- case {3,4,7,8} , abs_val = 0;
- end
- switch colorMode
- case {1,3,5,7} , cfs = wcodemat(cfs,nb_colors,'row',abs_val);
- end
- deb = [1];
- fin = [];
- for k = 1:nbtn
- fin(k) = deb(k)+sizes(1+depths(k))-1;
- deb(k+1) = fin(k)+1;
- end
- nbrows = (2.^(dmax-depths));
- NBrowtot = sum(nbrows);
- NBcoltot = sizes(1);
- matcfs = zeros(NBrowtot,NBcoltot);
- ypos = zeros(nbtn,1);
- if nbtn>1
- for k = 1:nbtn
- ypos(ord(k)) = sum(nbrows(ord([1:k-1])));
- end
- end
- ypos = NBrowtot+1-ypos-nbrows;
- ymin = (ypos-1)/NBrowtot;
- ymax = (ypos-1+nbrows)/NBrowtot;
- ylim = [0 1];
- alfa = 1/(2*NBrowtot);
- ydata = [(1-alfa)*ylim(1)+alfa*ylim(2) (1-alfa)*ylim(2)+alfa*ylim(1)];
- if NBrowtot==1
- ydata(1) = 1/2; ydata(2) = 1;
- end
- axes(axe);
- tag_axe = get(axe,'Tag');
- xlim = get(axe,'Xlim');
- set(axe,'Xlim',xlim,'Ylim',ylim,'Nextplot','replace');
- imgcfs = image(...
- 'Parent',axe, ...
- 'Xdata',[1:NBcoltot], ...
- 'Ydata',ydata, ...
- 'Cdata',matcfs, ...
- 'Tag',tag_img_cfs, ...
- 'Userdata',[depths posis ymin ymax] ...
- );
- NBdraw = 0;
- for k = 1:nbtn
- d = depths(k);
- z = cfs(deb(k):fin(k));
- z = z(ones(1,2^d),:);
- z = wkeep1(z(:)',NBcoltot);
- if find(colorMode==[2 4 6 8])
- z = wcodemat(z,nb_colors,'row',abs_val);
- end
- r1 = ypos(k);
- r2 = ypos(k)+nbrows(k)-1;
- matcfs(r1:r2,:) = z(ones(1,nbrows(k)),:);
- if dmax<=flg_line & nbtn~=1
- line(...
- 'Parent',axe, ...
- 'Xdata',[0.5 NBcoltot+0.5], ...
- 'Ydata',[ymin(k) ymin(k)], ...
- 'Linewidth',lwidth ...
- );
- end
- NBdraw = NBdraw+1;
- if NBdraw==10 | k==nbtn
- set(imgcfs,'Xdata',[1:NBcoltot], ...
- 'Ydata',ydata,'Cdata',matcfs);
- NBdraw = 0;
- end
- end
- set(axe,'Ydir','reverse',...
- 'layer','top', ...
- 'Xlim',xlim,'Ylim',ylim, ...
- 'Ytick',[],'YtickLabel',[], ...
- 'Box','on','Tag',tag_axe);
- case 4
- % Image Coding Value.
- %-------------------
- codemat_v = wimgcode('get',get(axe,'Parent'));
- Xbase = 1; Ybase = 1;
- hight = (Ybase./(2.^depths))';
- axes(axe);
- tag_axe = get(axe,'Tag');
- set(axe,'Xlim',[0 Xbase],'Ylim',[0 Ybase],'Nextplot','add');
- % Trunc images.
- %--------------
- newsiz = sizes;
- s0 = newsiz(:,1);
- for k=2:size(newsiz,2)
- sx = s0;
- for p = 1:k , sx = floor((sx+1)/2); end
- newsiz(:,k) = sx;
- end
- newsiz = newsiz(:,depths+1);
- size0 = sizes(:,1)';
- sizes = sizes(:,depths+1);
- deb = [1];
- fin = [];
- for k = 1:nbtn
- fin(k) = deb(k)+prod(sizes(:,k))-1;
- deb(k+1) = fin(k)+1;
- end
- large = (Xbase./(2.^depths))';
- xpos = zeros(1,nbtn);
- ypos = zeros(1,nbtn);
- pos = rem(tab,order);
- vdiv = 1./(2.^[1:dmax]);
- for k = 1:nbtn
- d = depths(k);
- dx = 0; dy = 0;
- for j = 1:d
- l = d-j+1;
- if pos(l,k)==1
- dy = dy + vdiv(j);
- elseif pos(l,k)==2
- dx = dx + vdiv(j);
- dy = dy + vdiv(j);
- elseif pos(l,k)==3
- elseif pos(l,k)==0
- dx = dx + vdiv(j);
- end
- end
- xpos(k) = Xbase*dx;
- ypos(k) = Ybase*dy;
- end
-
- % newsiz = sizes;
- xmin = xpos+(0.5*large)./(newsiz(2,:));
- if sizes(2,:)>1
- xmax = xpos+((newsiz(2,:)-0.5).*large)./newsiz(2,:);
- else
- xmax = xpos+large;
- end
- ymin = ypos+(0.5*hight)./newsiz(1,:);
- if sizes(1,:)>1
- ymax = ypos+((newsiz(1,:)-0.5).*hight)./newsiz(1,:);
- else
- ymax = ypos+hight;
- end
- NBdraw = 1;
- if nbtn==1 , flg_code = 0; else flg_code = 1; end
- for k = 1:nbtn
- tmp = zeros(sizes(:,k)');
- tmp(:) = cfs(deb(k):fin(k));
- img = image([xmin(k) xmax(k)],[ymin(k) ymax(k)], ...
- flipud(wimgcode('cod',flg_code,tmp,nb_colors,codemat_v,[depths(k) size0])),...
- 'Tag',tag_img_cfs,...
- 'Userdata',[depths(k) posis(k) xpos(k) xmax(k) ypos(k) ymax(k)],...
- 'Parent',axe);
- xx1 = xpos(k);
- xx2 = xx1+hight(k);
- yy1 = ypos(k);
- yy2 = yy1+hight(k);
- if dmax<=flg_line
- lin = line(...
- 'Parent',axe, ...
- 'Xdata',[xx1 xx2 xx2 xx1 xx1], ...
- 'Ydata',[yy1 yy1 yy2 yy2 yy1], ...
- 'Linewidth',lwidth ...
- );
- end
- if NBdraw==4 , NBdraw = 1; else , NBdraw = NBdraw+1; end
- end
- lin = line(...
- 'Parent',axe, ...
- 'Xdata',[0 1 1 0 0],...
- 'Ydata',[0 0 1 1 0],...
- 'Linewidth',2 ...
- );
- set(axe,'Tag',tag_axe);
- drawnow
- end