mfbox_pre_varthresholdg.m

toolbox for spm 5 for data, model free analysis

function varargout=mfbox_pre_varthresholdg(varargin)
% MFBOX_PRE_VARTHRESHOLDG M-file for mfbox_pre_varthresholdg.fig
%
% See also: GUIDE, GUIDATA, GUIHANDLES
%
% Copyright by Peter Gruber and Fabian J. Theis
% Signal Processing & Information Theory group
% Institute of Biophysics, University of Regensburg, Germany
% Homepage: http://research.fabian.theis.name
%           http://www-aglang.uni-regensburg.de
%
% This file is free software, subject to the 
% GNU GENERAL PUBLIC LICENSE, see gpl.txt

% Begin initialization code - DO NOT EDIT
gui_Singleton = 1;
gui_State = struct('gui_Name',mfilename, ...
    'gui_Singleton',gui_Singleton, ...
    'gui_OpeningFcn',@mfbox_pre_varthresholdg_OpeningFcn, ...
    'gui_OutputFcn',@mfbox_pre_varthresholdg_OutputFcn, ...
    'gui_LayoutFcn',[], ...
    'gui_Callback',[]);
if (nargin&&(ischar(varargin{1})))
    gui_State.gui_Callback = str2func(varargin{1});
end

if (nargout)
    [varargout{1:nargout}] = gui_mainfcn(gui_State,varargin{:});
else
    gui_mainfcn(gui_State,varargin{:});
end
% End initialization code - DO NOT EDIT


% --- Executes just before mfbox_pre_varthresholdg is made visible.
function mfbox_pre_varthresholdg_OpeningFcn(hObject,eventdata,handles,varargin)
% varargin{1} contains data X
% varargin{2} contains data mask
% varargin{3} contains data grid
% varargin{4} contains timeline
% varargin{5} (optional) contains the parameters
% varargin{6} (optional) contains runflag (0: param edit only, 1: param+run, 2: run only)
% varargin{7} (optional) contains private struct for mfbox_databackend('redraw_mask',mask,private)

handles.output = hObject;

handles.X = varargin{1};
if (isstruct(handles.X)), s = [handles.X.dim,handles.X.timesteps];
else, s = size(handles.X);
end
handles.dim = s(1:(end-1));
if (length(handles.dim)==1), handles.dim = [handles.dim,1]; end
handles.n = s(end);

if (nargin>4), handles.mask = varargin{2}; 
else handles.mask = logical(ones(handles.dim));
end
handles.dim = size(handles.mask);

if (nargin>5), handles.grid = varargin{3};
else handles.grid = mfbox_mkgrid(handles.dim)';
end
if (nargin>6), handles.timeline = varargin{4};
else handles.timeline = 0:(handles.n-1);
end
if (nargin>7), handles.params = varargin{5};
else handles.params = [];
end
if (nargin>8), handles.runflag = varargin{6};
else handles.runflag = 2;
end
if (nargin>9), handles.private = varargin{7}; end

ms = size(handles.mask);
ng = mat2cell(mfbox_mkgrid(ms),prod(ms),ones(1,length(ms)));
g = {};
for i=1:size(handles.grid,1), g{i} = reshape(handles.grid(i,:),handles.dim); end
if (length(handles.dim)>1)
    if (isstruct(handles.X))
        es = zeros(handles.dim); e = zeros(handles.dim);
        for i=1:handles.n
            v = reshape(mfbox_databackend('getdata',handles.X,i),handles.dim);
            es = es+v.^2; e = e+v;
        end
        handles.var = sqrt((es/handles.n)-(e/handles.n).^2);
    else
        es = zeros(handles.dim); e = zeros(handles.dim);
        handles.X = reshape(handles.X,[],handles.n);
        for i=1:handles.n
            v = handles.X(:,i);
            v = reshape(v,handles.dim);
            es = es+v.^2; e = e+v;
        end
        handles.var = sqrt((es/handles.n)-(e/handles.n).^2);
    end
else
    handles.var = [];
end
if (~isempty(handles.var))
    handles.var = reshape(interpn(g{:},handles.var,ng{:}),ms);
end
handles.maxvar = max(handles.var(:));

set(handles.priority,'String',mat2str(handles.params.priority));
set(handles.varmin,'Value',handles.params.varmin,'Min',0,'Max',1);
set(handles.varmax,'Value',handles.params.varmax,'Min',0,'Max',1);

settext(handles);

if (handles.runflag<1)
    if (isfield(handles,'private')), redrawmask(handles); end
    set(handles.run,'String','Set ...');
end

guidata(hObject,handles);

uiwait(handles.figure1); 


% --- Outputs from this function are returned to the command line.
function varargout=mfbox_pre_varthresholdg_OutputFcn(hObject,eventdata,handles) 
% varargout{1} contains the (edited) parameters
% varargout{2} (optional) contains (changed) private handle

if (isstruct(handles))
    close(handles.figure1);
    varargout{1} = handles.params;
    if (nargout>1), varargout{2} = handles.private; end
else
    varargout{1} = [];
end


% --- Executes on button press in run.
function run_Callback(hObject,eventdata,handles)

handles.params.priority = str2double(get(handles.priority,'String'));
handles.params.varmin = get(handles.varmin,'Value');
handles.params.varmax = get(handles.varmax,'Value');

guidata(hObject,handles);

uiresume(gcf);


function priority_Callback(hObject,eventdata,handles)

set(hObject,'String',mat2str(floor(str2double(get(hObject,'String')))));


% --- Executes during object creation, after setting all properties.
function priority_CreateFcn(hObject,eventdata,handles)

if (ispc&&(isequal(get(hObject,'BackgroundColor'),get(0,'defaultUicontrolBackgroundColor'))))
    set(hObject,'BackgroundColor','white');
end


% --- Executes on slider movement.
function varmin_Callback(hObject,eventdata,handles)

settext(handles);

if (isfield(handles,'private')), redrawmask(handles); end

guidata(hObject,handles);


% --- Executes during object creation, after setting all properties.
function varmin_CreateFcn(hObject,eventdata,handles)

if (isequal(get(hObject,'BackgroundColor'),get(0,'defaultUicontrolBackgroundColor')))
    set(hObject,'BackgroundColor',[.9,.9,.9]);
end


% --- Executes on slider movement.
function varmax_Callback(hObject,eventdata,handles)

settext(handles);

if (isfield(handles,'private')), redrawmask(handles); end

guidata(hObject,handles);


% --- Executes during object creation, after setting all properties.
function varmax_CreateFcn(hObject,eventdata,handles)

if (isequal(get(hObject,'BackgroundColor'),get(0,'defaultUicontrolBackgroundColor')))
    set(hObject,'BackgroundColor',[.9,.9,.9]);
end


% ---
function settext(handles)
v = get(handles.varmin,'Value'); w = get(handles.varmax,'Value');
if (v<w), set(handles.varlabel,'String',sprintf('%1.2f<=var<=%1.2f',v*handles.maxvar,w*handles.maxvar));
else set(handles.varlabel,'String',sprintf('var<%1.2f,%1.2f<var',w*handles.maxvar,v*handles.maxvar));
end


function redrawmask(handles)
v = get(handles.varmin,'Value'); w = get(handles.varmax,'Value');
if (v<w), mask = (handles.var>=v*handles.maxvar)&(handles.var<=w*handles.maxvar);
else mask = (handles.var<w*handles.maxvar)|(handles.var>v*handles.maxvar);
end
mfbox_databackend('redraw_mask',mask,handles.private);