my_fft2.m

来自「两个滤波器的代码」· M 代码 · 共 17 行

clear all;
close all;
load slice10.mat;
koko=size(T2_10);
intRows=koko(1);
intCols=koko(2);
Spectrum = fftshift(fft2(T2_10));
HighpassFilter=Fspecial('laplacian');
HighpassFilterSpectrum=fftshift(fft2(HighpassFilter,intRows,intCols));
figure('name','Highpass Filter','NumberTitle','off','MenuBar','none');
colormap(gray);
subplot(1,2,1);
imagesc(log10(1+abs(HighpassFilterSpectrum)));
title('Log Magnitude');
subplot(1,2,2);
imagesc(angle(HighpassFilterSpectrum));
title('Phase');