function Radii = ClusterAreas(FilePrefix,N)
%   ================================================
%   Lambert-Liao
%   2008 March 7
%   Cluster Areas
%   ================================================
%
%   Radii = ClusterAreas(FilePrefix,N)


RadiusData=zeros(N,2);
tic;

RadiusData=zeros(N,2);

for i = 1:N
    CurrentColorImage = imread([FilePrefix num2str(i-1,'%06g') '.png']);
    CurrentColorImage = CurrentColorImage(:,91:170,:);
    
    [a_r,b_r]=CorrelationRadius2(CurrentColorImage(:,:,1),0.8,5);
    RadiusData(i,1)=a_r;

    [a_g,b_g]=CorrelationRadius2(CurrentColorImage(:,:,2),0.8,5);
    RadiusData(i,2)=a_g;    

    ElapsedTime = toc;
    
%    [a2,b2]=length_phi_test(CurrentColorImage(:,:,2),0.1,0.5);  
%    [a1,b1]=length_phi_test(CurrentColorImage(:,:,1),0.1,0.5);     

 %   Radii(i,1)=b1;
 %   Radii(i,2)=b2;
    RemainingTime = (N - i)*ElapsedTime/i;
    fprintf(['Remaining time: ' ,num2str(RemainingTime), ' s.  Elapsed time: ' ,num2str(ElapsedTime), ' s.  \n']);
 %   fprintf(['Red Correlation Radius: ' ,num2str(b1), '. Green Correlation Radius: ' ,num2str(b2), ' \n']);     
    
end

Radii=RadiusData;