|

楼主 |
发表于 2025-4-15 20:41:48
|
显示全部楼层
第5步
I = imread('mohu.jpg');
%灰度化
I_gray = rgb2gray(I);
subplot(2,3,1);imshow(I_gray);title('模糊');
%傅里叶变换
F = fft2(I_gray);
Fc = fftshift(F);
subplot(2,3,2);imshow(log(abs(Fc)+1),[]);
title('模糊频谱');
[m,n]=size(I_gray);
core = zeros(m,n);
for i = 1:m
for j = 1:n
core(i,j) = (abs(i-m/2)+abs(j-n/2))/40+1;
end
end
subplot(2,3,3);imshow(core,[]);
Fc1 = Fc.*core;
subplot(2,3,4);imshow(log(abs(Fc1)+1),[]);
title('新的频谱');
F1 = ifftshift(Fc1);
I1_gray = abs(ifft2(F1));
subplot(2,3,5);imshow(uint8(I1_gray));title('清晰图像');
|
|