用matlab实现图像灰度变换课程设计报告

1、- .-. -可修编- . z.- . -可修编- .课程设计报告册课程名称： MATLAB课程设计 课题名称： 灰度变换增强 专业班级：姓 名： Bob Wang 学 号： 15164 课程设计主要场所： 信息楼220 时间：指导教师：成 绩：前言数字图像处理技术是20世界60年代开展起来的一门新兴学科，随着图像处理理论和方法的进一步完善，使得数字图像处理技术在各个领域得到了广泛应用，并显示出广阔的应用前景。MATLAB既是一种直观、高效的计算机语言，同时又是一个科学计算平台。它为数据分析和数据可视化、算法和应用程序开发提供了最核心的数学和高级图形工具。根据它提供的500多个数学和工程函数，

2、工程技术人员和科学工作者可以在它的集成环境互或变成以完成各自的计算。MATLAB中集成了功能强大的图像处理工具箱。由于MATLAB语言的语法特征与C语言极为相似，而且更加简单，更加符合科技人员对数学表达式的书写格式，而且这种语言可移植性好、可扩展性强，再加上其中有丰富的图像处理函数，所以MATLAB在图像处理的应用中具有很大的优势。MATLAB是一种以矩阵运算为根底的交互式程序语言，能够满足科学、工程计算和绘图的要求，与其它计算机语言相比，其特点是简洁和智能化，适应科技专业人员的思维方式和书写习惯，使得编程和调试效率大大提高。我们学习掌握MATLAB，也可以说是在科学工具上与国际接轨。目录课程

3、设计目的2设计任务及容2课题设计实验条件3涉及知识3具体设计过程及调试45.1、图像的读入和显示5.1.1、翻开图像45.1.2、显示原图像55.1.3、图像灰度处理75.1.4、显示灰阶后图像8 5.2、直方图均衡化 5.2.1、生成直方图.10 5.2.2、直方图均衡化.125.3、灰度变换5.3.1、线性变换.95.3.2、分段线性变换.95.3.3、非线性变换 .9六、心得体会.17七、参考文献.18八、程序清单.19-. z.设计目的本次课程设计旨在提升学生的动手能力，加强并稳固对专业理论知识的理解和实际运用。通过使用MATLAB软件，根据所分配的课题，锻炼我们的自学能力，提高我们实

4、践能力，为毕业设计做铺垫。1、熟悉并掌握应用MATLAB语言相关知识和程序设计方法。2、掌握MATLAB GUI程序设计3、熟悉MATLAB图像处理工具箱4、学会运用MATLAB工具箱对图像进展处理和分析设计任务及容利用MATLAB GUI设计实现图像处理的图形用户界面，利用MATLAB图像处理工具箱实现以下的图像处理功能：采用matlab Guide设计图像界面；可以读入并显示图像，设计实现以下功能的按钮，点击相应的按钮可以实现相应的功能。自选图像一幅，进展线性变换、分段线性变换和非线性变换，比照源图像与变换后的图像特征变换，体会该变换的特点。自选图像一幅，生成并绘制图像的直方图。直方图均衡

5、化方法研究，比拟变换前后图像的直方图。题目研究的实验条件计算机、相关书籍、MATLAB7以上版本课题所涉及知识面数字图像处理和MATLAB软件知识具体设计过程及调试、图像的读入和显示采用matlab Guide设计图像界面；可以读入并显示图像，设计实现以下功能的按钮，点击相应的按钮可以实现相应的功能。图形界面5.1.1、 翻开图像从计算机中查找要翻开的图片 5.1.2、显示原图像 把选中的图片显示在第一显示区 5.1.3、图像灰度处理 利用灰度处理函数rgb2gray()对原真彩图像灰阶化 5.1.4、灰阶后的图像显示 把灰阶后的图像显示在第二显示区，对两幅图片进展比拟，观察两者特点。显示原图

6、像和灰阶后的图像结论：对图像灰阶处理就是对原图像黑白化5.2、直方图均衡化直方图均衡化是一种常用的灰度增强算法，是将原图的直方图经过变换函数修正为均匀的直方图，然后按照均衡后直方图修整原图像。5.2.1、生成直方图图像的直方图是表示数字图像中每一灰度级与该灰度级出现的频数，函数为imhist(*,map)显示原图像的灰度直方图5.2.2、直方图均衡化均衡化函数为histeq(I,n),其中n默认为64直方图均衡化比拟5.3、灰度变换*些情况下，需要对图像的灰度级整个围或者其中的*一段扩展或压缩到记录器件输入灰度级动态围之。灰度变换可分为线性变换、分段线性变换、非线性变换和其他的灰度变换。5.3

7、.1、线性变换具体实现形式为：g(*,y)=(d-c)/(b-a)f(*,y)-a+c线性变换结果比照作用：使曝光不充分图像中黑的更黑，白得更白，从而提高图像比照度5.3.2、分段线性变换具体实现形式：g(*,y)=c/df(*,y) 0=f(*,y)=ag(*,y)= (d-c)/(b-a)f(*,y)-a+c a=f(*,y)=bg(*,y)= (f-b)/(e-a)f(*,y)-b+d b=f(*,y)=c分段线性变换结果比照5.3.3、非线性变换具体实现形式：g(*,y)=a+lnf(*,y)+1/b*lnc非线性变换结果比照六、心得体会这次MATLAB课程设计让我学到了很多，不仅仅是

8、知识，还学到了许多学习方法。MATLAB软件在这之前接触的很少，仅仅是在数字信号处理的实验课上有过一点了解。为了做好这次的课程设计，我查阅了大量的资料，并上网搜索了许多与此相关的知识，这个过程使我受益匪浅，那就认识了很多关于MATLAB的论坛，论坛上面提供了很多关于MATLAB学习的方法及其自学的方法，那些论坛对于将来学习MATLAB及其它方面的知识也有很大的帮助。虽然在规定的时间完成了教师所布置的课题，但是还是觉得不免有些缺乏。由于时间局限，对许多MATLAB所运用到的知识仅仅是有所了解，而并没有完全地掌握，甚至可以说是只知其然，而不知其所以然。这使我意识到对于MATLAB的学习任务还很艰巨

9、，还有许多知识需要去了解，去深入研究。通过这次MATLAB课设，学到了MATLAB GUI程序设计，并通过MATLAB所生成的图像调试结果， 以到达系统的要求。MATLAB作为第四代计算机语言，具有相当强大的矩阵运算和操作功能，其程序编写也比其它语言要简单，是一种对学习非常有帮助的软件，在之后的时间，我依然会对它进展学习。七、参考文献郭仕剑、邱志模、陆静芳.MATLAB入门与实践.：人民邮电，2008智星.MATLAB程序设计及应用.:清华大学,2002肖伟、忠.MATLAB程序设计与应用M.：清华大学，2005瞿亮.基于MATLAB的控制系统计算机仿真.：清华大学，交通大学.2005八、程序

10、清单一 翻开图片和灰阶化global imglobal *global yglobal z*=0.002;y=0.02;z=0.04;filename,pathname=. uigetfile(); str=pathname filename; im=imread(str); a*es(handles.a*es1); imshow(im); title(原图); im = rgb2gray(im); a*es(handles.a*es2); imshow(im);title(); 二 线性变换global im;global J;J=imadjust(im,0.3,0.7,);a*es(han

11、dles.a*es1);imshow(im);title(原图);a*es(handles.a*es2);imhist(im);title(灰度直方图);a*es(handles.a*es3);imshow(J);title(线性变换后图像);a*es(handles.a*es4);imhist(J);title(线性变换后灰度直方图);三 分段线性变换global im;global H;H=double(im); M,N=size(H);%for i=1:Mfor j=1:Nif H(i,j)=30 H(i,j)=H(i,j);elseif im(i,j)=150 H(i,j)=(200-

12、30)/(150-30)*(H(i,j)-30)+30;else H(i,j)=(255-200)/(255-150)*(H(i,j)-150)+200;endendend%a*es(handles.a*es1);imshow(im);title();a*es(handles.a*es2);imhist(im);title();a*es(handles.a*es3);imshow(uint8(H);title();a*es(handles.a*es4);imhist(uint8(H);title();四 非线性变换global im;global J;global H;J=double (im

13、) ;H=(log(J+1)/10;a*es(handles.a*es1);imshow(im);title();a*es(handles.a*es2);imhist(im);title();a*es(handles.a*es3);imshow(H);title();a*es(handles.a*es4);imhist(H);title();五 生成灰度直方图global im;a*es(handles.a*es1);imshow(im);title(原图);a*es(handles.a*es2);imhist(im);title(灰度直方图);六 直方图均衡化global im;global

14、 J;J=histeq(im);a*es(handles.a*es1);imshow(im);title(原图);a*es(handles.a*es2);imshow(J);title(直方图均衡化后图像);a*es(handles.a*es3);imhist(im);title(原图像直方图);a*es(handles.a*es4);imhist(J);title(直方图均衡化后图像直方图);平滑处理用3*3屏蔽窗口的8近邻均值进展滤波for(int j=1;jheight-1;j+)for(int i=1;i127.5)p_tempj*wide+i=averg;利用巴特沃斯Butterworth低通滤波器对受噪声干扰的图像进展平滑处理I=imread(aaa.jpg);imshow(I);J1=imnoise(I,salt & pepper); % 叠加椒盐噪声figure,imshow(J1);f=double(J1); % 数据类型转换，MATLAB不支持图像的无符号整型的计算g=fft2(f); % 傅立叶变换g=fftshift(g); % 转换数据矩阵M,N=size(g);nn=2; % 二阶巴特沃斯(Butterworth)低通滤波器d0=