数字图象处理实验一

深圳大学实验报告

课程名称：数字图像处理

实验项目名称：数字图像的基本操作与代数运算

学院：信息工程学院

专业：电子信息工程

指导教师：**

报告人：许儒满学号：********** 班级： 1 实验时间：2016/3/7

实验报告提交时间：

教务部制

%实验一

clc;

clear all;

close all;

X1=imread('eight.tif');

subplot(2,2,1);

imshow(X1);

title('eight');

c1=imfinfo('eight.tif')

X2=imread('D:\matleb\bin\exp1_img\flowers.tif'); subplot(2,2,2);

imshow(X2);

title('flowers');

c2=imfinfo('D:\matleb\bin\exp1_img\flowers.tif')

X3=imread('D:\matleb\bin\exp1_img\RGBimageA.tiff'); subplot(2,2,3);

imshow(X3);

title('RGBimageA');

c3=imfinfo('D:\matleb\bin\exp1_img\RGBimageA.tiff')

imshow(X(:,:,1));

title('Rͼ');

subplot(2,2,3)

imshow(X(:,:,2));

title('Gͼ');

subplot(2,2,4)

imshow(X(:,:,3));

title('Bͼ');

3%实验三

clc;

clear all;

close all;

X=imread('D:\matleb\bin\exp1_img\RGBimageA.tiff'); X2=imcrop(X,[100,100,400,425]);

imwrite(X2,'RGBimageA2.tiff')

subplot(1,2,1)

imshow(X);

title('Ô-ͼ');

imshow(X2);

title('¼ô²Ãͼ');

4

%实验四

clc;

clear all;

close all;

X=imread('D:\matleb\bin\exp1_img\RGBimageB.tiff'); imwrite(X,'RGBimageB.jpg','jpg');

5

%实验五

clc;

clear all;

close all;

X=imread('D:\matleb\bin\exp1_img\eight.png');

J=imnoise(X, 'gaussian',0,0.02);

subplot(2,2,1); imshow(X); title('原图像');

subplot(2,2,2); imshow(J); title('加高斯噪声'); [m,n] = size(X);

K = zeros(m, n);

K1=K;K2=K;K3=K;

for i = 1 : 5

J = imnoise(X, 'gaussian', 0, 0.02); J1 = im2double(J);

K1 = K1 + J1;

end

K1= K1/5 ;

for i = 1 : 50

J = imnoise(X, 'gaussian', 0, 0.02); J2 = im2double(J);

K2 = K2 + J2;

end

K2= K2/50 ;

for i = 1 : 500

J = imnoise(X, 'gaussian', 0, 0.02); J3 = im2double(J);

K3= K3 + J3;

end

K3= K3/5 ; %求图像的平均

figure;

subplot(2,2,1)

imshow(K1)

title('5副相加求平均的图像');

subplot(2,2,2)

imshow(K2)

title('50副相加求平均的图像');

subplot(2,2,3)

imshow(K3)

title('500副相加求平均的图像');

由图可知，50副加随机噪声的图相加求平均结果最佳。

6

%实验六

clc;

clear all;

close all;

X=imread('D:\matleb\bin\exp1_img\cat1.jpg'); J=imread('D:\matleb\bin\exp1_img\cat2.jpg'); J1=imresize(J,[240,220]);

K1=imadd(X,J1);

K2=imsubtract(X,J1);

K3=immultiply(X,J1);

subplot(2,3,1)

imshow(X)

title('cat1原图');

subplot(2,3,2)

imshow(J)

title('cat2原图');

subplot(2,3,3)

imshow(J1)

title('cat2缩放');

subplot(2,3,4)

imshow(K1)