clc
clear all
close all

% imsint1=zeros(100);
% 
% imsint1 = imnoise(imsint1,'gaussian',0,0.1);
%  for i=1:2:10 
%      imsint1((40+i):(60-i),(40+i):(60-i)) = 250;
%      imsint1((40+(i+1)):(60-(i+1)),(40+(i+1)):(60-(i+1))) = 0;
%  end

% imsint1(25:30,25:30)= 1;

%im = imsint1;
a=1

%while(a<2)
%imo = aquisicaoim;
imo = imread('10.03.09 - Darkfield\P800\100309 L AGUA RL P800 c.jpg');

im = imo(200:1:899,550:1:1249);
figure(3),imshow(imo);
figure(4),imshow(im);

PQ = paddedsize(size(im));
 
S = fftshift(fft2(im));
S = abs(S);
S = mat2gray(log(1+S));

%figure, imshow(S);
[srad, sang, Si]=specxture(im);
%figure(2), imshow(im)
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 
% % filtered image
% subplot(2, 2, 4);
% imagesc(L_filtered);
% colorbar;
% axis square;
% set(gca, 'TickDir', 'out');
% title('filtered image');
% xlabel('x');
% ylabel('y');
