数字图像处理上机报告.doc

数字图像处理上机报告 DIGITAL IMAGE PROCESSING 指导老师： 院 系： 班 级： 姓 名： 学 号： 一、DIP第一次上机实验 本次上机实验内容为灰度图象空域处理。 题目及要求： 0. 显示灰度图象p1~p6及直方图； 6. 对灰度图像p6的偏黑部分实施灰度展宽，对偏亮的部分实施灰度压缩； 编程： 0. 显示灰度图像p1~p6 方法1： I=imread(p1.bmp); subplot(3,2,1); imshow(I); I=imread(p2.bmp); subplot(3,2,2); imshow(I); I=imread(p3.bmp); subplot(3,2,3); imshow(I); I=imread(p4.bmp); subplot(3,2,4); imshow(I); I=imread(p5.bmp); subplot(3,2,5); imshow(I); I=imread(p6.bmp); subplot(3,2,6); imshow(I); 方法2： subplot(231);imshow(p1.bmp); subplot(232);imshow(p2.bmp); subplot(233);imshow(p3.bmp); subplot(234);imshow(p4.bmp); subplot(235);imshow(p5.bmp); subplot(236);imshow(p6.bmp); 方法1与方法2效果一样，均是读取p1~p6这6副图像并显示，Matlab运行结果如下图1.1所示： 图1.1 读取6副图像并显示 显示灰度图像直方图 方法1： I=imread(p1.bmp); whos; I1=rgb2gray(I);%灰度转换 subplot(2,2,1); title(原始图像); imshow(I1); subplot(2,2,2); title(图象的直方图,64); imhist(I1,64); subplot(2,2,3); title(图象的直方图,128); imhist(I1,128); subplot(2,2,4); title(图象的直方图,256); imhist(I1); 图1.2读取RGB图像并显示其直方图 方法2： I=imread(p1.bmp); subplot(2,2,1); title(原始图像); imshow(I); subplot(2,2,2); title(图象的直方图,64); imhist(I,64); subplot(2,2,3); title(图象的直方图,128); imhist(I,128); subplot(2,2,4); title(图象的直方图,256); imhist(I); 方法1为读取RGB图像并显示其直方图，方法2均是读取灰度图像并显示其直方图，一开始,绘制直方图报错,Function IMHIST expected its first input, I or X, to be two-dimensional.原来是RGB的,实际是三个二维矩阵了，而直方图函数只接受一个二维矩阵，所以,需要首先将RGB转成GRAY,单色灰度图像。 Matlab运行结果如下图1.2及图1.3所示： 图1.3 读取灰度图像并显示其直方图 6. 对灰度图像p6的偏黑部分实施灰度展宽，对偏亮的部分实施灰度压缩； 方法1：利用倍数灰度延展及压缩方法实现 a=imread(p6.bmp); a=double(a); b=size(a); for i=1:b(1) for j=1:b(2) if a(i,j)127 a1(i,j)=2.5*a(i,j); else a1(i,j)=0.8*a(i,j); end end end a=uint8(a); a1=uint8(a1); subplot(1,2,1); imshow(a); title(original photo); subplot(1,2,2); imshow(a1); title(photo processed); 效果图如下图1.3： 图1.4 方法2：以对数函数进行灰度值压缩扩展（减小了图像动态范围，效果太差） a=imread(p6.bmp); a=double(a); b=size(a); for i=1:b(1) for j=1:b(2) if a(i