% clc
% clear all
% close all

a=1;
max_cfa=1;
anterior_cfa=0;

%while(a<2)
%for t=1:1:7
%imo = aquisicaoim;
im = imread('C:\Programas\MATLAB71\work\Ensaios Mestrado - 05.05.08\Imagens 2ºparte Trabalho\13.04.09 - Teste FFT\Teste Clareira 1.jpg');
%imwrite(imo,'C:\Programas\MATLAB71\work\Ensaios Mestrado - 05.05.08\Imagens 2ºparte Trabalho\07.04.09 - Teste FFT\P280_ok.jpg');
im = im(200:1:899,550:1:1249);
%figure(3),imshow(imo);
figure(4),imshow(im);
hold on;

[srad, sang, S]=specxture(im);

c_lim = 0.2
actual_cfa = sum(srad((size(srad,2)*c_lim:size(srad,2))))
if anterior_cfa < actual_cfa 
    if max_cfa < actual_cfa
        max_cfa = actual_cfa
    end
    if max_cfa == actual_cfa
        text(20,100,'Imagem Focada!','Color','g');
    else
        text(20,100,'A focar imagem....','Color','g');
    end
else
   text(20,100,'A desfocar imagem....','Color','r');
end
anterior_cfa = actual_cfa;
text(20,130,int2str(actual_cfa),'Color','g');
text(20,150,int2str(max_cfa),'Color','b');

fst_h = sang(1:(size(sang,2)/2));
snd_h = sang((size(sang,2)/2)+1:(size(sang,2)));

r = fst_h(1,:)-snd_h(1,:);
mtrcs_fft = sqrt(sum(r(:,1)^2));


text(20,200,'Regularidade Textura:','Color','y');
text(20,230,int2str(mtrcs_fft),'Color','y');

num_pix = sum(srad);
rad_circ = 0;
sum_pix = 0;
j=1;
while (sum_pix < (num_pix*0.67))
    sum_pix = sum_pix + srad(j);
    j=j+1;
    rad_circ=j;
end

text(20,270,'Raio p/ 67% pixeis:','Color','y');
text(20,290,int2str(rad_circ),'Color','y');

% figure(1);
% clf reset;
% colormap gray;
% % luminance image
% subplot(2, 2, 1);
% imshow(im);
% title('imagem original');
% xlabel('x');
% ylabel('y');
% 
% % spectrum
% subplot(2, 2, 2);
% imshow(S);
% axis xy;
% axis square;
% title('representação spectral da frequência');
% xlabel('fx');
% ylabel('fy');
% 
% subplot(2,2,3)
%     plot(srad),
%     title('destribuição; frequência/raio')
% 
% subplot(2,2,4)
%     plot(sang), 
%     title('destribuição; frequência/ângulo')


%end
