数字图像处理91280new

PLC课程报告学号:200711710专业：电子信息工程日期：2010.12.07基于MATLAB的图像处理MATLAB中集成了功能强大的图像处理工具箱。由于MATLAB语言的语法特征与C语言极为相似，而且更加简单，更加符合科技人员对数学表达式的书写格式，而且这种语言可移植性好、可扩展性强，再加上其中有丰富的图像处理函数，所以MATLAB在图像处理的应用中具有很大的优势。

灰度由于RGB图像是三维图像，所以图像数据是一个三维数组，为了显示灰度图像，把三维图像降为二维，可以只取其中的二维数据，实现方法程序为：I=imread(‘zw1.jpg');[d1,d2,d3]=size(I);if(d3>1)I=rgb2gray(I);%如果是灰度图就不用先变换endI=double(I)/255;I1=uint8(255*I*0.5+0.5);imshow(I1);imwrite(I,'test.jpg')二值图像threshold=graythresh(GRAY);（二值图像）BW=im2bw(GRAY,threshold);figure,imshow(BW),title(‘BinaryImage’);%二值图像利用imshow函数显示图像。源代码：I1=imread('zw.bmp');I2=imread('sky.jpg');I3=imread('zw1.jpg');I4=imread('zw2.jpg');subplot(2,2,1),imshow(I1);subplot(2,2,2),imshow(I2);subplot(2,2,3),imshow(I3);subplot(2,2,4),imshow(I4);图像平移的试验程序和显示结果

I=imread('zw1.jpg');

figure,imshow(I);

I=double(I);

I_movesult=zeros(size(I));

H=size(I);

Move_x=50;

Move_y=50;

I_movesult(Move_x+1:H(1),Move_y+1:H(2),1:H(3))=I(1:H(1)-Move_x,1:H(2)-Move_y,1:H(3));

imshow(uint8(I_movesult));垂直、水平和对角镜像的程序如下：其中，I1为原始图像；I2为垂直图像I3为水平镜像；I4为对角镜像I1=imread('zw1.jpg');I1=double(I1);figure(1),imshow(uint8(I1));H=size(I1);figure(2),I2(1:H(1),1:H(2),1:H(3))=I1(H(1):-1:1,1:H(2),1:H(3));imshow(uint8(I2));figure(3),I3(1:H(1),1:H(2),1:H(3))=I1(1:H(1),H(2):-1:1,1:H(3));imshow(uint8(I3));figure(4),I4(1:H(1),1:H(2),1:H(3))=I1(H(1):-1:1,H(2):-1:1,1:H(3));imshow(uint8(I4));各种噪声中的图像变换

（1）对给定的图像添加均值为0，方差分别为200、400的高斯噪声，以及概率分别为0.1、0.2的椒盐噪声，显示添加噪声后的图像。

源代码：

a)加入噪声强度为d=0.2的椒盐噪声

I=imread('zw1.jpg');

J=imnoise(I,'salt&pepper',0.2);%给图像加入椒盐噪声

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

subplot(1,2,2);imshow(J);title('加入椒盐噪声的图像')

图像显示：b)加入噪声强度为d=0.1的椒盐噪声：

I=imread('zw1.jpg');

J=imnoise(I,'salt&pepper',0.1);%给图像加入椒盐噪声

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

subplot(1,2,2);imshow(J);title('加入椒盐噪声的图像')

图像显示：c)加入均值m=0,方差v=200的高斯噪声：I=imread('zw1.jpg');J=imnoise(I,'gaussian',0,200);%给图像加入高斯噪声subplot(1,2,1);imshow(I);title('原始图像');subplot(1,2,2);imshow(J);title('加入高斯噪声图像')图像显示：d)加入均值m=0,方差v=400的高斯噪声：

I=imread('zw1.jpg');

J=imnoise(I,'gaussian',0,400);%给图像加入高斯噪声

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

subplot(1,2,2);imshow(J);title('加入高斯噪声图像')

图像显示：）对加噪图像分别进行3×3、5×5、9×9的均值滤波，显示并比较滤波结果。

源代码：

I=imread('zw1.jpg');

J=imnoise(I,'salt&pepper',0.2);%给图像加入椒盐噪声

h=fspecial('average',3);

s=fspecial('average',5);

p=fspecial('average',9);

gd=imfilter(J,h);

qw=imfilter(J,s);

er=imfilter(J,p);

subplot(2,2,1);imshow(J);title('加入椒盐噪声的图像');

su