www.gusucode.com > wavelet工具箱matlab源码程序 > wavelet/wavelet/wmoreres.m
function varargout = wmoreres(option,varargin)
%WMORERES "More information" on wavelet residuals tool.
% VARARGOUT = WMORERES(OPTION,VARARGIN)
%
%------------------------------------------------------------
% Internal options:
% OPTION = 'create', 'select', 'update_bins', 'close'
% M. Misiti, Y. Misiti, G. Oppenheim, J.M. Poggi 01-Jul-98.
% Last Revision: 04-Jul-2013.
% Copyright 1995-2013 The MathWorks, Inc.
% $Revision: 1.8.4.17 $ $Date: 2013/08/23 23:45:27 $
% MemBloc1 of stored values.
%---------------------------
n_membloc1 = [mfilename '_MB1'];
ind_loc_struct = 1;
% Tag property.
%--------------
tag_figRes = 'Fig_Residual';
% Parameters initialization.
%---------------------------
default_bins = 50;
curr_color = wtbutils('colors','res');
nb_bins = default_bins;
if ~isequal(option,'create') , fig = varargin{1}; end
switch option
case 'create'
%**************************************************************%
%** OPTION = 'create' - Create the current figure **%
%**************************************************************%
caller = varargin{1};
oldFig = wfindobj('figure','Tag',tag_figRes);
existFig = ~isempty(oldFig);
if existFig
existFig = 0;
for k = 1:length(oldFig)
ls = wmemtool('rmb',oldFig(k),n_membloc1,ind_loc_struct);
if isequal(ls.caller,caller)
existFig = oldFig(k);
break
end
end
end
if existFig ~= 0
try delete(existFig); catch , end %#ok<*CTCH>
return
end
% Get Globals.
%-------------
[Def_Txt_Height,Def_Btn_Height,Def_Btn_Width,Pop_Min_Width, ...
X_Spacing,Y_Spacing,Def_EdiBkColor, ...
Def_FraBkColor,Def_ShadowColor] = ...
mextglob('get',...
'Def_Txt_Height','Def_Btn_Height','Def_Btn_Width', ...
'Pop_Min_Width','X_Spacing','Y_Spacing', ...
'Def_EdiBkColor','Def_FraBkColor','Def_ShadowColor' ...
);
% Window initialization.
%-----------------------
CallerTitle = get(caller,'Name');
strTitle = getWavMSG('Wavelet:commongui:Lab_MoreOnRes',CallerTitle);
[fig,pos_win,win_units,str_numwin,...
pos_frame0,Pos_Graphic_Area] = ...
wfigmngr('create',strTitle,[],'ExtFig_Tool_3',mfilename,1,1,0);
varargout{1} = fig;
set(fig,'Tag',tag_figRes);
% Add Help for Tool.
%------------------
% See Below, after detection of typeDATA
% Menu construction for current figure.
%--------------------------------------
m_files = wfigmngr('getmenus',fig,'file');
uimenu(m_files,...
'Label',getWavMSG('Wavelet:commongui:Lab_SaveRes'), ...
'Position',1, ...
'Enable','On', ...
'Callback', ...
[mfilename '(''save'',' str_numwin ');'] ...
);
% Begin waiting.
%---------------
set(wfindobj('figure'),'Pointer','watch');
% Setting local structure.
%-------------------------
ls = struct('caller',caller,'toggle',varargin{2},...
'handleRES',varargin{3},'handleORI',[],'handleTHR',[]);
Nb_Hdls_Inputs = length(varargin)-3;
blocPAR = 0;
if Nb_Hdls_Inputs>0
ls.handleORI = varargin{4};
if Nb_Hdls_Inputs>1
ls.handleTHR = varargin{5};
if Nb_Hdls_Inputs>2 , blocPAR = 1; end
end
end
if ishandle(ls.handleRES)
[typeDATA,propDATA] = getDataType(ls.handleRES);
resVal = wd2uiorui2d('uint2d',get(ls.handleRES,propDATA));
elseif ishandle(ls.handleORI)
[typeDATA,propDATA] = getDataType(ls.handleORI);
resVal = wd2uiorui2d('uint2d',get(ls.handleORI,propDATA));
if ishandle(ls.handleTHR)
if isequal(propDATA,'YData');
xdata = get(ls.handleORI,'XData');
xTHR = get(ls.handleTHR,'XData');
if ~isequal(xdata,xTHR)
typeDATA = 'regress';
yTHR = get(ls.handleTHR,'YData');
resVal = resVal-interp1(xTHR,yTHR,xdata);
else
resVal = resVal-...
wd2uiorui2d('uint2d',get(ls.handleTHR,propDATA));
end
else
resVal = resVal-...
wd2uiorui2d('uint2d',get(ls.handleTHR,propDATA));
end
else
resVal = zeros(size(resVal));
end
else
error(message('Wavelet:FunctionArgVal:Invalid_Handle'))
end
wmemtool('wmb',fig,n_membloc1,ind_loc_struct,ls);
% Add Help for Tool.
%-------------------
switch typeDATA
case {'line','regress'}
helpName = 'RESI1D_GUI';
helpItem = getWavMSG('Wavelet:commongui:HLP_1D_Res');
case 'image'
helpName = 'RESI2D_GUI';
helpItem = getWavMSG('Wavelet:commongui:HLP_2D_Res');
otherwise
helpName = 'Error';
end
wfighelp('addHelpTool',fig,helpItem,helpName);
% Creating the Command part of the window.
%=========================================
% Building the "Data, Wavelet and Level" block.
%----------------------------------------------
if blocPAR , utanapar('create_copyB',caller,fig); end
% General parameters initialization.
%-----------------------------------
dx = X_Spacing; dx2 = 2*dx;
dy = Y_Spacing; dy2 = 2*dy;
d_txt = (Def_Btn_Height-Def_Txt_Height);
x_frame0 = pos_frame0(1);
xlocINI = pos_frame0([1 3]);
cmd_width = pos_frame0(3);
push_width = (cmd_width-4*dx)/2;
pop_width = Pop_Min_Width;
ybottomINI = pos_win(4)-3.5*Def_Btn_Height-dy2;
% Position property of objects.
%------------------------------
x_fra = x_frame0+3*dx/2;
y_fra = ybottomINI-13*Def_Btn_Height;
w_fra = cmd_width-3*dx;
h_fra = 10*Def_Btn_Height;
pos_fra = [x_fra y_fra w_fra h_fra];
y_low = ybottomINI-11*Def_Btn_Height/4;
w_uic = 15*push_width/8;
px = x_frame0+(cmd_width-w_uic)/2;
pos_hist_chist = [px, y_low-3*Def_Btn_Height, ...
w_uic, 3*Def_Btn_Height/2];
y_low = pos_hist_chist(2)-3*Def_Btn_Height;
pos_corr_spec = [px, y_low , 1.04*w_uic, 3*Def_Btn_Height/2];
y_low = pos_corr_spec(2)-3*Def_Btn_Height;
pos_info_stat = [px, y_low , w_uic, 3*Def_Btn_Height/2];
y_low = pos_fra(2)-2*Def_Btn_Height;
pos_txt_bin = [px, y_low+d_txt/2, 2*pop_width, ...
Def_Txt_Height];
pos_edi_bin = [px+2*pop_width+dx, y_low, pop_width, ...
Def_Btn_Height];
y_low = ybottomINI-3*Def_Btn_Height-Def_Btn_Height/2;
w_uic = 1.2*(Def_Btn_Width+dx2);
x_uic = x_fra+(w_fra-w_uic)/2;
pos_txt_axe = [x_uic, y_low, w_uic, Def_Txt_Height];
% String property of objects.
%----------------------------
str_hist_chist = getWavMSG('Wavelet:commongui:Str_His_CHis');
str_txt_bin = getWavMSG('Wavelet:commongui:Str_NbBins');
str_edi_bin = sprintf('%.0f',default_bins);
str_corr_spec = getWavMSG('Wavelet:commongui:Str_AutoSpec');
str_info_stat = getWavMSG('Wavelet:commongui:Str_DesStat');
str_txt_axe = getWavMSG('Wavelet:commongui:Str_SelAxes');
% Callbacks property of objects.
%-------------------------------
cba_chk = [mfilename '(''select'',' str_numwin ');'];
cba_edi_bin = [mfilename '(''update_bins'',' str_numwin ');'];
% Uicontrols definition.
%-----------------------
commonProp = {'Parent',fig,'Units',win_units};
comChkProp = [commonProp,'Style','Checkbox'];
comFraProp = [commonProp, ...
'BackgroundColor',Def_FraBkColor, ...
'ForegroundColor',Def_ShadowColor, ...
'Style','frame' ...
];
comPropTxtCENT = [commonProp, ...
'Style','text', ...
'HorizontalAlignment','center', ...
'BackgroundColor',Def_FraBkColor ...
];
fra_utl = uicontrol(comFraProp{:}, ...
'Style','frame', ...
'Position',pos_fra, ...
'Visible','off' ...
);
txt_axe = uicontrol(comPropTxtCENT{:}, ...
'Position',pos_txt_axe, ...
'String',str_txt_axe, ...
'Visible','off' ...
);
chk_hist_chist = uicontrol(comChkProp{:}, ...
'Position',pos_hist_chist, ...
'String',str_hist_chist, ...
'TooltipString', ...
getWavMSG('Wavelet:commongui:Str_HandCH'), ...
'UserData',1, ...
'Value',1, ...
'Tag','Hist_CHist', ...
'Callback',cba_chk, ...
'Visible','off' ...
);
txt_bin = uicontrol(commonProp{:}, ...
'Style','text', ...
'Position',pos_txt_bin, ...
'String',str_txt_bin, ...
'BackgroundColor',Def_FraBkColor, ...
'Tag','Txt_NumOfBins', ...
'Visible','off' ...
);
edi_bin = uicontrol(commonProp{:}, ...
'Style','Edit', ...
'Position',pos_edi_bin, ...
'String',str_edi_bin, ...
'BackgroundColor',Def_EdiBkColor, ...
'Tag','Edi_NumOfBins', ...
'Callback',cba_edi_bin, ...
'Visible','off' ...
);
chk_corr_spec = uicontrol(comChkProp{:}, ...
'Position',pos_corr_spec, ...
'String',str_corr_spec, ...
'UserData',1, ...
'Value',1, ...
'Tag','Autocorr', ...
'Callback',cba_chk, ...
'Visible','off' ...
);
chk_info_stat = uicontrol(comChkProp{:}, ...
'Position',pos_info_stat, ...
'String',str_info_stat, ...
'UserData',1, ...
'Value',1, ...
'Tag','DesStat', ...
'Callback',cba_chk, ...
'Visible','off' ...
);
% End of Uicontrols definitions depending on signal or image case.
%-----------------------------------------------------------------
switch typeDATA
case 'line'
set([chk_info_stat;chk_corr_spec;chk_hist_chist;fra_utl;txt_axe]...
,'Visible','On');
case 'image'
% Adding colormap GUI.
%---------------------
pos_txt_bin(2) = pos_txt_bin(2)+10*Def_Btn_Height;
pos_edi_bin(2) = pos_edi_bin(2)+10*Def_Btn_Height;
set(txt_bin,'Position',pos_txt_bin);
set(edi_bin,'Position',pos_edi_bin);
try
pop_pal_caller = cbcolmap('get',caller,'pop_pal');
ok = true;
catch
ok = false;
end
if ok
prop_pal = get(pop_pal_caller,{'String','Value','UserData'});
utcolmap('create',fig, ...
'xloc',xlocINI, ...
'bkcolor',Def_FraBkColor, ...
'Enable','on');
pop_pal_loc = cbcolmap('get',fig,'pop_pal');
set(pop_pal_loc,'String',prop_pal{1},'Value',prop_pal{2}, ...
'UserData',prop_pal{3});
set(fig,'Colormap',get(caller,'Colormap'));
end
case 'regress'
pos_txt_bin(2) = pos_txt_bin(2)+10*Def_Btn_Height;
pos_edi_bin(2) = pos_edi_bin(2)+10*Def_Btn_Height;
set(txt_bin,'Position',pos_txt_bin);
set(edi_bin,'Position',pos_edi_bin);
end
set([txt_bin;edi_bin],'Visible','On');
% Creating the Graphical part of the window.
%===========================================
% Frame Stats. construction.
%---------------------------
infos_hdls = utstats('create',fig,...
'xloc',Pos_Graphic_Area([1,3]), ...
'bottom',Pos_Graphic_Area(2)+dy2);
fra_sta = infos_hdls(1);
pos_fra_sta = get(fra_sta,'Position');
% Axes Positions.
%----------------
gra_width = Pos_Graphic_Area(3);
xspace = gra_width/12;
ecy_up = 0.06*pos_win(4);
ecy_down = 0.06*pos_win(4);
ecy_mid = 0.08*pos_win(4);
y_rem = pos_frame0(4)-(pos_fra_sta(2)+pos_fra_sta(4)) ...
-ecy_up-ecy_down;
switch typeDATA
case 'line' , y_rem = y_rem-2*ecy_mid;
case 'image' , y_rem = y_rem-ecy_mid;
case 'regress' , y_rem = y_rem-ecy_mid;
end
axe_width = gra_width-2*xspace;
half_width = axe_width/2-xspace/2;
switch typeDATA
case 'line'
axe_height = y_rem/3;
y_low = pos_fra_sta(2)+pos_fra_sta(4)+ecy_down;
pos_ax_corr = [xspace y_low half_width axe_height];
pos_ax_spec = [2*xspace+half_width y_low half_width axe_height];
y_low = y_low+axe_height+ecy_mid;
pos_ax_hist = [xspace y_low half_width axe_height];
pos_ax_cumhist = [2*xspace+half_width y_low half_width axe_height];
y_low = y_low+axe_height+ecy_mid;
pos_ax_signal = [xspace y_low axe_width axe_height];
case 'image'
axe_height = y_rem/3;
y_low = pos_fra_sta(2)+pos_fra_sta(4)+ecy_down;
pos_ax_hist = [xspace y_low half_width axe_height];
pos_ax_cumhist = [2*xspace+half_width y_low half_width axe_height];
y_low = y_low+axe_height+ecy_mid;
axe_height = 2*y_rem/3;
cx = xspace+axe_width/2;
cy = y_low+axe_height/2;
imageSize = size(resVal);
imageSize = imageSize([2 1]);
[w_use,h_use] = wpropimg(imageSize,axe_width,axe_height,'pixels');
pos_ax_signal = [cx-w_use/2,cy-h_use/2,w_use,h_use];
case 'regress'
axe_height = y_rem/2;
y_low = pos_fra_sta(2)+pos_fra_sta(4)+ecy_down;
pos_ax_hist = [xspace y_low half_width axe_height];
pos_ax_cumhist = [2*xspace+half_width y_low half_width axe_height];
y_low = y_low+axe_height+ecy_mid;
pos_ax_signal = [xspace y_low axe_width axe_height];
end
% Axes Definitions.
%------------------
commonProp = {...
'Parent',fig,...
'Units',win_units,...
'Visible','Off',...
'box','on',...
'NextPlot','Replace' ...
};
axe_signal = axes(commonProp{:},'Position',pos_ax_signal);
axe_hist = axes(commonProp{:},'Position',pos_ax_hist);
axe_cumhist = axes(commonProp{:},'Position',pos_ax_cumhist);
switch typeDATA
case 'line'
axe_corr = axes(commonProp{:},'Position',pos_ax_corr);
axe_spec = axes(commonProp{:},'Position',pos_ax_spec);
end
% Displaying signal or image.
%----------------------------
switch typeDATA
case {'line','regress'}
len = length(resVal);
if ishandle(ls.handleORI)
xdata = get(ls.handleORI,'XData');
elseif ishandle(ls.handleTHR)
xdata = get(ls.handleTHR,'XData');
else
xdata = (1:len);
end
hdl_Sig = line('XData',xdata,'YData',resVal,...
'Color',curr_color,'Visible','Off','Parent',axe_signal);
set(axe_signal,'XLim',[xdata(1),xdata(end)])
case 'image'
X = wd2uiorui2d('d2uint',wcodemat(resVal,128,'mat',0));
hdl_Sig = image(X,'Parent',axe_signal);
end
wtitle(getWavMSG('Wavelet:commongui:Str_Residuals'),'Parent',axe_signal);
% Setting units to normalized.
%-----------------------------
wfigmngr('normalize',fig);
% Displaying histogram.
%----------------------
his = wgethist(resVal(:),nb_bins);
[xx,imod] = max(his(2,:)); %#ok<ASGLU>
mode_val = (his(1,imod)+his(1,imod+1))/2;
his(2,:) = his(2,:)/length(resVal(:));
wplothis(axe_hist,his,curr_color);
wtitle(getWavMSG('Wavelet:commongui:Str_HIST'),'Parent',axe_hist);
% Displaying cumulated histogram.
%--------------------------------
for i=6:4:length(his(2,:));
his(2,i) = his(2,i)+his(2,i-4);
his(2,i+1) = his(2,i);
end
wplothis(axe_cumhist,[his(1,:);his(2,:)],curr_color);
wtitle(getWavMSG('Wavelet:commongui:Str_CumHist'),'Parent',axe_cumhist);
switch typeDATA
case 'line'
% Displaying Autocorrelations.
%-----------------------------
[corr,lags] = wautocor(resVal);
lenLagsPos = (length(lags)-1)/2;
lenKeep = min(200,lenLagsPos);
first = lenLagsPos+1-lenKeep;
last = lenLagsPos+1+lenKeep;
Xval = lags(first:last);
Yval = corr(first:last);
line('XData',Xval,'YData',Yval,'Color',curr_color,...
'Parent',axe_corr);
set(axe_corr,'XLim',[Xval(1) Xval(end)],...
'YLim',[min(0,1.1*min(Yval)) 1]);
wtitle(getWavMSG('Wavelet:commongui:Str_AutoCor'),'Parent',axe_corr);
% Displaying Spectrum.
%---------------------
[sp,f] = wspecfft(resVal);
line('XData',f,'YData',sp,'Color',curr_color,'Parent',axe_spec);
set(axe_spec,'XLim',[min(f) max(f)]);
wtitle(getWavMSG('Wavelet:commongui:Str_Spectrum'),'Parent',axe_spec);
xlabel(getWavMSG('Wavelet:commongui:Str_Frequency'),'Parent',axe_spec);
ylabel(getWavMSG('Wavelet:commongui:Str_Energy'),'Parent',axe_spec);
end
% Displaying statistics.
%-----------------------
switch typeDATA
case {'line','regress'}
errtol = 1.0E-12;
mean_val = mean(resVal);
max_val = max(resVal);
min_val = min(resVal);
range_val = max_val-min_val;
std_val = std(resVal);
med_val = median(resVal);
case 'image'
errtol = 1.0E-12;
mean_val = mean(resVal(:));
if abs(mean_val)<errtol , mean_val = 0; end
max_val = max(resVal(:));
if abs(max_val)<errtol , max_val = 0; end
min_val = min(resVal(:));
if abs(min_val)<errtol , min_val = 0; end
range_val = max_val-min_val;
if abs(range_val)<errtol , range_val = 0; end
std_val = std(resVal(:));
if abs(std_val)<errtol , std_val = 0; end
med_val = median(resVal(:));
if abs(med_val)<errtol , med_val = 0; end
end
medDev_val = median(abs(resVal(:)-med_val));
if abs(medDev_val)<errtol , medDev_val = 0; end
meanDev_val = mean(abs(resVal(:)-mean_val));
if abs(meanDev_val)<errtol , meanDev_val = 0; end
L1_val = norm(resVal(:),1);
L2_val = norm(resVal(:),2);
LM_val = norm(resVal(:),Inf);
utstats('display',fig, ...
[mean_val; med_val ; mode_val; ...
max_val ; min_val ; range_val; ...
std_val ; medDev_val; meanDev_val; ...
L1_val ; L2_val ; LM_val]);
% Set axes visible.
%------------------
set([infos_hdls;axe_signal;axe_hist;axe_cumhist;hdl_Sig],...
'Visible','on');
switch typeDATA
case 'line' , set([axe_corr;axe_spec],'Visible','on');
end
% Store values.
%--------------
group_hdls = struct(...
'chk_hist_chist', chk_hist_chist, ...
'txt_bin', txt_bin, ...
'edi_bin', edi_bin, ...
'chk_corr_spec', chk_corr_spec, ...
'chk_info_stat', chk_info_stat ...
);
set(axe_signal,'UserData',get(axe_signal,'Position'));
set(axe_hist,'UserData',get(axe_hist,'Position'));
set(axe_cumhist,'UserData',get(axe_cumhist,'Position'));
Glb_Infos = struct(...
'typeDATA', typeDATA, ...
'resVal', resVal, ...
'infos_hdls', infos_hdls, ...
'axe_signal', axe_signal, ...
'axe_hist', axe_hist, ...
'axe_cumhist', axe_cumhist ...
);
switch typeDATA
case 'line'
set(axe_corr,'UserData',get(axe_corr,'Position'));
set(axe_spec,'UserData',get(axe_spec,'Position'));
Glb_Infos.axe_corr = axe_corr;
Glb_Infos.axe_spec = axe_spec;
case {'image','regress'}
end
wfigmngr('storeValue',fig,'group_hdls',group_hdls);
wfigmngr('storeValue',fig,'Glb_Infos',Glb_Infos);
wfigmngr('storeValue',fig,'Old_nb_bins',nb_bins);
% Dynvtool Attachment.
%---------------------
switch typeDATA
case 'line' , axe_ind = [axe_signal axe_corr axe_spec];
case {'image','regress'} , axe_ind = axe_signal;
end
axe_cmd = [];
axe_act = [];
dynvtool('attach',fig,axe_ind,axe_cmd,axe_act,[1 0],'','','');
dynvtool('put',fig);
% End waiting.
%-------------
set(fig,'Visible','On');
wwaiting('off',fig);
case 'select'
%**************************************************************%
%** OPTION = 'select' - GROUP AXES SELECTION **%
%**************************************************************%
sel_chk_btn = gcbo;
% Get stored structure.
%----------------------
group_hdls = wfigmngr('getValue',fig,'group_hdls');
chk_hist_chist = group_hdls.chk_hist_chist;
txt_bin = group_hdls.txt_bin;
edi_bin = group_hdls.edi_bin;
chk_corr_spec = group_hdls.chk_corr_spec;
chk_info_stat = group_hdls.chk_info_stat;
Glb_Infos = wfigmngr('getValue',fig,'Glb_Infos');
infos_hdls = Glb_Infos.infos_hdls;
axe_signal = Glb_Infos.axe_signal;
axe_hist = Glb_Infos.axe_hist;
axe_cumhist = Glb_Infos.axe_cumhist;
axe_corr = Glb_Infos.axe_corr;
axe_spec = Glb_Infos.axe_spec;
fra_sta = infos_hdls(1);
% Get the current selection.
%---------------------------
hist_chist = (get(chk_hist_chist,'UserData')~=0);
corr_spec = (get(chk_corr_spec,'UserData')~=0);
info_stat = (get(chk_info_stat,'UserData')~=0);
% Get the axes original positions.
%---------------------------------
pos_axe_signal = get(axe_signal,'UserData');
pos_axe_hist = get(axe_hist,'UserData');
pos_axe_cumhist = get(axe_cumhist,'UserData');
pos_axe_corr = get(axe_corr,'UserData');
pos_axe_spec = get(axe_spec,'UserData');
pos_fra_sta = get(fra_sta,'Position');
% Graphical parameters.
%----------------------
yspace = 2*pos_fra_sta(2);
heigh_max = pos_axe_signal(2)+pos_axe_signal(4)+yspace;
% Redraw depending on the current selection.
%-------------------------------------------
switch sel_chk_btn
case chk_hist_chist
if hist_chist
set(chk_hist_chist,'Value',0,'UserData',0);
if info_stat
y_rem = heigh_max-(pos_fra_sta(2)+ ...
pos_fra_sta(4))-7*yspace/4;
y_low = pos_fra_sta(2)+pos_fra_sta(4)+yspace/3;
if corr_spec
axe_height = y_rem/2;
pos_axe_corr(2) = y_low;
pos_axe_corr(4) = axe_height;
pos_axe_spec(2) = y_low;
pos_axe_spec(4) = axe_height;
y_low = y_low+axe_height+yspace/3;
pos_axe_signal(2) = y_low;
pos_axe_signal(4) = axe_height;
else
axe_height = y_rem;
pos_axe_signal(2) = y_low+yspace/6;
pos_axe_signal(4) = axe_height;
end
else
y_rem = heigh_max-9*yspace/4;
y_low = pos_fra_sta(2)+yspace/3;
if corr_spec
axe_height = y_rem/2;
pos_axe_corr(2) = y_low;
pos_axe_corr(4) = axe_height;
pos_axe_spec(2) = y_low;
pos_axe_spec(4) = axe_height;
y_low = y_low+axe_height+yspace/3;
pos_axe_signal(2) = y_low;
pos_axe_signal(4) = axe_height;
else
axe_height = y_rem;
pos_axe_signal(2) = y_low+yspace/6;
pos_axe_signal(4) = axe_height;
end
end
else
set(chk_hist_chist,'Value',1,'UserData',1);
if info_stat
y_rem = heigh_max-(pos_fra_sta(2)+ ...
pos_fra_sta(4))-8*yspace/4;
y_low = pos_fra_sta(2)+pos_fra_sta(4)+yspace/3;
if corr_spec
axe_height = y_rem/3;
pos_axe_corr(2) = y_low;
pos_axe_corr(4) = axe_height;
pos_axe_spec(2) = y_low;
pos_axe_spec(4) = axe_height;
y_low = y_low+axe_height+yspace/3;
pos_axe_hist(2) = y_low;
pos_axe_hist(4) = axe_height;
pos_axe_cumhist(2) = y_low;
pos_axe_cumhist(4) = axe_height;
y_low = y_low+axe_height+yspace/3;
pos_axe_signal(2) = y_low;
pos_axe_signal(4) = axe_height;
else
axe_height = y_rem/2;
pos_axe_hist(2) = y_low;
pos_axe_hist(4) = axe_height;
pos_axe_cumhist(2) = y_low;
pos_axe_cumhist(4) = axe_height;
y_low = y_low+axe_height+yspace/3;
pos_axe_signal(2) = y_low;
pos_axe_signal(4) = axe_height;
end
else
y_rem = heigh_max-10*yspace/4;
y_low = pos_fra_sta(2)+yspace/3;
if corr_spec
axe_height = y_rem/3;
pos_axe_corr(2) = y_low;
pos_axe_corr(4) = axe_height;
pos_axe_spec(2) = y_low;
pos_axe_spec(4) = axe_height;
y_low = y_low+axe_height+yspace/3;
pos_axe_hist(2) = y_low;
pos_axe_hist(4) = axe_height;
pos_axe_cumhist(2) = y_low;
pos_axe_cumhist(4) = axe_height;
y_low = y_low+axe_height+yspace/3;
pos_axe_signal(2) = y_low;
pos_axe_signal(4) = axe_height;
else
axe_height = y_rem/2;
pos_axe_hist(2) = y_low;
pos_axe_hist(4) = axe_height;
pos_axe_cumhist(2) = y_low;
pos_axe_cumhist(4) = axe_height;
y_low = y_low+axe_height+yspace/3;
pos_axe_signal(2) = y_low;
pos_axe_signal(4) = axe_height;
end
end
end
case chk_corr_spec
if corr_spec
set(chk_corr_spec,'Value',0,'UserData',0);
if info_stat
y_rem = heigh_max-(pos_fra_sta(2)+ ...
pos_fra_sta(4))-7*yspace/4;
y_low = pos_fra_sta(2)+pos_fra_sta(4)+yspace/3;
if hist_chist
axe_height = y_rem/2;
pos_axe_hist(2) = y_low;
pos_axe_hist(4) = axe_height;
pos_axe_cumhist(2) = y_low;
pos_axe_cumhist(4) = axe_height;
y_low = y_low+axe_height+yspace/3;
pos_axe_signal(2) = y_low;
pos_axe_signal(4) = axe_height;
else
axe_height = y_rem;
pos_axe_signal(2) = y_low+yspace/6;
pos_axe_signal(4) = axe_height;
end
else
y_rem = heigh_max-9*yspace/4;
y_low = pos_fra_sta(2)+yspace/3;
if hist_chist
axe_height = y_rem/2;
pos_axe_hist(2) = y_low;
pos_axe_hist(4) = axe_height;
pos_axe_cumhist(2) = y_low;
pos_axe_cumhist(4) = axe_height;
y_low = y_low+axe_height+yspace/3;
pos_axe_signal(2) = y_low;
pos_axe_signal(4) = axe_height;
else
axe_height = y_rem;
pos_axe_signal(2) = y_low+yspace/6;
pos_axe_signal(4) = axe_height;
end
end
else
set(chk_corr_spec,'Value',1,'UserData',1);
if info_stat
y_rem = heigh_max-(pos_fra_sta(2)+ ...
pos_fra_sta(4))-8*yspace/4;
y_low = pos_fra_sta(2)+pos_fra_sta(4)+yspace/3;
if hist_chist
axe_height = y_rem/3;
pos_axe_corr(2) = y_low;
pos_axe_corr(4) = axe_height;
pos_axe_spec(2) = y_low;
pos_axe_spec(4) = axe_height;
y_low = y_low+axe_height+yspace/3;
pos_axe_hist(2) = y_low;
pos_axe_hist(4) = axe_height;
pos_axe_cumhist(2) = y_low;
pos_axe_cumhist(4) = axe_height;
y_low = y_low+axe_height+yspace/3;
pos_axe_signal(2) = y_low;
pos_axe_signal(4) = axe_height;
else
axe_height = y_rem/2;
pos_axe_corr(2) = y_low;
pos_axe_corr(4) = axe_height;
pos_axe_spec(2) = y_low;
pos_axe_spec(4) = axe_height;
y_low = y_low+axe_height+yspace/3;
pos_axe_signal(2) = y_low;
pos_axe_signal(4) = axe_height;
end
else
y_rem = heigh_max-10*yspace/4;
y_low = pos_fra_sta(2)+yspace/3;
if hist_chist
axe_height = y_rem/3;
pos_axe_corr(2) = y_low;
pos_axe_corr(4) = axe_height;
pos_axe_spec(2) = y_low;
pos_axe_spec(4) = axe_height;
y_low = y_low+axe_height+yspace/3;
pos_axe_hist(2) = y_low;
pos_axe_hist(4) = axe_height;
pos_axe_cumhist(2) = y_low;
pos_axe_cumhist(4) = axe_height;
y_low = y_low+axe_height+yspace/3;
pos_axe_signal(2) = y_low;
pos_axe_signal(4) = axe_height;
else
axe_height = y_rem/2;
pos_axe_corr(2) = y_low;
pos_axe_corr(4) = axe_height;
pos_axe_spec(2) = y_low;
pos_axe_spec(4) = axe_height;
y_low = y_low+axe_height+yspace/3;
pos_axe_signal(2) = y_low;
pos_axe_signal(4) = axe_height;
end
end
end
case chk_info_stat
if info_stat
set(chk_info_stat,'Value',0,'UserData',0);
y_low = pos_fra_sta(2)+yspace/3;
if corr_spec
y_rem = heigh_max-10*yspace/4;
if hist_chist
axe_height = y_rem/3;
pos_axe_corr(2) = y_low;
pos_axe_corr(4) = axe_height;
pos_axe_spec(2) = y_low;
pos_axe_spec(4) = axe_height;
y_low = y_low+axe_height+yspace/3;
pos_axe_hist(2) = y_low;
pos_axe_hist(4) = axe_height;
pos_axe_cumhist(2) = y_low;
pos_axe_cumhist(4) = axe_height;
y_low = y_low+axe_height+yspace/3;
pos_axe_signal(2) = y_low;
pos_axe_signal(4) = axe_height;
else
axe_height = y_rem/2;
pos_axe_corr(2) = y_low;
pos_axe_corr(4) = axe_height;
pos_axe_spec(2) = y_low;
pos_axe_spec(4) = axe_height;
y_low = y_low+axe_height+yspace/3;
pos_axe_signal(2) = y_low+yspace/6;
pos_axe_signal(4) = axe_height;
end
else
y_rem = heigh_max-9*yspace/4;
if hist_chist
axe_height = y_rem/2;
pos_axe_hist(2) = y_low;
pos_axe_hist(4) = axe_height;
pos_axe_cumhist(2) = y_low;
pos_axe_cumhist(4) = axe_height;
y_low = y_low+axe_height+yspace/3;
pos_axe_signal(2) = y_low;
pos_axe_signal(4) = axe_height;
else
axe_height = y_rem;
pos_axe_signal(2) = y_low+yspace/6;
pos_axe_signal(4) = axe_height;
end
end
else
set(chk_info_stat,'Value',1,'UserData',1);
y_low = pos_fra_sta(2)+pos_fra_sta(4)+yspace/3;
if corr_spec
if hist_chist
y_rem = heigh_max-(pos_fra_sta(2)+ ...
pos_fra_sta(4))-8*yspace/4;
axe_height = y_rem/3;
pos_axe_corr(2) = y_low;
pos_axe_corr(4) = axe_height;
pos_axe_spec(2) = y_low;
pos_axe_spec(4) = axe_height;
y_low = y_low+axe_height+yspace/3;
pos_axe_hist(2) = y_low;
pos_axe_hist(4) = axe_height;
pos_axe_cumhist(2) = y_low;
pos_axe_cumhist(4) = axe_height;
y_low = y_low+axe_height+yspace/3;
pos_axe_signal(2) = y_low;
pos_axe_signal(4) = axe_height;
else
y_rem = heigh_max-(pos_fra_sta(2)+ ...
pos_fra_sta(4))-7*yspace/4;
axe_height = y_rem/2;
pos_axe_corr(2) = y_low;
pos_axe_corr(4) = axe_height;
pos_axe_spec(2) = y_low;
pos_axe_spec(4) = axe_height;
y_low = y_low+axe_height+yspace/3;
pos_axe_signal(2) = y_low;
pos_axe_signal(4) = axe_height;
end
else
y_rem = heigh_max-(pos_fra_sta(2)+ ...
pos_fra_sta(4))-7*yspace/4;
if hist_chist
axe_height = y_rem/2;
pos_axe_hist(2) = y_low;
pos_axe_hist(4) = axe_height;
pos_axe_cumhist(2) = y_low;
pos_axe_cumhist(4) = axe_height;
y_low = y_low+axe_height+yspace/3;
pos_axe_signal(2) = y_low;
pos_axe_signal(4) = axe_height;
else
y_rem = heigh_max-(pos_fra_sta(2)+ ...
pos_fra_sta(4))-6*yspace/4;
axe_height = y_rem;
pos_axe_signal(2) = y_low;
pos_axe_signal(4) = axe_height;
end
end
end
end
% Set the current positions.
%---------------------------
set(axe_signal,'Position',pos_axe_signal);
set(axe_hist,'Position',pos_axe_hist);
set(axe_cumhist,'Position',pos_axe_cumhist);
set(axe_corr,'Position',pos_axe_corr);
set(axe_spec,'Position',pos_axe_spec);
% Set enability and visibility for each controlable axes or frames.
%------------------------------------------------------------------
set(axe_signal,'Visible','on');
set(get(axe_signal,'Children'),'Visible','on');
if get(chk_hist_chist,'Value')
set(txt_bin,'Visible','on');
set(edi_bin,'Visible','on');
set(axe_hist,'Visible','on');
set(get(axe_hist,'Children'),'Visible','on');
set(axe_cumhist,'Visible','on');
set(get(axe_cumhist,'Children'),'Visible','on');
else
set(txt_bin,'Visible','off');
set(edi_bin,'Visible','off');
set(axe_hist,'Visible','off');
set(get(axe_hist,'Children'),'Visible','off');
set(axe_cumhist,'Visible','off');
set(get(axe_cumhist,'Children'),'Visible','off');
end
if get(chk_corr_spec,'Value')
set(axe_corr,'Visible','on');
set(get(axe_corr,'Children'),'Visible','on');
set(axe_spec,'Visible','on');
set(get(axe_spec,'Children'),'Visible','on');
else
set(axe_corr,'Visible','off');
set(get(axe_corr,'Children'),'Visible','off');
set(axe_spec,'Visible','off');
set(get(axe_spec,'Children'),'Visible','off');
end
if get(chk_info_stat,'Value')
set(infos_hdls,'Visible','on');
else
set(infos_hdls,'Visible','off');
end
case 'update_bins'
%**************************************************************%
%** OPTION = 'update_bins' - UPDATE HISTOGRAMS WITH NEW BINS **%
%**************************************************************%
edi_bin = gcbo;
% Get stored structure.
%----------------------
Old_nb_bins = wfigmngr('getValue',fig,'Old_nb_bins');
Glb_Infos = wfigmngr('getValue',fig,'Glb_Infos');
resVal = Glb_Infos.resVal;
axe_hist = Glb_Infos.axe_hist;
axe_cumhist = Glb_Infos.axe_cumhist;
% Return if no current display.
%------------------------------
vis = get(axe_hist,'Visible');
if isequal(vis(1:2),'of') , return, end
% Check the bins number.
%-----------------------
if ~isempty(Old_nb_bins)
default_bins = Old_nb_bins;
end
nb_bins = wstr2num(get(edi_bin,'String'));
if isempty(nb_bins) || (nb_bins<2)
nb_bins = default_bins;
set(edi_bin,'String',sprintf('%.0f',default_bins))
end
if default_bins==nb_bins , return; end
% Waiting message.
%-----------------
wwaiting('msg',fig,getWavMSG('Wavelet:commongui:WaitCompute'));
% Save Userdata.
%---------------
pos_axe_hist = get(axe_hist,'UserData');
pos_axe_cumhist = get(axe_cumhist,'UserData');
% Updating histograms.
%---------------------
if ~isempty(resVal)
Old_nb_bins = nb_bins;
wfigmngr('storeValue',fig,'Old_nb_bins',Old_nb_bins);
his = wgethist(resVal(:),nb_bins);
his(2,:) = his(2,:)/length(resVal(:));
wplothis(axe_hist,his,curr_color);
wtitle(getWavMSG('Wavelet:commongui:Str_HIST'),'Parent',axe_hist);
for i=6:4:length(his(2,:));
his(2,i) = his(2,i)+his(2,i-4);
his(2,i+1) = his(2,i);
end
wplothis(axe_cumhist,[his(1,:);his(2,:)],curr_color);
wtitle(getWavMSG('Wavelet:commongui:Str_CumHist'),'Parent',axe_cumhist);
end
% Restore Userdata.
%------------------
set(axe_hist,'UserData',pos_axe_hist);
set(axe_cumhist,'UserData',pos_axe_cumhist);
% End waiting.
%-------------
wwaiting('off',fig);
case 'save'
%---------------------------------%
% Option: SAVE : Saving residuals %
%---------------------------------%
% Get arguments.
%---------------
fig = varargin{1};
% Testing file.
%--------------
[filename,pathname,ok] = utguidiv('test_save',fig, ...
'*.mat',getWavMSG('Wavelet:commongui:Lab_SaveRes'));
if ~ok, return; end
% Begin waiting.
%--------------
wwaiting('msg',fig,getWavMSG('Wavelet:commongui:WaitSave'));
% Get the residuals to save.
%---------------------------
Glb_Infos = wfigmngr('getValue',fig,'Glb_Infos');
resVal = Glb_Infos.resVal;
typeDATA = Glb_Infos.typeDATA;
% Saving transformed Signal.
%---------------------------
[name,ext] = strtok(filename,'.');
if isempty(ext) || isequal(ext,'.')
ext = '.mat'; filename = [name ext];
end
switch typeDATA
case 'image'
X = resVal; %#ok<NASGU>
map = get(fig,'Colormap'); %#ok<NASGU>
saveStr = {'X','map'};
case 'regress'
ls = wmemtool('rmb',fig,n_membloc1,ind_loc_struct);
xdata = get(ls.handleTHR,'XData');
ydata = interp1(resVal,xdata); %#ok<NASGU>
saveStr = {'xdata','ydata'};
otherwise
try
eval([name ' = resVal;']);
saveStr = {name};
catch
x = resVal; %#ok<NASGU>
saveStr = {'x'};
end
end
try
save([pathname filename],saveStr{:});
catch
errargt(mfilename,getWavMSG('Wavelet:commongui:SaveFail'),'msg');
end
% End waiting.
%-------------
wwaiting('off',fig);
case 'close'
%**************************************************************%
%** OPTION = 'close' - Close the current figure **%
%**************************************************************%
fig = varargin{1};
ls = wmemtool('rmb',fig,n_membloc1,ind_loc_struct);
toggle = ls.toggle;
if ishandle(toggle) , set(toggle,'Value',0); end
otherwise
errargt(mfilename,getWavMSG('Wavelet:moreMSGRF:Unknown_Opt'),'msg');
error(message('Wavelet:FunctionArgVal:Invalid_ArgVal'));
end
%-------------------------------------------------------------------------%
function [type,prop] = getDataType(handle)
type = get(handle,'Type');
switch type
case 'line' , prop = 'YData';
case 'image' , prop = 'CData';
end
%-------------------------------------------------------------------------%
%-------------------------------------------------------------------------%
function [sp,f] = wspecfft(signal)
%WSPECFFT FFT spectrum of a signal.
%
% f is the frequency
% sp is the energy, the square of the FFT transform
% The input signal is empty.
%---------------------------
if isempty(signal)
sp = [];f =[];return
end
% Compute the spectrum.
%----------------------
n = length(signal);
XTF = fft(fftshift(signal));
m = ceil(n/2) + 1;
% Compute the output values.
%---------------------------
f = linspace(0,0.5,m);
sp = (abs(XTF(1:m))).^2;
%-------------------------------------------------------------------------%
%-------------------------------------------------------------------------%
function [c,lags] = wautocor(a,maxlag)
%WAUTOCOR Auto-correlation function estimates.
% [C,LAGS] = WAUTOCOR(A,MAXLAG) computes the
% autocorrelation function c of a one dimensional
% signal a, for lags = [-maxlag:maxlag].
% The autocorrelation c(maxlag+1) = 1.
% If nargin==1, by default, maxlag = length(a)-1.
if nargin == 1, maxlag = size(a,2)-1;end
lags = -maxlag:maxlag;
if isempty(a) , c = []; return; end
epsi = sqrt(eps);
a = a(:);
a = a - mean(a);
nr = length(a);
if std(a)>epsi
% Test of the variance.
%----------------------
mr = 2 * maxlag + 1;
nfft = 2^nextpow2(mr);
nsects = ceil(2*nr/nfft);
if nsects>4 && nfft<64
nfft = min(4096,max(64,2^nextpow2(nr/4)));
end
c = zeros(nfft,1);
minus1 = (-1).^(0:nfft-1)';
af_old = zeros(nfft,1);
n1 = 1;
nfft2 = nfft/2;
while (n1<nr)
n2 = min( n1+nfft2-1, nr );
af = fft(a(n1:n2,:), nfft);
c = c + af.* conj( af + af_old);
n1 = n1 + nfft2;
af_old = minus1.*af;
end
if n1==nr
af = ones(nfft,1)*a(nr,:);
c = c + af.* conj( af + af_old );
end
mxlp1 = maxlag+1;
c = real(ifft(c));
c = [ c(mxlp1:-1:2,:); c(1:mxlp1,1) ];
% Compute the autocorrelation function.
%--------------------------------------
cdiv = c(mxlp1,1);
c = c / cdiv;
else
% If the variance is too small.
%-------------------------------
c = ones(size(lags));
end
%-------------------------------------------------------------------------%