波变换

开发平台：
Matlab

compwav.m：源码内容
							function varargout = compwav(varargin)
%COMPWAV Pattern detection tool using an adapted wavelet versus a well known one.
%   VARARGOUT = COMPWAV(VARARGIN)
% COMPWAV M-file for compwav.fig
%      COMPWAV, by itself, creates a new COMPWAV or raises the existing
%      singleton*.
%
%      H = COMPWAV returns the handle to a new COMPWAV or the handle to
%      the existing singleton*.
%
%      COMPWAV('Property','Value',...) creates a new COMPWAV using the
%      given property value pairs. Unrecognized properties are passed via
%      varargin to compwav_OpeningFcn.  This calling syntax produces a
%      warning when there is an existing singleton*.
%
%      COMPWAV('CALLBACK') and COMPWAV('CALLBACK',hObject,...) call the
%      local function named CALLBACK in COMPWAV.M with the given input
%      arguments.
%
%      *See GUI Options on GUIDE's Tools menu.  Choose "GUI allows only one
%      instance to run (singleton)".
%
% See also: GUIDE, GUIDATA, GUIHANDLES
% Last Modified by GUIDE v2.5 08-Jul-2003 15:49:15
%
%   M. Misiti, Y. Misiti, G. Oppenheim, J.M. Poggi 01-Feb-2003.
%   Last Revision: 24-Feb-2004.
%   Copyright 1995-2004 The MathWorks, Inc.
%   $Revision: 1.1.6.2 $ $Date: 2004/03/15 22:39:57 $ 
%*************************************************************************%
%                BEGIN initialization code - DO NOT EDIT                  %
%                ----------------------------------------                 %
%*************************************************************************%
gui_Singleton = 0;
gui_State = struct('gui_Name',       mfilename, ...
                   'gui_Singleton',  gui_Singleton, ...
                   'gui_OpeningFcn', @compwav_OpeningFcn, ...
                   'gui_OutputFcn',  @compwav_OutputFcn, ...
                   'gui_LayoutFcn',  [], ...
                   'gui_Callback',   []);
if nargin & isstr(varargin{1})
    gui_State.gui_Callback = str2func(varargin{1});
end
if nargout
    [varargout{1:nargout}] = gui_mainfcn(gui_State, varargin{:});
else
    gui_mainfcn(gui_State, varargin{:});
end
%*************************************************************************%
%                END initialization code - DO NOT EDIT                    %
%*************************************************************************%
%*************************************************************************%
%                BEGIN Opening Function                                   %
%                ----------------------                                   %
% --- Executes just before compwav is made visible.                      %
%*************************************************************************%
function compwav_OpeningFcn(hObject, eventdata, handles, varargin)
% This function has no output args, see OutputFcn.
% hObject    handle to figure
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)
% varargin   unrecognized PropertyName/PropertyValue pairs from the
%            command line (see VARARGIN)
% Choose default command line output for compwav
handles.output = hObject;
% Update handles structure
guidata(hObject, handles);
% UIWAIT makes compwav wait for user response (see UIRESUME)
% uiwait(handles.Fig_PatCwtComp);
%%!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!%
% TOOL INITIALISATION Introduced manualy in the automatic generated code  %
%%!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!%
Init_Tool(hObject,eventdata,handles);
%*************************************************************************%
%                END Opening Function                                     %
%*************************************************************************%
%*************************************************************************%
%                BEGIN Output Function                                    %
%                ---------------------                                    %
% --- Outputs from this function are returned to the command line.        %
%*************************************************************************%
function varargout = compwav_OutputFcn(hObject, eventdata, handles)
% varargout  cell array for returning output args (see VARARGOUT);
% hObject    handle to figure
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)
% Get default command line output from handles structure
varargout{1} = handles.output;
%*************************************************************************%
%                END Output Function                                      %
%*************************************************************************%
%=========================================================================%
%                BEGIN Create Functions                                   %
%                ----------------------                                   %
% --- Executes during object creation, after setting all properties.      %
%=========================================================================%
function EdiPop_CreateFcn(hObject,eventdata,handles)
% hObject    handle to Edi_Object or Pop_Object (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    empty - handles not created until after all CreateFcns called
% Hint: edit controls usually have a white background on Windows.
% Hint: popupmenu controls usually have a white background on Windows.
% See ISPC and COMPUTER.
if ispc
    % set(hObject,'BackgroundColor','white');
    set(hObject,'BackgroundColor',get(0,'defaultUicontrolBackgroundColor'));
else
    set(hObject,'BackgroundColor',get(0,'defaultUicontrolBackgroundColor'));
end
%--------------------------------------------------------------------------
function Sli_CreateFcn(hObject, eventdata, handles)
% hObject    handle to Sli_Object (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    empty - handles not created until after all CreateFcns called
% Hint: slider controls usually have a light gray background, change
%       'usewhitebg' to 0 to use default.  See ISPC and COMPUTER.
usewhitebg = 1;
if usewhitebg
    set(hObject,'BackgroundColor',[.9 .9 .9]);
else
    set(hObject,'BackgroundColor',get(0,'defaultUicontrolBackgroundColor'));
end
%--------------------------------------------------------------------------
%=========================================================================%
%                END Create Functions                                     %
%=========================================================================%
%=========================================================================%
%                BEGIN Callback Functions                                 %
%                ------------------------                                 %
%=========================================================================%
% --- Executes on button press in Rad_Trans.
function Rad_TwoPatterns_Callback(hObject, eventdata, handles)
% hObject    handle to Rad_Trans (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)
set(hObject,'Value',1);
set(handles.Rad_Test_Sig,'Value',0);
% set([handles.Rad_Super,handles.Rad_Trans],'Enable','on')
set([handles.Pop_Test_Sig],'Enable','off')
%--------------------------------------------------------------------------
% --- Executes on button press in Rad_Trans.
function Rad_Trans_Callback(hObject, eventdata, handles)
% hObject    handle to Rad_Trans (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)
set(hObject,'Value',1);
set(handles.Rad_Super,'Value',0);
%--------------------------------------------------------------------------
% --- Executes on button press in Rad_Super.
function Rad_Super_Callback(hObject, eventdata, handles)
% hObject    handle to Rad_Super (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)
set(hObject,'Value',1);
set(handles.Rad_Trans,'Value',0);
%--------------------------------------------------------------------------
% --- Executes on button press in Rad_Test_Sig.
function Rad_Test_Sig_Callback(hObject, eventdata, handles)
% hObject    handle to Rad_Test_Sig (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)
set(hObject,'Value',1);
set(handles.Rad_TwoPatterns,'Value',0);
% set([handles.Rad_Super,handles.Rad_Trans],'Enable','off')
set([handles.Rad_Test_Sig,handles.Pop_Test_Sig],'Enable','on')
%--------------------------------------------------------------------------
% --- Executes on selection change in Pop_Test_Sig.
function Pop_Test_Sig_Callback(hObject, eventdata, handles)
% hObject    handle to Pop_Test_Sig (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)
%--------------------------------------------------------------------------
% --- Executes on button press in Chk_Noise.
function Chk_Noise_Callback(hObject, eventdata, handles)
% hObject    handle to Chk_Noise (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)
%--------------------------------------------------------------------------
% --- Executes on button press in Chk_Triangle.
function Chk_Triangle_Callback(hObject, eventdata, handles)
% hObject    handle to Chk_Triangle (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)
%--------------------------------------------------------------------------
% --- Executes on selection change in Pop_Wav_Fam.
function Pop_Wav_Fam_Callback(hObject, eventdata, handles)
% hObject    handle to Pop_Wav_Fam (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)
cbanapar('cba_fam',gcbf,[],guidata(gcbo));
%--------------------------------------------------------------------------
% --- Executes on selection change in Pop_Wav_Num.
function Pop_Wav_Num_Callback(hObject, eventdata, handles)
% hObject    handle to Pop_Wav_Num (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)
cbanapar('cba_num',gcbf,[],guidata(gcbo));
%--------------------------------------------------------------------------
% --- Executes on button press in Rad_Contours.
function Rad_Contours_Callback(hObject, eventdata, handles)
% hObject    handle to Rad_Contours (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)
set(hObject,'Value',1);
set(handles.Rad_Images,'Value',0);
%--------------------------------------------------------------------------
% --- Executes on slider movement.
function Sli_Thresh_Callback(hObject, eventdata, handles)
% hObject    handle to Sli_Thresh (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)
Tresh = get(hObject,'Value');
set(handles.Edi_Thresh,'string',Tresh);
% --- Executes on button press in Edi_Thresh.
function Edi_Thresh_Callback(hObject, eventdata, handles)
% hObject    handle to Edi_Thresh (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)
FieldDefault ='0.7';
Val_Thresh = str2double(get(hObject,'String'));
if Val_Thresh < 0.1 || Val_Thresh > 0.9 || isnan(Val_Thresh)
    set(hObject,'String',FieldDefault);
    set(handles.Sli_Thresh,'Value',str2double(FieldDefault));
else
    set(handles.Sli_Thresh,'Value',Val_Thresh);
end
%--------------------------------------------------------------------------
% --- Executes on button press in Rad_Images.
function Rad_Images_Callback(hObject, eventdata, handles)
% hObject    handle to Rad_Images (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)
set(hObject,'Value',1);
set(handles.Rad_Contours,'Value',0);
%--------------------------------------------------------------------------
function Pus_Run_Callback(hObject, eventdata, handles)
% hObject    handle to Pus_Run (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)
% Get figure handle.
%-------------------
hFig = handles.output;
% Cleaning.
%-----------
wwaiting('msg',hFig,'Wait ... cleaning');
cleanTOOL([                                 ...
            handles.Axe_RunSigAdapWave ,    ...
            handles.Axe_RunSigWave ,        ...
            handles.Axe_AdapWaveDetect ,    ...
            handles.Axe_WaveDetect          ...
            ]);
% Get parameter values.
%----------------------
tool_PARAMS = wtbxappdata('get',hObject,'tool_PARAMS');
X           = tool_PARAMS.X;
Y           = tool_PARAMS.Y;
PSI         = tool_PARAMS.PSI;
NC_PSI      = tool_PARAMS.NC_PSI;
% Get uicontrol values.
%----------------------
RadTrans    = get(handles.Rad_Trans,'Value');
RadSuper    = get(handles.Rad_Super,'Value');
ChkNoise    = get(handles.Chk_Noise,'Value');
ChkTriangle = get(handles.Chk_Triangle,'Value');
RadContours = get(handles.Rad_Contours,'Value');
RadImages   = get(handles.Rad_Images,'Value');
Str_Thresh  = get(handles.Edi_Thresh,'String');
Thresh      = str2double(Str_Thresh);
% Running signal construction.
%-----------------------------
titleSTR{1} = 'Running Signal';
radSIG = handles.Rad_Test_Sig;
radVAL = get(radSIG,'Value');
if radVAL==1
    popSIG = handles.Pop_Test_Sig;
    popVAL = get(popSIG,'Value');
    popSTR = get(popSIG,'String');
    sigNAM = deblank(popSTR(popVAL,:));
    switch sigNAM{1}
        case 'testsig0'
            load pthaar
        case 'testsig1'
            load ptpssin1
        case 'testsig2'
            load ptpssin2
        case 'testsig3'
            load ptsine
        case 'testsig4'
            load ptsumsin
        case 'testsig5'
            load ptsinpol
        case 'testsig6'
            load ptodtri
        case 'testsig7'
            load ptodlin
        case 'testsig8'
            load ptodpoly
    end
end
if RadTrans
    nbFormeBASE = 8; scaleBASE = 8; rapport = 2; PONDERATIONS =  [1,1];
    signalBASE = getSignal('translate',Y,rapport,PONDERATIONS,nbFormeBASE);
    lenSIG  = length(signalBASE);
    scaleMIN = 1;
    scaleMAX = 2*scaleBASE;
    NoiseLevel = 1.25;
    titleSTR{2} = 'F((t-20)/8) + sqrt(2)timesF((t-40)/4)';
else
    nbFormeBASE = 3; scaleBASE = 32; rapport = 4; PONDERATIONS =  [1,1];
    signalBASE = getSignal('superpose',Y,rapport,PONDERATIONS,nbFormeBASE);
    lenSIG  = length(signalBASE);
    scaleMIN = 1;
    scaleMAX = 2*scaleBASE;
    NoiseLevel = 0.75;   
    titleSTR{2} = 'F((t-40)/32) + 2timesF((t-40)/8)';
end
if radVAL==1, titleSTR{2} = sigNAM{1}; end
signal = signalBASE;
if ChkNoise
    NoiseLevel = NoiseLevel*std(signal);
    Noise = NoiseLevel*randn(1,lenSIG);
    signal = signal + Noise;
    titleSTR{2} = [titleSTR{2} ' + N'];
end
if ChkTriangle
    L = lenSIG/2;
    Triangle = [1:floor(L) , ceil(L):-1:1]/ceil(L);
    signal = signal + Triangle;
    titleSTR{2} = [titleSTR{2} ' + T'];
end
intervalleSIG = [0,nbFormeBASE*scaleBASE];
stepSIG = (intervalleSIG(2)-intervalleSIG(1))/(lenSIG-1);
xvalSIG = linspace(intervalleSIG(1),intervalleSIG(2),lenSIG);
stepScales = 1;
scales  = [scaleMIN:stepScales:scaleMAX];
    
% Running signal (adapted wavelet column) axes display.
%------------------------------------------------------
lw = 1;
axes(handles.Axe_RunSigAdapWave);
line(xvalSIG,signal,'color','b','linewidth',lw);
Xlim = intervalleSIG;
ext = abs(max(signal) - min(signal)) / 100;
Ylim = [min(signal)-ext max(signal)+ext];
set(handles.Axe_RunSigAdapWave,'Xlim',Xlim,'Ylim',Ylim);
if radVAL,
    [Tit1 Tit2] = deal(titleSTR{:});
    titleSTR = [Tit1 ': ' Tit2];
end
setAxesTitle(handles.Axe_RunSigAdapWave,titleSTR);
% Set the axes visible.
%----------------------
set(handles.Axe_RunSigAdapWave,'Visible','on');
% Running signal (wavelet column) axes display.
%----------------------------------------------
lw = 1;
axes(handles.Axe_RunSigWave);
line(xvalSIG,signal,'color','b','linewidth',lw);
Xlim = intervalleSIG;
ext = abs(max(signal) - min(signal)) / 100;
Ylim = [min(signal)-ext max(signal)+ext];
set(handles.Axe_RunSigWave,'Xlim',Xlim,'Ylim',Ylim);
setAxesTitle(handles.Axe_RunSigWave,titleSTR);
% Set the axes visible.
%----------------------
set(handles.Axe_RunSigWave,'Visible','on');
% Get the wavelet used.
%----------------------
Wav = cbanapar('get',hFig,'wav');                           
% Compute the cwt.
%-----------------
C_psi = cwt({signal,stepSIG},scales,PSI);
C_wav = cwt({signal,stepSIG},scales,Wav);
% Find the coordinates of max coefs.
%-----------------------------------
% [ROW_psi,COL_psi] = maxCoefs(C_psi);
% [ROW_wav,COL_wav] = maxCoefs(C_wav);
% Coefficients representation.
%-----------------------------
if RadContours
    plotCONTOUR(C_psi,scales,xvalSIG,...
        handles.Axe_AdapWaveDetect,handles.Axe_ColBar,'psi',Thresh);
%     if ~radVAL, LocalGrid(handles.Axe_AdapWaveDetect,RadTrans); end;
    LocalGrid(handles.Axe_AdapWaveDetect,RadTrans,radVAL);
    plotCONTOUR(C_wav,scales,xvalSIG,...
        handles.Axe_WaveDetect,handles.Axe_ColBar,Wav,Thresh);
%     if ~radVAL, LocalGrid(handles.Axe_WaveDetect,RadTrans); end;
    LocalGrid(handles.Axe_WaveDetect,RadTrans,radVAL);
    set([handles.Fra_ColPar,handles.Txt_PAL,handles.Txt_NBC, ...
         handles.Txt_BRI,handles.Pop_PAL,handles.Sli_NBC, ...
         handles.Edi_NBC,handles.Pus_BRI_M,handles.Pus_BRI_P],'Enable', 'Off');
elseif RadImages   
    plotCOEFS(C_psi,scales,xvalSIG,handles,...
        handles.Axe_AdapWaveDetect,handles.Axe_ColBar,'psi');
    plotCOEFS(C_wav,scales,xvalSIG,handles,...
        handles.Axe_WaveDetect,handles.Axe_ColBar,Wav);
    set([handles.Fra_ColPar,handles.Txt_PAL,handles.Txt_NBC, ...
         handles.Txt_BRI,handles.Pop_PAL,handles.Sli_NBC, ...
         handles.Edi_NBC,handles.Pus_BRI_M,handles.Pus_BRI_P],'Enable', 'On');
end
% Set the axes visible.
%----------------------
set(handles.Axe_AdapWaveDetect,'Visible','on');
set(handles.Axe_WaveDetect,'Visible','on');
% Init DynVTool.
%---------------
axe_IND = [];
axe_CMD = [...
        handles.Axe_RunSigAdapWave , ...
        handles.Axe_RunSigWave , ...
        handles.Axe_AdapWaveDetect , ...
        handles.Axe_WaveDetect ...
        ];
axe_ACT = [];
dynvtool('init',hFig,axe_IND,axe_CMD,axe_ACT,[1 0],'','','');
% End waiting.
%-------------
wwaiting('off',hFig);
%--------------------------------------------------------------------------
% --- Executes on button press in Pus_CloseWin.
function Pus_CloseWin_Callback(hObject, eventdata, handles)
% hObject    handle to Pus_CloseWin (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)
% Get nwavetool figure parameters.
%---------------------------------
tool_PARAMS = wtbxappdata('get',hObject,'tool_PARAMS');
% Enable the compare button on in nwavetool figure.
%--------------------------------------------------
set(tool_PARAMS.Pus_Compare,'Enable','On');
% Close current figure.
%----------------------
close(gcbf)
%--------------------------------------------------------------------------
%=========================================================================%
%                END Callback Functions                                   %
%=========================================================================%
%=========================================================================%
%                BEGIN Special Callback Functions                         %
%                --------------------------------                         %
%=========================================================================%
% --- Executes on button press in Compare from nwavtool figure.
function NwavtoolCall_Callback(hObject, eventdata, handles,Nwavetool_PARAMS)
% hObject    handle to Pus_Compare (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)
% Save nwavtool information in 'Userdata' field of compwav figure.
%-----------------------------------------------------------------
if isfield(Nwavetool_PARAMS,'demoMODE')
    F1 = allchild(0);
end
compwav('userdata',Nwavetool_PARAMS);
if isfield(Nwavetool_PARAMS,'demoMODE')
    F2 = allchild(0);
    hFig = setdiff(F2,F1);
    set(hFig,'HandleVisibility','On')
end
%=========================================================================%
%                END Special Callback Functions                           %
%=========================================================================%
%=========================================================================%
%                BEGIN Tool Initialization                                %
%                -------------------------                                %
%=========================================================================%
function Init_Tool(hObject, eventdata, handles)
% WTBX -- Install DynVTool.
%--------------------------
dynvtool('Install_V3',hObject,handles);
% WTBX -- Install MENUS.
%-----------------------
wfigmngr('extfig',hObject,'ExtFig_GUIDE');
wfigmngr('attach_close',hObject);
set(hObject,'HandleVisibility','On');
hdl_Menus = Install_MENUS(hObject);
% WTBX -- Install ANAPAR FRAME.
%------------------------------
wnameDEF  = 'db1';
utanapar('Install_V3_CB',hObject,'wtype','cwt');
set(hObject,'Visible','On'); pause(0.01) %%% MiMi : BUG MATLAB %%%
cbanapar('set',hObject,'wav',wnameDEF);
% WTBX -- Install COLORMAP FRAME
%-------------------------------
utcolmap('Install_V3',hObject,'enable','On');
default_nbcolors = 240;
cbcolmap('set',hObject,'pal',{'jet',default_nbcolors});
% Get nwavetool figure parameters.
%---------------------------------
Nwavetool_PARAMS = get(hObject,'Userdata');
% Initialize uicontrols.
%-----------------------
set(handles.Rad_Trans,'Value',Nwavetool_PARAMS.RadTrans);
set(handles.Rad_Super,'Value',Nwavetool_PARAMS.RadSuper);
set(handles.Chk_Noise,'Value',Nwavetool_PARAMS.ChkNoise);
set(handles.Chk_Triangle,'Value',Nwavetool_PARAMS.ChkTriangle);
set(Nwavetool_PARAMS.Pus_Compare,'Enable','Off');
% Set current tool figure parameters.
%------------------------------------
tool_PARAMS.X       = Nwavetool_PARAMS.X;
tool_PARAMS.Y       = Nwavetool_PARAMS.Y;
tool_PARAMS.PSI     = Nwavetool_PARAMS.PSI;
tool_PARAMS.NC_PSI  = Nwavetool_PARAMS.NC_PSI;
tool_PARAMS.Pus_Compare = Nwavetool_PARAMS.Pus_Compare;
wtbxappdata('set',hObject,'tool_PARAMS',tool_PARAMS);
% Set Title in the RunSigAdapWave axes.
%--------------------------------------
titleSTR = 'Running Signal';
setAxesTitle(handles.Axe_RunSigAdapWave,titleSTR);
% Set Title in the RunSigWave axes.
%----------------------------------
titleSTR = 'Running Signal';
setAxesTitle(handles.Axe_RunSigWave,titleSTR);
% Set Title in the AdapWaveDetect signal axes.
%---------------------------------------------
titleSTR = 'Adapted Wavelet';
setAxesTitle(handles.Axe_AdapWaveDetect,titleSTR);
% Set Title in the pattern detection axes.
%-----------------------------------------
titleSTR = 'Wavelet';
setAxesTitle(handles.Axe_WaveDetect,titleSTR);
% Resize Figure.
%---------------------
redimfig('On',hObject);
% Initialize colorbar.
%---------------------
cmap = get(hObject,'Colormap');
axes(handles.Axe_ColBar);
image([1:size(cmap,1)])
set(handles.Axe_ColBar,...
    'Xtick',[],'Xticklabel',[],'Ytick',[],'Yticklabel',[] ...    
    );
dynvzaxe('exclude',hObject,handles.Axe_ColBar)
ColorBarVisibility(handles.Axe_ColBar,'Off');
% First compare display.
%-----------------------
compwav('Pus_Run_Callback',hObject, eventdata, handles);
ColorBarVisibility(handles.Axe_ColBar,'On');
% End Of initialization.
%-----------------------
set(hObject,'HandleVisibility','Callback');
%=========================================================================%
%                END Tool Initialization                                  %
%=========================================================================%
%=========================================================================%
%                BEGIN CleanTOOL function                                 %
%                ------------------------                                 %
%=========================================================================%
function cleanTOOL(handles)
hLINES  = findobj(handles,'type','line');
delete(hLINES);
hLINES  = findobj(handles,'type','patch');
delete(hLINES);
hIMAGES = findobj(handles,'type','Image');
delete(hIMAGES);
%--------------------------------------------------------------------------
function ColorBarVisibility(Axe_ColBar,status)
hIMAGES = findobj(Axe_ColBar,'type','image');
set(Axe_ColBar,'Visible',status);
set(hIMAGES,'Visible',status);
%-------------------------------------------------------------------------
%=========================================================================%
%                END CleanTOOL function                                   %
%=========================================================================%
%=========================================================================%
%                BEGIN Internal Functions                                 %
%                ------------------------                                 %
%=========================================================================%
function hdl_Menus = Install_MENUS(hFig)
% Add Help for Tool.
%------------------
% wfighelp('addHelpTool',hFig,'&Continuous Analysis','CW1D_GUI');
% Add Help Item.
%----------------
% wfighelp('addHelpItem',hFig,'Continuous Transform','CW_TRANSFORM');
hdl_Menus = [];
%-------------------------------------------------------------------------
function signalBASE = getSignal(typeSIG,FormeBASE,rapport,PONDERATIONS,nbFormeBASE)
Forme_1 = FormeBASE;
len_F1  = length(Forme_1);
x_new   = linspace(0,1,len_F1/rapport);
x_old   = linspace(0,1,len_F1);        
Forme_2 = (rapport^0.5)*interp1(x_old,Forme_1,x_new);
len_F2  = length(Forme_2);
Forme_1 = PONDERATIONS(1)*Forme_1;
Forme_2 = PONDERATIONS(2)*Forme_2;
signalBASE = zeros(1,nbFormeBASE*len_F1);
switch typeSIG
	case 'superpose'
        deb = floor(((nbFormeBASE-1)/2-1/4)*len_F1) + 1; fin = deb + len_F1-1;
        signalBASE(deb:fin) = Forme_1;
        deb = deb + floor((len_F1-len_F2)/2); fin = deb + len_F2-1;
        signalBASE(deb:fin) = signalBASE(deb:fin) + Forme_2;
	case 'translate'
        deb = 2*len_F1; fin = deb + len_F1-1;
        signalBASE(deb:fin) = Forme_1;
        deb = floor(5*len_F1-len_F2/2); fin = deb + len_F2-1;        
        signalBASE(deb:fin) = signalBASE(deb:fin) + Forme_2;
end
%-------------------------------------------------------------------------
function plotCONTOUR(C,scales,xvalSIG,Axe,Axe_ColBar,Wav,Thresh)
% Compute and plot contours.
%---------------------------
axes(Axe);
maxi = max(max(abs(C))); 
D = abs(C)/maxi;
% xval = [0.7:0.025:1];
xval = [Thresh:0.025:1];
contour(xvalSIG,scales,D,xval);
% Set Title in the pattern detection axes.
%-----------------------------------------
if isequal(Wav,'psi')
    titleSTR = 'Adapted Wavelet';
else
    titleSTR = [Wav,' Wavelet'];
end
setAxesTitle(Axe,titleSTR);
setAxesATTRB(scales,xvalSIG,Axe);
% Set Xtick Values.
%------------------
NBC = 64;
xV = [xval(1),xval(end)]; 
setXTicks(Axe_ColBar,'real',NBC,xV);
%-------------------------------------------------------------------------
function plotCOEFS(C,scales,xvalSIG,handles,Axe,Axe_ColBar,Wav)
axes(Axe);
NBC = 240;
abs_mode = 1;
C = wcodemat(C,NBC,'mat',abs_mode);
cbcolmap('set',handles.output,'pal',{'jet',NBC});
xL = get(handles.Axe_RunSigWave,'Xlim');
yL = [scales(1) scales(end)];
img  = image(xL,yL,C);
% Set Title in the pattern detection axes.
%-----------------------------------------
if isequal(Wav,'psi')
    titleSTR = 'Adapted Wavelet';
else
    titleSTR = [Wav,' Wavelet'];
end
setAxesTitle(Axe,titleSTR);
setAxesATTRB(scales,xvalSIG,Axe);
setXTicks(Axe_ColBar,'real',NBC,[0 1],5);
%-------------------------------------------------------------------------
function plotSURF(C,scales,xvalSIG,Axe,Axe_ColBar,Wav)
axes(Axe);
NBC = 240;
colormap(jet(NBC));
maxi = max(max(abs(C)));
D = abs(C)/maxi;
s = surf(xvalSIG,scales,D); 
shading interp
set(Axe,'View',[-18 72],'Xlim',[xvalSIG(1),xvalSIG(end)],'Ylim',[1,size(C,1)]);
% Set Title in the pattern detection axes.
%-----------------------------------------
if isequal(Wav,'psi')
    titleSTR = 'Adapted Wavelet';
else
    titleSTR = [Wav,' Wavelet'];
end
setAxesTitle(Axe,titleSTR);
setAxesATTRB(scales,xvalSIG,Axe);
setXTicks(Axe_ColBar,'int',NBC,[0 NBC]);
%-------------------------------------------------------------------------
function setXTicks(Axe_ColBar,type,NBC,xV,nbTICS)
axes(Axe_ColBar);
colormap(jet(NBC))
image([1:NBC])
if nargin<5 , nbTICS = 7; end
xlim = get(Axe_ColBar,'Xlim');
alfa = (xlim(2)-xlim(1))/(xV(2)-xV(1));
beta = xlim(1)-alfa*xV(1);
xCOL  = linspace(xV(1),xV(2),nbTICS);
xtics = alfa*xCOL + beta;
switch type
    case 'int'  , xlabs = int2str(xCOL(:));
    case 'real' , xlabs = num2str(xCOL(:),2);
end
set(Axe_ColBar, ...
        'XTick',xtics,'XTickLabel',xlabs, ...
        'YTick',[],'YTickLabel',[] ...
        );
%-------------------------------------------------------------------------
function setAxesATTRB(scales,xvalSIG,axe)
nb_SCALES = length(scales);
nb    = min(5,nb_SCALES);
level = '';
for k=1:nb , level = [level ' '  num2str(scales(k))]; end
if nb<nb_SCALES , level = [level ' ...']; end
nb    = ceil(nb_SCALES/20);
ytics = nb:nb:nb_SCALES;
tmp   = scales(nb:nb:nb*length(ytics));
ylabs = num2str(tmp(:));
set(axe, ...
        'YTick',ytics,          ...
        'YTickLabel',ylabs,     ...
        'YDir','normal',        ...
        'Box','On'              ...
        );
%-------------------------------------------------------------------------
function [ROW,COL] = maxCoefs(C)
[nb_ROW,nb_COL] = size(C);
D = abs(C(:));
[dummy,idx] = max(D);
col_INF = fix(idx/nb_ROW);
ROW = idx - nb_ROW*col_INF;
if ROW~=0
    COL = col_INF+1;
else
    COL = col_INF; ROW = nb_ROW;
end
%-------------------------------------------------------------------------
function LocalGrid(Axe_Detect,RadTrans,radVAL)
axes(Axe_Detect)
% if RadTrans
%     if radVAL
%         line([10 10],[0 8],'color','k','LineStyle','--');
%         line([0 10],[8 8],'color','k','LineStyle','--');
%         line([30 30],[0 4],'color','k','LineStyle','--');
%         line([0 30],[4 4],'color','k','LineStyle','--');
%         line([50 50],[0 4],'color','k','LineStyle','--');
%         line([0 50],[4 4],'color','k','LineStyle','--');
%         line(10,8,'color','k','Marker','*');
%         line(30,4,'color','k','Marker','*');
%         line(50,4,'color','k','Marker','*');
%     else
%         line([20 20],[0 8],'color','k','LineStyle','--');
%         line([0 20],[8 8],'color','k','LineStyle','--');
%         line([40 40],[0 4],'color','k','LineStyle','--');
%         line([0 40],[4 4],'color','k','LineStyle','--');
%         line(20,8,'color','k','Marker','*');
%         line(40,4,'color','k','Marker','*');
%     end
% else
%         line([40 40],[0 32],'color','k','LineStyle','--');
%         line([0 40],[32 32],'color','k','LineStyle','--');
%         line([40 40],[0 8],'color','k','LineStyle','--');
%         line([0 40],[8 8],'color','k','LineStyle','--');
%         line(40,32,'color','k','Marker','*');
%         line(40,8,'color','k','Marker','*');
% end
% if ~radVAL
    if RadTrans
        line([20 20],[0 8],'color','k','LineStyle','--');
        line([0 20],[8 8],'color','k','LineStyle','--');
        line([40 40],[0 4],'color','k','LineStyle','--');
        line([0 40],[4 4],'color','k','LineStyle','--');
        line(20,8,'color','k','Marker','*');
        line(40,4,'color','k','Marker','*');
    else
        line([40 40],[0 32],'color','k','LineStyle','--');
        line([0 40],[32 32],'color','k','LineStyle','--');
        line([40 40],[0 8],'color','k','LineStyle','--');
        line([0 40],[8 8],'color','k','LineStyle','--');
        line(40,32,'color','k','Marker','*');
        line(40,8,'color','k','Marker','*');
    end
% else
%         line([20 20],[0 8],'color','k','LineStyle','--');
%         line([0 20],[8 8],'color','k','LineStyle','--');
%         line([40 40],[0 4],'color','k','LineStyle','--');
%         line([0 40],[4 4],'color','k','LineStyle','--');
%         line(20,8,'color','k','Marker','*');
%         line(40,4,'color','k','Marker','*');
%     if strcmp(sigNAM(end),'t')
%         line([10 10],[0 8],'color','k','LineStyle','--');
%         line([0 10],[8 8],'color','k','LineStyle','--');
%         line([30 30],[0 4],'color','k','LineStyle','--');
%         line([0 30],[4 4],'color','k','LineStyle','--');
%         line([50 50],[0 4],'color','k','LineStyle','--');
%         line([0 50],[4 4],'color','k','LineStyle','--');
%         line(10,8,'color','k','Marker','*');
%         line(30,4,'color','k','Marker','*');
%         line(50,4,'color','k','Marker','*');
        %    elseif strcmp(sigNAM(end),'s')
%     else
%         line([40 40],[0 32],'color','k','LineStyle','--');
%         line([0 40],[32 32],'color','k','LineStyle','--');
%         line([40 40],[0 8],'color','k','LineStyle','--');
%         line([0 40],[8 8],'color','k','LineStyle','--');
%         line(40,32,'color','k','Marker','*');
%         line(40,8,'color','k','Marker','*');
%     end
% end
%-------------------------------------------------------------------------
%=========================================================================%
%                END Internal Functions                                   %
%=========================================================================%
%=========================================================================%
%                BEGIN General Utilities                                  %
%                -----------------------                                  %
%=========================================================================%
function setAxesTitle(axe,label)
axes(axe); 
title(label,'Color','k','FontWeight','demi','Units','normalized','Fontsize',9);
%--------------------------------------------------------------------------
function setAxesXlabel(axe,label)
axes(axe); 
xlabel(label,'Color','k','FontWeight','demi','Units','normalized','Fontsize',9);
%--------------------------------------------------------------------------
function setAxesYlabel(axe,label)
axes(axe); 
ylabel(label,'Color','k','FontWeight','demi','Units','normalized','Fontsize',9);
%--------------------------------------------------------------------------
%=========================================================================%
%                END Tool General Utilities                               %
%=========================================================================%
%=========================================================================%
%                      BEGIN Demo Utilities                               %
%                      ---------------------                              %
%=========================================================================%
function closeDEMO(hFig,eventdata,handles,varargin)
close(hFig);
%----------------------------------------------------------
function demoPROC(hFig,eventdata,handles,varargin)
waveDEM = varargin{1};
handles = guidata(hFig);
cbanapar('set',hFig,'wav',waveDEM);
Pus_Run = handles.Pus_Run; 
Pus_Run_Callback(Pus_Run,[],handles);
%=========================================================================%
%                   END Tool Demo Utilities                               %
%=========================================================================%

						