A Figura 1. (c) mostra a imagem de teste reconstruída apenas do espectro MAGNITUDE. Podemos dizer que os valores de intensidade dos pixels de baixa frequência são comparativamente mais altos que os pixels de alta frequência.
A Figura 1. (d) mostra a imagem de teste reconstruída apenas do espectro PHASE. Podemos dizer que os valores de intensidade dos pixels de alta frequência (bordas, linhas) são comparativamente mais altos que os pixels de baixa frequência.
Por que essa contradição mágica da mudança (ou troca) de intensidade está presente entre a imagem de teste reconstruída apenas no espectro MAGNITUDE e a imagem de teste reconstruída apenas no espectro PHASE, que quando combinadas formam a imagem de teste original?
clc;
clear all;
close all;
i1=imread('C:\Users\Admin\Desktop\rough\Capture1.png');
i1=rgb2gray(i1);
f1=fftn(i1);
mag1=abs(f1);
s=log(1+fftshift(f1));
phase1=angle(f1);
r1=ifftshift(ifftn(mag1));
r2=ifftn(exp(1i*phase1));
figure,imshow(i1);
figure,imshow(s,[]);
figure,imshow(uint8(r1));
figure,imshow(r2,[]);
r2=histeq(r2);
r3=histeq(uint8(r2));
figure,imshow(r2);
figure,imshow(r3);
Na sua linha,
mag1=abs(f1);
você mantém a intensidade total da imagem inalterada (teste isso somando as intensidades de todos os pixels). Rejeitar as informações da fase no espaço de Fourier apenas leva a uma redistribuição espacial da intensidade no espaço real, de modo que r1 terá a mesma integridade total que i1.Na sua linha,
phase1=angle(f1);
você está normalizando as amplitudes de cada pixel (no espaço de Fourier) para 1, para que a intensidade total da imagem seja alterada. Como a fase carrega grande parte das informações espaciais da imagem, os principais recursos da imagem são preservados.fonte