Changeset 81 for aedes_readfid.m


Ignore:
Timestamp:
Feb 23, 2009, 12:24:36 PM (10 years ago)
Author:
tjniskan
Message:
  • Added (again) some "chewing gum and ironwire" into aedes_readfid.m.

The handling and sorting of VNMR data should be rewritten in the near
future.

  • Added a function in aedes_readfid.m for creating fse/fsems

phasetables. It is, however, not in use at the moment so phasetables
should still be in the tablib directory.

M aedes_readfid.m
M aedes_revision.m

File:
1 edited

Legend:

Unmodified
Added
Removed
  • aedes_readfid.m

    r80 r81  
    12651265  %if ~strcmpi(seqcon(2:3),'cc') && ( isempty(Dat.phasetable) | ...
    12661266  %        isfield(procpar,'flash_converted') )
    1267   if ( isempty(Dat.phasetable) && ~strcmpi(seqcon(2:3),'sc') ) || isfield(procpar,'flash_converted')
     1267  if (( isempty(Dat.phasetable) && ~strcmpi(seqcon(2:3),'sc') ) && ...
     1268      ~(strcmpi(seqcon(2:3),'cc') && ArrayLength==size(kspace,3))) || ...
     1269      isfield(procpar,'flash_converted')
    12681270                                 
    12691271    % Sort uncompressed arrayed data
     
    12881290  else
    12891291    %% Sort RARE type data
    1290     if UsePhaseTable
     1292    if UsePhaseTable && ~isempty(Dat.phasetable)
    12911293      Dat.phasetable = Dat.phasetable';
    12921294      kspace(:,Dat.phasetable(:),:) = kspace;
     
    15581560end
    15591561
     1562%========================================================
     1563% A function for creating phasetables for fse and fsems
     1564%========================================================
     1565function phasetable=l_CreateFsemsPhaseTable(procpar)
     1566
     1567phasetable = [];
     1568if ~isfield(procpar,'etl') || ~isfield(procpar,'kzero')
     1569  return
     1570end
     1571etl = procpar.etl;
     1572kzero = procpar.kzero;
     1573nv = procpar.nv;
     1574
     1575t = (-nv/2+1):(nv/2);
     1576t = t(:);
     1577tt = flipud(t);
     1578tt = tt(:);
     1579
     1580phasetable = [reshape(t(1:nv/2),etl*2,[]);...
     1581  flipud(reshape(tt(1:nv/2),etl*2,[]))];
     1582phasetable = circshift(fliplr(phasetable),[0 kzero-1]);
    15601583
    15611584
Note: See TracChangeset for help on using the changeset viewer.

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