实验二 图像增强

实验二图像增强

一、实验目的

1、掌握图像增强的空域点处理方法原理，会针对不同的图像采用合适的方法进

行增强处理。

2、掌握线性滤波和中值滤波的原理，分析两种空域滤波方法的优缺点。

3、掌握图像增强的频域处理方法，分别适用于处理什么类型的图像。

二、实验内容

1、对一幅图像求反，分析图像中暗色区域的白色或灰色细节有何变化。

2、运行课本例4.1、例4.3的程序，理解灰度变换对图像的影响。

3、对“circuit.tif”图像进行直方图均衡化和规定化处理，调整图像的对比度。

4、对“eight.tif”图像添加不同的噪声，分别使用不同的平滑滤波器和中值滤

波器进行去噪处理，比较两种方法的去噪效果。

5、运行课本例4.14、4.16的程序，理解低通滤波器和同态滤波的原理。

三、实验步骤

1

I=imread('cameraman.tif');

imshow(I);

I=double(I);

I=256-1-I;

I=uint8(I);

figure

imshow(I);

2

(1)

I=imread('pout.tif');

imshow(I);

I=double(I);

[M,N]=size(I);

for i=1:M

for j=1:N

if I(i,j)<=30

I(i,j)=I(i,j);

elseif I(i,j)<=150

I(i,j)=(200-30)/(150-30)*(I(i,j)-30)+30;

else

I(i,j)=(255-200)/(255-150)*(I(i,j)-150)+200; end

end

end

figure(2);

imshow(uint8(I));

(2)

I=imread('concordorthophoto.png'); figure;

imshow(I);

I=double(I);

[M,N]=size(I);

for i=1:M

for j=1:N

if I(i,j)<=50

I(i,j)=40;

elseif I(i,j)<=180

I(i,j)=220;

else

I(i,j)=40;

end

end

end

I=uint8(I);

figure;

imshow(I);

3

（1）

I=imread('circuit.tif'); figure

subplot(221);

imshow(I);

subplot(222);

imhist(I)

I1=histeq(I);

figure;

subplot(221);

imshow(I1)

subplot(222);

imhist(I1)

（2）

I=imread('circuit.tif');

[M,N]=size(I);

for i=1:8:257

counts(i)=i;

end

Q=imread('circuit.tif');

N=histeq(Q,counts);

figure

subplot(221);

imshow(N);

subplot(222);

imhist(N);

axis([0 260 0 5000]);

4

（1）

I=imread('eight.tif');

J=imnoise(I,'salt & pepper',0.02);

subplot(231),imshow(I); title('原图像');

subplot(232),imshow(J); title('添加椒盐噪声图像'); k1=filter2(fspecial('average',3),J);

k2=filter2(fspecial('average',5),J);

k3=filter2(fspecial('average',7),J);

k4=filter2(fspecial('average',9),J);

subplot(233); imshow(uint8(k1)); title('3*3模板平滑滤波'); subplot(234); imshow(uint8(k2)); title('5*5模板平滑滤波'); subplot(235); imshow(uint8(k3)); title('7*7模板平滑滤波'); subplot(236); imshow(uint8(k4)); title('9*9模板平滑滤波');

（2）

I=imread('eight.tif');

J=imnoise(I,'salt & pepper',0.02);

subplot(231),imshow(I); title('原图像');

subplot(232),imshow(J); title('添加椒盐噪声图像');

k1=medfilt2(J);

k2=medfilt2(J,[5 5]);

k3=medfilt2(J,[7 7]);

k4=medfilt2(J,[9 9]);

subplot(233); imshow(uint8(k1)); title('3*3模板中值滤波'); subplot(234); imshow(uint8(k2)); title('5*5模板中值滤波'); subplot(235); imshow(uint8(k3)); title('7*7模板中值滤波'); subplot(236); imshow(uint8(k4)); title('9*9模板中值滤波');