1、1目录一、 课程设计目的.3二、 课程设计要求.3三、 课程设计的内容.3四、 题目分析.3五、 总体设计.4六、 具体设计.51、文件.51.1、打开.51.2、保存.51.3、退出.52、编辑56.2.1、灰度.56.2.2、亮度.66.2.3、截图.76.2.4、缩放.73、旋转96.3.1、上下翻转.96.3.2、左右翻转.96.3.3任意角度翻转.96.4、噪声106.5、滤波106.6、直方图统计116.7、频谱分析126.7.1、频谱图 .126.7.2、通过高通滤波器 .126.7.3、通过低通滤波器 .136.8、灰度图像处理 . .146.8.1、二值图像 .146.8.2
2、、创建索引图像 .146.9、颜色模型转换 146.10、操作界面设计15七、 程序调试及结果分析.15八、 心得体会.16九、 参考文献.17十、 附录.182基于 MATLAB 的图像处理的课程设计摘要： 数字图像处理技术是 20世纪 60年代发展起来的一门新兴学科，随着图像处理理论和方法的进一步完善，使得数字图像处理技术在各个领域得到了广泛应用，并显示出广阔的应用前景。MATLAB 既是一种直观、高效的计算机语言,同时又是一个科学计算平台。它为数据分析和数据可视化、算法和应用程序开发提供了最核心的数学和高级图形工具。根据它提供的 500多个数学和工程函数，工程技术人员和科学工作者可以在它
3、的集成环境中交互或编程以完成各自的计算。MATLAB中集成了功能强大的图像处理工具箱。由于 MATLAB语言的语法特征与C语言极为相似，而且更加简单，更加符合科技人员对数学表达式的书写格式，而且这种语言可移植性好、可扩展性强，再加上其中有丰富的图像处理函数，所以 MATLAB在图像处理的应用中具有很大的优势。本文从预处理、边缘检测、车牌定位、字符分割、字符识别五个方面，具体介绍了车牌自动识别的原理。并用 MATLAB软件编程来实现每一个部分,最后识别出汽车牌照。关键词：MATLAB，数字图像处理，车牌识别系统3一、课程设计目的Matlab技术课程设计是利用 Matlab应用软件的基础知识和基本
4、程序设计方法，解决选修专业课程中数字信号处理的实际应用问题，从而加深对理论知识的掌握，并把所学的知识系统、高效的贯穿到实践中来，避免理论与实践的脱离。同时提高独立编程水平，并在实践中不断完善理论基础，有助于培养综合能力。二、课程设计的内容学习 MATLAB程序设计，利用 MATLAB图像处理工具箱，设计和实现通过设计一个车牌识别系统三、课程设计的要求：按照软件工程方法，根据需求进行程序的功能分析和界面设计，给出设计详细说明。然后按照自己拟定的功能要求进行程序设计和调试。 以下几点是程序必须实现的功能。1）设计图像用户界面，让用户能够对图像进行调整，显示和对比变换前后的图像。2）图像直方图统计和
5、直方图均衡，要求显示直方图统计，比较直方图均衡后的效果。3）显示识别结果四、题目分析（一）研究背景：随着我国经济的发展，综合实力的增强和城市化进程的加快，国内各大城市交通管理能力将面临重大考验。近年来，各主要城市都将城市智能交通系统（ITS）的建设作为改善城市交通状况的重点，以缓解城市交通管理的压力，车牌识别系统作为智能交通系统的主要技术之一也得到了很大的发展。（二）现状目前，国内的 ITS 产业仍然处于起步阶段，在产业链中，电子地图商、导航终端商、嵌入式操作系统均已经具备产业规模，而所欠缺的环节则是对实时交通信息的采集和处理系统。车牌识别系统是城市智能交通管理系统的一个重要组成部分，广泛用于
6、是针对公路行驶的机动车辆进行实时监控的智能化交通系统，电子收费，交通违规管理，安全停车管理等重要领域。 4五、总体设计车牌识是一辆汽车独一无二的信息，因此，对车辆牌照的识别技术可以作为辨识一样车最为有效地方法。车牌识别系统包括摄取的汽车图像，车牌号码的识别，车牌图像的采集和预处理，牌照区域的定位和提取，牌照字符的分割和识别等几个部分组成，如下图所示。图1 车牌识别系统流程图六、具体设计界面效果图输入图像 预处理 车牌定位字符识别输出结果 归一化字符分割5（1）灰度校正由于牌照图象在拍摄时受到种种条件的限制和干扰，图象的灰度值往往与实际景物不完全匹配，这将直接影响到图象的后续处理。如果造成这种影
7、响的原因主要是由于被摄物体的远近不同，使得图象中央区域和边缘区域的灰度失衡，或是由于摄像头在扫描时各点的灵敏度有较大的差异而产生图象灰度失真，或是由于曝光不足而使得图像的灰度变化范围很窄。这时就可以采用灰度校正的方法来处理，增强灰度的变化范围、丰富灰度层次，以达到增强图象的对比6度和分辨率。我们发现车辆牌照图象的灰度取值范围大多局限在 r=(50,200)间，而且总体上灰度偏低，图象较暗图6 灰度增强后的图像（2）平滑处理对于受噪声干扰严重的图象，由于噪声点多在频域中映射为高频分量，因此可以在通过低通滤波器来滤除噪声，但实际中为了简化算法，也可以直接在空域中用求邻域平均值的方法来削弱噪声的影响
8、，这种方法称为图象平滑处理。例如，某一象素点的邻域S 有两种表示方法：8邻域和4邻域分别对应的邻域平均值为， 。其中，M 为邻域中除中心象素点f(i,j) sjijfjig),(,1),(之外包括的其它象素总数，对于4邻M=4，8 邻域M=8。然而，邻域平均值的平滑处理会使得图象灰度急剧变化的地方，尤其是物体边缘区域和字符轮廓等部分产生模糊作用。为了克服这种平均化引起的图象模糊现象，我们给中心点象素值与其邻域平均值的差值设置一固定的阈值，只有大于该阈值的点才能替换为邻域平均值，而差值不大于阈值时，仍保留原来的值，从而减少由于平均化引起的图象模糊。边缘提取是较经典的算法，此处边缘的提取采用的是R
9、oberts算子。综上所述，结合 MATLAB实验过程，得出不是每一种图像处理之初都适合滤波和边界增强。本次汽车车牌的识别，为了保存更多的有用信息，经过多次比较，选择图 9作为后期处理的依据。7（6）操作界面设计为了使整个操作界面更加的美观，对背景颜色和字体颜色进行设置。二、程序调试及结果分析在程序设计过程中，碰到很多的问题。a) 一个函数只能对灰度图像处理，不能对 RGB图像处理，那么如何才能对 RGB图像处理呢？这主要是对 MATLAB函数的不够清楚，用到的很多函数是针对二维数据的，而 RGB图像的数据是一个三维矩阵，所以处理要与灰度图像不同，在开始的时候，我认为应该找一个能够应用于三维矩
10、阵的函数，结果却没找到，后来想到可以把三维数据进行降维处理，同样使用二维的函数，只要是同样处理三次。比如，彩色图像的滤波处理，直方图均衡等。b) 同一个操作对灰度图像可以使用，当用户选择的是彩色图像时，该操作就会出错？这是由于设计程序时，本身程序是有针对性的，有些程序只能对灰度图像有效，有些对彩色图像有效，但是用户并不清楚这些，所以在设计的过程中就要考虑全面，要分开设计。在本次课程设计中，我均对每个程序的开始时，用 if isrgb（x）进行判断。c) 在没有加入噪声的情况下，点击“中值滤波”或是其它滤波，会提示错误？这是由于在滤波程序设计的过程中，开始用到的变量是handles.noise_
11、img，而这个变量是在加入噪声时候才定义的，所以在没有加入噪声的情况下，点击各个滤波就会弹出变量没有定义的错误，解决方法就是在文件打开的时候就给定义handles.noise_img=x。d) 本次设计存在一个比较大的问题，就是每次操作都是独立的，比如：要对图像加入噪声，然后在此基础上进行亮度调整，截图，频谱分析等等就不行。我认为应该是整个程序设计过程中对变量的设置没有做好，应该每次操作后，把处理后的数据保存在一个全局变量，这样还要对处理后的数据进行在处理时，只要把这个全局变量作为原始数据带入就可。三， 心得体会1、在这次课程设计过程中，感触很深，由于对MATLAB图像处理的函数不熟悉，导致自
12、己走了很多的弯路，比如在设计图像左转90度和右转90度时，由于开始并不知道imrotate函数，只知道上下翻转flipud和左右翻转fliplr函数，想着要怎么用这两个函数来实现左转和右转呢，如果当当只用矩阵转置的话，并不能达到要求，后来想可以结合flipud和fliplr函数来实现，下面程序是实现左转90度：axes(handles.axes2);8x=(handles.img);if isrgb(handles.img) a=x(:,:,1);b=x(:,:,2);c=x(:,:,3);e=a;g=b;f=c;y(:,:,1)=e;y(:,:,2)=g;y(:,:,3)=f;axes(ha
13、ndles.axes2);for k=1:3 f(:,:,k)=flipud(y(:,:,k); endimshow(f); elsem=x;y=flipud(m);imshow(y);end同理，也可以通过转置和fliplr函数实现右转90度。但是后来发现其实只要用 imrotate函数就可以解决问题。通过这次经历后，后来在设计其它程序时，尽量找 MATLAB自带的现成函数，而不是一碰到问题就自己想算法，这样可以节省较多的时间。2、在理工科的专业应用背景下，用 matlab 进行相关计算与仿真编程的优势非常突出。特定的问题处理算法，我们通常都以 M文件的文本形式给定最终的解决方案，自己设计的
14、程序是在 MATLAB环境下，用 MATLAB语言编写的，这对于有安装 MATLAB软件的计算机上运行并不存在什么问题，关键是一般的计算机很少有去安装 MATLAB软件的，那么要在这些机子上运行该程序要怎么办呢，难道要先安装 MATLAB软件，这显然太麻烦了，对于一个通用的，比较成熟的解决方案，我们当然期望它能应用到更多的场合，而.exe （可执行）文件可运行于所有的通用 WINDOWS操作系统，为此，将 M文件转换成.exe 文件倒是个不错的想法。通过查找资料，了解到在生成可执行性文件之前，需要进行一系列编译环境方面的配置：要将用 Matlab语言编写的函数文件编译成可独立执行的*.exe
15、文件(即可脱离 Matalab环境的执行程序)，首先要安装和配置好 Matlab Compiler，一般来说，在安装 Matlab时就已经安装了相应版本的 Matlab Compiler。只是不同版本的 Matlab，其编译器的使用方法有一定的差异，这一点要引起一定的注意。在确定安装好 Matlab Compiler后，还需要对 Compiler进行适当的配置，方法是在 Matlab命令窗口输入：Mbuild setup9然后根据提示执行相应的操作，使用者可根据自己计算机中现有编译器的情况选择合适的编译器，如 VC+ 6.0、VC+7.0、Bland C 的编译器等，目前Matlab好象还不支
16、持 VC+8.0(我计算机安装的就是 VC+2005，Matlab 就无法识别)。当然，如果你的计算机里根本就没有安装其他任何语言的编译器，也可选择 Matlab自带的 Lcc编译器，其实这个编译器对大多数用户已经够用了(我就是选择的 Matlab自带的 Lcc编译器)。配置好编译器后，自然就是对自己编写的 M文件进行编译了。将 M文件编译为独立可执行文件的语法是：mcc m fun1.m fun2.m其中 fun1就是最后的可执行文件的名称。3、通过本次课程设计，使自己对 MATLAB GUI设计流程有了比较深刻的体会，同时也了解了一般软件设计的过程。在设计过程中碰到了很多的问题，通过这些问
17、题，使自己分析问题，解决问题的能力得到了较大的提高。四， 参考文献【1】 郑阿奇，曹戈，赵阳.MATLAB 实用教程M.北京：电子工业出版社【2】 程卫国，冯峰，姚东，徐听.MATLAB5.3 应用指南M.北京：人民邮电出版社【3】 陈杨.MATLAB 6.X图像编程与图像处理M.西安：西安电子科技大学出版社五， 附录function varargout = two(varargin)% Begin initialization code - DO NOT EDITgui_Singleton = 1;gui_State = struct(gui_Name, mfilename, .gui_Si
18、ngleton, gui_Singleton, .gui_OpeningFcn, two_OpeningFcn, .10gui_OutputFcn, two_OutputFcn, .gui_LayoutFcn, , .gui_Callback, );if nargin endif nargoutvarargout1:nargout = gui_mainfcn(gui_State, varargin:);elsegui_mainfcn(gui_State, varargin:);end% End initialization code - DO NOT EDIT% - Executes just
19、 before erzhi is made visible.function two_OpeningFcn(hObject, eventdata, handles, varargin)% This function has no output args, see OutputFcn.-%接下来是菜单的程序-function file_Callback(hObject, eventdata, handles)-%文件打开程序function open_Callback(hObject, eventdata, handles)name,path=uigetfile(*.bmp,载入图像);if i
20、sequal(name,0)|isequal(path,0)errordlg(没有选中文件,出错);return;elsex=imread(path,name);axes(handles.axes1);imshow(x);handles.img=x;handles.noise_img=x;guidata(hObject,handles)end-function save_Callback(hObject, eventdata, handles)%文件保存filename,pathname = uiputfile(*.bmp,图片保存为);if isequal(filename,pathname
21、,0,0)errordlg(没有保存,出错);return;elsefile=strcat(pathname,filename);(handles.axes2);i=getimage(gca);imwrite(i,file);end-function exit_Callback(hObject, eventdata, handles)clc;close all;close(gcf);11-%编辑菜单的程序-function edit_Callback(hObject, eventdata, handles)-function huidu_Callback(hObject, eventdata,
22、 handles)%灰度处理axes(handles.axes2);if isrgb(handles.img) y=rgb2gray(handles.img); %RGBimshow(y);elsemsgbox(这已经是灰度图像,转换失败);end -function liangdu_Callback(hObject, eventdata, handles)%亮度处理prompt=输入参数 1,输入参数 2,输入 gamma;defans=0 0.7,0 1,1;p=inputdlg(prompt,输入参数,1,defans);p1=str2num(p1);p2=str2num(p2);p3=
23、str2num(p3);gamma=p3;x=(handles.img);y=imadjust(x,p1,p2,gamma);axes(handles.axes2);imshow(y);-function jietu_Callback(hObject, eventdata, handles)%截图set(handles.axes2,HandleVisibility,ON);axes(handles.axes2);y=imcrop(handles.img);imshow(y); handles.Timage=y;-function fangda_Callback(hObject, eventda
24、ta, handles)-function lingjz_Callback(hObject, eventdata, handles)%邻近插值放大axes(handles.axes2);prompt=输入放大倍数:;defans=2;p=inputdlg(prompt,输入放大倍数,1,defans);p1=str2num(p1);y=imresize(handles.img,p1,nearest); %最近邻插值法放大imshow(y);-function shuangxianfa_Callback(hObject, eventdata, handles)%双线放大axes(handles.
25、axes2);prompt=输入放大倍数:;defans=2;p=inputdlg(prompt,输入放大倍数,1,defans);p1=str2num(p1);12y=imresize(handles.img,p1,bilinear); %最近邻插值法放大imshow(y);-function shuangsanci_Callback(hObject, eventdata, handles)%双三次放大axes(handles.axes2);prompt=输入放大倍数:;defans=2;p=inputdlg(prompt,输入放大倍数,1,defans);p1=str2num(p1);y=
26、imresize(handles.img,p1,bicubic); imshow(y);-function suoxiao_Callback(hObject, eventdata, handles)%邻近插值缩小axes(handles.axes2);prompt=输入放大倍数:;defans=0.2;p=inputdlg(prompt,输入放大倍数,1,defans);p1=str2num(p1);y=imresize(handles.img,p1,nearest); %最近邻插值法缩小imshow(y);-function shuangxian_Callback(hObject, even
27、tdata, handles)%双线性缩小axes(handles.axes2);prompt=输入放大倍数:;defans=0.2;p=inputdlg(prompt,输入放大倍数,1,defans);p1=str2num(p1);y=imresize(handles.img,p1,bilinear); %最近邻插值法缩小imshow(y); -function shuangsancisuox_Callback(hObject, eventdata, handles)%双三次缩小axes(handles.axes2);prompt=输入放大倍数:;defans=0.2;p=inputdlg(
28、prompt,输入放大倍数,1,defans);p1=str2num(p1);y=imresize(handles.img,p1,bicubic); imshow(y); -function xuanzhuan_Callback(hObject, eventdata, handles)-function updown_Callback(hObject, eventdata, handles)%上下翻转axes(handles.axes2);x=(handles.img);if isrgb(handles.img) for k=1:3 y(:,:,k)=flipud(x(:,:,k);%上下翻转
29、函数endimshow(y);elsex=(handles.img);y=flipud(x); 13imshow(y);end-function leftright_Callback(hObject, eventdata, handles)%左右翻转axes(handles.axes2);if isrgb(handles.img) x=(handles.img);for k=1:3 y(:,:,k)=fliplr(x(:,:,k);%左右翻转函数endimshow(y);elsex=(handles.img);y=fliplr(x);imshow(y);end-function zeft90_
30、Callback(hObject, eventdata, handles)%左转 90度axes(handles.axes2);x=(handles.img);y=imrotate(x,90); imshow(y);-function right90_Callback(hObject, eventdata, handles)%右转 90度axes(handles.axes2);x=(handles.img);y=imrotate(x,-90); imshow(y);-function other_Callback(hObject, eventdata, handles)%任意角度旋转axes(
31、handles.axes2);prompt=输入参数 1:;defans=30;p=inputdlg(prompt,输入参数,1,defans);p1=str2num(p1);y=imrotate(handles.img,p1); imshow(y);-function gs_Callback(hObject, eventdata, handles)%加入高斯噪声axes(handles.axes2);prompt=输入参数 1:,输入参数 2;defans=0,0.02;p=inputdlg(prompt,输入参数,1,defans);p1=str2num(p1);p2=str2num(p2
32、); y=imnoise(handles.img,gaussian,p1,p2);imshow(y);handles.noise_img=y;guidata(hObject,handles);-function jy_Callback(hObject, eventdata, handles)%加入椒盐噪声prompt=输入参数 1:; %对话框的设置，用户输入的是字符串defans=0.02; %缺省值p=inputdlg(prompt,输入参数,1,defans); p1=str2num(p1); %字符串转化为数值axes(handles.axes2);x=(handles.img);y=
33、imnoise(x,salt imshow(y);handles.noise_img=y;guidata(hObject,handles);-function cx_Callback(hObject, eventdata, handles)%加入乘性噪声axes(handles.axes2);prompt=输入参数 1:;defans=0.02;p=inputdlg(prompt,输入参数,1,defans);p1=str2num(p1);y=imnoise(handles.img,speckle,p1);imshow(y);handles.noise_img=y;guidata(hObjec
34、t,handles);-function zhifangtutongji_Callback(hObject, eventdata, handles)-function red_Callback(hObject, eventdata, handles)%R直方图set(handles.axes2,HandleVisibility,ON);axes(handles.axes2);x=imhist(handles.img(:,:,1); %直方图统计x1=x(1:10:256);horz=1:10:256;bar(horz,x1);set(handles.axes2,xtick,0:50:255);
35、-function gray_Callback(hObject, eventdata, handles)%G直方图set(handles.axes2,HandleVisibility,ON);axes(handles.axes2);if isrgb(handles.img) x=imhist(handles.img(:,:,2); %直方图统计x1=x(1:10:256);horz=1:10:256;bar(horz,x1);set(handles.axes2,xtick,0:50:255);elsemsgbox(这是灰度图像,旋转失败);end-function blue_Callback(
36、hObject, eventdata, handles)%B直方图set(handles.axes2,HandleVisibility,ON);axes(handles.axes2);if isrgb(handles.img) x=imhist(handles.img(:,:,3); %直方图统计x1=x(1:10:256);horz=1:10:256;bar(horz,x1);%axis(0 255 0 150000);set(handles.axes2,xtick,0:50:255);%set(handles.axes2,ytick,0:2000:15000);elsemsgbox(这是灰
37、度图像,旋转失败);end-15function junheng_Callback(hObject, eventdata, handles)%直方图均衡set(handles.axes2,HandleVisibility,ON);axes(handles.axes2);if isrgb(handles.img) a=histeq(handles.img(:,:,1);b=histeq(handles.img(:,:,2);c=histeq(handles.img(:,:,3);k(:,:,1)=a;k(:,:,2)=b;k(:,:,3)=c;imshow(k);elseh=histeq(han
38、dles.img); %直方图均衡imshow(h);end-%频谱分析-function pinpu_Callback(hObject, eventdata, handles)-function pinputu_Callback(hObject, eventdata, handles)%显示频谱图axes(handles.axes2);x=(handles.img);if isrgb(x)m=fft2(x(:,:,1);y=fftshift(m);imshow(log(abs(y),);elsem=fft2(x);y=fftshift(m);imshow(log(abs(y),);end-f
39、unction frequency_Callback(hObject, eventdata, handles)%低通滤波器axes(handles.axes2);x=(handles.img);if isrgb(x)msgbox(这是彩色图像,不能通过低通滤波器,失败);elsey1=imnoise(x,salt % 叠加椒盐噪声f=double(y1); % 数据类型转换，MATLAB 不支持图像的无符号整型的计算g=fft2(f); % 傅立叶变换g=fftshift(g); % 转换数据矩阵M,N=size(g);nn=2; % 二阶巴特沃斯(Butterworth)低通滤波器d0=10; %截止频率为 10m=fix(M/2); n=fix(N/2);for i=1:Mfor j=1:Nd=sqrt(i-m)2+(j-n)2);h=1/(1+0.414*(d/d0)(2*nn); % 计算低通滤波器传递函数result(i,j)=h*g(i,j);endendresult=ifftshift(result);y2=ifft2(result);y3=uint8(real(y2);16imshow(y3); % 显示滤波处理后的图像end-function gaoto
