gusucode.com > SOM算法,matlab程序,是经典的教科书的例子程序 > somtoolbox/preprocess.m
function preprocess(sData,arg2)
%PREPROCESS A GUI for data preprocessing.
%
% preprocess(sData)
%
% preprocess(sData)
%
% Launches a preprocessing GUI. The optional input argument can be
% either a data struct or a struct array of such. However, primarily
% the processed data sets are loaded to the application using the
% tools in the GUI. Also, the only way to get the preprocessed data
% sets back into the workspace is to use the tools in the GUI (press
% the button DATA SET MANAGEMENT).
%
% For a more throughout description, see online documentation.
% See also SOM_GUI.
%%%%%%%%%%%% DETAILED DESCRIPTION %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
% IN FILES: preprocess.html,preproc.jpg,sDman.jpg,clip.jpg,delay.jpg,window.jpg,selVect.jpg
% Contributed to SOM Toolbox vs2, February 2nd, 2000 by Juha Parhankangas
% Copyright (c) by Juha Parhankangas and the SOM Toolbox team
% http://www.cis.hut.fi/projects/somtoolbox/
% Juha Parhankangas 050100
global no_of_sc % every Nth component in 'relative values' is drawn stronger.
no_of_sc=5;
if nargin < 1 | nargin > 2
error('Invalid number of input arguments');
return;
end
if nargin == 1, arg2=[]; end
if ~isstr(sData) %%% Preprocess is started...
data.LOG{1}='% Starting the ''Preprocess'' -window...';
data.LOG{2}=cat(2,'preprocess(',...
sprintf('%s);',inputname(1)));
pre_h=findobj(get(0,'Children'),'Tag','Preprocess');
if ~isempty(pre_h)
figure(pre_h);
msgbox('''Preprocess''-figure already exists.');
return;
end
h0 = figure('Color',[0.8 0.8 0.8], ...
'PaperPosition',[18 180 576 432], ...
'PaperUnits','points', ...
'Position',[595 216 600 775], ...
'Tag','Preprocess');
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[0.752941176470588 0.752941176470588 0.752941176470588], ...
'FontWeight','demi', ...
'HorizontalAlignment','left', ...
'ListboxTop',0, ...
'Position',[0.015 0.06064516129032258 0.9550000000000001 0.1458064516129032], ...
'Style','text', ...
'Tag','StaticText1');
data.results_h = h1;
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[0.752941176470588 0.752941176470588 0.752941176470588], ...
'Callback','preprocess close', ...
'FontWeight','demi', ...
'ListboxTop',0, ...
'Position',[0.8067 0.0142 0.1667 0.0348],...
'String','CLOSE', ...
'Tag','Pushbutton1');
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[0.8 0.8 0.8], ...
'FontWeight','demi', ...
'HorizontalAlignment','left', ...
'ListboxTop',0, ...
'Position',[0.01833333333333333 0.2141935483870968 0.07000000000000001 0.01806451612903226], ...
'String','LOG', ...
'Style','text', ...
'Tag','StaticText2');
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[0.752941176470588 0.752941176470588 0.752941176470588], ...
'Callback','preprocess sel_comp',...
'FontWeight','demi', ...
'ListboxTop',0, ...
'Position',[0.7983333333333333 0.2090322580645161 0.1666666666666667 0.03483870967741935], ...
'String',' ', ...
'Style','popupmenu', ...
'Tag','sel_comp_h', ...
'Value',1);
data.sel_comp_h=h1;
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[1 1 1], ...
'ListboxTop',0, ...
'Position',[0.0183 0.2568 0.2133 0.1290], ...
'Style','text', ...
'Tag','StaticText3');
data.sel_cdata_h=h1;
h1 = axes('Parent',h0, ...
'CameraUpVector',[0 1 0], ...
'CameraUpVectorMode','manual', ...
'Color',[1 1 1], ...
'Position',[0.2583 0.2568 0.2133 0.1290], ...
'Tag','Axes1', ...
'XColor',[0 0 0], ...
'XTickLabel',['0 ';'0.5';'1 '], ...
'XTickLabelMode','manual', ...
'XTickMode','manual', ...
'YColor',[0 0 0], ...
'YTickMode','manual', ...
'ZColor',[0 0 0]);
data.sel_chist_h=h1;
h2 = text('Parent',h1, ...
'Color',[0 0 0], ...
'HandleVisibility','off', ...
'HorizontalAlignment','center', ...
'Position',[0.4960629921259843 -0.08080808080808044 9.160254037844386], ...
'Tag','Axes1Text4', ...
'VerticalAlignment','cap');
set(get(h2,'Parent'),'XLabel',h2);
h2 = text('Parent',h1, ...
'Color',[0 0 0], ...
'HandleVisibility','off', ...
'HorizontalAlignment','center', ...
'Position',[-0.0551181102362206 0.4848484848484853 9.160254037844386], ...
'Rotation',90, ...
'Tag','Axes1Text3', ...
'VerticalAlignment','baseline');
set(get(h2,'Parent'),'YLabel',h2);
h2 = text('Parent',h1, ...
'Color',[0 0 0], ...
'HandleVisibility','off', ...
'HorizontalAlignment','right', ...
'Position',[-1.2283 5.7980 9.1603], ...
'Tag','Axes1Text2', ...
'Visible','off');
set(get(h2,'Parent'),'ZLabel',h2);
h2 = text('Parent',h1, ...
'Color',[0 0 0], ...
'HandleVisibility','off', ...
'HorizontalAlignment','center', ...
'Position',[0.4960629921259843 1.070707070707071 9.160254037844386], ...
'Tag','Axes1Text1', ...
'VerticalAlignment','bottom');
set(get(h2,'Parent'),'Title',h2);
h1 = axes('Parent',h0, ...
'CameraUpVector',[0 1 0], ...
'CameraUpVectorMode','manual', ...
'Color',[0.7529 0.7529 0.7529], ...
'Position',[0.4950000000000001 0.2567741935483871 0.4766666666666667 0.1290322580645161], ...
'Tag','Axes2', ...
'XColor',[0 0 0], ...
'XTickMode','manual', ...
'YColor',[0 0 0], ...
'YTick',[0 0.5 1], ...
'YTickMode','manual', ...
'ZColor',[0 0 0]);
data.vector_h=h1;
h2 = text('Parent',h1, ...
'Color',[0 0 0], ...
'HandleVisibility','off', ...
'HorizontalAlignment','center', ...
'Position',[0.4982456140350879 -0.08080808080808044 9.160254037844386], ...
'Tag','Axes2Text4', ...
'VerticalAlignment','cap');
set(get(h2,'Parent'),'XLabel',h2);
h2 = text('Parent',h1, ...
'Color',[0 0 0], ...
'HandleVisibility','off', ...
'HorizontalAlignment','center', ...
'Position',[-0.1018 0.4848 9.1603], ...
'Rotation',90, ...
'Tag','Axes2Text3', ...
'VerticalAlignment','baseline');
set(get(h2,'Parent'),'YLabel',h2);
h2 = text('Parent',h1, ...
'Color',[0 0 0], ...
'HandleVisibility','off', ...
'HorizontalAlignment','right', ...
'Position',[-1.045614035087719 5.797979797979799 9.160254037844386], ...
'Tag','Axes2Text2', ...
'Visible','off');
set(get(h2,'Parent'),'ZLabel',h2);
h2 = text('Parent',h1, ...
'Color',[0 0 0], ...
'HandleVisibility','off', ...
'HorizontalAlignment','center', ...
'Position',[0.4982456140350879 1.070707070707071 9.160254037844386], ...
'Tag','Axes2Text1', ...
'VerticalAlignment','bottom');
set(get(h2,'Parent'),'Title',h2);
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[0.8 0.8 0.8], ...
'FontWeight','demi', ...
'HorizontalAlignment','left', ...
'ListboxTop',0, ...
'Position',[0.01833333333333333 0.3922580645161291 0.17 0.01806451612903226], ...
'String','STATISTICS', ...
'Style','text', ...
'Tag','StaticText4');
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[0.8 0.8 0.8], ...
'FontWeight','demi', ...
'HorizontalAlignment','left', ...
'ListboxTop',0, ...
'Position',[0.2583333333333334 0.3922580645161291 0.1633333333333333 0.01806451612903226], ...
'String','HISTOGRAM', ...
'Style','text', ...
'Tag','StaticText5');
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[0.8 0.8 0.8], ...
'FontWeight','demi',...
'FontSize',6,...
'HorizontalAlignment','left',...
'String',{'LEFT: NEW SELECTION';'RIGHT: ADD TO SELECTION'}, ...
'ListboxTop',0, ...
'Position',[0.5016666666666667 0.38 0.235 0.03741935483870968], ...
'Style','text', ...
'Tag','StaticText6', ...
'UserData','[ ]');
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[0.752941176470588 0.752941176470588 0.752941176470588], ...
'Callback','preprocess selall', ...
'FontWeight','demi', ...
'ListboxTop',0, ...
'Position',[0.8066666666666668 0.3922580645161291 0.1666666666666667 0.03483870967741935], ...
'String','SELECT ALL', ...
'Tag','Pushbutton2', ...
'UserData','[ ]');
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[0.7529 0.7529 0.7529], ...
'Position',[0.01833333333333333 0.4503225806451613 0.23 0.3225806451612903], ...
'String',' ', ...
'Style','listbox', ...
'Tag','Listbox1', ...
'Value',1);
data.comp_names_h=h1;
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[1 1 1], ...
'Position',[0.4950000000000001 0.4503225806451613 0.2333333333333333 0.3225806451612903], ...
'String',' ', ...
'Style','listbox', ...
'Tag','Listbox2', ...
'Value',1);
data.vect_mean_h = h1;
h1 = axes('Parent',h0, ...
'CameraUpVector',[0 1 0], ...
'CameraUpVectorMode','manual', ...
'Color',[1 1 1], ...
'Position',[0.7383333333333334 0.4503225806451613 0.2333333333333333 0.3225806451612903], ...
'Tag','Axes3', ...
'XColor',[0 0 0], ...
'XTickMode','manual', ...
'YColor',[0 0 0], ...
'YTickMode','manual', ...
'ZColor',[0 0 0]);
data.sel_cplot_h = h1;
h2 = text('Parent',h1, ...
'Color',[0 0 0], ...
'HandleVisibility','off', ...
'HorizontalAlignment','center', ...
'Position',[0.4964028776978418 -0.03212851405622486 9.160254037844386], ...
'Tag','Axes3Text4', ...
'VerticalAlignment','cap');
set(get(h2,'Parent'),'XLabel',h2);
h2 = text('Parent',h1, ...
'Color',[0 0 0], ...
'HandleVisibility','off', ...
'HorizontalAlignment','center', ...
'Position',[-0.05035971223021596 0.493975903614458 9.160254037844386], ...
'Rotation',90, ...
'Tag','Axes3Text3', ...
'VerticalAlignment','baseline');
set(get(h2,'Parent'),'YLabel',h2);
h2 = text('Parent',h1, ...
'Color',[0 0 0], ...
'HandleVisibility','off', ...
'HorizontalAlignment','right', ...
'Position',[-3.1942 1.7028 9.1603], ...
'Tag','Axes3Text2', ...
'Visible','off');
set(get(h2,'Parent'),'ZLabel',h2);
h2 = text('Parent',h1, ...
'Color',[0 0 0], ...
'HandleVisibility','off', ...
'HorizontalAlignment','center', ...
'Position',[0.4964028776978418 1.028112449799197 9.160254037844386], ...
'Tag','Axes3Text1', ...
'VerticalAlignment','bottom');
set(get(h2,'Parent'),'Title',h2);
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[0.752941176470588 0.752941176470588 0.752941176470588], ...
'Callback','preprocess plxy', ...
'FontWeight','demi', ...
'ListboxTop',0, ...
'Position',[0.265 0.4683870967741936 0.125 0.03483870967741935], ...
'String','XY-PLOT', ...
'Tag','Pushbutton3');
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[0.752941176470588 0.752941176470588 0.752941176470588], ...
'Callback','preprocess hist', ...
'FontWeight','demi', ...
'ListboxTop',0, ...
'Position',[0.265 0.5303225806451613 0.125 0.03483870967741935], ...
'String','HISTOGRAM', ...
'Tag','Pushbutton4');
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[0.752941176470588 0.752941176470588 0.752941176470588], ...
'Callback','preprocess bplo', ...
'FontWeight','demi', ...
'ListboxTop',0, ...
'Position',[0.265 0.5922580645161291 0.125 0.03483870967741935], ...
'String','BOX PLOT', ...
'Tag','Pushbutton5');
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[0.752941176470588 0.752941176470588 0.752941176470588], ...
'Callback','preprocess plot', ...
'FontWeight','demi', ...
'ListboxTop',0, ...
'Position',[0.265 0.654195483870968 0.125 0.03483870967741935], ...
'String','PLOT', ...
'Tag','Pushbutton6');
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[1 1 1], ...
'ListboxTop',0, ...
'Position',[0.4088888888888889 0.5333333333333333 0.06 0.03268817204301075], ...
'String','30', ...
'Style','edit', ...
'Tag','EditText1');
data.no_of_bins_h = h1;
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[0.8 0.8 0.8], ...
'FontWeight','demi', ...
'HorizontalAlignment','left', ...
'ListboxTop',0, ...
'Position',[0.01833333333333333 0.775483870967742 0.2016666666666667 0.01806451612903226], ...
'String','COMPONENT LIST', ...
'Style','text', ...
'Tag','StaticText7');
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[0.8 0.8 0.8], ...
'FontWeight','demi', ...
'HorizontalAlignment','left', ...
'ListboxTop',0, ...
'Position',[0.4950000000000001 0.775483870967742 0.1966666666666667 0.01806451612903226], ...
'String','AVERAGE', ...
'Style','text', ...
'Tag','StaticText8');
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[0.8 0.8 0.8], ...
'FontWeight','demi', ...
'HorizontalAlignment','left', ...
'ListboxTop',0, ...
'Position',[0.7383333333333334 0.775483870967742 0.225 0.01806451612903226], ...
'String','RELATIVE VALUES', ...
'Style','text', ...
'Tag','StaticText9');
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[0.8 0.8 0.8], ...
'FontSize',10, ...
'FontWeight','demi', ...
'HorizontalAlignment','left', ...
'ListboxTop',0, ...
'Position',[0.01833333333333333 0.8154838709677419 0.2033333333333333 0.0232258064516129], ...
'String','COMPONENTS', ...
'Style','text', ...
'Tag','StaticText10');
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[0.8 0.8 0.8], ...
'FontSize',10, ...
'FontWeight','demi', ...
'HorizontalAlignment','left', ...
'ListboxTop',0, ...
'Position',[0.4950000000000001 0.8154838709677419 0.2 0.0232258064516129], ...
'String','VECTORS', ...
'Style','text', ...
'Tag','StaticText11');
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[0.752941176470588 0.752941176470588 0.752941176470588], ...
'Callback','preprocess sD_management', ...
'FontSize',5, ...
'FontWeight','demi', ...
'ListboxTop',0, ...
'Position',[0.01833333333333333 0.8503225806451613 0.1666666666666667 0.03483870967741935], ...
'String','DATA SET MANAGEMENT', ...
'Tag','Pushbutton7');
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'Callback','preprocess sel_sD', ...
'ListboxTop',0, ...
'Position',[0.01833333333333333 0.8890322580645161 0.1666666666666667 0.03483870967741935], ...
'String',' ', ...
'Style','popupmenu', ...
'Tag','PopupMenu2', ...
'Value',1);
data.sD_set_h = h1;
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[1 1 1], ...
'ListboxTop',0, ...
'Position',[0.2516666666666667 0.8503225806451613 0.7216666666666667 0.07354838709677419], ...
'Style','text', ...
'Tag','StaticText12');
data.sD_name_h = h1;
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[0.8 0.8 0.8], ...
'FontSize',10, ...
'FontWeight','demi', ...
'HorizontalAlignment','left', ...
'ListboxTop',0, ...
'Position',[0.01833333333333333 0.9341935483870968 0.1616666666666667 0.02064516129032258], ...
'String','DATA SETS', ...
'Style','text', ...
'Tag','StaticText13');
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[0.8 0.8 0.8], ...
'FontSize',10, ...
'FontWeight','demi', ...
'HorizontalAlignment','left', ...
'ListboxTop',0, ...
'Position',[0.2516666666666667 0.9341935483870968 0.2833333333333333 0.02064516129032258], ...
'String','SELECTED DATA SET', ...
'Style','text', ...
'Tag','StaticText14');
if ~isstruct(sData),
sData=som_data_struct(sData);
end
ui_h=uimenu('Label','&Normalization');
uimenu(ui_h,'Label','Scale [0,1]','Callback','preprocess zscale');
uimenu(ui_h,'Label','Scale var=1','Callback','preprocess vscale');
uimenu(ui_h,'Label','HistD','Callback','preprocess histeq');
uimenu(ui_h,'Label','HistC','Callback','preprocess histeq2');
uimenu(ui_h,'Label','Log','Callback','preprocess log');
uimenu(ui_h,'Label','Eval (1-comp)','Callback','preprocess eval1');
ui_h=uimenu('Label','&Components');
uimenu(ui_h,'Label','Move component','Callback','preprocess move');
uimenu(ui_h,'Label','Copy component','Callback','preprocess copy');
uimenu(ui_h,'Label','Add: N binary types','Callback','preprocess oneo');
uimenu(ui_h,'Label','Add: zeros','Callback','preprocess zero');
uimenu(ui_h,'Label','Remove component','Callback','preprocess remove');
uimenu(ui_h,'Label','Remove selected vectors',...
'Callback','preprocess remove_vects');
uimenu(ui_h,'Label','Select all components',...
'Callback','preprocess sel_all_comps');
ui_h=uimenu('Label','&Misc');
ui_h1=uimenu(ui_h,'Label','Calculate');
ui_h2=uimenu(ui_h,'Label','Process');
uimenu(ui_h,'Label','Get LOG-file','Callback','preprocess LOG');
uimenu(ui_h,'Label','Indices of the selected vectors',...
'Callback','preprocess get_inds');
uimenu(ui_h,'Label','Undo','Callback','preprocess undo');
uimenu(ui_h1,'Label','Number of values','Callback','preprocess noof');
uimenu(ui_h1,'Label','Number of selected vectors',...
'Callback','preprocess no_of_sel');
uimenu(ui_h1,'Label','Correlation','Callback','preprocess corr');
uimenu(ui_h2,'Label','Unit length','Callback','preprocess unit');
uimenu(ui_h2,'Label','Eval','Callback','preprocess eval2');
uimenu(ui_h2,'Label','Clipping','Callback','preprocess clipping');
uimenu(ui_h2,'Label','Delay','Callback','preprocess delay');
uimenu(ui_h2,'Label','Windowed','Callback','preprocess window');
uimenu(ui_h2,'Label','Select vectors','Callback','preprocess select');
len=getfield(size(sData(1).data),{1});
data.selected_vects=find(ones(1,len));
data.sD_set=sData;
set(data.vector_h,'ButtonDownFcn','preprocess(''vector_bdf'',''down'')');
set(gcf,'UserData',data);
if ~set_sD_stats;
return;
end
sel_sD;
return; %%% Preprocess-window is ready.
else
arg=sData;
if strcmp(arg,'rename')
rename(arg2);
elseif strcmp(arg,'sel_sD')
if isempty(arg2)
sel_sD;
else
sel_sD(arg2);
end
elseif strcmp(arg,'zscale')
if isempty(arg2)
zero2one_scale;
else
zero2one_scale(arg2);
end
elseif strcmp(arg,'vscale');
if isempty(arg2)
var_scale;
else
var_scale(arg2);
end
elseif strcmp(arg,'histeq2')
if isempty(arg2)
hist_eq2;
else
hist_eq2(arg2);
end
elseif strcmp(arg,'log')
if isempty(arg2)
logarithm;
else
logarithm(arg2);
end
elseif strcmp(arg,'eval1')
if isempty(arg2)
eval1;
else
eval1(arg2);
end
elseif strcmp(arg,'eval2')
if isempty(arg2)
eval2;
else
eval2(arg2);
end
elseif strcmp(arg,'histeq');
if isempty(arg2)
hist_eq;
else
hist_eq(arg2);
end
elseif strcmp(arg,'selall')
if isempty(arg2)
select_all;
else
select_all(arg2);
end
elseif strcmp(arg,'sel_button');
if isempty(arg2)
sel_button;
else
sel_button(arg2);
end
elseif strcmp(arg,'clear_button')
if isempty(arg2)
clear_button;
else
clear_button(arg2)
end
elseif(strcmp(arg,'move'))
if isempty(arg2)
move_component;
else
move_component(arg2);
end
elseif(strcmp(arg,'copy'))
if isempty(arg2)
copy_component;
else
copy_component(arg2);
end
elseif strcmp(arg,'oneo')
if isempty(arg2)
one_of_n;
else
one_of_n(arg2);
end
elseif strcmp(arg,'zero')
if isempty(arg2)
add_zeros;
else
add_zeros(arg2);
end
elseif strcmp(arg,'remove')
if isempty(arg2)
remove_component;
else
remove_component(arg2);
end
elseif strcmp(arg,'remove_vects')
if isempty(arg2)
remove_vects;
else
remove_vects(arg2);
end
elseif strcmp(arg,'noof')
if isempty(arg2)
no_of_values;
else
no_of_values(arg2);
end
elseif strcmp(arg,'corr');
if isempty(arg2)
correlation;
else
correlation(arg2);
end
elseif strcmp(arg,'unit')
if isempty(arg2)
unit_length;
else
unit_length(arg2);
end
elseif strcmp(arg,'clip_data')
clip_data(arg2);
elseif strcmp(arg,'copy_delete')
copy_delete(arg2);
elseif strcmp(arg,'and_or_cb')
and_or_cb(arg2);
elseif strcmp(arg,'all_sel_cb')
all_sel_cb(arg2);
elseif strcmp(arg,'clip_exp_cb')
clip_exp_cb(arg2);
elseif strcmp(arg,'window_cb')
window_cb(arg2);
elseif strcmp(arg,'set_state_vals')
set_state_vals(arg2);
elseif strcmp(arg,'vector_bdf')
vector_bdf(arg2);
elseif strcmp(arg,'sD_management');
if isempty(arg2)
sD_management;
else
sD_management(arg2);
end
elseif strcmp(arg,'clipping')
if isempty(arg2)
clipping;
else
clipping(arg2);
end
elseif strcmp(arg,'delay')
if isempty(arg2)
delay;
else
delay(arg2);
end
elseif strcmp(arg,'window');
if isempty(arg2)
window;
else
window(arg2);
end
elseif strcmp(arg,'select');
if isempty(arg2)
select;
else
select(arg2);
end
elseif strcmp(arg,'import')
if isempty(arg2)
import;
else
import(arg2);
end
elseif strcmp(arg,'export')
if isempty(arg2)
export;
else
export(arg2);
end
elseif strcmp(arg,'undo');
if isempty(arg2)
undo;
else
undo(arg2);
end
elseif strcmp(arg,'delay_data')
if isempty(arg2)
delay_data;
else
delay_data(arg2);
end
elseif strcmp(arg,'eval_windowed')
if isempty(arg2)
eval_windowed;
else
eval_windowed(arg2);
end
elseif strcmp(arg,'get_inds')
if isempty(arg2)
get_selected_inds;
else
get_selected_inds(arg2);
end
elseif strcmp(arg,'no_of_sel')
if isempty(arg2)
no_of_selected;
else
no_of_selected(arg2);
end
elseif strcmp(arg,'sel_comp');
if isempty(arg2)
sel_comp;
else
sel_comp(arg2);
end
elseif strcmp(arg,'sel_all_comps')
if isempty(arg2)
select_all_comps;
else
select_all_comps(arg2);
end
elseif strcmp(arg,'refresh')
set_var_names;
elseif any(strcmp(arg,{'close_c','close_d','close_s','close_w','close_sD'}))
if isempty(arg2)
close_func(arg)
else
close_func(arg,arg2);
end
end
switch arg
case 'sD_stats'
sD_stats;
case 'LOG'
log_file;
otherwise
pro_tools(arg);
end
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function set_compnames(sData,h)
%SET_COMPNAMES
%
% set_compnames(sData,h)
%
% ARGUMENTS
%
% sData (struct) som_data_struct
% h (scalar) handle to a list box object
%
%
% This function sets the component names of sData to the list box
% indicated by 'h'.
%
pre_h=findobj(get(0,'Children'),'Tag','Preprocess');
if isempty(pre_h)
error('Figure ''Preprocess'' does not exist. Closing program...');
close_preprocess;
end
udata=get(pre_h,'UserData');
set(h,'Value',[]);
for i=1:length(sData.comp_names)
tmp=sprintf('#%d: ',i);
names{i,1}=cat(2,tmp, sData.comp_names{i});
end
set(h,'String',names,'Max',2);
set(udata.sel_comp_h,'String',names);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function draw_vectors(vectors,h)
%DRAW_VECTORS
%
% draw_vectors(vectors,h)
%
% ARGUMENTS
%
% vectors (vector) vector of 0's and 1's
% h (scalar) handle to an axis object
%
%
% This function draws an horizontal bar of 'vectors' in the axis
% indicated by 'h'.
%
%
pre_h=findobj(get(0,'Children'),'Tag','Preprocess');
subplot(h);
hold off;
cla;
set(h,'YLim',[0 1]);
set(h,'YTick',[]);
set(h,'XLim',[0 length(vectors)+1]);
hold on;
comp_no=get(getfield(get(pre_h,'UserData'),'sel_comp_h'),'Value');
comp=getfield(get(pre_h,'UserData'),'sData');
comp=comp.data(:,comp_no);
Max = max(comp);
Min = min(comp);
lims=get(gca,'YLim');
lims(1)=Min;
h=abs(0.1*Max);
lims(2)=Max;
if Max - Min <= eps
tmp=Max;
lims(1)=tmp-1;
lims(2)=tmp+1;
end
lims(2)=lims(2)+h;
if ~all(isnan(lims))
set(gca,'YLim',lims);
end
h=(lims(2)-lims(1))/4;
set(gca,'YTickMode','auto');
t=1:length(vectors);
h=plot(t,comp);
set(h,'ButtonDownFcn','preprocess(''vector_bdf'',''down'')');
indices =find(vectors);
vectors(indices)=0.1*(getfield(get(gca,'YLim'),...
{2})-getfield(get(gca,'YLim'),{1}));
plot(indices,vectors(indices)+getfield(get(gca,'YLim'),{1}),...
'ored','MarkerSize',4);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function vect_means(sData,handle,indices)
%VECT_MEANS
%
% vect_means(sData,handle,indices)
%
% ARGUMENTS
%
% sData (struct) som_data_struct
% handle (scalar) handle to the static text box object
% indices (vector) indices of selected vectors
%
%
% This function calculates means of selected vectors' components
% and writes them in the static text box indicated by 'handle'.
%
%
sData= sData.data(indices,:);
for i=1:length(sData(1,:))
names{i}=sprintf('#%d: ',i);
end
for i=1:length(sData(1,:))
tmp=sData(:,i);
tmp=cat(2,names{i},sprintf('%-10.3g',mean(tmp(find(~isnan(tmp))))));
string{i}=tmp;
end
set(handle,'String',string);
set(handle,'HorizontalAlignment','left');
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function vector_bdf(arg)
%VECTOR_BDF A button down function.
%
% vector_bdf(arg)
%
% ARGUMENTS
%
% arg (string) 'down' or 'up', tells the mouse button's state.
%
%
% This function selects vectors in the vector-window and plots maxima,
% minima and means of the selected vectors. It also writes means of the
% selected vectors' components in a static text box and takes care of
% changes of the chosen component's data.
%
% See also VECTOR_MEANS, SEL_COMP
%
%
arg2=arg(6:length(arg));
if ~isempty(arg2)
LOG=1;
else
LOG=0;
end
arg=arg(1:4);
%%% arg's first "word" is 4 letters long and it can be:
%%%
%%% 'key '
%%% 'down'
%%% 'drag'
%%% 'up '
if strcmp(arg,'key ') %string is 'key' + 1 space!!!
if ~LOG
key=get(gcf,'CurrentCharacter');
else
key=arg2
end
if ~strcmp(key,'<') & ~strcmp(key,'>')
return;
end
data=get(gcf,'UserData');
sel=data.selected_vects;
if length(sel) == 1
if strcmp(key,'<') & sel ~= 1
data.selected_vects=sel-1;
set(gcf,'UserData',data);
elseif strcmp(key,'>') & sel ~= length(data.sData.data(:,1))
data.selected_vects = sel + 1;
set(gcf,'UserData',data);
end
else
if strcmp(key,'<') & sel(1) ~= 1
data.selected_vects=cat(2,sel(1)-1,sel);
set(gcf,'UserData',data);
elseif strcmp(key,'>') & sel(length(sel)) ~= length(sel)
data.selected_vects=cat(2,sel,sel(length(sel))+1);
set(gcf,'UserData',data);
end
end
cplot_mimema;
pro_tools('plot_hist');
pro_tools('c_stat');
vects=zeros(1,length(data.sData.data(:,1)));
vects(data.selected_vects)=1;
draw_vectors(vects,data.vector_h);
if ~LOG
data=get(gcf,'UserData');
data.LOG{length(data.LOG)+1}=...
sprintf('preprocess(''vector_bdf'',''key %s'');',key);
%string is 'key'+2spaces+%s
set(gcf,'UserData',data);
end
return;
end
switch arg
case 'down'
set(gcf,'WindowButtonUpFcn','preprocess(''vector_bdf'',''up '')');
set(gcf,'WindowButtonMotionFcn','preprocess(''vector_bdf'',''drag'')');
switch get(gcf,'SelectionType')
case 'normal'
data.lims1=round(getfield(get(gca,'CurrentPoint'),{1,1}));
data.lims2=[];
case 'alt'
tmp=round(getfield(get(gca,'CurrentPoint'),{1,1}));
if isempty(get(gca,'UserData'))
data.lims1=tmp;
data.lims2=[];
else
data.lims1=cat(2,getfield(get(gca,'UserData'),'lims1'),tmp);
data.lims2=getfield(get(gca,'UserData'),'lims2');
end
end
coords=get(gca,'CurrentPoint');
h=line([coords(1),coords(1)],get(gca,'YLim'),'EraseMode','xor');
set(h,'Color','red');
h2=line([coords(1),coords(1)],get(gca,'YLim'),'EraseMode','xor');
set(h2,'Color','red');
data.h=h;
data.h2=h2;
set(gca,'UserData',data);
case 'drag'
coords=get(gca,'CurrentPoint');
lim=get(gca,'XLim');
h2=getfield(get(gca,'UserData'),'h2');
if lim(1) >= coords(1)
set(h2,'XData',[lim(1) lim(1)]);
elseif lim(2) <= coords(2)
set(h2,'XData',[lim(2) lim(2)]);
else
set(h2,'XData',[coords(1) coords(1)]);
end
case 'up ' % string is 'up' + 2 spaces!!!
set(gcf,'WindowButtonUpFcn','');
set(gcf,'WindowButtonMotionFcn','');
if ~LOG
data=get(gca,'UserData');
delete(data.h);
delete(data.h2);
tmp=round(getfield(get(gca,'CurrentPoint'),{1,1}));
data.lims2=cat(2,data.lims2,tmp);
tmp_data=sort(cat(1,data.lims1,data.lims2));
high=getfield(get(gca,'XLim'),{2})-1;
vectors=zeros(1,high);
tmp_data(find(tmp_data<1))=1;
tmp_data(find(tmp_data>high))=high;
for i=1:getfield(size(tmp_data),{2})
vectors(tmp_data(1,i):tmp_data(2,i))=1;
end
selected_vects=find(vectors);
else
pre_h=findobj(get(0,'Children'),'Tag','Preprocess');
len=size(getfield(getfield(get(pre_h,'UserData'),'sData'),'data'));
vectors=zeros(1,len(1));
i=1;
while i <= length(arg2) & (isspace(arg2(i)) | ~isletter(arg2(i)))
i=i+1;
end
arg3=arg2(i:length(arg2));
selected_vects=str2num(arg2(1:i-1));
if ~isempty(arg3) & ~all(isspace(arg3))
selected_vects=unique(cat(2,selected_vects,...
getfield(get(pre_h,'UserData'),'selected_vects')));
end
vectors(selected_vects)=1;
set(pre_h,'CurrentAxes',getfield(get(pre_h,'UserData'),'vector_h'));
set(0,'CurrentFigure',pre_h);
end
draw_vectors(vectors,gca);
sData=getfield(get(gcf,'UserData'),'sData');
h=getfield(get(gcf,'UserData'),'vect_mean_h');
vect_means(sData,h,selected_vects);
if ~LOG
set(gca,'UserData',data);
end
data=get(gcf,'UserData');
data.undo.sData=data.sData;
data.undo.selected=data.selected_vects;
data.selected_vects=selected_vects;
if ~LOG
data.LOG{length(data.LOG)+1}='% Vector selection by using the mouse...';
tmp=sprintf('preprocess(''vector_bdf'',''up %s'');',...
num2str(data.selected_vects));
if length(tmp) > 500
tmp=textwrap({tmp},500);
data.LOG{length(data.LOG)+1}=cat(2,tmp{1},''');');
for i=2:length(tmp)-1
data.LOG{length(data.LOG)+1}=...
cat(2,sprintf('preprocess(''vector_bdf'',''up %s',...
tmp{i}),'add'');');
end
data.LOG{length(data.LOG)+1}=...
cat(2,sprintf('preprocess(''vector_bdf'',''up %s',...
tmp{length(tmp)}(1:length(tmp{length(tmp)})-3)),' add'');');
else
data.LOG{length(data.LOG)+1}=tmp;
end
end
set(gcf,'UserData',data);
cplot_mimema;
sel_comp;
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function sel_button(varargin)
%SEL_BUTTON A Callback function. It performs the operations needed
% when vector components are selected.
%
% See also SEL_COMP
%
if nargin == 1
LOG=1;
pre_h=findobj(get(0,'Children'),'Tag','Preprocess');
string=getfield(get(pre_h,'UserData'),'comp_names_h');
string=getfield(get(string,'String'),{str2num(varargin{1})});
set(0,'CurrentFigure',pre_h);
else
LOG=0;
val=get(getfield(get(gcf,'UserData'),'comp_names_h'),'Value');
end
sel_button_h=getfield(get(gcf,'UserData'),'sel_button_h');
sel_comps_h=getfield(get(gcf,'UserData'),'sel_comps_h');
comp_names_h=getfield(get(gcf,'UserData'),'comp_names_h');
if ~LOG
string=getfield(get(comp_names_h,'String'),{get(comp_names_h,'Value')});
end
tmp_string=get(sel_comps_h,'String');
if iscell(tmp_string)
for i=1:length(string)
if ~any(strcmp(string{i},tmp_string))
tmp_string=cat(1,tmp_string,string(i));
end
end
string=tmp_string;
end
set(sel_comps_h,'String',string);
set(comp_names_h,'Value',[]);
sel_comp;
if ~LOG
data=get(gcf,'UserData');
data.LOG{length(data.LOG)+1}='% Select components';
data.LOG{length(data.LOG)+1}=sprintf('preprocess(''sel_button'',''%s'');',...
num2str(val));
set(gcf,'UserData',data);
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function clear_button(varargin)
%CLEAR_BUTTON Function callback evaluated when a 'Clear'-button is
% pressed. It removes texts from the 'selected components'
% -window and the 'selected component data' -window and
% clears the 'histogram' -axis.
%
%
if nargin==1
LOG=1;
pre_h=findobj(get(0,'Children'),'Tag','Preprocess');
set(0,'CurrentFigure',pre_h);
else
LOG=0;
end
sel_comp_h=getfield(get(gcf,'UserData'),'sel_comp_h');
sel_cdata_h=getfield(get(gcf,'UserData'),'sel_cdata_h');
sel_cplot_h=getfield(get(gcf,'UserData'),'sel_cplot_h');
sel_chist_h=getfield(get(gcf,'UserData'),'sel_chist_h');
vector_h=getfield(get(gcf,'UserData'),'vector_h');
set(sel_comp_h,'Value',1);
set(sel_cdata_h,'String',' ');
subplot(sel_chist_h);
hold off;
cla;
selected=getfield(get(gcf,'UserData'),'selected_vects');
dims=size(getfield(getfield(get(gcf,'UserData'),'sData'),'data'));
vectors=zeros(1,dims(1));
vectors(selected)=1;
subplot(vector_h);
draw_vectors(vectors,vector_h);
if ~LOG
data=get(gcf,'UserData');
data.LOG{length(data.LOG)+1}='% Remove components from the selected list.';
data.LOG{length(data.LOG)+1}='preprocess(''clear_button'',''foo'');';
set(gcf,'UserData',data);
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function sel_comp(varargin)
%SEL_COMP performs the operations needed when vector components are
% chosen. It writes maxima, minima, mean and standard deviation
% of the chosen component to a text box window and draws a
% histogram of the chosen component of selected vectors'
%
%
pre_h=findobj(get(0,'Children'),'Tag','Preprocess');
set(0,'CurrentFigure',pre_h);
sel_comp_h=getfield(get(pre_h,'UserData'),'sel_comp_h');
if nargin == 1
set(sel_comp_h,'Value',str2num(varargin{1}));
elseif ~isempty(gcbo)
no=get(sel_comp_h,'Value');
data=get(gcf,'UserData');
data.LOG{length(data.LOG)+1}='% Select one component';
data.LOG{length(data.LOG)+1}=cat(2,'preprocess(''sel_comp'',''',...
num2str(no),''');');
set(gcf,'UserData',data);
end
pro_tools('c_stat');
pro_tools('plot_hist');
data=get(gcf,'UserData');
sData=data.sData;
vector_h=data.vector_h;
len=length(sData.data(:,1));
vects=zeros(1,len);
vects(data.selected_vects)=1;
draw_vectors(vects,vector_h);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function cplot_mimema
global no_of_sc
sData=getfield(get(gcf,'UserData'),'sData');
sel_cplot_h=getfield(get(gcf,'UserData'),'sel_cplot_h');
selected=getfield(get(gcf,'UserData'),'selected_vects');
set(sel_cplot_h,'YLim',[0 length(sData.data(1,:))+1]);
subplot(sel_cplot_h);
hold off;
cla;
hold on;
for i=1:length(sData.data(1,:))
Max=max(sData.data(:,i));
Min=min(sData.data(:,i));
tmp=sData.data(selected,i);
selMax=max(tmp);
selMin=min(tmp);
Mean=abs(mean(tmp(find(~isnan(tmp)))));
Median=abs(median(tmp(find(~isnan(tmp)))));
if Max ~= Min & ~all(isnan(sData.data(:,i)))
if rem(i,no_of_sc) % no_of_sc is defined in the beginning of this file...
line([abs(selMin-Min)/(Max-Min) (selMax-Min)/(Max-Min)],...
[i i],'Color','black');
plot(abs(Mean-Min)/(Max-Min),i,'oblack');
plot(abs(Median-Min)/(Max-Min),i,'xblack');
else
line([abs(selMin-Min)/(Max-Min) (selMax-Min)/(Max-Min)],...
[i i],'Color','black','LineWidth',2);
plot(abs(Mean-Min)/(Max-Min),i,'oblack','LineWidth',2);
plot(abs(Median-Min)/(Max-Min),i,'xblack','LineWidth',2);
end
else
if rem(i,no_of_sc) % N is defined in the beginning of this file.
plot(mean(get(gca,'XLim')),i,'oblack');
plot(mean(get(gca,'XLim')),i,'xblack');
else
plot(mean(get(gca,'XLim')),i,'oblack','LineWidth',2);
plot(mean(get(gca,'XLim')),i,'xblack','LineWidth',2);
end
end
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function bool=set_sD_stats
%SET_SD_STATS Writes the data set names to popup menu.
%
%
bool=1;
data=get(gcf,'UserData');
for i=1:length(data.sD_set)
% if ~isvalid_var_name({data.sD_set(i).name})
% close_preprocess;
% bool=0;
% return;
% end
string{i}=cat(2,sprintf('#%d: ',i),data.sD_set(i).name);
end
set(data.sD_set_h,'String',string);
data.sData=data.sD_set(get(data.sD_set_h,'Value'));
data.sData.MODIFIED=0;
data.sData.INDEX=1;
set(gcf,'UserData',data);
write_sD_stats;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function write_sD_stats
%WRITE_SD_STATS writes data's name, length and dimension to text box.
%
%
pre_h=findobj(get(0,'Children'),'Tag','Preprocess');
sD_name_h=getfield(get(pre_h,'UserData'),'sD_name_h');
sData=getfield(get(pre_h,'UserData'),'sData');
dims=size(sData.data);
string{1}=cat(2,'Name: ',sData.name);
string{2}=cat(2,'Length: ',sprintf('%d',dims(1)));
string{3}=cat(2,'Dim: ',sprintf('%d',dims(2)));
set(sD_name_h,'String',string);
set(sD_name_h,'HorizontalAlignment','left');
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function sel_sD(varargin)
%SEL_SD sets new data to UserData's 'sData'.
%
%
if nargin==1
LOG=1;
index=str2num(varargin{1});
pre_h=findobj(get(0,'Children'),'Tag','Preprocess');
set(0,'CurrentFigure',pre_h);
else
LOG=0;
end
sD_set_h=getfield(get(gcf,'UserData'),'sD_set_h');
comp_names_h=getfield(get(gcf,'UserData'),'comp_names_h');
vector_h=getfield(get(gcf,'UserData'),'vector_h');
vect_mean_h=getfield(get(gcf,'UserData'),'vect_mean_h');
if ~LOG
index=get(sD_set_h,'Value');
end
data=get(gcf,'UserData');
data.undo = [];
INDEX=data.sData.INDEX;
data.sData=rmfield(data.sData,'MODIFIED');
data.sData=rmfield(data.sData,'INDEX');
tmp=data.sD_set(index);
tmp.MODIFIED=0;
tmp.INDEX=index;
data.sD_set(INDEX)=data.sData;
data.sData=tmp;
len=getfield(size(tmp.data),{1});
data.selected_vects=find(ones(1,len));
if ~LOG
data.LOG{length(data.LOG)+1}='% Select a new data set.';
data.LOG{length(data.LOG)+1}=sprintf('preprocess(''sel_sD'',''%s'');',...
num2str(index));
end
set(gcf,'UserData',data);
write_sD_stats;
set_compnames(tmp,comp_names_h);
draw_vectors(ones(1,len),vector_h);
vect_means(tmp,vect_mean_h,data.selected_vects);
clear_button;
sel_comp;
cplot_mimema;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function indices=get_indices
pre_h=findobj(get(0,'Children'),'Tag','Preprocess');
comp_names_h=getfield(get(pre_h,'UserData'),'comp_names_h');
indices = get(comp_names_h,'Value');
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function sD_management(varargin)
if nargin ~= 1
pre_h=findobj(get(0,'Children'),'Tag','Preprocess');
preh_udata=get(pre_h,'UserData');
preh_udata.LOG{length(preh_udata.LOG)+1}=...
'% Starting the ''Data Set Management'' -window...';
preh_udata.LOG{length(preh_udata.LOG)+1}=...
'preprocess(''sD_management'',''foo'');';
set(pre_h,'UserData',preh_udata);
end
man_h=findobj(get(0,'Children'),'Tag','Management');
if ~isempty(man_h)
figure(man_h);
return;
end
h0 = figure('BackingStore','off', ...
'Color',[0.8 0.8 0.8], ...
'Name','Data Set Management', ...
'PaperPosition',[18 180 576 432], ...
'PaperUnits','points', ...
'Position',[753 523 324 470], ...
'RendererMode','manual', ...
'Tag','Management');
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[1 1 1], ...
'FontWeight','demi', ...
'HorizontalAlignment','left', ...
'ListboxTop',0, ...
'Max',2, ...
'Position',[0.02777777777777778 0.0723404255319149 0.7716049382716049 0.1914893617021277], ...
'String',' ', ...
'Style','edit', ...
'Tag','EditText1');
data.new_c_name_h = h1;
h1 = uicontrol('Parent',h0, ...
'Callback','preprocess rename comp',...
'Units','normalized', ...
'FontSize',6, ...
'FontWeight','demi', ...
'ListboxTop',0, ...
'Position',[0.8240740740740741 0.2106382978723404 0.154320987654321 0.05319148936170213], ...
'String','RENAME', ...
'Tag','Pushbutton1');
h1 = uicontrol('Parent',h0, ...
'Callback','preprocess close_sD',...
'Units','normalized', ...
'FontWeight','demi', ...
'ListboxTop',0, ...
'Position',[0.8240740740740741 0.01914893617021277 0.154320987654321 0.05319148936170213], ...
'String','CLOSE', ...
'Tag','Pushbutton2');
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[0.8 0.8 0.8], ...
'FontWeight','demi', ...
'HorizontalAlignment','left', ...
'ListboxTop',0, ...
'Position',[0.02777777777777778 0.2680851063829787 0.345679012345679 0.02978723404255319], ...
'String','COMPONENTS:', ...
'Style','text', ...
'Tag','StaticText1');
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[1 1 1], ...
'HorizontalAlignment','left', ...
'Position',[0.02777777777777778 0.3170212765957447 0.3549382716049382 0.5319148936170213], ...
'String',' ', ...
'Style','listbox', ...
'Tag','Listbox1', ...
'Value',1);
data.sets_h=h1;
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[1 1 1], ...
'HorizontalAlignment','left', ...
'Position',[0.6234567901234568 0.3170212765957447 0.3549382716049382 0.5319148936170213], ...
'String',' ', ...
'Style','listbox', ...
'Tag','Listbox2', ...
'Value',1);
data.variables_h = h1;
h1 = uicontrol('Parent',h0, ...
'Callback','preprocess export',...
'Units','normalized', ...
'FontWeight','demi', ...
'ListboxTop',0, ...
'Position',[0.4259259259259259 0.551063829787234 0.154320987654321 0.0425531914893617], ...
'String','->', ...
'Tag','Pushbutton4');
h1 = uicontrol('Parent',h0, ...
'Callback','preprocess import',...
'Units','normalized', ...
'FontWeight','demi', ...
'ListboxTop',0, ...
'Position',[0.4259259259259259 0.625531914893617 0.154320987654321 0.0425531914893617], ...
'String','<-', ...
'Tag','Pushbutton3');
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[0.8 0.8 0.8], ...
'FontWeight','demi', ...
'HorizontalAlignment','left', ...
'ListboxTop',0, ...
'Position',[0.02777777777777778 0.8531914893617022 0.2993827160493827 0.02978723404255319], ...
'String','DATA SETS', ...
'Style','text', ...
'Tag','StaticText2');
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[0.8 0.8 0.8], ...
'FontWeight','demi', ...
'HorizontalAlignment','left', ...
'ListboxTop',0, ...
'Position',[0.6234567901234568 0.8531914893617022 0.2561728395061728 0.02978723404255319], ...
'String','WORKSPACE', ...
'Style','text', ...
'Tag','StaticText3');
h1 = uicontrol('Parent',h0, ...
'Callback','preprocess rename set',...
'Units','normalized', ...
'BackgroundColor',[1 1 1], ...
'HorizontalAlignment','left', ...
'ListboxTop',0, ...
'Position',[0.1820987654320987 0.9127659574468086 0.7808641975308641 0.0425531914893617], ...
'Style','edit', ...
'Tag','EditText2');
data.new_name_h = h1;
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[0.8 0.8 0.8], ...
'FontWeight','demi', ...
'HorizontalAlignment','left', ...
'ListboxTop',0, ...
'Position',[0.02777777777777778 0.9127659574468086 0.1388888888888889 0.02978723404255319], ...
'String','NAME:', ...
'Style','text', ...
'Tag','StaticText4');
ui_h=uimenu('Label','&Tools');
uimenu(ui_h,'Label','Copy','Callback','preprocess copy_delete copy');
uimenu(ui_h,'Label','Delete','Callback','preprocess copy_delete delete');
uimenu(ui_h,'Label','Refresh','Callback','preprocess refresh');
set(gcf,'UserData',data);
set_var_names;
sD_names;
sD_stats;
%%% Subfunction: set_var_names %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function set_var_names
variables_h=getfield(get(gcf,'UserData'),'variables_h');
value=get(variables_h,'Value');
len=evalin('base','length(who)');
names=cell(len,1);
for i=1:len
string=cat(2,'getfield(who,{',num2str(i),'})');
names(i)=evalin('base',string);
end
set(variables_h,'String',names);
if(value > length(names))
set(variables_h,'Value',1);
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%% Subfunction: sD_names %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function sD_names
sets_h=getfield(get(gcf,'UserData'),'sets_h');
pre_h=findobj(get(0,'Children'),'Tag','Preprocess');
sD_set = getfield(get(pre_h,'UserData'),'sD_set');
for i=1:length(sD_set)
names{i,1}=cat(2,sprintf('#%d: ',i),sD_set(i).name);
end
set(sets_h,'String',names);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%% Subfunction: sD_stats %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function sD_stats
man_h=findobj(get(0,'Children'),'Tag','Management');
c_names_h=getfield(get(man_h,'UserData'),'new_c_name_h');
sD_name_h=getfield(get(man_h,'UserData'),'new_name_h');
pre_h=findobj(get(0,'Children'),'Tag','Preprocess');
INDEX=getfield(getfield(get(pre_h,'UserData'),'sData'),'INDEX');
MODIFIED=getfield(getfield(get(pre_h,'UserData'),'sData'),'MODIFIED');
value=get(getfield(get(man_h,'UserData'),'sets_h'),'Value');
if value==INDEX
data=get(pre_h,'UserData');
sData=rmfield(data.sData,[{'INDEX'};{'MODIFIED'}]);
data.sD_set(INDEX)=sData;
data.sData.MODIFIED=0;
set(pre_h,'UserData',data);
end
sData=getfield(getfield(get(pre_h,'UserData'),'sD_set'),{value});
string1=[{sData.name}];
set(sD_name_h,'String',string1);
set(c_names_h,'String',sData.comp_names);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%% Subfunction: import %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function import(varargin)
if nargin==1
LOG=1;
man_h=findobj(get(0,'Children'),'Tag','Management');
set(0,'CurrentFigure',man_h);
name=varargin;
else
LOG=0;
end
variables_h=getfield(get(gcf,'UserData'),'variables_h');
pre_h=findobj(get(0,'Children'),'Tag','Preprocess');
if ~LOG
name=getfield(get(variables_h,'String'),{get(variables_h,'Value')});
end
errstr='Data to be imported must be real matrix or ''som_data_struct''.';
new_sD=evalin('base',name{1});
if isempty(pre_h)
errordlg('''Preprocess'' -figure does not exist. Terminating program...');
close_preprocess;
return;
end
if isstr(new_sD) | (~isstruct(new_sD) & ~isreal(new_sD))
errordlg(errstr);
return;
elseif isstruct(new_sD) & length(new_sD) > 1
errordlg(errstr)
return;
elseif ~isstruct(new_sD)
new_sD=som_data_struct(new_sD);
new_sD.name=name{1};
end
new_sD_names=fieldnames(new_sD);
right_names=fieldnames(som_data_struct(1));
for i=1:length(new_sD_names)
if ~any(strcmp(new_sD_names(i),right_names));
errordlg(errstr);
return;
end
end
data=get(pre_h,'UserData');
data.sD_set(length(data.sD_set) + 1)=new_sD;
if ~LOG
data.LOG{length(data.LOG)+1}='% Import a data set from the workspace.';
data.LOG{length(data.LOG)+1}=cat(2,'preprocess(''import'',''',...
name{1},''');');
end
set(pre_h,'UserData',data);
sD_names;
sD_stats;
old =gcf;
set(0,'CurrentFigure',pre_h);
set_sD_stats;
set(0,'CurrentFigure',old);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%% Subfunction: export %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function export(varargin)
if nargin == 1
LOG=1;
man_h=findobj(get(0,'Children'),'Tag','Management');
set(0,'CurrentFigure',man_h);
index=str2num(varargin{1});
else
LOG=0;
end
sets_h=getfield(get(gcf,'UserData'),'sets_h');
pre_h=findobj(get(0,'Children'),'Tag','Preprocess');
if ~LOG
index=get(sets_h,'Value');
end
if isempty(pre_h)
errordlg('''Preprocess''-figure does not exist. Terminating program...');
close(findobj(get(0,'Children'),'Tag','Management'));
close(findobj(get(0,'Children'),'Tag','PlotWin'));
return;
end
sData=getfield(getfield(get(pre_h,'UserData'),'sD_set'),{index});
if ~isvalid_var_name({sData.name})
return;
end
assignin('base',sData.name,sData);
disp(sprintf('Data set ''%s'' is set to the workspace.',sData.name));
if ~LOG
data=get(pre_h,'UserData');
data.LOG{length(data.LOG)+1}='% Export a data set to the workspace.';
data.LOG{length(data.LOG)+1}=cat(2,'preprocess(''export'',''',...
num2str(index),''');');
set(pre_h,'UserData',data);
end
set_var_names;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%% Subfunction: rename %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function rename(arg)
i=1;
while i <= length(arg) & arg(i) ~= ' '
i=i+1;
end
arg2=arg(i+1:length(arg));
arg=arg(1:i-1);
if ~isempty(arg2)
LOG=1;
i=1;
if arg2(1) ~= '{'
while i <= length(arg2) & arg2(i) ~= ' '
i=i+1;
end
index=str2num(arg2(i+1:length(arg2)));
arg2=arg2(1:i-1);
else
while i <= length(arg2) & arg2(i) ~= '}'
i=i+1;
end
index=str2num(arg2(i+1:length(arg2)));
arg2=arg2(1:i);
end
else
LOG=0;
end
new_name_h=getfield(get(gcf,'UserData'),'new_name_h');
new_c_name_h=getfield(get(gcf,'UserData'),'new_c_name_h');
pre_h=findobj(get(0,'Children'),'Tag','Preprocess');
if isempty(pre_h)
errordlg('''Preprocess'' -figure does not exist. Terminating program...');
close_preprocess;
return;
end
switch arg
case 'set'
if LOG
name={arg2};
else
name=get(new_name_h,'String');
end
if ~isempty(name{1}) & ~any(isspace(name{1}))
if ~isvalid_var_name(name)
sD_stats;
return;
end
if ~LOG
index=get(getfield(get(gcf,'UserData'),'sets_h'),'Value');
end
data=get(pre_h,'UserData');
tmp_set.name=name{1};
data.sD_set(index).name=name{1};
if data.sData.INDEX == index
data.sData.name=name{1};
end
if ~LOG
data.LOG{length(data.LOG)+1}='% Rename a data set.';
data.LOG{length(data.LOG)+1}=cat(2,'preprocess(''rename'',''set ',...
name{1},' ',...
num2str(index),...
''');');
end
set(pre_h,'UserData',data);
sD_names;
string=get(data.sD_set_h,'String');
string{index}=cat(2,sprintf('#%d: ',index),name{1});
set(data.sD_set_h,'String',string);
string=get(data.sD_name_h,'String');
string{1}=cat(2,'Name: ',name{1});
if index==data.sData.INDEX
set(data.sD_name_h,'String',string);
end
else
sD_stats;
end
case 'comp'
if ~LOG
names=get(new_c_name_h,'String');
index=get(getfield(get(gcf,'UserData'),'sets_h'),'Value');
else
names=eval(arg2);
end
if check_cell_names(names)
data=get(pre_h,'UserData');
sData=data.sD_set(index);
if length(sData.comp_names)==length(names)
data.sD_set(index).comp_names=names;
if index == data.sData.INDEX
for i=1:length(names)
names{i}=cat(2,sprintf('#%d: ',i),names{i});
end
set(data.comp_names_h,'String',names);
set(data.sel_comp_h,'String',names);
end
if ~LOG
data.LOG{length(data.LOG)+1}='% Rename components.';
str='preprocess(''rename'',''comp {';
for i=1:length(names)-1
str=cat(2,str,'''''',names{i},''''',');
end
str=cat(2,str,'''''',names{length(names)},'''''} ',...
num2str(index),''');');
data.LOG{length(data.LOG)+1}=str;
else
set(new_c_name_h,'String',names);
end
set(pre_h,'UserData',data);
else
errordlg('There are less components in data.');
sD_stats;
return;
end
else
sD_stats;
end
end
%%% Subfunction: check_cell_names %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function bool=check_cell_names(names)
bool = 1;
if isempty(names)
bool= 0;
return;
end
for i=1:length(names)
if isempty(names{i}) | isspace(names{i})
bool = 0;
return;
end
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%% Subfunction: isvalid_var_name %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function bool=isvalid_var_name(name)
bool=1;
tmp=name{1};
if ~((tmp(1)>='a' & tmp(1)<='z') | (tmp(1)>='A' & tmp(1)<='Z'))
errordlg('Invalid name.');
bool=0;
return;
end
for j=1:length(tmp)
if ~((tmp(j)>='a' & tmp(j)<='z') | ...
(tmp(j)>='A' & tmp(j)<='Z') | ...
(j>1 & tmp(j) == '_') | ...
(tmp(j)>='0' & tmp(j) <= '9')) | tmp(j) == '.'
errordlg('Invalid name.');
bool=0;
return;
end
if j == length(tmp) & tmp(j) == '_'
errordlg('Invalid name.');
bool=0;
return;
end
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%% Subfunction: copy_delete %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function copy_delete(arg)
i=1;
while i <= length(arg) & arg(i) ~= ' '
i=i+1;
end
arg2=arg(i+1:length(arg));
arg=arg(1:i-1);
if ~isempty(arg2)
index=str2num(arg2);
LOG=1;
else
LOG=0;
end
sets_h=getfield(get(gcf,'UserData'),'sets_h');
if ~LOG
index=get(sets_h,'Value');
end
pre_h=findobj(get(0,'Children'),'Tag','Preprocess');
if isempty(pre_h)
errordlg('''Preprocess'' -figure does not exist. Terminating program.');
close_preprocess;
return;
end
switch arg
case 'copy'
data=get(pre_h,'UserData');
data.sD_set(length(data.sD_set)+1)=data.sD_set(index);
if ~LOG
data.LOG{length(data.LOG)+1}='% Copy a data set.';
data.LOG{length(data.LOG)+1}=cat(2,'preprocess(''copy_delete'',''',...
'copy ',num2str(index),''');');
end
set(pre_h,'UserData',data);
sD_names;
old=gcf;
set(0,'CurrentFigure',pre_h);
set_sD_stats;
set(0,'CurrentFigure',old);
case 'delete'
if length(get(sets_h,'String')) == 1
msgbox('No data left. Closing program...')
close_preprocess;
return;
end
data=get(pre_h,'UserData');
if ~isempty(data.undo) & any(strcmp('index',fieldnames(data.undo)))
if data.undo.index > index
data.undo.index = data.undo.index-1;
elseif data.undo.index==index;
data.undo=[];
end
end
set1=data.sD_set(1:index-1);
set2=data.sD_set(index+1:length(data.sD_set));
if ~isempty(set1)
data.sD_set=[set1 set2];
else
data.sD_set=set2;
end
if ~LOG
data.LOG{length(data.LOG)+1}='% Delete a data set.';
data.LOG{length(data.LOG)+1}=cat(2,'preprocess(''copy_delete'',''',...
'delete ',num2str(index),''');');
end
set(pre_h,'UserData',data);
set(sets_h,'Value',1);
sD_names;
sD_stats;
old = gcf;
set(0,'CurrentFigure',pre_h);
for i=1:length(data.sD_set)
string{i}=cat(2,sprintf('#%d: ',i),data.sD_set(i).name);
end
set(data.sD_set_h,'String',string);
data.sData=data.sD_set(get(data.sD_set_h,'Value'));
data.sData.MODIFIED=0;
data.sData.INDEX=1;
set(gcf,'UserData',data);
write_sD_stats;
sData=getfield(get(gcf,'UserData'),'sData');
if sData.INDEX > index
value=get(getfield(get(gcf,'UserData'),'sD_set_h'),'Value');
set(getfield(get(gcf,'UserData'),'sD_set_h'),'Value',value-1);
sData.INDEX = sData.INDEX -1;
elseif sData.INDEX == index
set(getfield(get(gcf,'UserData'),'sD_set_h'),'Value',1);
end
sel_sD;
set(0,'CurrentFigure',old);
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function clipping(varargin)
if nargin ~= 1
pre_h=findobj(get(0,'Children'),'Tag','Preprocess');
preh_udata=get(pre_h,'UserData');
preh_udata.LOG{length(preh_udata.LOG)+1}=...
'% Starting the ''Clipping'' -window...';
preh_udata.LOG{length(preh_udata.LOG)+1}='preprocess(''clipping'',''foo'');';
set(pre_h,'UserData',preh_udata);
end
clip_h=findobj(get(0,'Children'),'Tag','Clipping');
if ~isempty(clip_h)
figure(clip_h);
return;
end
h0 = figure('Color',[0.8 0.8 0.8], ...
'PaperPosition',[18 180 575 432], ...
'PaperUnits','points', ...
'Position',[718 389 300 249], ...
'Tag','Clipping');
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[0.752941176470588 0.752941176470588 0.752941176470588], ...
'ListboxTop',0, ...
'Position',[0.03 0.03614457831325301 0.4666666666666667 0.9236947791164658], ...
'Style','frame', ...
'Tag','Frame1');
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[0.8 0.8 0.8], ...
'ListboxTop',0, ...
'Position',[0.05333333333333334 0.5983935742971887 0.42 0.3333333333333333], ...
'Style','frame', ...
'Tag','Frame2');
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[0.8 0.8 0.8], ...
'ListboxTop',0, ...
'Style','frame', ...
'Position',[0.05333333333333334 0.33 0.42 0.24], ...
'Tag','Frame3');
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[0.8 0.8 0.8], ...
'ListboxTop',0, ...
'Style','frame', ...
'Position',[0.05333333333333334 0.06 0.42 0.24],...
'Tag','Frame4');
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'ListboxTop',0, ...
'Position',[0.5133333333333334 0.6385542168674698 0.4666666666666667 0.321285140562249], ...
'Style','frame', ...
'Tag','Frame5');
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[0.8 0.8 0.8], ...
'ListboxTop',0, ...
'Position',[0.5366666666666667 0.6666666666666666 0.42 0.2650602409638554], ...
'Style','frame', ...
'Tag','Frame6');
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[1 1 1], ...
'ListboxTop',0, ...
'Position',[0.31 0.823293172690763 0.15 0.09638554216867469], ...
'Style','edit', ...
'Tag','EditText1');
data.big_val_h = h1;
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[1 1 1], ...
'ListboxTop',0, ...
'Position',[0.31 0.7148594377510039 0.15 0.09638554216867469], ...
'Style','edit', ...
'Tag','EditText2');
data.small_val_h = h1;
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[1 1 1], ...
'ListboxTop',0, ...
'Position',[0.31 0.606425702811245 0.15 0.09638554216867469], ...
'Style','edit', ...
'Tag','EditText3');
data.equal_val_h=h1;
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[0.8 0.8 0.8], ...
'FontSize',6, ...
'FontWeight','demi', ...
'HorizontalAlignment','left', ...
'ListboxTop',0, ...
'Position',[0.06000000000000001 0.8473895582329316 0.22 0.05622489959839357], ...
'String','Bigger than', ...
'Style','text', ...
'Tag','StaticText1');
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[0.8 0.8 0.8], ...
'FontSize',6, ...
'FontWeight','demi', ...
'HorizontalAlignment','left', ...
'ListboxTop',0, ...
'Position',[0.06000000000000001 0.7389558232931727 0.24 0.04819277108433735], ...
'String','Smaller than', ...
'Style','text', ...
'Tag','StaticText2');
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[0.8 0.8 0.8], ...
'FontSize',6, ...
'FontWeight','demi', ...
'HorizontalAlignment','left', ...
'ListboxTop',0, ...
'Position',[0.06000000000000001 0.610441767068273 0.22 0.07228915662650602], ...
'String','Equal to', ...
'Style','text', ...
'Tag','StaticText3');
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[0.8 0.8 0.8], ...
'ListboxTop',0, ...
'Position',[0.07000000000000001 0.465863453815261 0.06333333333333334 0.07228915662650602], ...
'Style','radiobutton', ...
'Value',1,...
'Tag','Radiobutton1');
data.and_button_h=h1;
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[0.8 0.8 0.8], ...
'ListboxTop',0, ...
'Position',[0.07000000000000001 0.3734939759036144 0.06333333333333334 0.07228915662650602], ...
'Style','radiobutton', ...
'Tag','Radiobutton2');
data.or_button_h=h1;
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[0.8 0.8 0.8], ...
'FontSize',6, ...
'FontWeight','demi', ...
'HorizontalAlignment','left', ...
'Position',[0.1466666666666667 0.45 0.2333333333333333 0.07228915662650602], ...
'String','AND', ...
'Style','text', ...
'Tag','StaticText4');
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[0.8 0.8 0.8], ...
'FontSize',6, ...
'FontWeight','demi', ...
'HorizontalAlignment','left', ...
'ListboxTop',0, ...
'String','OR', ...
'Position',[0.1466666666666667 0.35 0.2333333333333333 0.07228915662650602], ...
'Style','text', ...
'Tag','StaticText5');
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[0.8 0.8 0.8], ...
'ListboxTop',0, ...
'Position',[0.07000000000000001 0.1967871485943775 0.06333333333333334 0.07228915662650602], ...
'Style','radiobutton', ...
'Value',1,...
'Tag','Radiobutton3');
data.all_button_h=h1;
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[0.8 0.8 0.8], ...
'ListboxTop',0, ...
'Position',[0.07000000000000001 0.09236947791164658 0.06333333333333334 0.07228915662650602], ...
'Style','radiobutton', ...
'Tag','Radiobutton4');
data.sel_vects_button_h=h1;
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[0.8 0.8 0.8], ...
'FontSize',6, ...
'FontWeight','demi', ...
'HorizontalAlignment','left', ...
'ListboxTop',0, ...
'Position',[0.1466666666666667 0.1927710843373494 0.2333333333333333 0.07228915662650602], ...
'String','All vectors', ...
'Style','text', ...
'Tag','StaticText6');
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[0.8 0.8 0.8], ...
'FontSize',6, ...
'FontWeight','demi', ...
'HorizontalAlignment','left', ...
'ListboxTop',0, ...
'Position',[0.1466666666666667 0.09638554216867469 0.3133333333333334 0.05622489959839357], ...
'String','Among selected', ...
'Style','text', ...
'Tag','StaticText7');
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[1 1 1], ...
'ListboxTop',0, ...
'Position',[0.7866666666666667 0.823293172690763 0.1366666666666667 0.09236947791164658], ...
'Style','edit', ...
'Tag','EditText4');
data.replace_val_h=h1;
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[0.8 0.8 0.8], ...
'FontSize',6, ...
'FontWeight','demi', ...
'HorizontalAlignment','left', ...
'ListboxTop',0, ...
'Position',[0.5633333333333334 0.8273092369477911 0.2066666666666667 0.07630522088353413], ...
'String','Replace', ...
'Style','text', ...
'Tag','StaticText8');
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'FontWeight','demi', ...
'ListboxTop',0, ...
'Position',[0.5700000000000001 0.6827309236947791 0.3566666666666667 0.08032128514056225], ...
'String','Replace', ...
'Tag','Pushbutton1');
data.OK_button_h=h1;
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[0.752941176470588 0.752941176470588 0.752941176470588], ...
'Callback','preprocess close_c',...
'FontWeight','demi', ...
'ListboxTop',0, ...
'Position',[0.6633333333333333 0.07228915662650602 0.2833333333333333 0.09638554216867469], ...
'String','Close', ...
'Tag','Pushbutton2');
data.state.and=1;
data.state.all=1;
data.state.big=[];
data.state.small=[];
data.state.equal=[];
data.state.replace=[];
set(data.or_button_h,'Callback','preprocess and_or_cb or');
set(data.and_button_h,'Callback','preprocess and_or_cb and');
set(data.and_button_h,'Value',1);
set(data.all_button_h,'Callback','preprocess all_sel_cb all');
set(data.sel_vects_button_h,'Callback','preprocess all_sel_cb sel');
set(data.big_val_h,'Callback','preprocess set_state_vals big');
set(data.small_val_h,'Callback','preprocess set_state_vals small');
set(data.equal_val_h,'Callback','preprocess set_state_vals equal');
set(data.replace_val_h,'Callback','preprocess set_state_vals replace');
set(data.OK_button_h,'Callback','preprocess clip_data clip');
set(h0,'UserData',data);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function select(varargin)
if nargin ~= 1
pre_h=findobj(get(0,'Children'),'Tag','Preprocess');
preh_udata=get(pre_h,'UserData');
preh_udata.LOG{length(preh_udata.LOG)+1}=...
'% Starting the ''Select'' -window...';
preh_udata.LOG{length(preh_udata.LOG)+1}='preprocess(''select'',''foo'');';
set(pre_h,'UserData',preh_udata);
end
sel_h=findobj(get(0,'Children'),'Tag','Select');
if ~isempty(sel_h)
figure(sel_h);
return;
end
h0 = figure('Color',[0.8 0.8 0.8], ...
'PaperPosition',[18 180 576 432], ...
'PaperUnits','points', ...
'Position',[750 431 168 365], ...
'Tag','Select');
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[0.752941176470588 0.752941176470588 0.752941176470588], ...
'ListboxTop',0, ...
'Position',[0.05357142857142857 0.2712328767123288 0.8333333333333333 0.6301369863013698], ...
'Style','frame', ...
'Tag','Frame1');
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[0.752941176470588 0.752941176470588 0.752941176470588], ...
'ListboxTop',0, ...
'Position',[0.05357142857142857 0.1041095890410959 0.8333333333333333 0.1397260273972603], ...
'Style','frame', ...
'Tag','Frame2');
h1 = uicontrol('Parent',h0,...
'Units','normalized', ...
'BackgroundColor',[0.8 0.8 0.8], ...
'ListboxTop',0, ...
'Position',[0.09523809523809523 0.6547945205479452 0.75 0.2273972602739726], ...
'Style','frame', ...
'Tag','Frame3');
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[0.8 0.8 0.8], ...
'ListboxTop',0, ...
'Position',[0.09523809523809523 0.4794520547945206 0.75 0.1506849315068493], ...
'Style','frame', ...
'Tag','Frame4');
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[0.8 0.8 0.8], ...
'ListboxTop',0, ...
'Position',[0.09523809523809523 0.2986301369863014 0.75 0.1506849315068493], ...
'Style','frame', ...
'Tag','Frame5');
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[1 1 1], ...
'ListboxTop',0, ...
'Position',[0.5535714285714285 0.8082191780821918 0.2678571428571429 0.06575342465753425], ...
'Style','edit', ...
'Tag','EditText1');
data.big_val_h=h1;
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[1 1 1], ...
'ListboxTop',0, ...
'Position',[0.5535714285714285 0.7342465753424657 0.2678571428571429 0.06575342465753425], ...
'Style','edit', ...
'Tag','EditText2');
data.small_val_h = h1;
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[1 1 1], ...
'ListboxTop',0, ...
'Position',[0.5535714285714285 0.6602739726027397 0.2678571428571429 0.06575342465753425], ...
'Style','edit', ...
'Tag','EditText3');
data.equal_val_h=h1;
h1 = uicontrol('Parent',h0, ...
'BackgroundColor',[0.8 0.8 0.8], ...
'Units','normalized', ...
'FontWeight','demi', ...
'FontSize',8,...
'HorizontalAlignment','left', ...
'ListboxTop',0, ...
'Position',[0.1071 0.8247 0.3929 0.0384], ...
'String','Bigger than', ...
'Style','text', ...
'Tag','StaticText1');
h1 = uicontrol('Parent',h0, ...
'BackgroundColor',[0.8 0.8 0.8], ...
'Units','normalized', ...
'FontWeight','demi', ...
'FontSize',8,...
'HorizontalAlignment','left', ...
'ListboxTop',0, ...
'Position',[0.1071 0.7507 0.4286 0.0329], ...
'String','Smaller than', ...
'Style','text', ...
'Tag','StaticText2');
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[0.8 0.8 0.8], ...
'FontWeight','demi', ...
'FontSize',8,...
'HorizontalAlignment','left', ...
'ListboxTop',0, ...
'Position',[0.1071 0.6630 0.3929 0.0493], ...
'String','Equal to', ...
'Style','text', ...
'Tag','StaticText3');
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[0.8 0.8 0.8], ...
'ListboxTop',0, ...
'Position',[0.125 0.5643835616438356 0.1130952380952381 0.04931506849315068], ...
'Style','radiobutton', ...
'Tag','Radiobutton1');
data.and_button_h = h1;
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[0.8 0.8 0.8], ...
'ListboxTop',0, ...
'Position',[0.125 0.5013698630136987 0.1130952380952381 0.04931506849315068], ...
'Style','radiobutton', ...
'Tag','Radiobutton2');
data.or_button_h = h1;
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[0.8 0.8 0.8], ...
'FontWeight','demi', ...
'FontSize',8,...
'HorizontalAlignment','left', ...
'ListboxTop',0, ...
'Position',[0.2619047619047619 0.5561643835616439 0.3809523809523809 0.05205479452054795], ...
'String','AND', ...
'Style','text', ...
'Tag','StaticText4');
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[0.8 0.8 0.8], ...
'FontWeight','demi', ...
'FontSize',8,...
'HorizontalAlignment','left', ...
'ListboxTop',0, ...
'Position',[0.2619047619047619 0.4986301369863014 0.3809523809523809 0.04657534246575343], ...
'String','OR', ...
'Style','text', ...
'Tag','StaticText5');
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[0.8 0.8 0.8], ...
'ListboxTop',0, ...
'Position',[0.125 0.3808219178082192 0.1130952380952381 0.04931506849315068], ...
'Style','radiobutton', ...
'Tag','Radiobutton3');
data.all_button_h = h1;
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[0.8 0.8 0.8], ...
'ListboxTop',0, ...
'Position',[0.125 0.3095890410958904 0.1130952380952381 0.04931506849315068], ...
'Style','radiobutton', ...
'Tag','Radiobutton4');
data.sel_vects_button_h = h1;
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[0.8 0.8 0.8], ...
'FontWeight','demi', ...
'FontSize',8,...
'HorizontalAlignment','left', ...
'ListboxTop',0, ...
'Position',[0.2619047619047619 0.3780821917808219 0.4166666666666666 0.04931506849315068], ...
'String','All vectors', ...
'Style','text', ...
'Tag','StaticText6');
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[0.8 0.8 0.8], ...
'FontWeight','demi', ...
'FontSize',8,...
'HorizontalAlignment','left', ...
'ListboxTop',0, ...
'Position',[0.2619047619047619 0.3123287671232877 0.5595238095238095 0.03835616438356165], ...
'String','Among selected', ...
'Style','text', ...
'Tag','StaticText7');
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[0.8 0.8 0.8], ...
'ListboxTop',0, ...
'Position',[0.0952 0.1178 0.7500 0.1068], ...
'Style','frame', ...
'Tag','Frame6');
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[1 1 1], ...
'ListboxTop',0, ...
'Position',[0.5298 0.1342 0.2738 0.0712], ...
'Style','edit', ...
'Tag','EditText4');
data.replace_val_h = h1;
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[0.8 0.8 0.8], ...
'FontSize',8,...
'FontWeight','demi', ...
'HorizontalAlignment','left', ...
'ListboxTop',0, ...
'Position',[0.1369047619047619 0.136986301369863 0.3214285714285714 0.06027397260273973], ...
'String','Vectors', ...
'Style','text', ...
'Tag','StaticText8');
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[0.752941176470588 0.752941176470588 0.752941176470588], ...
'FontWeight','demi', ...
'ListboxTop',0, ...
'Position',[0.05357142857142857 0.01917808219178082 0.3869047619047619 0.0684931506849315], ...
'String','OK', ...
'Tag','Pushbutton1');
data.OK_button_h = h1;
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[0.752941176470588 0.752941176470588 0.752941176470588], ...
'Callback','preprocess close_s',...
'FontWeight','demi', ...
'ListboxTop',0, ...
'Position',[0.5 0.01917808219178082 0.3869047619047619 0.0684931506849315], ...
'String','Close', ...
'Tag','Pushbutton2');
data.state.and=1;
data.state.all=1;
data.state.big=[];
data.state.small=[];
data.state.equal=[];
data.state.replace=[];
set(data.or_button_h,'Callback','preprocess and_or_cb or');
set(data.and_button_h,'Callback','preprocess and_or_cb and');
set(data.and_button_h,'Value',1);
set(data.all_button_h,'Callback','preprocess all_sel_cb all');
set(data.sel_vects_button_h,'Callback','preprocess all_sel_cb sel');
set(data.big_val_h,'Callback','preprocess set_state_vals big');
set(data.small_val_h,'Callback','preprocess set_state_vals small');
set(data.equal_val_h,'Callback','preprocess set_state_vals equal');
set(data.replace_val_h,'Callback','preprocess set_state_vals replace');
set(data.OK_button_h,'Callback','preprocess clip_data sel');
set(h0,'UserData',data);
%%% Subfunction: and_or_cb %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function and_or_cb(arg)
%AND_OR_CB A callback function. Checks that only one of the radiobox
% buttons 'AND' and 'OR' is pressed down.
%
%
and_button_h=getfield(get(gcf,'UserData'),'and_button_h');
or_button_h=getfield(get(gcf,'UserData'),'or_button_h');
data=get(gcf,'UserData');
switch arg
case 'or'
set(and_button_h,'Value',0);
set(or_button_h,'Value',1);
data.state.and=0;
case 'and'
set(or_button_h,'Value',0);
set(and_button_h,'Value',1);
data.state.and=1;
end
set(gcf,'UserData',data);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%% Subfunction: all_sel_cb %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function all_sel_cb(arg)
all_button_h=getfield(get(gcf,'UserData'),'all_button_h');
sel_vects_button_h=getfield(get(gcf,'UserData'),'sel_vects_button_h');
data=get(gcf,'UserData');
switch arg
case 'all'
set(sel_vects_button_h,'Value',0);
set(all_button_h,'Value',1);
data.state.all=1;
case 'sel'
set(all_button_h,'Value',0);
set(sel_vects_button_h,'Value',1);
data.state.all=0;
end
set(gcf,'UserData',data);
%%% Subfunction: set_state_vals %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function set_state_vals(arg)
%SET_STATE_VALS sets the values to the UserData's state-struct.
%
%
data=get(gcf,'UserData');
switch arg
case 'big'
big_val_h=getfield(get(gcf,'UserData'),'big_val_h');
val =str2num(get(big_val_h,'String'));
dims=size(val);
if dims(1) ~= 1 | dims(2) ~= 1
errordlg('Argument of the operation must be scalar.');
set(big_val_h,'String','');
return;
end
if isreal(val)
data.state.big=val;
else
errordlg('Limits of the operation must be real.');
set(big_val_h,'String','');
return;
end
case 'small'
small_val_h=getfield(get(gcf,'UserData'),'small_val_h');
val=str2num(get(small_val_h,'String'));
dims=size(val);
if dims(1) ~= 1 | dims(2) ~= 1
errordlg('Argument of the operation must be scalar.')
set(small_val_h,'String','');
return;
end
if isreal(val)
data.state.small=val;
else
errordlg('Limits of the operation must be real.');
set(small_val_h,'String','');
return;
end
case 'equal'
equal_val_h=getfield(get(gcf,'UserData'),'equal_val_h');
val = str2num(get(equal_val_h,'String'));
dims=size(val);
if dims(1) ~= 1 | dims(2) ~= 1
errordlg('Argument of the operation must be scalar.');
set(equal_val_h,'String','');
return;
end
if isreal(val)
data.state.equal=val;
else
errordlg('Limits of the operation must be real.');
set(equal_val_h,'String','');
return;
end
case 'replace'
replace_val_h=getfield(get(gcf,'UserData'),'replace_val_h');
val=str2num(get(replace_val_h,'String'));
dims=size(val);
if (dims(1) ~= 1 | dims(2) ~= 1) & ~strcmp(get(gcf,'Tag'),'Select')
errordlg('Argument of the operation must be scalar.');
set(replace_val_h,'String','');
return;
end
if isreal(val)
data.state.replace=val;
else
errordlg('Limits of the operation must be real.');
set(replace_val_h,'String','');
return;
end
end
set(gcf,'UserData',data);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%% Subfunction: clip_data %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function clip_data(arg)
%CLIP_DATA A callback function. Filters the data.
%
%
i=1;
while i <= length(arg) & arg(i) ~= ' '
i=i+1;
end
arg2=arg(i+1:length(arg));
arg=arg(1:i-1);
if ~isempty(arg2)
LOG=1;
if strcmp(arg,'sel')
c_h=findobj(get(0,'Children'),'Tag','Select');
else
c_h=findobj(get(0,'Children'),'Tag','Clipping');
end
set(0,'CurrentFigure',c_h);
i=1;
while i <= length(arg2) & arg2(i) ~= ' '
i=i+1;
end
BT=str2num(arg2(1:i-1));
i=i+1;
j=i;
while i <= length(arg2) & arg2(i) ~= ' '
i=i+1;
end
ST=str2num(arg2(j:i-1));
i=i+1;
j=i;
while i <= length(arg2) & arg2(i) ~= ' '
i=i+1;
end
EQ=str2num(arg2(j:i-1));
i=i+1;
j=i;
while i <= length(arg2) & arg2(i) ~= ' '
i=i+1;
end
AND_OR=str2num(arg2(j:i-1));
i=i+1;
j=i;
while i <= length(arg2) & arg2(i) ~= ' '
i=i+1;
end
ALL_AMONG=str2num(arg2(j:i-1));
i=i+1;
j=i;
while i <= length(arg2)
i=i+1;
end
VECT_REPL=str2num(arg2(j:i-1));
else
LOG=0;
end
if ~LOG
big_val_h=getfield(get(gcf,'UserData'),'big_val_h');
small_val_h=getfield(get(gcf,'UserData'),'small_val_h');
equal_val_h=getfield(get(gcf,'UserData'),'equal_val_h');
replace_val_h=getfield(get(gcf,'UserData'),'replace_val_h');
end
pre_h=findobj(get(0,'Children'),'Tag','Preprocess');
if isempty(pre_h)
errordlg('''Preprocess'' -figure does not exist. Terminating program...');
pro_tools('close');
return;
end
comp_names_h=getfield(get(pre_h,'UserData'),'comp_names_h');
selected=getfield(get(pre_h,'UserData'),'selected_vects');
sData=getfield(get(pre_h,'UserData'),'sData');
undo = sData;
state=getfield(get(gcf,'UserData'),'state');
if LOG
state.big=BT;
state.small=ST;
state.equal=EQ;
state.replace=VECT_REPL;
state.and=AND_OR;
state.all=ALL_AMONG;
end
if isempty(pre_h)
pro_tools('close');
end
if isempty(get(comp_names_h,'Value'))
clear_state_vals;
errordlg('There must be one component chosen for the operation.');
return;
end
n_th_comp=getfield(get_indices,{1});
if isempty(state.big) & isempty(state.small) & isempty(state.equal) & ...
strcmp(arg,'clip')
clear_state_vals;
errordlg('At least one limit must be chosen for the-operation.');
return;
end
if ~isempty(state.replace) & strcmp(arg,'sel')
if ~all(state.replace == round(state.replace)) | any(state.replace < 1)
errordlg('Indices of vectors must be positive integers.');
return;
elseif any(state.replace > length(sData.data(:,1)))
errordlg('Indices of the vectors to be selected are too big.');
return;
end
end
if isempty(state.replace) & strcmp(arg,'clip')
clear_state_vals;
errordlg('Replace value must be determined for Clipping-operation.');
return;
end
if isempty(state.big) & isempty(state.small) & isempty(state.equal) & ...
isempty(state.replace)
clear_state_vals;
return;
end
bt_indices=[];
lt_indices=[];
equal_indices=[];
if ~isempty(state.big)
if state.all
bt_indices=find(sData.data(:,n_th_comp) > state.big);
else
bt_indices=selected(find(sData.data(selected,n_th_comp) > state.big));
end
end
if ~isempty(state.small)
if state.all
lt_indices=find(sData.data(:,n_th_comp) < state.small);
else
lt_indices=selected(find(sData.data(selected,n_th_comp) < state.small));
end
end
if ~isempty(state.equal)
if isnan(state.equal)
if state.all
equal_indices=find(isnan(sData.data(:,n_th_comp)));
else
equal_indices=selected(find(isnan(sData.data(selected,n_th_comp))));
end
elseif state.all
equal_indices=find(sData.data(:,n_th_comp)==state.equal);
else
equal_indices=selected(find(sData.data(selected,n_th_comp)==state.equal));
end
end
if state.and
if ~isempty(bt_indices) | ~isempty(lt_indices) | ~isempty(equal_indices)...
| strcmp(arg,'sel')
if isempty(bt_indices) & isempty(lt_indices) & isempty(equal_indices) &...
isempty(state.replace)
clear_state_vals;
return;
end
if isempty(bt_indices)
if ~state.all
bt_indices=selected;
else
bt_indices=1:getfield(size(sData.data),{1});
end
end
if isempty(lt_indices)
if ~state.all
lt_indices=selected;
else
lt_indices=1:getfield(size(sData.data),{1});
end
end
if isempty(equal_indices)
if ~state.all
equal_indices=selected;
else
equal_indices=1:getfield(size(sData.data),{1});
end
end
indices=intersect(intersect(bt_indices,lt_indices),equal_indices);
if strcmp(arg,'sel')
if ~isempty(indices) | ~isempty(state.replace)
if isempty(state.replace)
NOTEMPTY=0;
if ~state.all
state.replace=selected;
else
state.replace=1:getfield(size(sData.data),{1});
end
else
NOTEMPTY=1;
end
if isempty(indices)
indices=selected;
end
indices=intersect(indices,state.replace);
if isempty(indices)
indices=selected;
end
data=get(pre_h,'UserData');
data.undo.sData=sData;
data.undo.selected=data.selected_vects;
data.selected_vects=indices;
if ~LOG
if ~NOTEMPTY
data.LOG{length(data.LOG)+1}='% Select vectors.';
data.LOG{length(data.LOG)+1}=cat(2,'preprocess(''clip_data'',''',...
arg,...
' ',num2str(state.big),...
' ',num2str(state.small),...
' ',num2str(state.equal),...
' ',num2str(state.and),...
' ',num2str(state.all),...
''');');
else
code=write_log_code(state.replace,...
arg,...
state.big,...
state.small,...
state.equal,...
state.and,...
state.all);
data.LOG(length(data.LOG)+1:length(data.LOG)+length(code))=code;
end
end
set(pre_h,'UserData',data);
old=gcf;
set(0,'CurrentFigure',pre_h);
sel_comp;
cplot_mimema;
vect_means(data.sData,data.vect_mean_h,data.selected_vects);
set(0,'CurrentFigure',old);
end
clear_state_vals;
return;
end
sData.data(indices,n_th_comp) = state.replace;
sData.MODIFIED=1;
end
else
indices=union(union(bt_indices,lt_indices),equal_indices);
if ~isempty(indices) | strcmp(arg,'sel')
if strcmp(arg,'sel')
if ~isempty(indices) | ~isempty(state.replace')
data=get(pre_h,'UserData');
data.undo.sData=sData;
data.undo.selected=data.selected_vects;
data.selected_vects=union(indices,state.replace);
if ~LOG
if isempty(state.replace);
data.LOG{length(data.LOG)+1}=cat(2,'preprocess(''clip_data'',''',...
arg,...
' ',num2str(state.big),...
' ',num2str(state.small),...
' ',num2str(state.equal),...
' ',num2str(state.and),...
' ',num2str(state.all),...
''');');
else
code=write_log_code(state.replace,...
arg,...
state.big,...
state.small,...
state.equal,...
state.and,...
state.all);
data.LOG(length(data.LOG)+1:length(data.LOG)+length(code))=code;
end
end
set(pre_h,'UserData',data);
old=gcf;
set(0,'CurrentFigure',pre_h);
sel_comp;
vect_means(data.sData,data.vect_mean_h,data.selected_vects);
cplot_mimema;
set(0,'CurrentFigure',old);
end
clear_state_vals;
return;
end
sData.data(indices,n_th_comp)=state.replace;
sData.MODIFIED=1;
end
end
if sData.MODIFIED
data=get(pre_h,'UserData');
data.sData=sData;
data.undo.sData=undo;
if ~LOG
if strcmp(arg,'sel')
data.LOG{length(data.LOG)+1}='% Select vectors';
else
data.LOG{length(data.LOG)+1}='% Clip values.';
end
if strcmp(arg,'clip') | isempty(state.replace)
data.LOG{length(data.LOG)+1}=cat(2,'preprocess(''clip_data'',''',arg,...
' ',num2str(state.big),...
' ',num2str(state.small),...
' ',num2str(state.equal),...
' ',num2str(state.and),...
' ',num2str(state.all),...
' ',num2str(state.replace),...
''');');
else
code=write_log_code(state.replace,...
arg,...
state.big,...
state.small,...
state.equal,...
state.and,...
state.all);
data.LOG(length(data.LOG)+1:length(data.LOG)+length(code))=code;
end
end
set(pre_h,'UserData',data);
old=gcf;
set(0,'CurrentFigure',pre_h)
vector_h=getfield(get(gcf,'UserData'),'vector_h');
vect_mean_h=getfield(get(gcf,'UserData'),'vect_mean_h');
set(gcf,'CurrentAxes',vector_h);
vect_means(sData,vect_mean_h,selected);
cplot_mimema;
sel_comp;
set(0,'CurrentFigure',old);
end
clear_state_vals;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%% Subfunction: clear_state_vals %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function clear_state_vals
%CLEAR_STATE_VALS Sets the fields of the UserData's state-struct empty.
%
%
data=get(gcf,'UserData');
set(data.big_val_h,'String','');
set(data.small_val_h,'String','');
set(data.equal_val_h,'String','');
set(data.replace_val_h,'String','');
data.state.big=[];
data.state.small=[];
data.state.equal=[];
data.state.replace=[];
set(gcf,'UserData',data);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function delay(varargin)
delay_h=findobj(get(0,'Children'),'Tag','Delay');
if nargin ~= 1
pre_h=findobj(get(0,'Children'),'Tag','Preprocess');
preh_udata=get(pre_h,'UserData');
preh_udata.LOG{length(preh_udata.LOG)+1}=...
'% Starting the ''Delay'' -window...';
preh_udata.LOG{length(preh_udata.LOG)+1}='preprocess(''delay'',''foo'');';
set(pre_h,'UserData',preh_udata);
end
if ~isempty(delay_h)
figure(delay_h);
return;
end
h0 = figure('Color',[0.8 0.8 0.8], ...
'PaperPosition',[18 180 576 432], ...
'PaperUnits','points', ...
'Position',[759 664 162 215], ...
'Tag','Delay');
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[0.752941176470588 0.752941176470588 0.752941176470588], ...
'ListboxTop',0, ...
'Position',[0.05555555555555555 0.2046511627906977 0.8950617283950617 0.7441860465116279], ...
'Style','frame', ...
'Tag','Frame1');
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[0.8 0.8 0.8], ...
'ListboxTop',0, ...
'Position',[0.08641975308641975 0.6976744186046512 0.8333333333333333 0.2232558139534884], ...
'Style','frame', ...
'Tag','Frame2');
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[0.8 0.8 0.8], ...
'ListboxTop',0, ...
'Position',[0.08641975308641975 0.227906976744186 0.8333333333333333 0.4418604651162791], ...
'Style','frame', ...
'Tag','Frame3');
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'Callback','preprocess delay_data',...
'FontWeight','demi', ...
'ListboxTop',0, ...
'Position',[0.0556 0.0326 0.4012 0.1163], ...
'String','OK', ...
'Tag','Pushbutton1');
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'Callback','preprocess close_d',...
'FontWeight','demi', ...
'ListboxTop',0, ...
'Position',[0.5494 0.0326 0.4012 0.1163], ...
'String','Close', ...
'Tag','Pushbutton2');
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[1 1 1], ...
'ListboxTop',0, ...
'Position',[0.4876543209876543 0.7534883720930232 0.3518518518518519 0.1255813953488372], ...
'Style','edit', ...
'Tag','EditText1');
data.delay_val_h = h1;
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[0.8 0.8 0.8], ...
'FontWeight','demi', ...
'HorizontalAlignment','left', ...
'ListboxTop',0, ...
'Position',[0.1173 0.7860 0.3086 0.0651], ...
'String','Delay', ...
'Style','text', ...
'Tag','StaticText1');
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'Callback','preprocess clip_exp_cb c_this',...
'BackgroundColor',[0.8 0.8 0.8], ...
'ListboxTop',0, ...
'Position',[0.1173 0.5349 0.1173 0.0837], ...
'Style','radiobutton', ...
'Tag','Radiobutton1');
data.c_this_button_h=h1;
data.mode='c_this';
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[0.8 0.8 0.8], ...
'Callback','preprocess clip_exp_cb c_all',...
'ListboxTop',0, ...
'Position',[0.1173 0.4047 0.1173 0.0837], ...
'Style','radiobutton', ...
'Tag','Radiobutton2');
data.c_all_button_h=h1;
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[0.8 0.8 0.8], ...
'Callback','preprocess clip_exp_cb e_all',...
'ListboxTop',0, ...
'Position',[0.1173 0.2651 0.1173 0.0837], ...
'Style','radiobutton', ...
'Tag','Radiobutton3');
data.e_all_button_h=h1;
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[0.8 0.8 0.8], ...
'FontWeight','demi', ...
'FontSize',8,...
'HorizontalAlignment','left', ...
'ListboxTop',0, ...
'Position',[0.26 0.5534883720930233 0.4135802469135802 0.06511627906976744], ...
'String','Clip this', ...
'Style','text', ...
'Tag','StaticText2');
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[0.8 0.8 0.8], ...
'FontWeight','demi', ...
'FontSize',8,...
'HorizontalAlignment','left', ...
'ListboxTop',0, ...
'Position',[0.26 0.413953488372093 0.3765432098765432 0.06511627906976744], ...
'String','Clip all', ...
'Style','text', ...
'Tag','StaticText3');
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[0.8 0.8 0.8], ...
'FontWeight','demi', ...
'FontSize',8,...
'HorizontalAlignment','left', ...
'ListboxTop',0, ...
'Position',[0.26 0.2744186046511628 0.4197530864197531 0.06511627906976744], ...
'String','Expand all', ...
'Style','text', ...
'Tag','StaticText4');
set(gcf,'UserData',data);
%%% Subfunction clip_exp_cb %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function clip_exp_cb(arg)
c_this_button_h=getfield(get(gcf,'UserData'),'c_this_button_h');
c_all_button_h=getfield(get(gcf,'UserData'),'c_all_button_h');
e_all_button_h=getfield(get(gcf,'UserData'),'e_all_button_h');
data=get(gcf,'UserData');
switch arg
case 'c_this'
set(c_all_button_h,'Value',0);
set(e_all_button_h,'Value',0);
set(c_this_button_h,'Value',1);
data.mode='c_this';
case 'c_all'
set(c_this_button_h,'Value',0);
set(e_all_button_h,'Value',0);
set(c_all_button_h,'Value',1);
data.mode='c_all';
case 'e_all'
set(c_this_button_h,'Value',0);
set(c_all_button_h,'Value',0);
set(e_all_button_h,'Value',1);
data.mode='e_all';
end
set(gcf,'UserData',data);
%%% Subfunction: delay_data %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function delay_data(varargin)
if nargin == 1
del_h=findobj(get(0,'Children'),'Tag','Delay');
set(0,'CurrentFigure',del_h);
LOG=1;
arg=varargin{1};
i=1;
while i <= length(arg) & arg(i) ~= ' '
i=i+1;
end
delay=str2num(arg(1:i-1));
no=str2num(arg(i+1:length(arg)));
else
LOG=0;
end
pre_h=findobj(get(0,'Children'),'Tag','Preprocess');
q='Delay operation is not evaluated.';
t='Warning';
if isempty(pre_h)
errordlg('''Preprocess'' -figure does not exist. Terminating program...');
pro_tools('close');
return;
end
sData=getfield(get(pre_h,'UserData'),'sData');
undo = sData;
data=get(gcf,'UserData');
if ~LOG
delay=str2num(get(data.delay_val_h,'String'));
if isempty(delay)
errordlg('Value of ''Delay'' must be defined.');
return
end
set(data.delay_val_h,'String','');
if round(delay) ~= delay
errordlg('Value of ''Delay'' must be integer.');
return;
end
end
comp_names_h=getfield(get(pre_h,'UserData'),'comp_names_h');
if isempty(get(comp_names_h,'Value'))
errordlg('There are not components chosen.');
return;
end
n_th_comp=getfield(get_indices,{1});
len=length(sData.data(:,1));
if LOG
switch no
case 1
data.mode='c_this';
preprocess('clip_exp_cb','c_this');
case 2
data.mode='c_all';
preprocess('clip_exp_cb','c_all');
case 3
data.mode='e_all';
preprocess('clip_exp_cb','e_all');
end
end
switch data.mode
case 'c_this'
MODE='1';
if delay > 0
sData.data(delay+1:len,n_th_comp)=sData.data(1:len-delay);
if delay >= len
errordlg(q,t);
return;
else
sData.data(1:delay,n_th_comp)=NaN;
end
elseif delay < 0
sData.data(1:len+delay,n_th_comp)=...
sData.data(abs(delay)+1:len,n_th_comp);
if abs(delay) >= len
errordlg(q,t);
return;
else
sData.data(len+delay+1:len,n_th_comp)=NaN;
end
end
if delay ~= 0
data=get(pre_h,'UserData');
sData.MODIFIED=1;
sData.comp_norm(n_th_comp)=[];
data.sData=sData;
data.undo.sData=undo;
set(pre_h,'UserData',data);
old = gcf;
set(0,'CurrentFigure',pre_h);
sel_comp;
cplot_mimema;
set(0,'CurrentFigure',old);
end
case 'c_all'
MODE='2';
if delay > 0
sData.data(delay+1:len,n_th_comp)=sData.data(1:len-delay,n_th_comp);
if delay >= len
errordlg(q,t);
return;
else
sData.data=sData.data(delay+1:len,:);
end
elseif delay < 0
sData.data(1:len+delay,n_th_comp)=sData.data(abs(delay)+1:len,n_th_comp);
if abs(delay) >= len
errordlg(q,t);
return;
else
sData.data=sData.data(1:len+delay,:);
end
end
if delay ~= 0
data=get(pre_h,'UserData');
sData.MODIFIED=1;
sData.comp_norm(:,:)={[]};
data.sData=sData;
data.undo.sData=undo;
data.undo.selected=data.selected_vects;
if delay > 0
data.selected_vects=...
data.selected_vects(find(data.selected_vects>delay));
data.selected_vects=data.selected_vects-delay;
elseif nargin == 1
data.selected_vects=...
data.selected_vects(find(data.selected_vects<=len-abs(delay)));
end
set(pre_h,'UserData',data);
old=gcf;
set(0,'CurrentFigure',pre_h);
vects=zeros(1,length(sData.data(:,1)));
vects(data.selected_vects)=1;
write_sD_stats;
draw_vectors(vects,data.vector_h);
sel_comp;
cplot_mimema;
set(0,'CurrentFigure',old);
end
case 'e_all'
MODE='3';
if delay > 0
sData.data(len+1:len+delay,:)=NaN;
sData.data(1+delay:delay+len,n_th_comp)=sData.data(1:len,n_th_comp);
sData.data(1:delay,n_th_comp)=NaN;
elseif delay < 0
delay=abs(delay);
sData.data(delay+1:len+delay,:)=sData.data;
sData.data(1:delay,:)=NaN;
sData.data(1:len,n_th_comp)=sData.data(delay+1:len+delay,n_th_comp);
sData.data(len+1:len+delay,n_th_comp)=NaN;
end
if delay ~= 0
data=get(pre_h,'UserData');
sData.MODIFIED=1;
sData.comp_norm(:,:)={[]};
data.sData=sData;
data.undo.sData=undo;
data.undo.selected=data.selected_vects;
set(pre_h,'UserData',data);
old=gcf;
set(0,'CurrentFigure',pre_h);
write_sD_stats;
pro_tools('selall');
set(0,'CurrentFigure',old);
end
end
if ~LOG
data=get(pre_h,'UserData');
data.LOG{length(data.LOG)+1}='% Delay a component.';
data.LOG{length(data.LOG)+1}=cat(2,'preprocess(''delay_data'',''',...
num2str(delay),' ',MODE,''');');
set(pre_h,'UserData',data);
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function window(varargin)
if nargin ~= 1
pre_h=findobj(get(0,'Children'),'Tag','Preprocess');
preh_udata=get(pre_h,'UserData');
preh_udata.LOG{length(preh_udata.LOG)+1}=...
'% Starting the ''Windowed'' -window...';
preh_udata.LOG{length(preh_udata.LOG)+1}='preprocess(''window'',''foo'');';
set(pre_h,'UserData',preh_udata);
end
win_h=findobj(get(0,'Children'),'Tag','Window');
if ~isempty(win_h)
figure(win_h);
return;
end
h0 = figure('Color',[0.8 0.8 0.8], ...
'PaperPosition',[18 180 576 432], ...
'PaperUnits','points', ...
'Position',[513 703 288 219], ...
'Tag','Window');
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[0.752941176470588 0.752941176470588 0.752941176470588], ...
'ListboxTop',0, ...
'Position',[0.03125 0.1552511415525114 0.9375 0.7990867579908676], ...
'Style','frame', ...
'Tag','Frame1');
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[0.8 0.8 0.8], ...
'ListboxTop',0, ...
'Position',[0.04861111111111111 0.7214611872146118 0.9027777777777777 0.2009132420091324], ...
'Style','frame', ...
'Tag','Frame2');
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[0.8 0.8 0.8], ...
'ListboxTop',0, ...
'Position',[0.04861111111111111 0.1780821917808219 0.2777777777777778 0.5251141552511416], ...
'Style','frame', ...
'Tag','Frame3');
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[0.8 0.8 0.8], ...
'ListboxTop',0, ...
'Position',[0.3611111111111111 0.1780821917808219 0.2777777777777778 0.5251141552511416], ...
'Style','frame', ...
'Tag','Frame4');
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[0.8 0.8 0.8], ...
'ListboxTop',0, ...
'Position',[0.6736111111111111 0.1780821917808219 0.2777777777777778 0.5251141552511416], ...
'Style','frame', ...
'Tag','Frame5');
h1 = uicontrol('Parent',h0, ...
'Callback','preprocess eval_windowed',...
'Units','normalized', ...
'FontWeight','demi', ...
'ListboxTop',0, ...
'Position',[0.03125 0.0319634703196347 0.2256944444444444 0.091324200913242], ...
'String','OK', ...
'Tag','Pushbutton1');
h1 = uicontrol('Parent',h0, ...
'Callback','preprocess close_w', ...
'Units','normalized', ...
'FontWeight','demi', ...
'ListboxTop',0, ...
'Position',[0.7430555555555555 0.0319634703196347 0.2256944444444444 0.091324200913242], ...
'String','Close', ...
'Tag','Pushbutton2');
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[1 1 1], ...
'ListboxTop',0, ...
'Position',[0.7083333333333333 0.7625570776255708 0.2083333333333333 0.1232876712328767], ...
'Style','edit', ...
'Tag','EditText1');
data.win_len_h=h1;
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[0.8 0.8 0.8], ...
'FontWeight','demi', ...
'HorizontalAlignment','left', ...
'ListboxTop',0, ...
'Position',[0.07638888888888888 0.8036529680365296 0.3784722222222222 0.0547945205479452], ...
'String','Window length', ...
'Style','text', ...
'Tag','StaticText1');
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[0.8 0.8 0.8], ...
'Callback','preprocess window_cb centered',...
'ListboxTop',0, ...
'Position',[0.06597222222222222 0.5616438356164384 0.06597222222222222 0.0821917808219178], ...
'Style','radiobutton', ...
'Tag','Radiobutton1');
data.centered_h=h1;
data.position='center';
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[0.8 0.8 0.8], ...
'Callback','preprocess window_cb previous',...
'ListboxTop',0, ...
'Position',[0.06597222222222222 0.4018264840182648 0.06597222222222222 0.0821917808219178], ...
'Style','radiobutton', ...
'Tag','Radiobutton2');
data.previous_h=h1;
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[0.8 0.8 0.8], ...
'Callback','preprocess window_cb next',...
'ListboxTop',0, ...
'Position',[0.06597222222222222 0.2465753424657534 0.06597222222222222 0.0821917808219178], ...
'Style','radiobutton', ...
'Tag','Radiobutton3');
data.next_h=h1;
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[0.8 0.8 0.8], ...
'Callback','preprocess window_cb mean',...
'ListboxTop',0, ...
'Position',[0.3784722222222222 0.5799086757990868 0.06597222222222222 0.0821917808219178], ...
'Style','radiobutton', ...
'Tag','Radiobutton4');
data.mean_h=h1;
data.mode='mean';
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[0.8 0.8 0.8], ...
'Callback','preprocess window_cb median',...
'ListboxTop',0, ...
'Position',[0.3784722222222222 0.4611872146118721 0.06597222222222222 0.0821917808219178], ...
'Style','radiobutton', ...
'Tag','Radiobutton5');
data.median_h=h1;
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[0.8 0.8 0.8], ...
'Callback','preprocess window_cb max',...
'ListboxTop',0, ...
'Position',[0.3784722222222222 0.3515981735159817 0.06597222222222222 0.0821917808219178], ...
'Style','radiobutton', ...
'Tag','Radiobutton6');
data.max_h=h1;
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'Callback','preprocess window_cb min',...
'BackgroundColor',[0.8 0.8 0.8], ...
'ListboxTop',0, ...
'Position',[0.3784722222222222 0.2374429223744292 0.06597222222222222 0.0821917808219178], ...
'Style','radiobutton', ...
'Tag','Radiobutton7');
data.min_h = h1;
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[0.8 0.8 0.8], ...
'Callback','preprocess window_cb clip',...
'ListboxTop',0, ...
'Position',[0.6909722222222222 0.5525114155251141 0.06597222222222222 0.0821917808219178], ...
'Style','radiobutton', ...
'Tag','Radiobutton8');
data.clip_h=h1;
data.eval_mode='clip';
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[0.8 0.8 0.8], ...
'Callback','preprocess window_cb expand',...
'ListboxTop',0, ...
'Position',[0.6909722222222222 0.2922374429223744 0.06597222222222222 0.0821917808219178], ...
'Style','radiobutton', ...
'Tag','Radiobutton9');
data.expand_h=h1;
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[0.8 0.8 0.8], ...
'FontWeight','demi', ...
'FontSize',8,...
'HorizontalAlignment','left', ...
'ListboxTop',0, ...
'Position',[0.132 0.5799 0.19 0.0548], ...
'String','Centered', ...
'Style','text', ...
'Tag','StaticText2');
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[0.8 0.8 0.8], ...
'FontWeight','demi', ...
'FontSize',8,...
'HorizontalAlignment','left', ...
'ListboxTop',0, ...
'Position',[0.132 0.4247 0.1667 0.0548], ...
'String','Previous', ...
'Style','text', ...
'Tag','StaticText3');
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[0.8 0.8 0.8], ...
'FontWeight','demi', ...
'FontSize',8,...
'HorizontalAlignment','left', ...
'ListboxTop',0, ...
'Position',[0.132 0.2648 0.1632 0.0548], ...
'String','Next', ...
'Style','text', ...
'Tag','StaticText4');
h1 = uicontrol('Parent',h0, ...,
'Units','normalized', ...
'BackgroundColor',[0.8 0.8 0.8], ...
'FontWeight','demi', ...
'FontSize',8,...
'HorizontalAlignment','left', ...
'ListboxTop',0, ...
'Position',[0.445 0.6027397260273972 0.19 0.0547945205479452], ...
'String','Mean', ...
'Style','text', ...
'Tag','StaticText5');
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[0.8 0.8 0.8], ...
'FontWeight','demi', ...
'FontSize',8,...
'HorizontalAlignment','left', ...
'ListboxTop',0, ...
'Position',[0.445 0.4795 0.1806 0.0548], ...
'String','Median', ...
'Style','text', ...
'Tag','StaticText6');
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[0.8 0.8 0.8], ...
'FontWeight','demi', ...
'FontSize',8,...
'HorizontalAlignment','left', ...
'ListboxTop',0, ...
'Position',[0.445 0.3699 0.1667 0.0548], ...
'String','Max', ...
'Style','text', ...
'Tag','StaticText7');
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[0.8 0.8 0.8], ...
'FontWeight','demi', ...
'FontSize',8,...
'HorizontalAlignment','left', ...
'ListboxTop',0, ...
'Position',[0.445 0.2557077625570776 0.1597222222222222 0.0547945205479452], ...
'String','Min', ...
'Style','text', ...
'Tag','StaticText8');
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[0.8 0.8 0.8], ...
'FontWeight','demi', ...
'FontSize',8,...
'HorizontalAlignment','left', ...
'ListboxTop',0, ...
'Position',[0.7535 0.5753 0.1354 0.054], ...
'String','Clip', ...
'Style','text', ...
'Tag','StaticText9');
h1 = uicontrol('Parent',h0, ...
'Units','normalized', ...
'BackgroundColor',[0.8 0.8 0.8], ...
'FontWeight','demi', ...
'FontSize',8,...
'HorizontalAlignment','left', ...
'ListboxTop',0, ...
'Position',[0.7534722222222222 0.3150684931506849 0.1527777777777778 0.0547945205479452], ...
'String','Expand', ...
'Style','text', ...
'Tag','StaticText10');
set(gcf,'UserData',data);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function window_cb(arg)
data=get(gcf,'UserData');
if any(strcmp(arg,[{'centered'},{'previous'},{'next'}]))
switch arg
case 'centered'
data.position='center';
set(data.previous_h,'Value',0);
set(data.next_h,'Value',0);
set(data.centered_h,'Value',1);
case 'previous'
data.position='previous';
set(data.centered_h,'Value',0);
set(data.next_h,'Value',0);
set(data.previous_h,'Value',1);
case 'next'
data.position='next';
set(data.centered_h,'Value',0);
set(data.previous_h,'Value',0);
set(data.next_h,'Value',1);
end
elseif any(strcmp(arg,[{'mean'},{'median'},{'min'},{'max'}]))
switch arg
case 'mean'
data.mode='mean';
set(data.median_h,'Value',0);
set(data.min_h,'Value',0);
set(data.max_h,'Value',0);
set(data.mean_h,'Value',1);
case 'median'
data.mode='median';
set(data.mean_h,'Value',0);
set(data.max_h,'Value',0);
set(data.min_h,'Value',0);
set(data.median_h,'Value',1);
case 'max'
data.mode='max';
set(data.mean_h,'Value',0);
set(data.median_h,'Value',0);
set(data.min_h,'Value',0);
set(data.max_h,'Value',1);
case 'min'
data.mode='min';
set(data.mean_h,'Value',0);
set(data.median_h,'Value',0);
set(data.max_h,'Value',0);
set(data.min_h,'Value',1);
end
elseif any(strcmp(arg,[{'clip','expand'}]))
switch arg
case 'clip'
data.eval_mode='clip';
set(data.expand_h,'Value',0);
set(data.clip_h,'Value',1);
case 'expand'
data.eval_mode='expand';
set(data.clip_h,'Value',0);
set(data.expand_h,'Value',1);
end
end
set(gcf,'UserData',data);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function eval_windowed(varargin)
if nargin == 1
LOG=1;
i=1;
arg=varargin{1};
while i <= length(arg) & arg(i) ~= ' '
i=i+1;
end
value=str2num(arg(1:i-1));
i=i+1;
j=i;
while i <= length(arg) & arg(i) ~= ' '
i=i+1;
end
position=arg(j:i-1);
i=i+1;
j=i;
while i <= length(arg) & arg(i) ~= ' '
i=i+1;
end
mode=arg(j:i-1);
i=i+1;
j=i;
while i <= length(arg) & arg(i) ~= ' '
i=i+1;
end
eval_mode=arg(j:i-1);
else
LOG=0;
end
data=get(gcf,'UserData');
if LOG
data.position=position;
data.eval_mode=eval_mode;
data.mode=mode;
end
pre_h=findobj(get(0,'Children'),'Tag','Preprocess');
if isempty(pre_h)
errordlg('''Preprocess''-window does not exist. Terminating program...');
pro_tools('close');
return;
end
comp_names_h=getfield(get(pre_h,'UserData'),'comp_names_h');
sData=getfield(get(pre_h,'UserData'),'sData');
undo=sData;
if isempty(get(comp_names_h,'Value'))
errordlg('There are not components chosen.');
return;
end
if ~LOG
if isempty(get(data.win_len_h,'String'))
errordlg('Window length must be defined');
return;
end
value=str2num(get(data.win_len_h,'String'));
end
set(data.win_len_h,'String','');
if ~LOG
if isempty(value) | value < 0 | value ~= round(value)
errordlg('Window length must be positive integer.');
return;
end
if value > length(sData.data(:,1))
errordlg('Length of window is too big.');
return;
end
end
index=getfield(get_indices,{1});
sData=eval_operation(sData,value,data.mode,data.eval_mode,data.position,index);
sData.comp_norm(index)={[]};
u_data=get(pre_h,'UserData');
u_data.sData=sData;
u_data.undo.sData=undo;
u_data.undo.selected=u_data.selected_vects;
if ~LOG
u_data.LOG{length(u_data.LOG)+1}=...
'% Evaluating the wanted ''windowed'' -operation.';
u_data.LOG{length(u_data.LOG)+1}=cat(2,'preprocess(''eval_windowed'',',...
'''',num2str(value),...
' ',data.position,' ',data.mode,...
' ',data.eval_mode,''');');
end
set(pre_h,'UserData',u_data);
old=gcf;
set(0,'CurrentFigure',pre_h);
if strcmp(data.eval_mode,'expand');
write_sD_stats;
pro_tools('selall');
else
sel_comp;
cplot_mimema;
end
set(0,'CurrentFigure',old);
%%% Subfunction: eval_operation %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function sData=eval_operation(sData,winlen,mode,evalmode,position,n)
len=length(sData.data(:,1));
dim=length(sData.data(1,:));
switch(position)
case 'center'
prev=round(winlen/2)-1;
next=winlen-round(winlen/2);
case 'previous'
prev=winlen-1;
next=0;
case 'next'
prev=0;
next=winlen-1;
end
switch(evalmode)
case 'clip'
for center=1:len
win=center-prev:center-prev+winlen-1;
win=win(find(win > 0 & win <= len));
str=cat(2,mode,'(sData.data(win(find(~isnan(sData.data(win,n)))),n))');
tmp(center)=eval(str);
end
sData.data(:,n)=tmp;
case 'expand'
for i=1:len+winlen-1
win=i-(winlen-1):i;
win=win(find(win > 0 & win <= len));
str=cat(2,mode,'(sData.data(win(find(~isnan(sData.data(win,n)))),n))');
tmp(i)=eval(str);
end
sData.data=cat(1,repmat(NaN,next,dim),sData.data,repmat(NaN,prev,dim));
sData.data(:,n)=tmp;
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function pro_tools(arg)
switch arg
case 'close'
close_preprocess;
case 'c_stat'
write_c_stats;
case 'plot_hist'
plot_hist;
case 'plot'
plot_button;
case 'plxy'
plxy_button;
case 'bplo'
bplo_button;
case 'hist'
hist_button;
end
%%% Subfunction close_preprocess %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function close_preprocess
pre_h=findobj(get(0,'Children'),'Tag','Preprocess');
man_h=findobj(get(0,'Children'),'Tag','Management');
clip_h=findobj(get(0,'Children'),'Tag','Clipping');
plot_h=findobj(get(0,'Children'),'Tag','PlotWin');
delay_h=findobj(get(0,'Children'),'Tag','Delay');
window_h=findobj(get(0,'Children'),'Tag','Window');
sel_h=findobj(get(0,'Children'),'Tag','Select');
if ~isempty(man_h)
close(man_h);
end
if ~isempty(clip_h)
close(clip_h);
end
if ~isempty(plot_h)
close(plot_h);
end
if ~isempty(delay_h)
close(delay_h);
end
if ~isempty(window_h)
close(window_h);
end
if ~isempty(sel_h)
close(sel_h);
end
if ~isempty(pre_h)
close(pre_h);
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%% Subfunction: undo %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function undo(varargin)
if nargin == 1
pre_h=findobj(get(0,'Children'),'Tag','Preprocess');
set(0,'CurrentFigure',pre_h);
LOG=1;
else
LOG=0;
end
data=get(gcf,'UserData');
if ~isempty(data.undo)
if any(strcmp('selected',fieldnames(data.undo)))
data.selected_vects=data.undo.selected;
end
if ~any(strcmp('index',fieldnames(data.undo)))
data.sData=data.undo.sData;
data.undo=[];
if ~LOG
data.LOG{length(data.LOG)+1}='% Undo the most recent operation.';
data.LOG{length(data.LOG)+1}='preprocess(''undo'',''foo'');';
end
set(gcf,'UserData',data);
set_compnames(data.sData,data.comp_names_h);
write_sD_stats;
vect_means(data.sData,data.vect_mean_h,data.selected_vects);
sel_comp;
cplot_mimema;
return;
end
% 'undo.sData' does not exist in sD_set - array
index=data.undo.index;
data.undo.sData=rmfield(data.undo.sData,[{'INDEX'};{'MODIFIED'}]);
if index<=length(data.sD_set)
rest=data.sD_set(index:length(data.sD_set));
else
rest=[];
end
data.sD_set=cat(2,data.sD_set(1:index-1),data.undo.sData,rest);
data.undo=[];
if ~LOG
data.LOG{length(data.LOG)+1}='% Undo the most recent operation.';
data.LOG{length(data.LOG)+1}='preprocess(''undo'',''foo'');';
end
set(gcf,'UserData',data);
set(getfield(get(gcf,'UserData'),'sD_set_h'),'Value',index);
set_sD_stats;
sel_sD;
else
msgbox('Can''t do...');
end
%%% Subfunction: write_c_stats %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function write_c_stats(varargin)
pre_h=findobj(get(0,'Children'),'Tag','Preprocess');
comp_names_h=getfield(get(pre_h,'UserData'),'comp_names_h');
sel_comp_h=getfield(get(pre_h,'UserData'),'sel_comp_h');
sel_chist_h=getfield(get(pre_h,'UserData'),'sel_chist_h');
if nargin==1
val1=varargin(1);
else
val1=get(sel_comp_h,'String');
end
if ~isempty(val1) & iscell(val1)
selected_vects=getfield(get(pre_h,'UserData'),'selected_vects');
sData=getfield(get(pre_h,'UserData'),'sData');
sel_cdata_h=getfield(get(pre_h,'UserData'),'sel_cdata_h');
name=getfield(get(sel_comp_h,'String'),{get(sel_comp_h,'Value')});
name=name{1};
i=2;
while ~isempty(str2num(name(i)))
value(i-1)=name(i);
i=i+1;
end
value=str2num(value);
data=sData.data(selected_vects,value);
string{1} = cat(2,'Min: ',sprintf('%-10.3g',min(data)));
string{2} = cat(2,'Mean: ',sprintf('%-10.3g',mean(data(find(~isnan(data))))));
string{3} = cat(2,'Max: ',sprintf('%-10.3g',max(data)));
string{4} = cat(2,'Std: ',sprintf('%-10.3g',std(data(find(~isnan(data))))));
string{5} = cat(2,'Number of NaNs: ',sprintf('%-10.3g',sum(isnan(data))));
string{6} = cat(2,'NaN (%):',...
sprintf('%-10.3g',100*sum(isnan(data))/length(data)));
string{7} = cat(2,'Number of values: ',sprintf('%-10.3g',...
length(find(~isnan(unique(data))))));
set(sel_cdata_h,'String',string);
set(sel_cdata_h,'HorizontalAlignment','left');
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%% Subfunction plot_hist %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function plot_hist
pre_h=findobj(get(0,'Children'),'Tag','Preprocess');
sel_chist_h=getfield(get(pre_h,'UserData'),'sel_chist_h');
sData=getfield(get(pre_h,'UserData'),'sData');
selected=getfield(get(pre_h,'UserData'),'selected_vects');
value=get(getfield(get(pre_h,'UserData'),'sel_comp_h'),'Value');
subplot(sel_chist_h);
hold off;
cla;
if all(isnan(sData.data(:,value)));
return;
end
hold on;
lim1=min(sData.data(:,value));
lim2=max(sData.data(:,value));
if lim2 - lim1 >= eps
x=lim1:(lim2-lim1)/(30-1):lim2;
set(sel_chist_h,'XLim',[lim1 lim2]);
elseif lim1 ~= 0
x=(lim1)/2:lim1/(30-1):lim1+(lim1)/2;
set(sel_chist_h,'Xlim',[lim1-abs(lim1/2) lim1+abs(lim1/2)]);
else
x=-1:2/(30-1):1;
set(sel_chist_h,'XLim',[-1 1]);
end
hist(sData.data(selected,value),x);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%% Subfunction: select_all %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function select_all(varargin)
if nargin == 1
LOG=1;
else
LOG=0;
end
data=get(gcf,'UserData');
data.selected_vects=(1:length(data.sData.data(:,1)));
if ~LOG
data.LOG{length(data.LOG)+1}='% Select all vectors.';
data.LOG{length(data.LOG)+1}='selall(''foo'');';
end
set(gcf,'UserData',data);
tmp=zeros(1,length(data.sData.data(:,1)));
tmp(data.selected_vects)=1;
draw_vectors(tmp,data.vector_h);
cplot_mimema;
vect_means(data.sData,data.vect_mean_h,data.selected_vects);
sel_comp;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%% Subfunction: plot_button %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function plot_button
%PLOT_BUTTON A callback function. Plots all the components and marks
% the chosen components.
%
%
sData=getfield(get(gcf,'UserData'),'sData');
selected=getfield(get(gcf,'UserData'),'selected_vects');
indices=get_indices;
if isempty(indices)
return;
end
h=findobj(get(0,'Children'),'Tag','PlotWin');
if isempty(h)
h= figure;
set(h,'Tag','PlotWin');
end
names=sData.comp_names(indices);
data=sData.data(:,indices);
set(0,'CurrentFigure',h);
hold off;
clf;
t=0:1/(getfield(size(data),{1})-1):1;
tmp=setdiff(1:length(data(:,1)),selected);
for i=1:length(names)
subplot(length(names),1,i);
hold on;
if max(data(:,i))- min(data(:,i)) <= eps
set(gca,'YLim',[max(data(:,i))-1 max(data(:,i))+1]);
end
plot(t,data(:,i));
if ~isempty(tmp);
data(tmp,i)=NaN;
end
plot(t,data(:,i),'red');
ylabel(names{i});
set(gca,'XTick',[]);
end
set(gcf,'Name','Plotted Data Components');
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%% Subfunction: plxy_button %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function plxy_button
%PLXY_BUTTON A callback function. XY-plots the first and the second
% components chosen.
%
%
sData=getfield(get(gcf,'UserData'),'sData');
selected=getfield(get(gcf,'UserData'),'selected_vects');
inds = get_indices;
if length(inds) < 2
errordlg('There must be two components chosen for XY-plot.');
return;
end
inds=inds(1:2);
names=getfield(sData,'comp_names',{inds});
h=findobj(get(0,'Children'),'Tag','PlotWin');
if isempty(h)
h= figure;
set(h,'Tag','PlotWin');
end
set(0,'CurrentFigure',h);
clf;
axes;
if max(sData.data(:,inds(1))) - min(sData.data(:,inds(1))) <= eps
set(gca,'XLim',[max(sData.data(:,inds(1)))-1 max(sData.data(:,inds(1)))+1]);
end
if max(sData.data(:,inds(2))) - min(sData.data(:,inds(2))) <= eps
set(gca,'YLim',[max(sData.data(:,inds(2)))-1 max(sData.data(:,inds(2)))+1]);
end
hold on;
plot(sData.data(:,inds(1)),sData.data(:,inds(2)),'o');
x=sData.data(selected,inds(1));
y=sData.data(selected,inds(2));
plot(x,y,'ored','MarkerSize',4);
xlabel(names(1));
ylabel(names(2));
set(h,'Name','Plotted Data Components');
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%% Sub_function: bplo_button %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function bplo_button
%BPLO_BUTTON A callback function. Box-plots the first component chosen.
sData=getfield(get(gcf,'UserData'),'sData');
selected=getfield(get(gcf,'UserData'),'selected_vects');
if length(selected) == 1
errordlg('There are too few vectors chosen for box-plotting.');
else
indices=get_indices;
if isempty(indices)
return;
end
for i=1:length(indices)
if length(unique(sData.data(selected,indices(i))))==1
errordlg('All the values are the same. Operation can''t be evaluated.');
return;
end
end
names=getfield(sData,'comp_names',{indices});
h= findobj(get(0,'Children'),'Tag','PlotWin');
if isempty(h)
h= figure;
set(h,'Tag','PlotWin');
end
data=sData.data(selected,indices);
set(0,'CurrentFigure',h);
hold off;
clf;
hold on;
for i=1:getfield(size(data),{2})
subplot(getfield(size(data),{2}),1,i);
if ~all(isnan(data(:,i)))
boxplot(data(:,i));
end
name=names{i};
tmp=get(get(gca,'YLabel'),'String');
ylabel(cat(2,sprintf('[%s] ',name),tmp));
end
set(h,'Name','Box-plot');
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%% Subfunction: hist_button %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function hist_button
no_of_bins_h=getfield(get(gcf,'UserData'),'no_of_bins_h');
selected=getfield(get(gcf,'UserData'),'selected_vects');
sData=getfield(get(gcf,'UserData'),'sData');
n=str2num(get(no_of_bins_h,'String'));
s1='Invalid number of bins.';
s2=sprintf('\nSet new value to the box under the ''Histogram''-button.');
if isempty(n)
errordlg(cat(2,s1,s2));
else
indices=get_indices;
if isempty(indices)
return;
end
n=round(n);
if n < 1
errordlg('Number of bins must be positive integer.');
else
h= findobj(get(0,'Children'),'Tag','PlotWin');
if isempty(h)
h= figure;
set(h,'Tag','PlotWin');
end
set(0,'CurrentFigure',h);
hold off;
clf;
data=sData.data(selected,indices);
names=sData.comp_names(indices);
for i=1:length(names)
subplot(length(names),1,i);
hold on;
lim1=min(sData.data(:,indices(i)));
lim2=max(sData.data(:,indices(i)));
if n > 1
if lim2 - lim1 >= eps
x=lim1:(lim2-lim1)/(n-1):lim2;
set(gca,'XLim',[lim1 lim2]);
elseif lim1 ~= 0
x=lim1/2:lim1/(n-1):lim1/2+lim1;
if ~all(isnan([lim1 lim2]))
set(gca,'XLim',[lim1-abs(lim1/2) lim1+abs(lim1/2)]);
end
else
x=-1:2/(n-1):1;
set(gca,'XLim',[-1 1]);
end
else
x=1;
if lim2 ~= lim1
set(gca,'XLim',[lim1 lim2]);
else
set(gca,'XLim',[lim1/2 lim1/2+lim1]);
end
end
if ~all(isnan(data(:,i)))
hist(data(:,i),x);
end
name=names{i};
xlabel(name);
end
set(h,'Name','Histogram');
end
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%% Subfunction: no_of_values %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function no_of_values(varargin);
%NO_OF_VALUES A callback function. Calculates the number of different
% values of the chosen components.
%
%
if nargin==1;
LOG=1;
else
LOG=0;
end
pre_h=findobj(get(0,'Children'),'Tag','Preprocess');
results_h=getfield(get(pre_h,'UserData'),'results_h');
sData=getfield(get(pre_h,'UserData'),'sData');
selected=getfield(get(pre_h,'UserData'),'selected_vects');
str1='There must be one component chosen for ''Number of Values''-operation';
if ~LOG & isempty(get_indices)
errordlg(str1);
else
indices=get_indices;
data=sData.data(selected,indices);
string{1} = 'Number of different values:';
for i=1:getfield(size(data),{2})
tmp=data(:,i);
string{i+1}=cat(2,sprintf('#%d:',indices(i)),...
sprintf('%d',length(find(~isnan(unique(data(:,i)))))));
end
set(results_h,'String',string);
set(results_h,'HorizontalAlignment','left');
if ~LOG
data=get(pre_h,'UserData');
data.LOG{length(data.LOG)+1}='% Number of values';
data.LOG{length(data.LOG)+1}='preprocess(''noof'',''foo'');';
set(pre_h,'UserData',data);
end
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%% Subfunction: correlation %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function correlation(varargin)
if nargin == 1
LOG=1;
else
LOG=0;
end
pre_h=findobj(get(0,'Children'),'Tag','Preprocess');
results_h=getfield(get(pre_h,'UserData'),'results_h');
selected=getfield(get(pre_h,'UserData'),'selected_vects');
sData=getfield(get(pre_h,'UserData'),'sData');
if length(get_indices) < 2
errordlg('There must be two components chosen for Correlation');
else
indices=getfield(get_indices,{1:2});
data=sData.data(selected,indices);
inds=find(~isnan(data(:,1)) & ~isnan(data(:,2)));
value=getfield(corrcoef(data(inds,1),data(inds,2)),{1,2});
names=sData.comp_names(indices);
string{1}='Correlation between';
string{2}=cat(2,names{1},' and ',names{2},':');
string{3}=sprintf('%-10.3g',value);
set(results_h,'String',string);
set(results_h,'HorizontalAlignment','left');
if ~LOG
data=get(pre_h,'UserData');
data.LOG{length(data.LOG)+1}='% Correlation';
data.LOG{length(data.LOG)+1}='preprocess(''corr'',''foo'');';
set(pre_h,'UserData',data);
end
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%% Subfunction: unit_length %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function unit_length(varargin)
%UNIT_LENGTH A callback function Scales all the vectors to the unit
% length.
%
%
if nargin==1
LOG=1;
else
LOG=0;
end
vect_mean_h=getfield(get(gcf,'UserData'),'vect_mean_h');
sData=getfield(get(gcf,'UserData'),'sData');
sData.MODIFIED=1;
scaled=sData.data;
comp_names_h=getfield(get(gcf,'UserData'),'comp_names_h');
if ~LOG & isempty(get(comp_names_h,'Value'))
errordlg('There must be components chosen for the ''unit length''- operation');
return;
end
inds=get_indices;
for i=1:length(scaled(:,1));
x=find(~isnan(scaled(i,inds)));
scaled(i,inds(x))=(1/sqrt(sum(scaled(i,inds(x)).^2)))*scaled(i,inds(x));
end
data=get(gcf,'UserData');
data.undo.sData = sData;
data.sData.data=scaled;
for i=1:length(inds)
data.sData.comp_norm{inds(i)}=[];
end
if ~LOG
data.LOG{length(data.LOG)+1}='% Unit length';
data.LOG{length(data.LOG)+1}='preprocess(''unit'',''foo'');';
end
set(gcf,'UserData',data);
vects=zeros(1,length(sData.data(:,1)));
vects(data.selected_vects)=1;
draw_vectors(vects,data.vector_h);
vect_means(sData,vect_mean_h,data.selected_vects);
cplot_mimema;
plot_hist;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%% Subfunction: one_of_n %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function one_of_n(varargin)
if nargin==1
LOG=1;
else
LOG=0;
end
data=get(gcf,'UserData');
vector_h=getfield(get(gcf,'Userdata'),'vector_h');
comp_names_h=getfield(get(gcf,'Userdata'),'comp_names_h');
vect_mean_h=getfield(get(gcf,'UserData'),'vect_mean_h');
sData=data.sData;
undo=data.sData;
selected=getfield(get(gcf,'UserData'),'selected_vects');
msg='Creating over 10 new components. Stop operation?';
if ~LOG
if isempty(get(data.comp_names_h,'Value'))
errordlg('There must be one component chosen for ''Add: N binary types'' -operation');
return;
end
end
index=getfield(get_indices,{1});
tmp=unique(sData.data(:,index));
n=length(tmp);
if ~LOG
if n>10
answer=questdlg(msg,'Question','Yes','No','Yes');
if strcmp(answer,'Yes')
msgbox('Operation stopped.');
return;
end
end
end
dim1=getfield(size(sData.data),{1});
dim2=getfield(size(sData.data),{2});
sData.data=cat(2,sData.data,zeros(dim1,n));
dim=dim2+n;
for i=1:n
sData.data(:,dim-(n-i))=(sData.data(:,index) == tmp(i));
end
INDEX=sData.INDEX;
for i=1:n
sData.comp_names{dim2+i}=sprintf('%dNewVar',dim2+i);
end
tmp_norm=cat(1,sData.comp_norm,cell(n,1));
sData=som_data_struct(sData.data,...
'name',sData.name,...
'labels',sData.labels,...
'comp_names',sData.comp_names);
sData.MODIFIED=1;
sData.INDEX=INDEX;
sData.comp_norm=tmp_norm;
data.undo.sData=undo;
data.sData=sData;
data.selected_vects=1:length(sData.data(:,1));
if ~LOG
data.LOG{length(data.LOG)+1}='% Add: N binary types';
data.LOG{length(data.LOG)+1}='preprocess(''oneo'',''foo'');';
end
set(gcf,'UserData',data);
clear_button;
write_sD_stats;
set_compnames(sData,comp_names_h);
tmp=ones(1,length(sData.data(:,1)));
draw_vectors(tmp,vector_h);
vect_means(sData,vect_mean_h,1:length(sData.data(:,1)));
cplot_mimema;
sel_comp;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%% Subfunction: add_zeros %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function add_zeros(varargin)
if nargin == 1
LOG=1;
else
LOG=0;
end
data=get(gcf,'UserData');
vector_h=getfield(get(gcf,'Userdata'),'vector_h');
comp_names_h=getfield(get(gcf,'Userdata'),'comp_names_h');
vect_mean_h=getfield(get(gcf,'UserData'),'vect_mean_h');
sData=data.sData;
undo=sData;
dim1=getfield(size(sData.data),{1});
dim2=getfield(size(sData.data),{2});
sData.data=cat(2,sData.data,zeros(dim1,1));
INDEX=sData.INDEX;
sData.comp_names{dim2+1}=sprintf('%dNewVar',dim2+1);
tmp_norm=cat(1,sData.comp_norm,cell(1,1));
sData=som_data_struct(sData.data,...
'name',sData.name,...
'labels',sData.labels,...
'comp_names',sData.comp_names);
sData.MODIFIED=1;
sData.INDEX=INDEX;
sData.comp_norm=tmp_norm;
data.sData=sData;
data.undo.sData=undo;
data.selected_vects=1:length(sData.data(:,1));
if ~LOG
data.LOG{length(data.LOG)+1}='% Add: zeros';
data.LOG{length(data.LOG)+1}='preprocess(''zero'',''foo'');';
end
set(gcf,'UserData',data);
clear_button;
write_sD_stats;
set_compnames(sData,comp_names_h);
tmp=ones(1,length(sData.data(:,1)));
draw_vectors(tmp,vector_h);
vect_means(sData,vect_mean_h,1:length(sData.data(:,1)));
cplot_mimema;
sel_comp;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%% Subfunction: move_component %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function move_component(varargin)
%MOVE_COMPONENT A callback function. Moves one component of vectors to
% the position wanted.
%
%
if nargin == 1
LOG=1;
i=1;
while varargin{1}(i) ~= ' '
value(i)=varargin{1}(i);
i=i+1;
end
value=str2num(value); % the new place
index=str2num(varargin{1}(i:length(varargin{1}))); % index of the chosen
% component
else
LOG=0;
end
data=get(gcf,'UserData');
sData=data.sData;
undo=sData;
prompt='Enter the number of the new component place:';
if isempty(get(data.comp_names_h,'Value'))
errordlg('There must be one component chosen for ''Move Component''-operation');
return;
end
if ~LOG
index=getfield(get_indices,{1});
answer=inputdlg(prompt);
if isempty(answer) | (iscell(answer) & isempty(answer{1}))
msgbox('No components moved');
return;
end
value=str2num(answer{1});
dims=size(value);
if dims(1) ~= 1 | dims(2) ~= 1 | ~isreal(value)
errordlg('The new component place must be positive integer.')
return;
end
if value <= 0 | round(value) ~= value
errordlg('The new component place must be positive integer.');
return;
end
if value > getfield(size(sData.data),{2})
errordlg('Too big value for the new component place.');
return;
end
end
sData.MODIFIED=1;
if index < value
indices1=setdiff(1:value,index);
indices2=setdiff(value+1:length(sData.data(1,:)),index);
elseif index > value
indices1=setdiff(1:value-1,index);
indices2=setdiff(value:length(sData.data(1,:)),index);
else
data.sData=sData;
data.undo.sData=undo;
set(gcf,'UserData',data);
return;
end
tmp1=sData.data(:,indices1);
tmp2=sData.data(:,indices2);
sData.data=cat(2,tmp1,sData.data(:,index),tmp2);
tmp1=sData.comp_names(indices1);
tmp2=sData.comp_names(indices2);
sData.comp_names=cat(1,tmp1,sData.comp_names(index),tmp2);
tmp1=sData.comp_norm(indices1);
tmp2=sData.comp_norm(indices2);
sData.comp_norm=cat(1,tmp1,sData.comp_norm(index),tmp2);
data.sData=sData;
data.undo.sData=undo;
if ~LOG
data.LOG{length(data.LOG)+1}='% Move component.';
data.LOG{length(data.LOG)+1}=sprintf('preprocess(''move'',''%s %s'');',...
num2str(value),num2str(index));
end
comp_names_h=getfield(get(gcf,'UserData'),'comp_names_h');
vect_mean_h=getfield(get(gcf,'UserData'),'vect_mean_h');
vector_h=getfield(get(gcf,'UserData'),'vector_h');
data.selected_vects=1:length(sData.data(:,1));
set(gcf,'UserData',data);
clear_button;
set_compnames(sData,comp_names_h);
draw_vectors(ones(1,length(sData.data(:,1))),vector_h);
vect_means(sData,vect_mean_h,data.selected_vects);
cplot_mimema;
sel_comp;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%% Subfunction: copy_component %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function copy_component(varargin)
%COPY_COMPONENT Copies one component of vectors to the position wanted.
%
%
if nargin == 1
LOG=1;
i=1;
while varargin{1}(i) ~= ' '
value(i)=varargin{1}(i);
i=i+1;
end
value=str2num(value); % the new place
index=str2num(varargin{1}(i:length(varargin{1}))); % index of the chosen
% component
else
LOG=0;
end
data=get(gcf,'UserData');
sData=data.sData;
undo=sData;
if ~LOG
prompt='Enter the number of the new component place:';
if isempty(get(data.comp_names_h,'Value'))
errordlg('There must be one component chosen for ''Copy Component''-operation');
return;
end
index=getfield(get_indices,{1});
answer=inputdlg(prompt);
if isempty(answer) | (iscell(answer) & isempty(answer{1}))
msgbox('No components moved');
return
end
value=str2num(answer{1});
dims=size(value);
if dims(1) ~= 1 | dims(2) ~= 1 | ~isreal(value)
errordlg('The new component place must be positive integer.')
return;
end
if value <= 0 | round(value) ~= value
errordlg('The new component place must be positive integer.');
return;
end
if value > getfield(size(sData.data),{2}) + 1
errordlg('Too big value for the new component place.');
return;
end
end
sData.MODIFIED=1;
indices1=1:value-1;
indices2=value:length(sData.data(1,:));
tmp1=sData.data(:,indices1);
tmp2=sData.data(:,indices2);
sData.data=cat(2,tmp1,sData.data(:,index),tmp2);
tmp1=sData.comp_names(indices1);
tmp2=sData.comp_names(indices2);
name=cell(1,1);
name{1}=cat(2,'Copied',sData.comp_names{index});
sData.comp_names=cat(1,tmp1,name,tmp2);
tmp1=sData.comp_norm(indices1);
tmp2=sData.comp_norm(indices2);
norm=cell(1,1);
norm{1}=sData.comp_norm{index};
sData.comp_norm=cat(1,tmp1,norm,tmp2);
data.sData=sData;
data.undo.sData=undo;
if ~LOG
data.LOG{length(data.LOG)+1}='% Copy component';
data.LOG{length(data.LOG)+1}=sprintf('preprocess(''copy'',''%s %s'');',...
num2str(value),num2str(index));
end
comp_names_h=getfield(get(gcf,'UserData'),'comp_names_h');
vect_mean_h=getfield(get(gcf,'UserData'),'vect_mean_h');
vector_h=getfield(get(gcf,'UserData'),'vector_h');
data.selected_vects=1:length(sData.data(:,1));
set(gcf,'UserData',data);
clear_button;
write_sD_stats;
set_compnames(sData,comp_names_h);
draw_vectors(ones(1,length(sData.data(:,1))),vector_h);
vect_means(sData,vect_mean_h,data.selected_vects);
cplot_mimema;
sel_comp;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%% Subfunction: remove_component %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function remove_component(varargin)
if nargin == 1
LOG=1;
value=str2num(varargin{1});
else
LOG=0;
end
data=get(gcf,'UserData');
vect_mean_h=getfield(get(gcf,'UserData'),'vect_mean_h');
vector_h=getfield(get(gcf,'UserData'),'vector_h');
comp_names_h=getfield(get(gcf,'UserData'),'comp_names_h');
sData=data.sData;
undo=sData;
prompt='Enter the number of component to be removed.';
dim=length(sData.data(1,:));
if ~LOG
answer=inputdlg(prompt);
if isempty(answer) | (iscell(answer) & isempty(answer{1}))
msgbox('Components not removed.');
return;
end
value=str2num(answer{1});
dims=size(value);
if dims(1) ~= 1 | dims(2) ~= 1 | ~isreal(value)
errordlg('Number of the component to be removed must be positive integer.')
return;
end
if value <= 0 | round(value) ~= value
errordlg('Number of the component to be removed must be positive integer.');
return;
end
if value > getfield(size(sData.data),{2})
errordlg('There are less components.');
return;
end
end
sD_set_h=getfield(get(gcf,'UserData'),'sD_set_h');
index=get(sD_set_h,'Value');
if value == 1 & getfield(size(sData.data),{2}) == 1
if length(get(sD_set_h,'String')) == 1
msgbox('No data left. Closing program...')
pro_tools('close');
return;
end
set1=data.sD_set(1:index-1);
set2=data.sD_set(index+1:length(data.sD_set));
data.sD_set=[set1 set2];
set(gcf,'UserData',data);
set_sD_stats;
sel_sD;
data=get(gcf,'UserData');
data.undo.sData=undo;
data.undo.index=index;
set(gcf,'UserData',data);
return;
end
dims=size(sData.data);
tmp_data=cat(2,sData.data(:,1:value-1),sData.data(:,value+1:dims(2)));
tmp_norm=cat(1,sData.comp_norm(1:value-1),sData.comp_norm(value+1:dims(2)));
names=cat(1,sData.comp_names(1:value-1),sData.comp_names(value+1:dims(2)));
INDEX=sData.INDEX;
comp_norm=sData.comp_norm;
sData=som_data_struct(tmp_data,...
'name',sData.name,...
'labels',sData.labels,...
'comp_names',names);
sData.comp_norm=tmp_norm;
sData.MODIFIED=1;
sData.INDEX=INDEX;
data=get(gcf,'UserData');
data.sData=sData;
data.undo.sData=undo;
data.selected_vects=1:length(sData.data(:,1));
if ~LOG
data.LOG{length(data.LOG)+1}='% Remove component';
data.LOG{length(data.LOG)+1}=sprintf('preprocess(''remove'',''%s'');',...
answer{1});
end
set(gcf,'UserData',data);
clear_button;
write_sD_stats;
set_compnames(sData,comp_names_h);
tmp=ones(1,length(sData.data(:,1)));
draw_vectors(tmp,vector_h);
vect_means(sData,vect_mean_h,1:length(sData.data(:,1)));
cplot_mimema;
sel_comp;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%% Subfunction: remove_vects %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function remove_vects(varargin)
if nargin==1
LOG=1;
tmp_str=varargin{1};
else
LOG=0;
tmp_str='_foo';
end
data=get(gcf,'UserData');
vect_mean_h=data.vect_mean_h;
vector_h=data.vector_h;
sData=data.sData;
undo=sData;
if length(data.selected_vects) == getfield(size(sData.data),{1})
if LOG
answer='Yes';
else
answer=questdlg('Do you want to delete this data set?');
end
if strcmp(answer,'No')
return;
else
index=get(data.sD_set_h,'Value');
if length(get(data.sD_set_h,'String')) == 1
msgbox('No data left. Closing program...')
pro_tools('close');
return;
end
set1=data.sD_set(1:index-1);
set2=data.sD_set(index+1:length(data.sD_set));
data.sD_set=[set1 set2];
set(gcf,'UserData',data);
set(data.sD_set_h,'Value',1);
set_sD_stats;
sel_sD;
data=get(gcf,'UserData');
data.undo.sData=undo;
data.undo.index=index;
if ~LOG
data.LOG{length(data.LOG)+1}='% Remove selected vectors';
data.LOG{length(data.LOG)+1}=cat(2,'preprocess(''remove_vects'',''',...
tmp_str,''');');
end
set(gcf,'UserData',data);
return;
end
end
tmp=sData.data(data.selected_vects,:);
if ~LOG
answer=questdlg('Do you want to save removed values to workspace?');
else
if ~strcmp(tmp_str,'_foo')
answer='Yes';
else
answer='No';
end
end
old=gcf;
if strcmp(answer,'Yes')
if ~LOG
answer=inputdlg('Give the name of the output -variable.');
else
answer={tmp_str};
end
if isvalid_var_name(answer)
assignin('base',answer{1},tmp);
disp(sprintf('Removed values are set to workspace as''%s''.',answer{1}));
tmp_str=answer{1};
end
end
set(0,'CurrentFigure',old);
sData.data(data.selected_vects,:)=[];
sData.labels(data.selected_vects,:)=[];
sData.MODIFIED=1;
data.sData=sData;
data.selected=1:length(sData.data(:,1));
data.undo.sData=undo;
if ~LOG
data.LOG{length(data.LOG)}='% Remove selected vectors';
data.LOG{length(data.LOG)+1}=cat(2,'preprocess(''remove_vects'',''',...
tmp_str,''');');
end
set(gcf,'UserData',data);
draw_vectors(ones(1,length(data.selected)),data.vector_h);
write_sD_stats;
select_all('foo');
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%% Subfunction: eval1 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function eval1(varargin);
if nargin == 1
answer=varargin
LOG=1;
else
LOG=0;
end
pre_h=findobj(get(0,'Children'),'Tag','Preprocess');
if isempty(pre_h)
errordlg('''Preprocess''-figure does not exist. Terminating program...');
pro_tools('close');
return;
end
undo=getfield(get(pre_h,'UserData'),'sData');
if ~LOG
prompt={'Enter the expression to be evaluated.',...
'Enter the inverse normalization method (optional).'};
title='Single component eval.';
answer= inputdlg(prompt,title,1);
end
if ~isempty(answer)
tmp=[];
if ~isempty(answer{1})
[tmp,method]=build_expr(answer{1},'single');
if ~isstr(tmp)
sData=getfield(get(gcf,'UserData'),'sData');
tmp='Done.';
%if ~isempty(answer{2})
% sN=som_norm_struct('eval',{method,answer{2}});
%else
% sN=som_norm_struct('eval',{method});
%end
%sN=som_set(sN,'status','done');
params={answer{1};answer{2}};
ind=getfield(get_indices,{1});
x.type='';
x.method='eval';
x.params={answer{1};answer{2}};
x.status='';
sData.comp_norm{ind}=x;
data=get(gcf,'UserData');
data.undo.sData=undo;
data.sData=sData;
if ~LOG
data.LOG{length(data.LOG)+1}='% Eval (1-comp)';
data.LOG{length(data.LOG)+1}=cat(2,'preprocess eval1 ',...
sprintf('{''%s'' ''%s''};',answer{1},answer{2}));
end
set(pre_h,'UserData',data);
end
end
set(getfield(get(pre_h,'UserData'),'results_h'),'String',tmp);
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%% Subfunction: eval2
function eval2(varargin)
if nargin == 1
answer=varargin{1};
LOG=1;
else
LOG=0;
end
undo=getfield(get(gcf,'UserData'),'sData');
pre_h=findobj(get(0,'Children'),'Tag','Preprocess');
if isempty(pre_h)
errordlg('''Preprocess''-figure does not exist. Terminating program.');
pro_tools('close');
return;
end
if ~LOG
prompt='Enter the expression to be evaluated.';
title ='Eval';
answer=inputdlg(prompt,title,1);
end
if ~isempty(answer) & ~isempty(answer{1})
str=answer{1};
[answer,foo]=build_expr(answer{1},'multiple');
if ~isstr(answer)
answer='Done.';
data=get(gcf,'UserData');
data.undo.sData=undo;
if ~LOG
data.LOG{length(data.LOG)+1}='% Eval';
data.LOG{length(data.LOG)+1}=cat(2,'preprocess(''eval2'',',...
sprintf('{''%s''});',str));
end
set(gcf,'UserData',data);
end
end
set(getfield(get(pre_h,'UserData'),'results_h'),'String',answer);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%% Subfunction: zero2one_scale %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function zero2one_scale(varargin)
if nargin == 1
LOG=1;
else
LOG=0;
end
data=get(gcf,'UserData');
sData=data.sData;
undo=sData;
INDEX=sData.INDEX;
sData=rmfield(sData,[{'INDEX'};{'MODIFIED'}]);
if isempty(get(data.comp_names_h,'Value'))
errordlg('There must be components chosen for scaling.');
return;
end
sData=som_normalize(sData,'range',get_indices);
sData.MODIFIED=1;
sData.INDEX=INDEX;
data.sData=sData;
data.undo.sData=undo;
if ~LOG
data.LOG{length(data.LOG)+1}='% Scale [0,1]';
data.LOG{length(data.LOG)+1}='preprocess(''zscale'', ''foo'');';
end
set(gcf,'UserData',data);
vects=zeros(1,length(sData.data(:,1)));
vects(data.selected_vects)=1;
cplot_mimema;
plot_hist;
vect_means(sData,data.vect_mean_h,data.selected_vects);
draw_vectors(vects,data.vector_h);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%% Subfunction: var_scale %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function var_scale(varargin)
if nargin == 1
LOG=1;
else
LOG=0;
end
data=get(gcf,'UserData');
sData=data.sData;
undo=sData;
INDEX=sData.INDEX;
sData=rmfield(sData,[{'INDEX'};{'MODIFIED'}]);
if isempty(get(data.comp_names_h,'Value'))
errordlg('There must be components chosen for scaling.');
return;
end
sData=som_normalize(sData,'var',get_indices);
sData.INDEX=INDEX;
sData.MODIFIED=1;
data.sData=sData;
data.undo.sData=undo;
if ~LOG
data.LOG{length(data.LOG)+1}='% Scale var=1';
data.LOG{length(data.LOG)+1}='preprocess(''vscale'', ''foo'');';
end
set(gcf,'UserData',data);
vects=zeros(1,length(sData.data(:,1)));
vects(data.selected_vects)=1;
cplot_mimema;
plot_hist;
vect_means(sData,data.vect_mean_h,data.selected_vects);
draw_vectors(vects,data.vector_h);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%% Subfunction: hist_eq %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function hist_eq(varargin)
if nargin == 1
LOG=1;
else
LOG=0;
end
data=get(gcf,'UserData');
sData=data.sData;
undo=sData;
INDEX=sData.INDEX;
sData=rmfield(sData,[{'INDEX'},{'MODIFIED'}]);
if isempty(get(data.comp_names_h,'Value'))
errordlg('There must be components chosen for ''Histogram eq''.');
return;
end
sData=som_normalize(sData,'histD',get_indices);
sData.INDEX=INDEX;
sData.MODIFIED=1;
data.sData=sData;
data.undo.sData=undo;
if ~LOG
data.LOG{length(data.LOG)+1}='% Histogram eq';
data.LOG{length(data.LOG)+1}='preprocess(''histeq'', ''foo'');';
end
set(gcf,'UserData',data);
vects=zeros(1,length(sData.data(:,1)));
vects(data.selected_vects)=1;
cplot_mimema;
plot_hist;
vect_means(sData,data.vect_mean_h,data.selected_vects);
draw_vectors(vects,data.vector_h);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%% Subfunction: hist_eq2 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function hist_eq2(varargin)
if nargin == 1
LOG=1;
else
LOG=0;
end
data=get(gcf,'UserData');
sData=data.sData;
undo=sData;
INDEX=sData.INDEX;
sData=rmfield(sData,[{'INDEX'};{'MODIFIED'}]);
if isempty(get(data.comp_names_h,'Value'))
errordlg('There must be components chosen for ''Histogram eq2''.');
return;
end
inds=get_indices;
%%%[sData,ok]=som_normalize(sData,inds,'histC');
sData=som_normalize(sData,'histC',inds);
sData.INDEX=INDEX;
sData.MODIFIED=1;
data.sData=sData;
data.undo.sData=undo;
if ~LOG
data.LOG{length(data.LOG)+1}='% Histogram eq2';
data.LOG{length(data.LOG)+1}='preprocess(''histeq2'', ''foo'');';
end
set(gcf,'UserData',data);
vects=zeros(1,length(sData.data(:,1)));
vects(data.selected_vects)=1;
cplot_mimema;
plot_hist;
vect_means(sData,data.vect_mean_h,data.selected_vects);
draw_vectors(vects,data.vector_h);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%% Subfunction: logarithm %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function logarithm(varargin)
if nargin == 1
LOG=1;
else
LOG=0;
end
data=get(gcf,'UserData');
sData=data.sData;
undo=sData;
INDEX=sData.INDEX;
sData=rmfield(sData,[{'INDEX'},{'MODIFIED'}]);
if isempty(get(data.comp_names_h,'Value'))
errordlg('There must be components chosen for ''Log''.');
return;
end
Data=som_normalize(sData,'log',get_indices);
sData.INDEX=INDEX;
sData.MODIFIED=1;
data.sData=sData;
data.undo.sData=undo;
if ~LOG
data.LOG{length(data.LOG)+1}='% Log';
data.LOG{length(data.LOG)+1}='preprocess(''log'', ''foo'');';
end
set(gcf,'UserData',data);
vects=zeros(1,length(sData.data(:,1)));
vects(data.selected_vects)=1;
cplot_mimema;
plot_hist;
vect_means(sData,data.vect_mean_h,data.selected_vects);
draw_vectors(vects,data.vector_h);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function [answer,method]=build_expr(string,evaltype)
pre_h=findobj(get(0,'Children'),'Tag','Preprocess');
method=[];
if isempty(pre_h)
close_preprocess;
errordlg('''Preprocess'' -figure does not exist. Terminating program...');
return;
end
if isempty(string)
str = '[]';
return;
end
tmp=[];
[name,assign,skip]=check_assign(string,evaltype);
if ~strcmp(assign,'NOTASSIGN') & ~strcmp(assign,'error')
string=string(skip:length(string));
end
if ~strcmp(assign,'error')
if isempty(string)
answer='Illegal expression.';
return;
end
[str,skip]=check_token(string,evaltype);
method=string;
while ~strcmp(str,'error') & ~strcmp(tmp,'error') & skip < length(string)
if ~strcmp(tmp,')')
str=cat(2,str,tmp);
end
[tmp,skip2]=check_token(string(skip+1:length(string)),evaltype);
skip=skip+skip2;
end
if ~strcmp(tmp,')') & ~strcmp(tmp,'error')
str=cat(2,str,tmp);
elseif strcmp(tmp,'error')
str='error';
end
end
if ~strcmp(assign,'error') & ~strcmp(str,'error');
answer=evalin('caller',str,'lasterr');
else
answer='??? Illegal expression.';
end
data=get(pre_h,'UserData');
sData=data.sData;
if strcmp(assign,'NOTASSIGN') & strcmp(evaltype,'single') & ~isstr(answer)
if isempty(get(getfield(get(pre_h,'UserData'),'comp_names_h'),'Value'))
errordlg('There are not components chosen.');
answer='??? Illegal expression.';
return;
end
index=getfield(get_indices,{1});
if strcmp(assign,'NOTASSIGN')
if length(sData.data(:,index)) ~=length(answer) & ~isscalar(answer)
answer='??? Illegal assignment.';
else
sData.data(:,index)=answer;
sData.MODIFIED=1;
data.sData=sData;
set(pre_h,'UserData',data);
end
else
if length(sData.data(str2num(assign),index)) ~=length(answer) & ~isscalar(answer)
answer='??? Illegal assignment.';
else
sData.data(str2num(assign),index)=answer;
sData.MODIFIED=1;
data.sData=sData;
set(pre_h,'UserData',data);
end
end
elseif ~strcmp(assign,'error') & ~isstr(answer) & ~strcmp(assign,'NOTASSIGN')
switch name
case 'x'
if isempty(get(data.comp_names_h,'Value'))
return;
end
index = getfield(get_indices,{1});
if isempty(assign)
if length(sData.data(:,index)) ~= length(answer) & ~isscalar(answer)
answer='??? Illegal assignment.';
else
sData.data(:,index)=answer;
sData.MODIFIED=1;
data.sData=sData;
if strcmp(evaltype,'multiple')
data.sData.comp_norm(index)={[]};
end
set(pre_h,'UserData',data);
end
else
args=create_args(assign,'x');
if length(args) == 1
len=max(str2num(args{1}));
if ~isscalar(len)
answer='??? Illegal assignment.';
return;
elseif len > length(sData.data(:,1)) | min(str2num(args{1})) < 1
answer='??? Illegal assignment.';
return;
elseif ~all(size(sData.data(str2num(args{1}),index))) == size(answer) & ~isscalar(answer)
answer='??? Illegal assignment.';
return;
else
sData.data(str2num(args{1}),index)=answer;
sData.MODIFIED=1;
data.sData=sData;
if strcmp(evaltype,'multiple')
data.sData.comp_norm(index)={[]};
end
set(pre_h,'UserData',data);
end
else
len=max(str2num(args{1}));
dim=max(str2num(args{2}));
asize=size(answer);
msize=size(sData.data);
if ~isscalar(len) | ~isscalar(dim)
answer='??? Illegal assignment.';
return;
elseif len > length(sData.data(:,1)) | len < 1
answer='??? Illegal assignment.';
return;
elseif dim > 1 | dim > msize(2) | min(str2num(args{2})) < 1
answer='??? Illegal assignment.';
return;
end
len=length(str2num(args{1}));
dim=length(str2num(args{1}));
if ~all([len dim] == asize) & ~isscalar(answer)
answer='??? Illegal assignment.';
return;
else
tmp=sData.data(:,index);
tmp([str2num(args{1})],[str2num(args{2})])=answer;
sData.data(:,index)=tmp;
sData.MODIFIED=1;
data.sData=sData;
if strcmp(evaltype,'multiple')
data.sData.comp_norm(index)={[]};
end
set(pre_h,'UserData',data);
end
end
end
case 'xs'
if isempty(get(data.comp_names_h,'Value'))
return;
end
indices=get_indices;
if isempty(assign)
if ~all(size(answer) == size(sData.data(:,indices))) & ~isscalar(answer)
answer='??? Illegal assignment.';
else
sData.data(:,indices) = answer;
sData.MODIFIED=1;
data.sData=sData;
data.sData.comp_norm(indices)={[]};
set(pre_h,'UserData',data);
end
else
args=create_args(assign,'xs');
if length(args) == 1
len=max(str2num(args{1}));
if ~isscalar(len)
answer='??? Illegal assignment.';
return;
elseif len > length(sData.data(:,1)) | min(str2num(args{1})) < 1
answer='??? Illegal assignment.';
return;
end
if ~all(size(answer) == size(sData.data(str2num(args{1})))) &...
~isscalar(answer)
answer='??? Illegal assignment.';
return;
else
tmp=sData.data(:,indices);
tmp(str2num(args{1}))=answer;
sData.data(:,indices)=tmp;
sData.MODIFIED=1;
sData.comp_norm{indices}={[]};
data.sData=sData;
set(pre_h,'UserData',data);
end
else
len=max(str2num(args{1}));
dim=max(str2num(args{2}));
asize=size(answer);
msize=size(sData.data(:,indices));
if ~isscalar(len) | ~isscalar(dim)
answer='??? Illegal assignment.';
return;
elseif len > msize(1) | min(str2num(args{1})) < 1
answer='??? Illegal assignment.';
return;
elseif dim > msize(2) | min(str2num(args{2})) < 1
answer='??? Illegal assignment.';
return;
end
len=length(str2num(args{1}));
dim=length(str2num(args{2}));
if ~all([len dim] == asize) & ~isscalar(answer)
answer='??? Illegal assignment';
return;
else
tmp=sData.data(:,indices);
tmp([str2num(args{1})],[str2num(args{2})])=answer;
sData.MODIFIED=1;
sData.data(:,indices)=tmp;
data.sData=sData;
data.sData.comp_norm(indices)={[]};
set(pre_h,'UserData',data);
end
end
end
case 'D'
if isempty(assign)
if ~all(size(answer) == size(sData.data)) & ~isscalar(answer)
answer='??? Illegal assignment.';
else
if isscalar(answer)
sData.data(:,:)=answer;
else
sData.data=answer;
end
sData.MODIFIED=1;
data.sData=sData;
data.sData.comp_norm(1:length(sData.data(1,:)))={[]};
set(pre_h,'UserData',data);
end
else
args=create_args(assign,'D');
if length(args) == 1
len=max(str2num(args{1}));
if ~isscalar(len)
answer='??? Illegal assignment.';
return;
elseif len > length(sData.data(:,1)) | min(str2num(args{1})) < 1
answer='??? Illegal assignment.';
return;
end
if ~all(size(answer) == size(sData.data(str2num(args{1})))) &...
~isscalar(answer)
answer='??? Illegal assignment.';
else
sData.data(str2num(args{1}))=answer;
sData.MODIFIED=1;
data.sData=sData;
[i,j]=ind2sub(size(sData.data),str2num(args{1}));
data.sData.comp_norm(j)={[]};
set(pre_h,'UserData',data);
end
else
len=max(str2num(args{1}));
dim=max(str2num(args{2}));
asize=size(answer);
msize=size(sData.data);
if ~isscalar(len) | ~isscalar(dim)
answer='??? Illegal assignment.';
return;
elseif len > msize(1) | min(str2num(args{1})) < 1
answer='??? Illegal assignment.';
return;
elseif dim > msize(2) | min(str2num(args{2})) < 1
answer= '??? Illegal assignment.';
return;
end
len = length(str2num(args{1}));
dim = length(str2num(args{2}));
if ~all([len dim] == asize) & ~isscalar(answer)
answer='??? Illegal assignment.';
return;
else
sData.data([str2num(args{1})],[str2num(args{2})])=answer;
sData.MODIFIED=1;
data.sData=sData;
data.sData.comp_norm(str2num(args{2}))={[]};
set(pre_h,'UserData',data);
end
end
end
end
end
if sData.MODIFIED
selected=getfield(get(pre_h,'UserData'),'selected_vects');
vector_h=getfield(get(pre_h,'UserData'),'vector_h');
vect_mean_h=getfield(get(pre_h,'UserData'),'vect_mean_h');
vects=zeros(length(sData.data(:,1)));
vects(selected)=1;
draw_vectors(vects,vector_h);
vect_means(sData,vect_mean_h,selected);
pro_tools('plot_hist');
pro_tools('c_stat');
cplot_mimema;
end
%%% Subfunction: check_assign %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function [name,string,skip]=check_assign(string,evaltype)
reswords=[{'D'};{'x'};{'xs'}];
flag=0;
pre_h=findobj(get(0,'Children'),'Tag','Preprocess');
if isempty(pre_h)
man_h=findobj(get(0,'Children'),'Tag','Management');
clip_h=findobj(get(0,'Children'),'Tag','Clipping');
errordlg('''Preprocess'' -window does not exist. Terminating program.');
if ~isempty(man_h)
close man_h;
end
if ~isempty(clip_h)
close clip_h;
end
return;
end
EMPTY=isempty(get(getfield(get(pre_h,'UserData'),'comp_names_h'),'Value'));
[name,s]=give_token(string,evaltype);
skip=length(s);
if strcmp(evaltype,'single') & ~strcmp(name,'x')
string='NOTASSIGN';
return;
end
if strcmp(name,'other') & ~strcmp(s,'x')
string = 'error';
return;
end
if strcmp(name,[{'x'};{'xs'}])
comp_names_h=getfield(get(gcf,'UserData'),'comp_names_h');
if isempty(get(comp_names_h,'Value'))
errordlg('There are not components chosen.');
string='error';
return;
end
end
if skip == length(string) | ~strcmp(name,reswords)
string = 'NOTASSIGN';
return;
end
if (strcmp(name,'x') | strcmp(name,'xs')) & EMPTY
errordlg('There are not components chosen.');
string = 'error';
return;
end
[t,s]=give_token(string(length(name)+1),evaltype);
if strcmp(t,'(')
flag=1;
end
[foo,skip]=check_token(string,evaltype);
if length(name) ~= skip-1
skip=skip-1;
tmp=string(length(name)+1:skip);
else
tmp = [];
end
if flag & tmp(length(tmp)) ~= ')'
tmp(length(tmp)+1)=')';
end
if skip==length(string)
return;
end
skip=skip+1;
if length(string) ~= skip
[t,s]=give_token(string(skip+1:length(string)),evaltype);
else
string='NOTASSIGN';
return;
end
if ~strcmp(t,'=')
string = 'NOTASSIGN';
return;
end
string=tmp;
skip = skip+2;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%% Subfunction: isscalar %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function bool = isscalar(x)
m= size(x);
bool = m(1) == 1 & m(2) == 1;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%% Subfunction: create_args %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function args=create_args(string,type)
arg2='';
i=2;
j=1;
pre_h=findobj(get(0,'Children'),'Tag','Preprocess');
msize=size(getfield(getfield(get(pre_h,'UserData'),'sData'),'data'));
if string(i) == ':'
arg1=num2str(cat(2,'1:',num2str(msize(1))));
i=i+1;
j=j+length(arg1);
end
while string(i) ~=',' & string(i) ~=')'
arg1(j)=string(i);
i=i+1;
j=j+1;
end
if string(i) ==','
j=1;
i=i+1;
if string(i)==':'
switch type
case 'x'
arg2='1';
case 'cs'
arg2=num2str(get_indices);
case 'D'
arg2=num2str(cat(2,'1:',num2str(msize(2))));
end
i=i+1;
j=j+length(arg2);
end
while string(i) ~= ')'
arg2(j)=string(i);
j=j+1;
i=i+1;
end
end
args{1}=arg1;
if ~isempty(arg2)
args{2} = arg2;
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function [str,skip] = check_token(string,evaltype)
pre_h=findobj(get(0,'Children'),'Tag','Preprocess');
tmp_string=string;
[t,s]=give_token(tmp_string,evaltype);
skip=length(s);
if strcmp(t,'c')
if isempty(get(getfield(get(pre_h,'UserData'),'comp_names_h'),'Value'))
errordlg('There are no components chosen.');
str='error';
return;
end
index=getfield(get_indices,{1});
str=cat(2,'[',num2str(index),']');
if skip == length(tmp_string)
return;
end
tmp_string=tmp_string(skip+1:length(tmp_string));
[t,s] = give_token(tmp_string,evaltype);
if ~strcmp(t,'(')
return;
end
[args,skip2] = get_args(tmp_string(length(s)+1:length(tmp_string)),'c',...
evaltype);
skip=skip+skip2+2;
if strcmp(args,'error')
str = 'error'
return;
elseif ~strcmp(args,'all')
str=cat(2,'getfield(',str,',',args,')');
else
str=cat(2,'getfield(',str,',{[1]})');
end
elseif strcmp(t,'cs')
if isempty(get(getfield(get(pre_h,'UserData'),'comp_names_h'),'Value'))
errordlg('There are no components chosen.');
str='error';
return;
end
str =cat(2,'[',num2str(get_indices),']');
if length(s) == length(string)
return;
end
tmp_string=tmp_string(1+length(s):length(string));
[t,s]=give_token(tmp_string,evaltype);
if ~strcmp(t,'(')
return;
else
[args,skip2]=get_args(tmp_string(1+length(s):length(tmp_string)),'cs',...
evaltype);
skip=2+skip+skip2;
if strcmp(args,'error')
str='error';
return;
elseif ~strcmp(args,'all')
str = cat(2,'getfield(',str,',',args,')');
else
tmp_str=str;
str=cat(2,'[getfield(',str,',','{1})');
for i=2:length(get_indices)
str=cat(2,str,';getfield(',tmp_str,',',sprintf('{%d})',i));
end
str=cat(2,str,']');
end
end
elseif strcmp(t,'dim')
ind1=getfield(size(getfield(getfield(get(pre_h,'UserData'),'sData'),'data')),{2});
str=cat(2,'[',num2str(ind1),']');
if length(s)==length(string)
return;
end
tmp_string=string(1+length(s):length(string));
[t,s]=give_token(tmp_string,evaltype);
if ~strcmp(t,'(')
return;
end
skip=1+skip+length(s);
[args,skip2]=get_args(tmp_string(1+length(s):length(tmp_string)),'dim',...
evaltype);
if strcmp(args,'error')
str = 'error';
return;
else
skip=skip+skip2;
if ~strcmp(args,'all')
str=cat(2,'getfield(',str,',',args,')');
end
end
elseif strcmp(t,'dlen')
ind1=getfield(size(getfield(getfield(get(pre_h,'UserData'),'sData'),'data')),{1});
str=cat(2,'[',num2str(ind1),']');
if length(s)==length(string)
return;
end
tmp_string=string(1+length(s):length(string));
[t,s]=give_token(tmp_string,evaltype);
if ~strcmp(t,'(')
return;
end
skip=skip+length(s);
[args,skip2]=get_args(tmp_string(1+length(s):length(tmp_string)),'dlen',...
evaltype);
if strcmp(args,'error')
str='error';
return;
else
skip=1+skip+skip2;
if ~strcmp(args,'all')
str=cat(2,'getfield(',str,',',args,')');
end
end
elseif strcmp(t,'x')
if isempty(get(getfield(get(pre_h,'UserData'),'comp_names_h'),'Value'))
errordlg('There are not components chosen.');
str='error';
return;
end
len=getfield(size(getfield(getfield(get(pre_h,'UserData'),...
'sData'),'data')),{1});
index=num2str(getfield(get_indices,{1}));
h_str='findobj(get(0,''Children''),''Tag'',''Preprocess'')';
get_str=cat(2,'getfield(get(',h_str,',''UserData''),''sData'')');
get_str=cat(2,'getfield(',get_str,',''data'')');
str=cat(2,'getfield(',get_str,',{[1:',num2str(len),'],',index,'})');
if length(s) == length(string)
return;
end
tmp_string=string(1+length(s):length(string));
[t,s]=give_token(tmp_string,evaltype);
if ~strcmp(t,'(');
return;
end
skip=skip+length(s);
[args,skip2]=get_args(tmp_string(1+length(s):length(tmp_string)),'x',...
evaltype);
if strcmp(args,'error')
str = 'error';
return;
else
skip=1+skip+skip2;
if ~strcmp(args,'all')
str=cat(2,'getfield(',str,',',args,')');
end
end
elseif strcmp(t,'xs')
if isempty(get(getfield(get(pre_h,'UserData'),'comp_names_h'),'Value'))
errordlg('There are not components chosen.');
str='error';
return;
end
len=getfield(size(getfield(getfield(get(pre_h,'UserData'),...
'sData'),'data')),{1});
index=get_indices;
index=cat(2,'[',num2str(index),']');
h_str='findobj(get(0,''Children''),''Tag'',''Preprocess'')';
get_str=cat(2,'getfield(get(',h_str,',''UserData''),''sData'')');
get_str=cat(2,'getfield(',get_str,',''data'')');
str=cat(2,'getfield(',get_str,',{[1:',num2str(len),'],',index,'})');
if length(s) == length(string)
return;
end
tmp_string=string(1+length(s):length(string));
[t,s]=give_token(tmp_string,evaltype);
if ~strcmp(t,'(')
return;
end
skip=1+skip+length(s);
[args,skip2]=get_args(tmp_string(1+length(s):length(tmp_string)),'xs',...
evaltype);
if strcmp(args,'error')
str = 'error';
return;
elseif ~strcmp(args,'all')
str=cat(2,'getfield(',str,',',args,')');
skip=skip+skip2;
else
skip=skip+skip2;
[dlen,dim]=size(eval(str));
tmp_str=str;
str=cat(2,'[','getfield(',tmp_str,sprintf(',{1:%d,1})',dlen));
for i=2:dim
tmp=sprintf(',{1:%d,%d})',dlen,dim);
str=cat(2,str,';','getfield(',tmp_str,tmp);
end
str=cat(2,str,']');
end
elseif strcmp(t,'D')
get_h='findobj(get(0,''Children''),''Tag'',''Preprocess'')';
str=cat(2,'getfield(getfield(get(',get_h,',''UserData''),''sData''),''data'')');
if length(s) >= length(tmp_string)
return;
end
tmp_string=tmp_string(1+length(s):length(tmp_string));
[t,s]=give_token(tmp_string,evaltype);
if ~strcmp(t,'(')
return;
else
tmp_string=tmp_string(1+length(s):length(tmp_string));
skip = skip+length(s);
[args, skip2]=get_args(tmp_string,'D',evaltype);
if strcmp(args,'error')
str='error';
return;
elseif ~strcmp(args,'all')
str=cat(2,'getfield(',str,',',args,')');
skip=1+skip+skip2;
else
skip=1+skip+skip2;
[dlen,dim]=size(eval(str));
tmp_str=str;
str=cat(2,'[getfield(',str,sprintf(',{1:%d,1})',dlen));
for i=2:dim
tmp=sprintf(',{1:%d,%d}',dlen,i);
str=cat(2,str,';getfield(',tmp_str,tmp,')');
end
str=cat(2,str,']');
end
end
else
if strcmp(t,'(')
str = t;
str2='';
tmp_string=tmp_string(1+length(s):length(tmp_string));
while ~strcmp(str2,')') & ~isempty(tmp_string)
[str2,skip2]=check_token(tmp_string,evaltype);
if strcmp(str2,'error')
str='error';
return;
end
skip=skip+skip2;
tmp_string=tmp_string(skip2+1:length(tmp_string));
str=cat(2,str,str2);
end
if ~strcmp(str2,')')
str = 'error';
end
else
str = s;
end
end
%%% Subfunction: get_args %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function [str,skip] = get_args(string,flag,evaltype)
res_words=[{'D'};{'c'};{'cs'};{'dim'};{'dlen'};{'x'};{'xs'}];
NOTALL=1;
if isempty(string)
str='error'
skip=[];
return;
end
[t,s] = give_token(string,evaltype);
skip=length(s);
if any(strcmp(t,res_words));
[str,skip2] = check_token(string,evaltype);
string=string(1+length(s):length(string));
str=cat(2,'{[',str);
[t,s]=give_token(string,evaltype);
elseif t==')' | t==','
str = 'error';
return;
elseif strcmp(t,':');
if length(s) == length(string)
str='error';
return;
end
[t,s]=give_token(string(1+length(s):length(string)),evaltype);
if t == ')'
str = 'all';
return;
end
switch flag
case {'c','cs','dim','dlen'}
str= '{[1';
otherwise
str=cat(2,'{[',get_all('vect'));
end
NOTALL=0;
string=string(1+length(s):length(string));
[t,s]=give_token(string,evaltype);
skip=skip+1;
else
str = cat(2,'{[',s);
end
str2 =[];
if ~strcmp(t,',') & ~strcmp(t,')')
skip=skip-length(s);
end
while ~strcmp(t,',') & ~strcmp(t,')') & NOTALL;
str=cat(2,str,str2);
[t,s] = give_token(string,evaltype);
if length(s) == length(string)
str = 'error';
return;
end
string=string(1+length(s):length(string));
skip=skip+length(s);
[t,s]=give_token(string,evaltype);
if length(s) == length(string) & ~strcmp(t,')')
str = 'error';
return;
end
[str2,foo]=check_token(string,evaltype);
end
if NOTALL & ~strcmp(t,')')
skip=skip+1;
end
if strcmp(t,')')
str=cat(2,str,']}');
return
end
str=cat(2,str,']',',','[');
str2 = [];
[t,s] = give_token(string,evaltype);
if strcmp(t,')')
str = 'error'
return;
end
NOTALL=1;
string=string(1+length(s):length(string));
[t,s]=give_token(string,evaltype);
if strcmp(t,':');
switch flag
case {'c','dim','dlen','x'}
str=cat(2,str,'1');
case 'D'
str=cat(2,str,get_all('comp'));
case {'cs','xs'}
str=cat(2,str,'1:',num2str(length(get_indices)));
end
NOTALL=0;
if length(s) == length(string)
str='error';
return;
end
string=string(1+length(s):length(string));
[t,s]=give_token(string,evaltype);
end
if ~strcmp(t,')') & NOTALL
skip=skip-1;
end
while ~strcmp(t,')') & NOTALL
str=cat(2,str,str2);
skip=skip+length(s);
if length(s) == length(string) & ~strcmp(t,')')
str='error';
return;
end
[str2,foo]=check_token(string,evaltype);
string=string(1+length(s):length(string));
[t,s]=give_token(string,evaltype);
end
if ~strcmp(t,')')
str='error';
return;
end
str=cat(2,str,str2,']}');
skip=skip+length(s);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%% Subfunction: get_all %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function str=get_all(vect_or_comp)
pre_h=findobj(get(0,'Children'),'Tag','Preprocess');
switch vect_or_comp
case 'vect'
dim=getfield(size(getfield(getfield(get(pre_h,'UserData'),...
'sData'),'data')),{1});
str=cat(2,'1:',num2str(dim));
case 'comp'
dim=getfield(size(getfield(getfield(get(pre_h,'UserData'),...
'sData'),'data')),{2});
str=cat(2,'1:',num2str(dim));
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function [token,str]=give_token(string,evaltype)
n=length(string);
i=1;
char=string(i);
switch analyze_char(string(i));
case 'num'
token='num';
while i <= n & strcmp('num',analyze_char(string(i)))
str(i)=string(i);
i=i+1;
end
case 'other'
switch string(i)
case ':'
token = ':';
case ','
token = ',';
case '('
token = '(';
case ')'
token = ')';
case '='
token = '=';
otherwise
token='other';
end
str=string(i);
case 'alpha'
while i <= n & strcmp('alpha',analyze_char(string(i)))
str(i)=string(i);
i=i+1;
end
token = find_res_word(str,evaltype);
end
%%% Subfunction: analyze_char %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function type=analyze_char(char)
if ((char-0) >= ('0'-0) & (char-0) <= ('9'-0))
type='num';
elseif ((char-0) >= ('a'-0) & (char-0) <= ('z'-0)) ...
| ((char-0) >= ('A'-0) & (char-0) <= ('Z'-0))
type='alpha';
else
type='other';
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%% Subfunction: find_res_word %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function token = find_res_word(string,evaltype)
reswords=[{'D'};{'c'};{'cs'};{'dim'};{'dlen'};{'x'};{'xs'};{'other'}];
for i=1:length(reswords);
token=reswords{i};
if strcmp(string,reswords{i})
if strcmp(evaltype,'single') & ~strcmp(string,'x')
token = 'other';
end
return;
end
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function close_func(varargin)
switch varargin{1}
case 'close_c'
str='% Closing the ''Clipping'' -window...';
clip_h=findobj(get(0,'Children'),'Tag','Clipping');
close(clip_h);
case 'close_sD'
str='% Closing the ''Data Set Management'' -window...';
sD_h=findobj(get(0,'Children'),'Tag','Management');
close(sD_h);
case 'close_w'
str='% Closing the ''Windowed'' -window...';
win_h=findobj(get(0,'Children'),'Tag','Window');
close(win_h);
case 'close_s'
str='% Closing the ''Select'' -window...';
sel_h=findobj(get(0,'Children'),'Tag','Select');
close(sel_h);
case 'close_d'
str='% Closing the ''Delay'' -window...';
del_h=findobj(get(0,'Children'),'Tag','Delay');
close(del_h);
end
if nargin ~= 2
pre_h=findobj(get(0,'Children'),'Tag','Preprocess');
preh_udata=get(pre_h,'UserData');
str2=cat(2,'preprocess(''',varargin{1},''',''foo'');');
preh_udata.LOG{length(preh_udata.LOG)+1}=str;
preh_udata.LOG{length(preh_udata.LOG)+1}=str2;
set(pre_h,'UserData',preh_udata);
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function log_file
answer=inputdlg('Give the name of the outputfile:','LOG function',1,...
{'log_function'});
if isempty(answer)
return;
end
tmp=clock;
str =cat(2,'% Created: ',...
date,...
' ',sprintf('%d:%d\n%\n\n',tmp(4),tmp(5)));
pre_h=findobj(get(0,'Children'),'Tag','Preprocess');
LOG=getfield(get(pre_h,'UserData'),'LOG');
file=cat(2,pwd,'/',answer{1},'.m');
fid =fopen(file,'w');
arg=LOG{2}(12:length(LOG{2})-2);
fprintf(fid,'%s\n \n',cat(2,'function ',answer{1},'(',arg,')'));
fprintf(fid,'%s\n',str);
for i=1:length(LOG)
fprintf(fid,'%s\n',LOG{i});
end
fclose(fid);
disp(sprintf('LOG-file ''%s'' is done.',file));
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function get_selected_inds(varargin)
if nargin == 1
LOG=1;
answer = {varargin{1}};
else
LOG=0;
end
selected=getfield(get(gcf,'UserData'),'selected_vects');
if ~LOG
answer=inputdlg('Give the name of the output variable:',...
'',1,{'indices'});
end
if isempty(answer) | isempty(answer{1})
return;
else
assignin('base',answer{1},selected);
disp(cat(2,'Indices of the selected vectors are set to the workspace ',...
sprintf(' as ''%s''.',answer{1})));
if ~LOG
data=get(gcf,'UserData');
data.LOG{length(data.LOG)+1}=...
'% Saving indices of the selected vectors to the workspace.';
data.LOG{length(data.LOG)+1}=cat(2,'preprocess(''get_inds'',',...
'''',answer{1},''');');
set(gcf,'UserData',data);
end
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function no_of_selected(varargin)
if nargin == 1
pre_h=findobj(get(0,'Children'),'Tag','Preprocess');
set(0,'CurrentFigure',pre_h);
LOG = 1;
else
LOG = 0;
end
results_h=getfield(get(gcf,'UserData'),'results_h');
no=length(getfield(get(gcf,'UserData'),'selected_vects'));
str={sprintf('Number of selected vectors: %d\n', no)};
set(results_h,'String',str,'HorizontalAlignment','left');
if ~LOG
data=get(gcf,'UserData');
data.LOG{length(data.LOG)+1}='% Number of selected vectors';
data.LOG{length(data.LOG)+1}='preprocess(''no_of_sel'',''foo'');';
set(gcf,'UserData',data);
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function select_all_comps(varargin)
if nargin == 1
pre_h=findobj(get(0,'Children'),'Tag','Preprocess');
set(0,'CurrentFigure',pre_h);
LOG=1;
else
LOG=0;
end
comp_names_h=getfield(get(gcf,'UserData'),'comp_names_h');
set(comp_names_h,'Value',[1:length(get(comp_names_h,'String'))]);
sel_comp;
if ~LOG
data=get(gcf,'UserData');
data.LOG{length(data.LOG)+1}='% Select all components';
data.LOG{length(data.LOG)+1}='preprocess(''sel_all_comps'',''foo'');';
set(gcf,'UserData',data);
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function code=write_log_code(indices,arg1,arg2,arg3,arg4,arg5,arg6);
str=textwrap({num2str(indices)},500);
code{1}=sprintf('inds=[];');
for i=1:length(str);
code{i+1}=sprintf(' inds=cat(2,inds,[%s]);',str{i});
end
str=cat(2,'preprocess(''''clip_data'''',''''',arg1,' ',num2str(arg2),' ',...
num2str(arg3),' ',num2str(arg4),...
' ',num2str(arg5),' ',num2str(arg6),' ');
code{length(code)+1}=cat(2,'eval(cat(2,',...
'''',str,'''',...
',num2str(inds),'''''');''));');
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%