《图像校正及分割处理软件设计与实现模式识别与智能系统--本科毕设论文》由会员分享，可在线阅读，更多相关《图像校正及分割处理软件设计与实现模式识别与智能系统--本科毕设论文（24页珍藏版）》请在金锄头文库上搜索。
1、图像校正及分割处理软件设计与实现图像分析与模式识别课程期末大作业报告课程名称： 图像分析与模式识别 目录图像校正及分割处理软件设计与实现11 软件需求11.1 操作界面需求分析11.2 图像几何校正需求分析11.2 ISODATA聚类算法的图像分割需求分析22 算法原理及流程图22.1 Hough变换22.1.1 Hough变换原理22.1.2 图像几何校正流程32.2 ISODATA算法原理32.2.1 ISODATA算法原理步骤32.2.1 ISODATA算法流程图53 程序设计框图64 实验结果及分析74.1 图像几何74.1.1 支票图像几何校正结果74.1.1 支票图像几何校正结果分
2、析84.2 图像分割84.2.1图像分割结果84.2.1图像分割结果分析10附录：11附录一：图像几何校正代码11附录一：ISODATA聚类算法图像分割代码11附录一：MATLAB的GUI操作界面代码17 图像分析与模式识别课程期末大作业报告姓名：赵世瑜 学号：2013204067图像校正及分割处理软件设计与实现摘要：设计一种图像几何校正及图像分割处理软件，实现对倾斜支票图像校正为水平，及对已给的图像进行分割。图像几何校正主要通过边缘检测、Hough变换、求倾斜角、图像旋转等算法实现；图像分割采用ISODATA聚类算法实现。实验结果表明，实现了对倾斜支票图像几何校正和对图像的分割功能。关键词：
3、边沿检测；Hough变换、图像旋转、ISODATA算法；图像分割；1 软件需求需求分析是指对要解决的问题进行详细的分析，弄清楚问题的要求，包括整个系统所要实现的功能。根据题目，设计需要编写带有操作界面的图像处理软件，功能包含图像几何校和基于ISODATA聚类算法的图像分割。开发工具使用MATLAB。根据MATLAB的特点及需要实现的软件功能，软件需求分析如下：1.1 操作界面需求分析软件实现的功能是图像处理，因此操作界面中需要一个按钮，用于选择待处理的图像，称为“图像选择”按钮；选择的原始图像需要在操作界面中显示，因此需要一个用于显示原始图像的坐标控件；由于软件需要实现多功能，需要一个下拉框实
4、现功能的选择；一个开始功能处理的按钮，称为“开始处理”按钮；一个用于显示图像处理结果的坐标控件。最后需要一个退出操作界面的按钮，称为“退出”按钮。1.2 图像几何校正需求分析原始图像在拍照时出现了倾斜，需要将图像校正，消除倾斜。通过对图像的观察，图像是一张支票，支票上存在表格。因此想法是通过检测支票表格边框线的倾斜角，就是整张图像的倾斜程度，然后根据这个倾斜角旋转，便可得到校正的水平图像。要检测支票图像表格边框直线。需要对图像的边缘进行检测，图像的边缘检测方法很多，如：sobel算子边缘检测，prewitt算子边缘检测，roberts算子边缘检测，log算子边缘检测，canny算子边缘检测等，
5、由于原图的噪声并不是很严重，因此上面这些边缘检测办法都能胜任，这里使用log算子进行边缘检测，为了适应其它边缘检测方法，本设计需要尽量在改变边缘检测方法时尽可能少地修改代码。边缘检测后，采用Hough变换提取直线，然后在提取的直线上找两个点，计算出直线的倾斜角度。最后根据这个倾斜角度对图像进行旋转。1.2 ISODATA聚类算法的图像分割需求分析ISODATA聚类算法的图像分割功能，ISODATA比较复杂，因此主要是ISODATA算法的实现。设计的处理过程是，在操作界面选择了待处理的图像后，选择图像分割功能，点击开始处理按钮后，弹出图像，通过鼠标左键选择聚类中心点，然后按回车键或者单击鼠标右键
6、确认后，调用ISODATA对图像进行分割。2 算法原理及流程图设计的功能有图像几何和基于ISODATA聚类算法的图像分割。图像几何校正功能中，采用Hough变换对支票表格边框直线进行检测；ISODATA聚类算法的图像分割功能和核心是ISODATA聚类算法。对其原理及实现的流程图分析如下。2.1 Hough变换哈夫变换是利用图像全局特性而直接检测目标轮廓，就可将边缘像素连接起来组成区域封闭边界的一种方法。在预先知道区域形状的条件下，利用哈夫变换可以方便地得到边界曲线而将不连续的边缘像素点连接起来。2.1.1 Hough变换原理哈夫变换的基本思想是点-线对偶性。图像变换前在图像空间，变换后在参数空
7、间。在图像空间XY里，所有过点(x, y)的直线都满足方程: y=px+q写成： q=-px+y （1）可以认为代表参数空间PQ中过点(p,q)的一条直线。图1 图像空间 图2 参数空间图像空间的点(xi, yi)和(xj, yj)对应参数空间两条直线，它们的交点(p, q)对应图像空间XY中过(xi, yi)和(xj, yj)的直线，这样，图像空间XY中过(xi, yi)和(xj,yj)的直线上的所有点都对应参数空间PQ里的一条直线，这些直线相交于点(p, q) 。 在图像空间中共线的点对应在参数空间中相交的线，反过来，在参数空间中交于同一个点的所有直线在图像空间中都有共线的点相对应，这就是
8、点-线对偶性。根据点-线对偶性，当给定图像空间的一些边缘点，就可以通过哈夫变换确定连接这些点的直线方程，哈夫变换把图像空间中的直线检测问题转换到参数空间里点检测的问题。2.1.2 图像几何校正流程图像几何校正中使用的Hough变换算法采用MATLAB已经封装好的hough函数，这里不要详细阐述其流程，整个图像几何校正的流程如图3所示：图3 图像几何校正流程图2.2 ISODATA算法原理ISODATA算法的思想是，在每轮迭代过程中，样本重新调整类别之后计算类内及类间有关参数，并和设定的门限比较，确定是两类合并为一类还是一类分裂为两类，不断地“自组织”，以达到在各参数满足设计要求条件下，使各模式
9、到其类心的距离平方和最小。2.2.1 ISODATA算法原理步骤 预置 设定聚类分析控制参数：=预期的类数，=初始聚类中心个数(可以不等于c)，137=每一类中允许的最少模式数目，=类内各分量分布的距离标准差上界，(分裂用)=两类中心间的最小距离下界，(合并用)=在每次迭代中可以合并的类的最多对数，=允许的最多迭代次数。 将待分类的模式特征矢量读入； 选定初试聚类中心，可从待分类的模式特征矢量集中任选Nc个模式特征矢量作为初始聚类中心；（2）按最小距离原则将模式集中每一个模式分到某一类中；（3）依据每一类中允许的最少模式数目判断合并；（4）计算分类后的参数：分类中心，；类平均距离及中体平均距离
10、； 计算各类的中心 （2） 计算各类中模式到类心的平均距离 （3） 计算各个模式到其类内中心的总体平均距离 （4）（5）依据Ip, Nc判断停止、分类或合并；（6）计算各类类内距离的标准差矢量； （5）（7）对每一次聚类，求出类内距离标准差矢量中的最大分量；（8）在（7）所得结果中对任一最大值大于类内各分量分布的距离标准差上界则分裂；（9）计算各对聚类中心间的距离；（10）依据两类中心间的最小距离下界判断合并；（11）如果迭代次数Ip已达I次或过程收敛，则结束，否则，Ip=Ip+1，若需要调整参数，则转（1），若不改变参数则转（2）；2.2.1 ISODATA算法流程图图4 ISODATA聚类
11、算法流程图3 程序设计框图根据需求分析，对整个操作界面及其中的函数调用等，对整个程序进行分析。操作界面运行后，单击选择图像按钮选择待处理的图像。然后选择要对图形进行处理的功能，单击开始处理按钮，调用相应的处理方法进行处理，整个程序的设计框图如图5所示：图5 程序设计框图设计好的操作界面如下图所示：图6 操作界面设计结果4 实验结果及分析4.1 图像几何4.1.1 支票图像几何校正结果图像几何校正功能中，对支票图像的校正如下图所示：图7 图像结合校正处理中图8 图像结合校正处理中4.1.1 支票图像几何校正结果分析通过支票图像几何校正结果可知，设计达到了预计的结果。通过MATLAB运行提取计算斜
12、率的坐标为 (1200,522)和(1540, 511)。直线的斜率为slope = -0.0324，通过次斜率计算出的倾斜角度为angle = -1.8530。根据人眼观测结果显示支票的请倾斜角也比较小，而且根据校正后的支票图像结果显示已经水平。即达到了实验设计目的。4.2 图像分割4.2.1图像分割结果在运行图像分割功能时，需要通过鼠标选择预期聚类，其中选择之一如图9所示：图9 选择预期的聚类中心图10 图像分割处理中最终对图像的分割结果如下图所示：图11 图像分割结果（a）图11 图像分割结果（b）图13 图像分割结果（c）图14 图像分割结果（d）4.2.1图像分割结果分析通过多次选择
13、不同的预期聚类中心的实验结果表明，设计基本完成了所需的功能。同时也证明了设计的合理性，以及ISODATA聚类算法在图像分割中的运用。结果显示，把鸭子的眼睛分成了一类，把鸭子的嘴分成了一类，把鸭子的身体分成了一类，虽然有时分类结果不是非常好。但选择预期聚类中心时有一定的影响。总体来看，仍然达到了图像分割的目的。本学期在图像分析与模式识别课程中学到很多，学到的不仅仅是课程相关的知识。更重要的是学到了许多处理和解决问题的思想，相信对今后的学习会有很大的帮助。这都是刘辉老师教导的结果，在此对刘辉老师表示衷心感谢！附录：附录一：图像几何校正代码function jiaozheng = jiaozheng
14、(cheque)BW_laplace = edge(cheque,log);% log算子边缘检测BW2 = bwmorph(BW_laplace, thin,Inf);%细化H,T,R = hough(BW2); %Hough变换P = houghpeaks(H,5); %提取Hough变换后参数平面上的峰值点 lines = houghlines(BW2,T,R,P); %提取直线for k = 1:length(lines) xy = lines(k).point1; lines(k).point2; %提取坐标endslope=(xy(2,2)- xy(1,2)/(xy(2,1)-xy(1,1); %求斜率radian=atan(slope);%求角度angle=radian*180/pi; %转化成角度jiaozheng=imrotate(cheque,angle); %图像校正imshow(jiaozheng);close ;附录一：ISODATA聚类算法图像分割代码function isodata=isodata_zsy(I)
