%
% Use to display raw data images.
% 

file = 'bruce_156';
%file = '/working/fmri/bruce/run9/80.img';
%file='/home/mirl/jim/code/pr/process/DATA_FINAL_FLT0001';
%file = '/v/raid10/mirl/jim/diffusion/kirkman/diffA.BX';
%file = '/v/raid1a/jim/ryan/merged_n_swapped_6000710';

slices =156;
read=256;
phase=256;
size = read * phase;

f2=zeros(read,phase,slices);

fid=fopen(file,'r','native');
for i=1:slices
  f1 = fread(fid, size, 'int16');
  f2(:,:,i)=reshape(f1,read,phase);
end
fclose(fid);

for i=1:slices
     imagesc(abs(f2(:,:,i)));
     pause(1);
end