Exercício em Sala

Prévia do material em texto

Exercício em Sala
IMAGENS MÉDICAS 2
Profª. Ana Claudia Patrocínio
Maria Vitória Garcia - 11821EBI013
Uberlândia, 15 de Maio de 2023
http://www.feelt.ufu.br/pessoas/docentes/eduardo-lazaro-martins-naves
Código
clear all;
clc;
pkg load image;
img_in = imread('mamo.png');
subplot(2, 3, 1);
imshow(img_in);
title('Imagem Original');
% Filtro de média 3x3
mask3x3 = fspecial('average', 3);
img_out3x3 = imfilter(img_in, mask3x3);
subplot(2, 3, 2);
imshow(img_out3x3);
title('Filtro de média 3x3');
% Imagem Equalizada
img_equalizada = histeq(img_out3x3);
subplot(2, 3, 3);
imshow(img_equalizada);
title('Equalização do Histograma');
% Construção do filtro ideal
[M, N] = size(img_in);
fcorte = 50;
filtropbi = zeros(M, N);
cx = floor(M/2) + 1;
cy = floor(N/2) + 1;
for x = 1:M
for y = 1:N
D = sqrt((x - cx)^2 + (y - cy)^2);
if D < fcorte
filtropbi(x, y) = 1;
end
end
end
% Transformada de Fourier - espectro de frequências
DFT = fft2(img_in);
% Centralização do espectro de frequência
DFTC = fftshift(DFT);
% Multiplicação do filtro pelo espectro centralizado
GC = filtropbi .* DFTC;
% Descentralização do espectro resultante da operação acima
G = ifftshift(GC);
% Transformada inversa do espectro resultante descentralizado
img_filtrada = uint8(real(ifft2(G)));
subplot(2, 3, 4);
imshow(img_filtrada);
title('Passa Baixa Ideal');
% Imagem Equalizada
img_equalizadapbi = histeq(img_filtrada);
subplot(2, 3, 5);
imshow(img_equalizadapbi);
title('Equalização do Histograma P.B.I');
Resultados
O filtro de média 3x3 deixou a imagem mais realçada em alguns pontos do
tecido, o que pode ser favorável para visualizações de algumas áreas de interesse.
Posteriormente a primeira equalização e o passa baixa ideal a imagem já sofreu
mais alterações como nas bordas e na resolução, e a equalização do passa baixa
ideal (P.B.I) já não trouxe resultados favoráveis pois diminuiu contraste de cores
perdendo assim informações da imagem.