【实验名称】

图像锐化处理技术

?【实验目的】

1、通过本次实验掌握图像锐化技术；

2、通过本次实验掌握图像高频率通过滤波技术；

3、熟悉使用MATLAB库函数fft2, ifft2, fftshift, ifftshift, imfilter；

?【实验内容】

图1：

图2

?

【实验代码】

?题1

I=imread('peppers.png');
H=[1,1,1:1,-8,1:1,1,1];
w=0.8;
subplot(1,2,1);
imshow(I);
title('原图');
I1=imfilter(I,H);
I2=I-1*I1;
subplot(1,2,2);
imshow(I2);
title('锐化后');

题2

clc;clear all;close all;
%  读取原图像
Img1=imread('peppers.png');
figure
subplot(1,3,1);
imshow(Img1);
title('原图');
Img2=rgb2gray(Img1);
subplot(1,3,2);
imshow(Img2);
title('灰度图');
[m,n]=size(Img2);
mc=floor(m/2)+1;
nc=floor(n/2)+1;
[X,Y]=meshgrid(1:n,1:m);
D=sqrt((Y-mc).^2+(X-nc).^2);
D0=20;
Hhp=1-exp(-D.^2/(2*D0^2));
A=6;
H=(A-1)+Hhp;
% 傅里叶变换
G=fft2(Img2);
% 中心化
Gc=fftshift(G);
% 滤波处理（点乘）
GG=Gc.*H;
% 反中心化
GGG=ifftshift(GG);
% 逆傅里叶变换
ggg=ifft2(GGG);
% 取实部
gre=real(ggg);
% 规格化
gscale=255*gre/max(gre(:));
% 转换类型
g5=uint8(gscale);
subplot(1,3,3);
imshow(g5);
title('频域锐化处理');

?

?【实验结果】

题1

题2