Changeset 104


Ignore:
Timestamp:
Mar 11, 2010, 2:39:28 PM (9 years ago)
Author:
tjniskan
Message:
  • Added support for multiple receiver EPI data. (EXPERIMENTAL)

M aedes_readfid.m
M aedes_revision.m

Files:
2 edited

Legend:

Unmodified
Added
Removed
  • aedes_readfid.m

    r103 r104  
    12741274if isfield(procpar,'rcvrs') && ~isempty(procpar.rcvrs)
    12751275  % Multiple receivers used
    1276   nRcv = length(find(procpar.rcvrs{1}=='y'));
    1277   data = [];
    1278   kspace2 = [];
    1279   for ii=1:nRcv
    1280     tmp_kspace = l_ReconstructKspace(kspace(:,:,ii),procpar,Dat);
    1281     kspace2(:,:,:,ii)=tmp_kspace;
    1282     if Dat.ReturnFTData
    1283       tmp_data = l_CalculateFFT(tmp_kspace,procpar,Dat);
    1284       data(:,:,:,ii) = tmp_data;
    1285     end
    1286   end
    1287   kspace = kspace2;
    1288   kspace2=[];
     1276  if isfield(procpar,'apptype') && strcmp(procpar.apptype{1},'imEPI')
     1277    % Store multiple receiver data in EPI measurements in 5th dimension
     1278    % and calculate sum-of-squares image
     1279    nRcv = length(find(procpar.rcvrs{1}=='y'));
     1280    data = [];
     1281    kspace2 = [];
     1282    for ii=1:nRcv
     1283      tmp_kspace = l_ReconstructKspace(kspace(:,:,ii:nRcv:end),procpar,Dat);
     1284      kspace2(:,:,:,:,ii)=tmp_kspace;
     1285    end
     1286    kspace = kspace2;
     1287    data = fftshift(fftshift(abs(fft(fft(kspace,[],1),[],2)),1),2);
     1288    data = sqrt(mean(data.*conj(data),5));
     1289    kspace2=[];
     1290   
     1291    if Dat.OrientImages && ~isempty(procpar) && ...
     1292        isfield(procpar,'orient') && any(Dat.AcqType==[2 3 4])
     1293      orient = procpar.orient{1};
     1294      if any(strcmpi(orient,{'xyz','trans90','cor90','sag90'}))
     1295        data = flipdim(aedes_rot3d(data,1,3),2);
     1296      elseif strcmpi(orient,'oblique')
     1297        data = flipdim(flipdim(data,1),2);
     1298      else
     1299        data = flipdim(flipdim(data,1),2);
     1300      end
     1301    end
     1302
     1303   
     1304  else
     1305    nRcv = length(find(procpar.rcvrs{1}=='y'));
     1306    data = [];
     1307    kspace2 = [];
     1308    for ii=1:nRcv
     1309      tmp_kspace = l_ReconstructKspace(kspace(:,:,ii),procpar,Dat);
     1310      kspace2(:,:,:,ii)=tmp_kspace;
     1311      if Dat.ReturnFTData
     1312        tmp_data = l_CalculateFFT(tmp_kspace,procpar,Dat);
     1313        data(:,:,:,ii) = tmp_data;
     1314      end
     1315    end
     1316    kspace = kspace2;
     1317    kspace2=[];
     1318  end
    12891319else
    12901320  % Only one receiver
  • aedes_revision.m

    r103 r104  
    9393% bash-script every time it is called so that this file "aedes_revision.m" is
    9494% always in the list of committed files. DO NOT EDIT THE NEXT LINE!!!
    95 % - SVN Hook -
     95% - Svn Hook -
Note: See TracChangeset for help on using the changeset viewer.

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