Created
May 6, 2019 19:39
-
-
Save johnnyferreiradev/1d6a2f06a8f45361489971e4f7209682 to your computer and use it in GitHub Desktop.
Filtro Gaussiano Passa-Altas no domínio da Frequência
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
close all, | |
clear all, | |
clc | |
%Leitura da imagem | |
imgg=imread('lenna.jpg'); | |
%Transformar a imagem dergb para escala de cinza | |
img=rgb2gray(imgg); | |
%Transformar a imagem em double | |
img=im2double(img); | |
%Mostra a imagem | |
figure(1), imshow(img); | |
title('Imagem Original') | |
%Calcula a transformada de fourier da imagem | |
img_f=fft2(img); | |
%Reordena as frequencias (centraliza as frequencias baixas) | |
img_f_ord=fftshift(img_f); | |
%Parte real da transformada de fourier | |
real_img_f_ord=real(img_f_ord); | |
%Magnitude da transformada de fourier | |
abs_img_f_ord=abs(img_f_ord); | |
figure(2), imshow(real_img_f_ord); | |
title('Espectro da Imagem') | |
figure(3),imshow(log(abs_img_f_ord+1),[]); | |
title('Espectro da Imagem - Escala Logaritmica') | |
%Transformada inversa de fourier | |
img_if=ifft2(img_f); | |
figure(4), imshow(img_if); | |
title('Transformada inversa de fourier') | |
%Filtro passa-altas gaussiano | |
%dimensoes da imagem | |
i=1:465; | |
j=1:465; | |
sigma=1; % seção transversal onde ocorre a passagem da onda | |
D0=20; | |
for i=1:465, | |
for j=1:465, | |
w= i-256; | |
k= j-256; | |
D=sqrt(w*w + k*k); | |
filtroPA(i,j)=1-exp(-((D/D0)^2)/(2*sigma*2)); | |
end | |
end | |
figure(5), imshow(filtroPA); | |
title('Filtro passa-baixas') | |
%Multiplicação ponto-a-ponto | |
img_F_PA=filtroPA.*img_f_ord; | |
%Mostra o filtro utilizado | |
figure(6), imshow(real(img_F_PA)); | |
title('Filtro Passa-Altas com o espectro') | |
%Imagem filtrada - Passa-baixas | |
% ifft2 - retorna a transformada de Fourier | |
% inversa discreta bidimensional de uma matriz | |
% usando um algoritmo de transformada rápida de Fourier | |
img_if_PA=ifft2(img_F_PA); | |
figure(7),imshow(real(img_if_PA)); | |
title('Imagem Filtrada- Passa-Altas') | |
%Filtro | |
figure(8),surf(filtroPA); | |
title('Filtro') |
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment