Changeset 112 for aedes_readfid.m


Ignore:
Timestamp:
Mar 18, 2010, 1:16:34 PM (9 years ago)
Author:
tjniskan
Message:
  • Still trying to get EPI working

M aedes_readfid.m
M aedes_revision.m

File:
1 edited

Legend:

Unmodified
Added
Removed
  • aedes_readfid.m

    r111 r112  
    10561056  %% of block headers...
    10571057 
     1058  % Check the size of the FID-file
     1059  d=dir(fopen(file_fid));
     1060  file_sz = d.bytes/1024/1024; % File size in MB
     1061  nBlocks = ceil(file_sz/500); % Read data in 500 MB blocks
     1062 
    10581063  % Initialize waitbar
    10591064  if Dat.ShowWaitbar
     
    11071112  % Determine the number of values that will result from block header(s)
    11081113  nVals = (nbheaders*28)/prec;
     1114 
     1115 
    11091116 
    11101117  % Read the whole data including block headers etc...
     
    12911298    fprintf(1,['Processing block...%0',lnumstr,'d/%0',lnumstr,'d'],1,nBlocks);
    12921299    for ii=1:nBlocks
     1300      fprintf(1,repmat('\b',1,bsl));
     1301      fprintf(1,['Processing block...%0',lnumstr,'d/%0',lnumstr,'d'],ii,nBlocks);
    12931302      tmp_data = [];
    12941303      for kk=1:nRcv
     
    12981307          inds = inds(2:end);
    12991308        end
    1300           tmp_kspace = l_ReconstructKspace(kspace(:,:,inds),procpar,Dat);
    1301           tmp_data(:,:,:,:,kk) = fftshift(fftshift(fft(fft(tmp_kspace,[],1),[],2),1),2);
     1309        tmp_kspace = l_ReconstructKspace(kspace(:,:,inds),procpar,Dat);
     1310        tmp_data(:,:,:,:,kk) = fftshift(fftshift(fft(fft(tmp_kspace,[],1),[],2),1),2);
    13021311      end
    13031312      if ii==1
    13041313        data = sqrt(sum(tmp_data.*conj(tmp_data),5));
    13051314      elseif ii==nBlocks
    1306         data(:,:,:,((ii-1)*blksz+1):end) = sqrt(sum(tmp_data(:,:,:,2:end,:).*conj(tmp_data(:,:,:,2:end,:)),5));
     1315        data(:,:,:,((ii-1)*blksz+1):(nVols+1)) = sqrt(sum(tmp_data(:,:,:,2:end,:).*conj(tmp_data(:,:,:,2:end,:)),5));
    13071316      else
    13081317        data(:,:,:,((ii-1)*blksz+1):(ii*blksz)) = sqrt(sum(tmp_data(:,:,:,2:end,:).*conj(tmp_data(:,:,:,2:end,:)),5));
    13091318      end
    1310       fprintf(1,repmat('\b',1,bsl));
    1311       fprintf(1,['%0',lnumstr,'d/%0',lnumstr,'d'],ii,nBlocks);
    13121319    end
    13131320    fprintf(1,'\n')
Note: See TracChangeset for help on using the changeset viewer.

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