Changeset 106


Ignore:
Timestamp:
Mar 12, 2010, 2:03:50 PM (9 years ago)
Author:
tjniskan
Message:
  • Fiddling around with multireceiver EPI-data. It's much slower now

but should consume less memory.

M aedes_readfid.m
M aedes_revision.m

Files:
2 edited

Legend:

Unmodified
Added
Removed
  • aedes_readfid.m

    r105 r106  
    12791279    nRcv = length(find(procpar.rcvrs{1}=='y'));
    12801280    data = [];
     1281    tmp_data = [];
    12811282    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=[];
     1283    for kk=1:(size(kspace,3)/nRcv-1)
     1284      for ii=1:nRcv
     1285        tmp_kspace = l_ReconstructKspace(kspace(:,:,[ii kk*nRcv+ii]),procpar,Dat);
     1286        tmp_data(:,:,:,:,ii) = fftshift(fftshift(fft(fft(tmp_kspace,[],1),[],2),1),2);
     1287       
     1288        %tmp_kspace = l_ReconstructKspace(kspace(:,:,ii:nRcv:end),procpar,Dat);
     1289        %kspace2(:,:,:,:,ii)=tmp_kspace;
     1290      end
     1291      if kk==1
     1292        data = sqrt(mean(tmp_data.*conj(tmp_data),5));
     1293      else
     1294        data(:,:,:,kk+1) = sqrt(mean(tmp_data(:,:,:,2,:).*conj(tmp_data(:,:,:,2,:)),5));
     1295      end
     1296    end
     1297    %kspace = kspace2;
     1298    %data = fftshift(fftshift(fft(fft(kspace,[],1),[],2),1),2);
     1299    %data = sqrt(mean(data.*conj(data),5));
     1300    %kspace2=[];
    12901301   
    12911302    % Remove reference image if requested
  • aedes_revision.m

    r105 r106  
    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