Changeset 111 for aedes_readfid.m


Ignore:
Timestamp:
Mar 18, 2010, 10:53:54 AM (9 years ago)
Author:
tjniskan
Message:
  • Still trying to get multireceiver EPI-data to work...

M aedes_readfid.m
M aedes_revision.m

File:
1 edited

Legend:

Unmodified
Added
Removed
  • aedes_readfid.m

    r110 r111  
    11131113 
    11141114  % Remove block headers from the data
    1115   kspace=kspace(nVals+1:end,:);
     1115  kspace(1:nVals,:)=[];
     1116  %kspace=kspace(nVals+1:end,:);
    11161117 
    11171118  % Transform to complex values
     
    12791280    % and calculate sum-of-squares image
    12801281    nRcv = length(find(procpar.rcvrs{1}=='y'));
    1281     tmp_data = [];
    12821282    nVols = size(kspace,3)/nRcv-1;
    12831283    data = zeros(procpar.nv,procpar.np/2,procpar.ns,nVols+1,'single');
     
    12911291    fprintf(1,['Processing block...%0',lnumstr,'d/%0',lnumstr,'d'],1,nBlocks);
    12921292    for ii=1:nBlocks
     1293      tmp_data = [];
    12931294      for kk=1:nRcv
    12941295        inds = [kk ((ii-1)*blksz*nRcv+kk):nRcv:min((nRcv*ii*blksz),kssz(3))];
     
    12971298          inds = inds(2:end);
    12981299        end
    1299         tmp_kspace = l_ReconstructKspace(kspace(:,:,inds),procpar,Dat);
    1300         tmp_data(:,:,:,:,kk) = fftshift(fftshift(fft(fft(tmp_kspace,[],1),[],2),1),2);
     1300          tmp_kspace = l_ReconstructKspace(kspace(:,:,inds),procpar,Dat);
     1301          tmp_data(:,:,:,:,kk) = fftshift(fftshift(fft(fft(tmp_kspace,[],1),[],2),1),2);
    13011302      end
    13021303      if ii==1
    13031304        data = sqrt(sum(tmp_data.*conj(tmp_data),5));
    13041305      elseif ii==nBlocks
    1305         data(:,:,:,((ii-1)*blksz+1):end) = sqrt(sum(tmp_data(:,:,:,2,:).*conj(tmp_data(:,:,:,2,:)),5));
     1306        data(:,:,:,((ii-1)*blksz+1):end) = sqrt(sum(tmp_data(:,:,:,2:end,:).*conj(tmp_data(:,:,:,2:end,:)),5));
    13061307      else
    1307         data(:,:,:,((ii-1)*blksz+1):(ii*blksz)) = sqrt(sum(tmp_data(:,:,:,2,:).*conj(tmp_data(:,:,:,2,:)),5));
     1308        data(:,:,:,((ii-1)*blksz+1):(ii*blksz)) = sqrt(sum(tmp_data(:,:,:,2:end,:).*conj(tmp_data(:,:,:,2:end,:)),5));
    13081309      end
    13091310      fprintf(1,repmat('\b',1,bsl));
Note: See TracChangeset for help on using the changeset viewer.

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