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 EDITgui_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});endif (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);endhandles.dim = s(1:(end-1));if (length(handles.dim)==1), handles.dim = [handles.dim,1]; endhandles.n = s(end);if (nargin>4), handles.mask = varargin{2}; else handles.mask = logical(ones(handles.dim));endhandles.dim = size(handles.mask);if (nargin>5), handles.grid = varargin{3};else handles.grid = mfbox_mkgrid(handles.dim)';endif (nargin>6), handles.timeline = varargin{4};else handles.timeline = 0:(handles.n-1);endif (nargin>7), handles.params = varargin{5};else handles.params = [];endif (nargin>8), handles.runflag = varargin{6};else handles.runflag = 2;endif (nargin>9), handles.private = varargin{7}; endms = 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); endif (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);    endelse    handles.var = [];endif (~isempty(handles.var))    handles.var = reshape(interpn(g{:},handles.var,ng{:}),ms);endhandles.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 ...');endguidata(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 handleif (isstruct(handles))    close(handles.figure1);    varargout{1} = handles.params;    if (nargout>1), varargout{2} = handles.private; endelse    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); endguidata(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); endguidata(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));endfunction 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);endmfbox_databackend('redraw_mask',mask,handles.private);