实验5 锐化滤波器、傅立叶变换及简单的频域滤波操作

最新推荐文章于 2022-08-24 20:12:23 发布

遇见生活

最新推荐文章于 2022-08-24 20:12:23 发布

阅读量569

点赞数

分类专栏：数字图像处理

本文链接：https://blog.csdn.net/qq_53625091/article/details/115519855

版权

数字图像处理专栏收录该内容

8 篇文章 3 订阅

订阅专栏

本文通过实例展示了图像处理的一些基本操作，包括拉普拉斯滤波、傅立叶变换及其应用。首先，对图像进行拉普拉斯滤波以增强边缘；接着，对图像进行傅立叶变换，展示原图、傅立叶谱及移中后的谱；然后，进一步探讨傅立叶变换在相位谱和对数变换中的应用；最后，进行了逆傅立叶变换以观察其对图像的影响。此外，还介绍了使用空间滤波器和频率域滤波的方法。

摘要由CSDN通过智能技术生成

1. 读入moon.tif图像，进行拉普拉斯滤波操作，使用fspecial函数定义模板，要注意使用im2double函数与图像减操作，理解拉普拉斯滤波的作法。

f = imread('miaoa.jpg');
f2 = im2double(f) ;
w = fspecial('laplacian',0);
g1 = imfilter(f,w,'replicate');
g2 = imfilter(f2,w,'replicate');
subplot(131);
imshow(g1);
subplot(132);
imshow(g2);
g = f2 - g2;
subplot(133);
imshow(g);

2. 读入一幅自己拍摄的图像，转化为灰度图，进行傅立叶变换并显示，要求显示原图、傅立叶谱及移中后的傅立叶谱。

f = imread('miaoa.jpg');
f1 = rgb2gray(f);
subplot(1,3,1);
imshow(f);
title('原图');
F = fft2(f1);
S = abs(F);
subplot(1,3,2);
imshow(S,[]);
title('傅立叶谱');
Fc = fftshift(F);
subplot(1,3,3);
imshow(abs(Fc),[]);
title('变换后的傅立叶谱');

3.读入图像f1.tif，进行傅立叶变换并显示，要求显示原图、傅立叶谱、移中后的傅立叶谱及利用对数变换改进的移中后的傅立叶谱。

f = imread('miaoa.jpg');
f1 = rgb2gray(f);
subplot(2,2,1);
imshow(f);
title('原图');
F = fft2(f1);
S = abs(F);
subplot(2,2,2);
imshow(S,[]);
title('傅立叶谱');
Fc = fftshift(F);
subplot(2,2,3);
imshow(abs(Fc),[]);
title('变换后的傅立叶谱');
S2 = log(1 + abs(Fc));
subplot(2,2,4);
imshow(S2,[]);
title('对数变换后的傅立叶谱');

4.读入图像f1.tif，进行傅立叶变换并计处相位谱，要求显示原图、傅立叶谱及相位谱

f = imread('miaoa.jpg');
f1 = rgb2gray(f);
subplot(2,2,1);
imshow(f);
title('原图');
F = fft2(f1);
S = abs(F);
subplot(2,2,2);
imshow(S,[]);
title('傅立叶谱');
Fc = angle(F);
subplot(2,2,3);
imshow(abs(Fc),[]);
title('相位谱');

5.读入一幅自己拍摄的图像，转化为灰度图，进行傅立叶变换并显示，然后再进行逆傅立叶变换并显示，要求显示原图、傅立叶谱及逆傅立叶变换后的图，对比原图与逆傅立叶之后的图的区别，加深理解图像变换的意义

f = imread('miaoa.jpg');
f1 = rgb2gray(f);
subplot(2,2,1);
imshow(f);
title('原图');
F = fft2(f1);
F2 = ifft2(F);
S = abs(F);
subplot(2,2,2);
imshow(S,[]);
title('傅立叶谱');
subplot(2,2,3);
imshow(abs(F2),[]);
title('逆傅立叶');