波变换

开发平台：
Matlab

dw2dutil.m：源码内容
							function dw2dutil(option,win_dw2dtool,in3,in4)
%DW2DUTIL Discrete wavelet 2-D utilities.
%   DW2DUTIL(OPTION,WIN_DW2DTOOL,IN3,IN4)
%   M. Misiti, Y. Misiti, G. Oppenheim, J.M. Poggi 12-Mar-96.
%   Last Revision: 08-Jan-2001.
%   Copyright 1995-2002 The MathWorks, Inc.
%   $Revision: 1.14 $ $Date: 2002/06/17 12:18:47 $
% Get Globals.
%-------------
[Terminal_Prop,Def_AxeFontSize, Def_FraBkColor] = ...
    mextglob('get','Terminal_Prop','Def_AxeFontSize','Def_FraBkColor');
% Default values.
%----------------
max_lev_anal = 5;
default_nbcolors = 128;
% Tag property of objects.
%-------------------------
tag_m_savesyn  = 'Save_Syn';
tag_m_savecfs  = 'Save_Cfs';
tag_m_savedec  = 'Save_Dec';
tag_cmd_frame  = 'Cmd_Frame';
tag_pus_anal   = 'Pus_Anal';
tag_pus_deno   = 'Pus_Deno';
tag_pus_comp   = 'Pus_Comp';
tag_pus_hist   = 'Pus_Hist';
tag_pus_stat   = 'Pus_Stat';
tag_pop_declev = 'Pop_DecLev';
tag_pus_visu   = 'Pus_Visu';
tag_pus_big    = 'Pus_Big';
tag_pus_rec    = 'Pus_Rec';
tag_pop_viewm  = 'Pop_ViewM';
tag_pus_full   = ['Pus_Full.1';'Pus_Full.2';'Pus_Full.3';'Pus_Full.4'];
tag_btnaxeset  = 'Btn_Axe_Set';
tag_axefigutil = 'Axe_FigUtil';
tag_linetree   = 'Tree_lines';
tag_txttree    = 'Tree_txt';
tag_axeimgbig  = 'Axe_ImgBig';
tag_axeimgini  = 'Axe_ImgIni';
tag_axeimgvis  = 'Axe_ImgVis';
tag_axeimgsel  = 'Axe_ImgSel';
tag_axeimgdec  = 'Axe_ImgDec';
tag_axeimgsyn  = 'Axe_ImgSyn';
tag_axeimghdls = 'Img_Handles';
tag_imgdec     = 'Img_Dec';
% MemBloc1 of stored values.
%---------------------------
n_param_anal   = 'DWAn2d_Par_Anal';
ind_img_name   = 1;
ind_wav_name   = 2;
ind_lev_anal   = 3;
ind_img_t_name = 4;
ind_img_size   = 5;
ind_nbcolors   = 6;
ind_act_option = 7;
ind_simg_type  = 8;
ind_thr_val    = 9;
nb1_stored     = 9;
% MemBloc3 of stored values.
%---------------------------
n_miscella      = 'DWAn2d_Miscella';
ind_graph_area  =  1;
ind_pos_axebig  =  2;
ind_pos_axeini  =  3;
ind_pos_axevis  =  4;
ind_pos_axedec  =  5;
ind_pos_axesyn  =  6;
ind_pos_axesel  =  7;
ind_view_status =  8;
ind_save_status =  9;
ind_sel_funct   = 10;
nb3_stored      = 10;
% Miscellaneous values.
%----------------------
square_viewm = 1;
tree_viewm   = 2;
% View Status
%--------------------------------------------------------%
% 'none' : init
% 's_l*' : square        * = lev_dec (1 --> Level_Anal)
% 'f1l*' : full ini      * = lev_dec (1 --> Level_Anal)
% 'f2l*' : full syn      * = lev_dec (1 --> Level_Anal)
% 'f3l*' : full vis      * = lev_dec (1 --> Level_Anal)
% 'f4l*' : full dec      * = lev_dec (1 --> Level_Anal)
% 'b*l*' : big
%                first   * = index   (1 --> 4*Level_Anal)
%                second  * = lev_dec (1 --> Level_Anal)
% 't_l*' : tree          * = lev_dec (1 --> Level_Anal)
%--------------------------------------------------------%
% Handles of tagged objects.
%---------------------------
str_numwin  = sprintf('%.0f',win_dw2dtool);
children    = get(win_dw2dtool,'Children');
uic_handles = findobj(children,'flat','type','uicontrol');
axe_handles = findobj(children,'flat','type','axes');
txt_handles = findobj(uic_handles,'Style','text');
pop_handles = findobj(uic_handles,'Style','popupmenu');
pus_handles = findobj(uic_handles,'Style','pushbutton');
m_files     = wfigmngr('getmenus',win_dw2dtool,'file');
m_savesyn   = findobj(m_files,'Tag',tag_m_savesyn);
m_savecfs   = findobj(m_files,'Tag',tag_m_savecfs);
m_savedec   = findobj(m_files,'Tag',tag_m_savedec);
pus_anal    = findobj(pus_handles,'Tag',tag_pus_anal);
pus_deno    = findobj(pus_handles,'Tag',tag_pus_deno);
pus_comp    = findobj(pus_handles,'Tag',tag_pus_comp);
pus_hist    = findobj(pus_handles,'Tag',tag_pus_hist);
pus_stat    = findobj(pus_handles,'Tag',tag_pus_stat);
pop_declev  = findobj(pop_handles,'Tag',tag_pop_declev);
pus_visu    = findobj(pus_handles,'Tag',tag_pus_visu);
pus_big     = findobj(pus_handles,'Tag',tag_pus_big);
pus_rec     = findobj(pus_handles,'Tag',tag_pus_rec);
pop_viewm   = findobj(pop_handles,'Tag',tag_pop_viewm);
for k =1:size(tag_pus_full,1)
    pus_full(k) = (findobj(pus_handles,'Tag',tag_pus_full(k,:)))';
end
Axe_ImgBig = findobj(axe_handles,'flat','Tag',tag_axeimgbig);
Axe_ImgIni = findobj(axe_handles,'flat','Tag',tag_axeimgini);
Axe_ImgVis = findobj(axe_handles,'flat','Tag',tag_axeimgvis);
Axe_ImgSel = findobj(axe_handles,'flat','Tag',tag_axeimgsel);
Axe_ImgSyn = findobj(axe_handles,'flat','Tag',tag_axeimgsyn);
switch option
    case 'clean'
        % in3 = load option
        % in4 optional ('new_synt')
        %---------------------------
        if nargin==3 , init = 1; else , init = 0; end
        if init==1
            % in3 = type of loading.
            %-----------------------
            % 'load_img' , 'load_dec'
            % 'load_cfs' , 'demo'    
            %-----------------------
            strwin = sprintf('%.0f',win_dw2dtool);
            if strcmp(in3,'load_cfs')
                str_btn = 'Synthesize';
                cba_btn = ['dw2dmngr(''synthesize'',' strwin ');'];
            else
                str_btn = 'Analyze';
                cba_btn = ['dw2dmngr(''analyze'',' strwin ');'];
            end
            set(pus_anal,'String',xlate(str_btn),'Callback',cba_btn);
        end
        % Testing first use.
        %--------------------
        active_option = wmemtool('rmb',win_dw2dtool,n_param_anal,ind_act_option);
        if strcmp(active_option,'create') , first = 1; else first = 0; end
        % Cleaning stored values.
        %------------------------
        if init==1 , wmemtool('ini',win_dw2dtool,n_param_anal,nb1_stored); end
        % End of Cleaning when first is true.
        %------------------------------------
        if first , return; end
        % Cleaning axes.
        %---------------
        dw2dimgs('clean',win_dw2dtool);
        wmemtool('wmb',win_dw2dtool,tag_axeimghdls,1,[]);
        Axe_ImgDec = wmemtool('rmb',win_dw2dtool,tag_axeimgdec,1);
        cleanaxe([Axe_ImgIni,Axe_ImgVis,Axe_ImgSyn,...
                  Axe_ImgSel,Axe_ImgBig,Axe_ImgDec(:)']);
        set([Axe_ImgIni,Axe_ImgVis,Axe_ImgSel, ...
             Axe_ImgSyn,Axe_ImgDec,Axe_ImgBig],...
             'Visible','Off');
        bdy = 18;
        wboxtitl('set',Axe_ImgSel,'Image Selection',Def_AxeFontSize,...
                                        9,18,bdy,'off');
        axe_figutil = findobj(win_dw2dtool,              ...
                                    'type','axes',       ...
                                    'Tag',tag_axefigutil ...
                                    );
        t_lines     = findobj(axe_figutil,               ...
                                    'type','line',       ...
                                    'Tag',tag_linetree   ...
                                    );
        t_txt       = findobj(axe_figutil,               ...
                                    'type','text',       ...
                                    'Tag',tag_txttree    ...
                                    );
        set([t_lines' t_txt'],'Visible','off');
        dw2darro('clean',win_dw2dtool,'off');
        drawnow;
        % Setting enable property of objects.
        %------------------------------------
        set([m_savesyn, m_savecfs, m_savedec],'Enable','Off');
        cbanapar('enable',win_dw2dtool,'off');
        cbcolmap('enable',win_dw2dtool,'off');
        set([...
             pus_anal,   pus_deno,           ...
             pus_comp,   pus_hist, pus_stat, ...
             pop_declev, pus_big,            ...
             pus_visu,   pus_rec,  pus_big,  ...
             pus_full,   pop_viewm ],...
             'Enable','off'...
             );
        % Cleaning DynVTool.
        %-------------------
        dynvtool('stop',win_dw2dtool);
        if init==1
            % Setting string property of objects.
            %------------------------------------
            str_lev_data    = int2str([1:max_lev_anal]');
            cbanapar('set',win_dw2dtool,...
                'nam','',             ...
                'wav','haar',         ...
                'lev',{'String',str_lev_data,'Value',1})
            set(pop_declev,'String',str_lev_data,'Value',1);
        end
        % Cleaning buttons.
        %------------------
        val = get(pop_viewm,'Value');
        pos = zeros(4,4);
        pos(1,:) = get(pus_full(1),'Position');
        if val~=square_viewm
            pos(1,3) = (3*pos(1,3))/2;
            pos(2,:) = pos(1,:); pos(2,2) = pos(2,2)-pos(2,4);
            pos(3,:) = pos(1,:); pos(3,1) = pos(3,1)+pos(3,3);
            pos(4,:) = pos(3,:);
            pos(4,2) = pos(4,2)-pos(4,4);
            set(pus_full,'Visible','off'); drawnow;
            for num=1:4
                set(pus_full(num),...
                        'Position',pos(num,:),...
                        'Backgroundcolor',Def_FraBkColor,...
                        'String',sprintf('%.0f',num),...
                        'Userdata',0);
            end
            set(pus_full,'Visible','on'); drawnow;
        else
            for num=1:length(pus_full)
                set(pus_full(num),...
                        'Backgroundcolor',Def_FraBkColor,...
                        'String',sprintf('%.0f',num),...               
                        'Userdata',0);
            end
        end
        set(pop_viewm,'Value',square_viewm);
        cba_pus_big = ['dw2dmngr(''select'','   ...
                        str_numwin ','          ...
                        num2mstr(pus_big) ');'];
        set(pus_big,...
                'String',xlate('Full Size'),           ...
                'BackGroundColor',Def_FraBkColor,...
                'Callback',cba_pus_big);
        % Cleaning miscellaneous parameters.
        %-----------------------------------
        wmemtool('wmb',win_dw2dtool,n_miscella, ...
                       ind_view_status,'none',ind_save_status,'none');
        if init==1
            cbcolmap('set',win_dw2dtool,'pal',{'pink',default_nbcolors})
        end
    case 'clean2'
        % in3 = calling option
        %------------------
        call_opt = in3;
        % Cleaning axes.
        %---------------
        dw2dimgs('clean',win_dw2dtool);
        wmemtool('wmb',win_dw2dtool,tag_axeimghdls,1,[]);
        Axe_ImgDec = wmemtool('rmb',win_dw2dtool,tag_axeimgdec,1);
        cleanaxe([Axe_ImgVis,Axe_ImgSyn,Axe_ImgSel,Axe_ImgBig,Axe_ImgDec(:)']);
        set([Axe_ImgVis,Axe_ImgSel,Axe_ImgSyn,Axe_ImgDec,Axe_ImgBig],...
             'Visible','Off');
        bdy = 18;
        wboxtitl('set',Axe_ImgSel,'Image Selection',Def_AxeFontSize,...
                                        9,18,22,'off');
        axe_figutil = findobj(win_dw2dtool,              ...
                                    'type','axes',       ...
                                    'Tag',tag_axefigutil ...
                                    );
        t_lines     = findobj(axe_figutil,               ...
                                    'type','line',       ...
                                    'Tag',tag_linetree   ...
                                    );
        t_txt       = findobj(axe_figutil,               ...
                                    'type','text',       ...
                                    'Tag',tag_txttree    ...
                                    );
        set([t_lines' t_txt'],'Visible','off');
        dw2darro('clean',win_dw2dtool,'on');
        dw2dutil('pos_axe_init',win_dw2dtool,option);
        drawnow;
        % Cleaning buttons.
        %------------------
        val = get(pop_viewm,'Value');
        pos = zeros(4,4);
        pos(1,:) = get(pus_full(1),'Position');
        if val~=square_viewm
            pos(1,3) = (3*pos(1,3))/2;
            pos(2,:) = pos(1,:); pos(2,2) = pos(2,2)-pos(2,4);
            pos(3,:) = pos(1,:); pos(3,1) = pos(3,1)+pos(3,3);
            pos(4,:) = pos(3,:);
            pos(4,2) = pos(4,2)-pos(4,4);
            set(pus_full,'Visible','off'); drawnow;
            for num=1:4
                set(pus_full(num),...
                        'Position',pos(num,:),...
                        'Backgroundcolor',Def_FraBkColor,...
                        'String',sprintf('%.0f',num),...
                        'Userdata',0);
            end
            set(pus_full,'Visible','on'); drawnow;
        else
            for num=1:length(pus_full)
                set(pus_full(num),...
                        'Backgroundcolor',Def_FraBkColor,...
                        'String',sprintf('%.0f',num),...               
                        'Userdata',0);
            end
        end
        set(pop_viewm,'Value',square_viewm);
        % Cleaning miscellaneous parameters.
        %-----------------------------------
        wmemtool('wmb',win_dw2dtool,n_miscella, ...
                       ind_view_status,'none',ind_save_status,'none');
    case 'enable'
        % in3 = calling option.
        %----------------------
        switch in3
            case {'load_img','demo','load_cfs'}
                cbanapar('enable',win_dw2dtool,'on');
                cbcolmap('enable',win_dw2dtool,'on');
                set(pus_anal,'Enable','on');
                dw2darro('set_arrow',win_dw2dtool,in3);
            end
        switch in3
            case {'demo','analyze','synthesize'}
                set([pus_deno, pus_comp,   pus_stat,   ...
                     pus_hist, pop_declev,             ...
                     pus_visu, pus_rec,                ...
                     pus_big , pus_full,   pop_viewm], ...
                     'Enable','On'   ...
                     );
                set([m_savesyn, m_savecfs, m_savedec],'Enable','on');
            case 'load_dec'
                cbanapar('enable',win_dw2dtool,'on');
                cbcolmap('enable',win_dw2dtool,'on');
                set([pus_anal, pus_deno,   pus_comp, pus_stat, ...
                     pus_hist, pop_declev, pus_visu, pus_rec , ...
                     pus_big,  pop_viewm,  pus_full] ,...
                     'Enable','On'   ...
                     );
                dw2darro('set_arrow',win_dw2dtool,in3);
                set([m_savesyn, m_savecfs, m_savedec],'Enable','on');
            case {'comp','deno'}
                set([m_files , pus_anal , pus_deno , pus_comp],'Enable','off');
            case {'return_comp','return_deno'}
                set([m_files , pus_anal , pus_deno , pus_comp],'Enable','on');
        end
    case 'pos_axe_init'
        % in3 = calling option.
        %--------------------------
        % Getting Analysis parameters.
        %-----------------------------
        [Image_Size,Level_Anal] = wmemtool('rmb',win_dw2dtool,n_param_anal,...
                                                ind_img_size,ind_lev_anal);
        level       = cbanapar('get',win_dw2dtool,'lev');
        Axe_ImgDec  = wmemtool('rmb',win_dw2dtool,tag_axeimgdec,1);
        Img_Handles = wmemtool('rmb',win_dw2dtool,tag_axeimghdls,1);
        % Getting position parameters.
        %-----------------------------
        pos      = wmemtool('rmb',win_dw2dtool,n_miscella,ind_graph_area);
        term_dim = Terminal_Prop;
        % View boundary parameters.
        %--------------------------
        mxbig = 0.83; mybig = 0.83;
        mx    = 0.75; my    = 0.75;
        % Computing axes positions.
        %--------------------------
        cent_one = [ pos(1)+pos(3)/2 , pos(2)+pos(4)/2 ];
        w_theo  = pos(3);               h_theo  = pos(4);
        w_pos   = pos(3)*mxbig;         h_pos   = pos(4)*mybig;
        [w_one,h_one] = wpropimg(Image_Size,w_pos,h_pos);
        pos_axebig = [cent_one(1)-w_one/2 ,...
                      cent_one(2)-h_one/2 ,...
                      w_one               ,...
                      h_one];
        NBL = 2; NBC = 2;
        w_theo = pos(3)/NBC;  h_theo = pos(4)/NBL;
        w_pos   = w_theo*mx;  h_pos   = h_theo*my;
        X_cent = pos(1)+(w_theo/2)*[1:2:2*NBC-1];
        Y_cent = pos(2)+(h_theo/2)*[1:2:2*NBL-1];
        [w_used,h_used] = wpropimg(Image_Size,w_pos,h_pos);
        w_u2 = w_used/2;     h_u2 = h_used/2;
        pos_axeini = [X_cent(1)-w_u2,Y_cent(2)-h_u2,w_used,h_used];
        pos_axevis = [X_cent(2)-w_u2,Y_cent(2)-h_u2,w_used,h_used];
        pos_axesyn = [X_cent(1)-w_u2,Y_cent(1)-h_u2,w_used,h_used];
        pos_axesel = [X_cent(2)-w_u2,Y_cent(1)-h_u2,w_used,h_used];
        xl = pos_axesel(1);             yb = pos_axesel(2);
        la = pos_axesel(3)/2;           ha = pos_axesel(4)/2;
        ind = 1;
        for k = 1:max_lev_anal
            pos_axedec(ind:ind+3,1:4) = ...
                    [xl      ,yb     ,la     ,ha;
                     xl+la   ,yb     ,la     ,ha;
                     xl+la   ,yb+ha  ,la     ,ha;
                     xl      ,yb+ha  ,la     ,ha ...
                    ];
            ind = ind+4;
            yb = yb+ha;    la = la/2;     ha = ha/2;
        end
        % Axes visibility.
        %-----------------
        if isempty(Level_Anal) , Level_Anal = level; end 
        max_a = 4*Level_Anal;
        index = [1:max_a];
        rem_4 = rem(index,4);
        if strcmp(in3,'load_img')
            ind_On = [];    ind_Off = index;
        else
            ind_On  = [find(index<max_a & rem_4~=0) , max_a];
            ind_Off = [find(index<max_a & rem_4==0) , find(index>max_a)];
        end
        % Setting axes.
        %--------------
        xlim = [1 Image_Size(1)];
        ylim = [1 Image_Size(2)];
        set([Axe_ImgIni,Axe_ImgVis,Axe_ImgSyn,Axe_ImgSel],'Visible','On');
        set(Axe_ImgIni,...
                'Visible','On','Position',pos_axeini, ...
                'Xlim',xlim,'Ylim',ylim               ...
                );
        set(Axe_ImgVis,...
                'Visible','On','Position',pos_axevis, ...
                'Xlim',xlim,'Ylim',ylim               ...
                );
        set(Axe_ImgSyn,...
                'Visible','On','Position',pos_axesyn, ...
                'Xlim',xlim,'Ylim',ylim               ...
                );
        set(Axe_ImgSel,...
                'Visible','on','Position',pos_axesel, ...
                'Xlim',xlim,'Ylim',ylim...
                );
        for k = 1:4*max_lev_anal
            set(Axe_ImgDec(k),...
                    'Position',pos_axedec(k,:),...
                    'Xlim',xlim,'Ylim',ylim);
        end
        set(Axe_ImgDec(ind_On),'Visible','On');
        set(Axe_ImgBig,'Position',pos_axebig,'Xlim',xlim,'Ylim',ylim);
        % Axes Titles.
        %-------------
        axes(Axe_ImgIni);
        if strcmp(in3,'load_cfs')
            wtitle('Initial Image ... None','Parent',Axe_ImgIni);
        else
            wtitle('Original Image','Parent',Axe_ImgIni);
        end
        axes(Axe_ImgSyn);
        wtitle('Synthesized Image','Parent',Axe_ImgSyn);
        wsetxlab(Axe_ImgSel,'Decomposition at level ...');
        wboxtitl('pos',Axe_ImgSel,'on');
        % Storing axes positions.
        %------------------------
        wmemtool('wmb',win_dw2dtool,n_miscella,   ...
                       ind_pos_axeini,pos_axeini, ...
                       ind_pos_axevis,pos_axevis, ...
                       ind_pos_axesel,pos_axesel, ...
                       ind_pos_axesyn,pos_axesyn, ...
                       ind_pos_axebig,pos_axebig, ...
                       ind_pos_axedec,pos_axedec  ...
                       );
    otherwise
        errargt(mfilename,'Unknown Option','msg');
        error('*');
end

						