fmri_session_profile_ui.m

绝对经典,老外制作的功能强大的matlab实现PLS_TOOBOX

      set(h,'String',num2str(new_num_runs));
      setappdata(gcf,'SessionNumRuns',new_num_runs);
   end;
   
   return;						% EditRuns


%----------------------------------------------------------------------------
function EditDatamatPrefix

   datamat_prefix = deblank(strjust(get(gcbo,'String'),'left')); 
   set(gcbo,'String',datamat_prefix);

   if ~isempty(findstr(datamat_prefix,' '));
      msg = sprintf('ERROR: Datamat prefix cannot contain any space');
      set(findobj(gcf,'Tag','MessageLine'),'String',msg);
      return;
   end;

   setappdata(gcf,'SessionDatamatPrefix',datamat_prefix);

   return;						% EditDatamatPrefix


%----------------------------------------------------------------------------
function status = ChkModified()
% Output:
%          status = 0, error
% 	   status = 1, ok

   status = 1;

   [session_info,session_file,old_session_info] = GetSessionInfo;

   if (isequal(session_info,old_session_info) == 0),
      dlg_title = 'Session Information has been changed';
      msg = 'WARNING: The session information has been changed.  Do you want to save it?';
      response = questdlg(msg,dlg_title,'Yes','No','Cancel','Yes');

      switch response,
         case 'Yes'
 	      status = SaveSessionInfo;		
         case 'Cancel'
 	      status = 0;
         case 'No'
 	      status = 2;
      end; 
   end;

   return;						% ChkModified


%----------------------------------------------------------------------------
function SetSessionInfo(session_info,session_file)

   setappdata(gcf,'SessionDescription',session_info.description);
   setappdata(gcf,'SessionDatamatPrefix',session_info.datamat_prefix);
   setappdata(gcf,'SessionPLSDir',session_info.pls_data_path);
   setappdata(gcf,'SessionNumConditions',session_info.num_conditions);
   setappdata(gcf,'SessionConditions',session_info.condition);
   setappdata(gcf,'SessionNumRuns',session_info.num_runs);
   setappdata(gcf,'SessionRunInfo',session_info.run);

   if ~isfield(session_info,'across_run')
      session_info.across_run = 1;
   end

   setappdata(gcf,'SessionAcrossRun',session_info.across_run);

   if ~isfield(session_info,'condition_baseline')
      session_info.condition_baseline = {};
      for i=1:session_info.num_conditions
         session_info.condition_baseline=[session_info.condition_baseline, ...
                                                {[0 1]}];
      end
   end

   setappdata(gcf,'SessionConditionBaseline',session_info.condition_baseline);

   if ~isfield(session_info,'num_conditions0')
      session_info.num_conditions0 = session_info.num_conditions;
   end

   setappdata(gcf,'SessionNumConditions0',session_info.num_conditions0);

   if ~isfield(session_info,'condition0')
      session_info.condition0 = session_info.condition;
   end

   setappdata(gcf,'SessionConditions0',session_info.condition0);

   if ~isfield(session_info,'condition_baseline0')
      session_info.condition_baseline0 = session_info.condition_baseline;
   end

   setappdata(gcf,'SessionConditionBaseline0',session_info.condition_baseline0);

   if ~isfield(session_info,'behavname_all')
      session_info.behavname_all = {};
      session_info.behavdata_all = [];
   end
   setappdata(gcf,'SessionBehavNameAll',session_info.behavname_all);
   setappdata(gcf,'SessionBehavDataAll',session_info.behavdata_all);

   if ~isfield(session_info,'behavname_each')
      session_info.behavname_each = {};
      session_info.behavdata_each = [];
   end
   setappdata(gcf,'SessionBehavNameEach',session_info.behavname_each);
   setappdata(gcf,'SessionBehavDataEach',session_info.behavdata_each);

   if ~isfield(session_info,'behavname_all_single')
      session_info.behavname_all_single = {};
      session_info.behavdata_all_single = [];
   end
   setappdata(gcf,'SessionBehavNameAllSingle',session_info.behavname_all_single);
   setappdata(gcf,'SessionBehavDataAllSingle',session_info.behavdata_all_single);

   if ~isfield(session_info,'behavname_each_single')
      session_info.behavname_each_single = {};
      session_info.behavdata_each_single = [];
   end
   setappdata(gcf,'SessionBehavNameEachSingle',session_info.behavname_each_single);
   setappdata(gcf,'SessionBehavDataEachSingle',session_info.behavdata_each_single);

   setappdata(gcf,'OldSessionInfo',session_info);
   setappdata(gcf,'SessionFile',session_file);

   if isempty(session_file)
      set(gcf,'Name','New Event Related fMRI Session Information');
   else
      set(gcf,'Name',['Session File: ' session_file]);
   end;

   return;						% SetSessionInfo


%----------------------------------------------------------------------------
function [session_info,session_file,old_session_info] = GetSessionInfo(),

   session_info.description = getappdata(gcf,'SessionDescription');
   session_info.pls_data_path = getappdata(gcf,'SessionPLSDir');
   session_info.datamat_prefix = getappdata(gcf,'SessionDatamatPrefix');
   session_info.num_conditions = getappdata(gcf,'SessionNumConditions');
   session_info.condition = getappdata(gcf,'SessionConditions');
   session_info.condition_baseline = getappdata(gcf,'SessionConditionBaseline');
   session_info.num_conditions0 = getappdata(gcf,'SessionNumConditions0');
   session_info.condition0 = getappdata(gcf,'SessionConditions0');
   session_info.condition_baseline0 = getappdata(gcf,'SessionConditionBaseline0');
   session_info.num_runs = getappdata(gcf,'SessionNumRuns');
   session_info.run = getappdata(gcf,'SessionRunInfo');
   session_info.across_run = getappdata(gcf,'SessionAcrossRun');
   session_info.behavname_all = getappdata(gcf,'SessionBehavNameAll');
   session_info.behavdata_all = getappdata(gcf,'SessionBehavDataAll');
   session_info.behavname_each = getappdata(gcf,'SessionBehavNameEach');
   session_info.behavdata_each = getappdata(gcf,'SessionBehavDataEach');
   session_info.behavname_all_single = getappdata(gcf,'SessionBehavNameAllSingle');
   session_info.behavdata_all_single = getappdata(gcf,'SessionBehavDataAllSingle');
   session_info.behavname_each_single = getappdata(gcf,'SessionBehavNameEachSingle');
   session_info.behavdata_each_single = getappdata(gcf,'SessionBehavDataEachSingle');

   session_file = getappdata(gcf,'SessionFile');
   old_session_info = getappdata(gcf,'OldSessionInfo');

   return;						% GetSessionInfo


%----------------------------------------------------------------------------
function CreateSTDatamat(),

   %  make sure the session information has been saved
   %
   [session_info,session_file,old_session_info] = GetSessionInfo;

   if isempty(session_file) | isequal(session_info,old_session_info) == 0

       status = ChkModified;

       session_file = getappdata(gcf,'SessionFile');

       if isempty(session_file)
          % status = 2;
          return;
       end

       if (status == 0)
          return;
       elseif (status == 2)
          msg1 = 'ERROR: Please save the session information';
          msg =  [msg1 ' before creating the datamat.'];

          set(findobj(gcf,'Tag','MessageLine'),'String',msg);
          return;
       end
   end;

   all_onsets = [session_info.run(:).evt_onsets];
   all_onsets = reshape(all_onsets, [session_info.num_conditions0 length(all_onsets)/session_info.num_conditions0]);

   empty_cond_across_run = 0;
   any_empty_cond_within_run = 0;

   for i = 1:size(all_onsets,1)
      empty_cond_within_run = 0;

      for j = 1:size(all_onsets,2)
         if isequal(all_onsets(i,j),{-1})
            empty_cond_within_run = empty_cond_within_run + 1;
            any_empty_cond_within_run = 1;
         end
      end

      if isequal(size(all_onsets,2),empty_cond_within_run)
         empty_cond_across_run = 1;
         break;
      end
   end

   if empty_cond_across_run
      msg = 'ERROR: At least one condition has no onset for all the runs';

      set(findobj(gcf,'Tag','MessageLine'),'String',msg);
      return;
   end

   if ~session_info.across_run & any_empty_cond_within_run
      msg = 'ERROR: There is a condition without onset while merging data within each run';

      set(findobj(gcf,'Tag','MessageLine'),'String',msg);
      return;
   end

   num_runs = getappdata(gcf,'SessionNumRuns');

   data_path = session_info.run(1).data_path;
   img_file = session_info.run(1).data_files{1};
   try
      dims = rri_imginfo(fullfile(data_path,img_file));
      num_slices = dims(3); 
   catch
      msg = sprintf('ERROR: Cannot access the image file: %s',img_file);
      set(findobj(gcf,'Tag','MessageLine'),'String',msg);
      return;
   end;

   fmri_create_datamat_ui(session_file,num_runs,num_slices);

   return;						% CreateSTDatamat


%------------------------------------------------------------------------------
function ModifyBehav

   session_file = getappdata(gcbf,'SessionFile');

   if isempty(session_file)
      return;
   end

   load(session_file);

   if findstr('BfMRIsession.mat', session_file)
      datamat_file = fullfile(session_info.pls_data_path, ...
		[session_info.datamat_prefix, '_BfMRIdatamat.mat']);
   else
      datamat_file = fullfile(session_info.pls_data_path, ...
		[session_info.datamat_prefix, '_fMRIdatamat.mat']);
   end

   load(datamat_file);

   repeat = 1;

   while repeat
      [new_behavdata, new_behavname] = ...
         rri_edit_behav(num2str(behavdata),behavname,'Edit Behavior Data');

      if isempty(new_behavname) | isempty(new_behavdata)
         return;
      elseif size(st_datamat,1) == size(str2num(new_behavdata),1)
         repeat = 0;
      else
         msg = ['Number of rows should be equal to ' num2str(size(st_datamat,1))];
         uiwait(msgbox(msg, 'Error'));
         repeat = 1;
      end
   end

   if isequal(new_behavname, behavname) & isequal(str2num(new_behavdata), behavdata)
      return;
   end

   if ~exist('create_ver','var')
      create_ver = plsgui_vernum;
   end

   behavdata = str2num(new_behavdata);
   behavname = new_behavname;

   save(datamat_file,'behavdata','behavname','create_ver', ...
	'normalize_volume_mean','st_coords','st_datamat','st_dims', ...
	'st_evt_list','st_origin','st_sessionFile','st_voxel_size', ...
	'st_win_size');

   msg = 'Datamat has been modified with new behavior data';
   uiwait(msgbox(msg, 'Error'));

   return;						% ModifyBehav


%------------------------------------------------------------------------------
function delete_fig()
    if (CloseSessionInput == 1);
        uiresume
    end

    return


%----------------------------------------------------------------------------
function PathSessionInfo(init_flag)
%  init_flag = 0  for clear operation
%  init_flag = 1  for initialization
%

   if (init_flag == 0 & ChkModified == 0)
      return;			% error
   end;

   old_session_info = getappdata(gcf,'OldSessionInfo');
   session_info = old_session_info;
   session_info = rri_changepath_se(session_info);

   if isempty(session_info)
      session_info = old_session_info;
   end

   setappdata(gcf,'SessionPLSDir',session_info.pls_data_path);
   setappdata(gcf,'SessionRunInfo',session_info.run);

   return;						% PathSessionInfo


%----------------------------------------------------------------------------
function SelectMergeDataWithinRun()

if get(findobj(gcf,'Tag','MergeDataWithinRunButton'),'Value') == 0		% click itself
   set(findobj(gcf,'Tag','MergeDataWithinRunButton'),'Value',1);
else
   set(findobj(gcf,'Tag','MergeDataAcrossRunsButton'),'Value',0);
   setappdata(gcf,'SessionAcrossRun',0);
end

   return;					% SelectMergeDataWithinRun


%----------------------------------------------------------------------------
function SelectMergeDataAcrossRuns()

if get(findobj(gcf,'Tag','MergeDataAcrossRunsButton'),'Value') == 0		% click itself
   set(findobj(gcf,'Tag','MergeDataAcrossRunsButton'),'Value',1);
else
   set(findobj(gcf,'Tag','MergeDataWithinRunButton'),'Value',0);
   setappdata(gcf,'SessionAcrossRun',1);
end

   return;					% SelectMergeDataAcrossRuns