лаб6 / lab6_1

I = imread('candle.jpg');
if (size(I,3) == 3) I = rgb2gray(I); end
I = im2double(I);
S = fftshift(fft2(I)); % 2D-БПФ изображения, центрирован
A = abs(S); % амплитудный спектр
A_max = max(max(A));
eps = A_max*10^(-6); %добавка, чтобы не вычислить lg(0)
A_dB = 20*log10(A + eps); % перевод в дБ
figure('Name','Спектр изображения, в псевдоцвете');
imagesc(A_dB) % вывод двумерного спектра изображения
colorbar