beam_profile.m

mni = 0;
isum = 0
mnj = 0;
jsum = 0;
dir = 'one_C001S0001'

cd(['f:/Anurag/011206_1'])

c = 1

for n = 1001:2:1012
    n
    nstr = num2str(n);
    i = imread([dir,'000', nstr(2:4), '.tif']);
    nstr = num2str(n+1);    
    j = imread([dir,'000', nstr(2:4), '.tif']);
    
    
    [max loc1(c)] = max(mean(i(:,:), 2));
    [max loc2(c)] = max(mean(j(:,:), 2));
    
    isum = isum + double(i(:,:));
    mni = mni + mean(i(:,:), 2);
    jsum = jsum + double(j(:,:));
    mnj = mnj + mean(j(:,:), 2);
    c = c + 1;
    
end


i = imread([dir,'000001', '.tif']);
j = imread([dir,'000002', '.tif']);

%for n=1:1024
%    [y pix1] = max(i,n);
%end