Changeset 82 for aedes_readfid.m


Ignore:
Timestamp:
Jun 16, 2009, 10:54:45 AM (10 years ago)
Author:
tjniskan
Message:
  • Some minor modifications to aedes_readfid.m
  • Added BOLD strength estimation to

plugins/fmri_plugins/basic_fmri_analysis.m

  • Changed the image in about dialog
  • Added support for negative threshold in misclib/fmri_blob_overlay.m

M misclib/fmri_blob_overlay.m
M aedes_helpabout.m
M plugins/fmri_plugins/basic_fmri_analysis.m
M aedes_readfid.m
M aedes_revision.m

File:
1 edited

Legend:

Unmodified
Added
Removed
  • aedes_readfid.m

    r81 r82  
    492492  end
    493493 
     494  % If the sequence is a fast spin echo, try to construct the phasetable
     495  if isempty(phasetable) && isfield(procpar,'petable') && ...
     496      strncmpi(procpar.petable{1},'fse',3)
     497    phasetable = {'t1',l_CreateFsemsPhaseTable(procpar)};
     498  end
     499 
    494500  %% Abort and throw error if phasetable cannot be read and the FID-file
    495501  % has not been sorted
     
    13531359  % patches that hopefully work...
    13541360 
    1355  
    1356   if strcmp(seqcon,'nccnn') && length(size(kspace))==2 && ...
    1357           procpar.ns>=1 && AcqType~=1
     1361  if strcmp(seqcon,'nncnn')
     1362    kspace = permute(kspace,[1 3 2]);
     1363  elseif strcmp(seqcon,'nccnn') && length(size(kspace))==2 && ...
     1364      procpar.ns>=1 && AcqType~=1
    13581365    if ~isempty(Dat.phasetable)
    1359           kssz = size(kspace);
    1360           phsz = size(Dat.phasetable);
    1361           kspace=permute(reshape(kspace,procpar.np/2,...
    1362                 phsz(2),...
    1363                 kssz(2)/phsz(2)),[1 3 2]);
    1364           kspace = permute(reshape(kspace,...
    1365                 procpar.np/2,procpar.ns,kssz(2)/procpar.ns),[1 3 2]);
    1366           if Dat.Sorting
    1367                 kspace(:,Dat.phasetable(:),:)=kspace;
    1368           end
    1369         else
    1370           kspace = reshape(kspace,...
    1371                 [procpar.np/2 procpar.ns ...
    1372                 size(kspace,2)/procpar.ns]);
    1373           kspace=permute(kspace,[1 3 2]);
    1374         end
    1375        
     1366      kssz = size(kspace);
     1367      phsz = size(Dat.phasetable);
     1368      kspace=permute(reshape(kspace,procpar.np/2,...
     1369        phsz(2),...
     1370        kssz(2)/phsz(2)),[1 3 2]);
     1371      kspace = permute(reshape(kspace,...
     1372        procpar.np/2,procpar.ns,kssz(2)/procpar.ns),[1 3 2]);
     1373      if Dat.Sorting
     1374        kspace(:,Dat.phasetable(:),:)=kspace;
     1375      end
     1376    else
     1377      kspace = reshape(kspace,...
     1378        [procpar.np/2 procpar.ns ...
     1379        size(kspace,2)/procpar.ns]);
     1380      kspace=permute(kspace,[1 3 2]);
     1381    end
     1382   
    13761383  elseif strcmp(seqcon,'nscsn') && length(size(kspace))==3 && ...
    13771384          AcqType~=1
     
    15781585tt = tt(:);
    15791586
    1580 phasetable = [reshape(t(1:nv/2),etl*2,[]);...
    1581   flipud(reshape(tt(1:nv/2),etl*2,[]))];
     1587phasetable = [reshape(t(1:nv/2),[],etl);...
     1588  flipud(reshape(tt(1:nv/2),[],etl))];
    15821589phasetable = circshift(fliplr(phasetable),[0 kzero-1]);
    15831590
Note: See TracChangeset for help on using the changeset viewer.

Powered by Trac 1.0.9.Copyright © Juha-Pekka Niskanen 2008