河 北 大 学 数 学 与 计 算 机 学 院
课程名称：数据仓库与数据挖掘 成绩：
综合实验题目
基于正交最小二乘拟合算法的图像矢量化
课 程 号： 课 序 号：02 任课教师：邢红杰 班 级：2011软工 姓 名：连冰川 学 号：
填写日期：5月18日
河 北 大 学 数 学 与 计 算 机 学 院
基于正交最小二乘拟合算法的图像矢量化
2011软件工程 连冰川
1. 实验问题
形图形（或图像）在计算机里主要有两种存储和表示方法。矢量图是使用点、直线或多边形等基于数学方程的几何对象来描述图形，位图则使用像素来描述图像。一般来说，照片等相对杂乱的图像使用位图格式较为合适，矢量图则多用于工程制图、标志、字体等场合。矢量图可以任意放缩，图形不会有任何改变。而位图一旦放大后会产生较为明显的模糊，线条也会出现锯齿边缘等现象。 矢量图从本质上只是使用曲线方程对图形进行的精确描述，在以像素为基本显示单元的显示器或打印机上是无法直接表现的。将矢量图转换成以像素点阵来表示的信息，再加以显示或打印，这个过程称之为栅格化（Rasterization），栅格化的逆过程相对比较困难。假设有一个形状较为简单的图标，保存成一定分辨率的位图文件。我们希望将其矢量化，请你建立合理的数学模型，尽量准确地提取出图案的边界线条，并将其用方程表示出来。
问题分析与建模（或所使用方法的算法描述及算法流程）
2.1问题分析
题目要求我们针对形状较为简单的位图，通过建立合理的数学模型，将其矢量化，并准确地提取出图案边界线条，找出线条的曲线方程。
初步分析可知：由于位图放大后较为模糊，其线条将会出现明显地锯齿边缘，为了转换为精确、可以任意放缩的矢量图，我们需要通过对二值化图像进行矢量化和有效地拟合。
基于以上分析，为了得到有效地二值化图像，我们需要采取多种方法对图像进行分割，在多种分割的图像结果的比较中，选出作为完美的分割结果。
然而，图像有效地分割需要我们首先对扫描后的灰度图像进行处理，即进行多尺度边缘检测。在采取多种边缘检测算子的基础上，我们通过对灰度图像的膨胀与腐蚀的处理，得到较为完美的多尺度边缘检测结果，作为图像分割的基础。从而实现从位图到矢量图的完美转化。 2.2问题假设
（1）假设原位图是二维的、简单的图像； （2）假设原图像是单色图；
（3）假设所采用的边缘检测方法抽取到最佳线类型特征；
（4）假设细化后的图像不会造成原先连续图形的断裂，从而不破坏图形的连通性；
（5）假设细化后保持了原图形的基本特征，不会破坏图形间的拓扑关系； （6）假设运用的细化算法具有稳定性和收敛性； （7）假设经过预处理的二值图是确切的矢量图；
（8）假设基于模糊识别的非线性系统的参数收敛性最佳； （9）假设基于正交最小二乘的非线性系统模糊识别得到的矢量化方程误差最
第1页
小。
2.3模型建立与求解
（一）基于图像腐蚀、膨胀的多尺度边缘检测
我们在扫描位图得到的灰度图像的基础上，通过采取图像腐蚀、膨胀的方法
对灰度图像进行多尺度边缘检测的初步处理。 1.边缘检测算子的选取
针对本文选取简单的位图进行图像分割、目标区域识别、区域形状提取等处理，从而实现边缘检测算子的提取，将选取的算子作为图像腐蚀、膨胀的基础。
在进行边缘检测时，求出梯度f(x,y)梯度的局部最大值和方
?f?f?x?f?y?????fxcos??fysin?向:
?r?x?r?y?r
??f?????f??r??0达到最大值，扫描得到的灰度图
由导数的算术意义可知当??时，?r像的方程设为fxsing??fycos?g?0，根据导数中最大值的求法可知
?fx??g?arcta?n???f??y?或者表示为
??f?g?????r??maxfx2?fy2?g??。由此，得到梯度模
对,梯度模算子具有位移不变性和各方向的良好性质，
?f?于灰度变化的方向，即边界的方向，我们可以通过?g?arctan?x?得到。
?f??y?（1）边缘算子图像处理
为了加强图像处理的效果，我们通过六种边缘算子处理的比较，选出最优的表现结果。
1）Roberts边缘算子
G?f?x,y??? (1)
??f(x,y)?f(x?1,y?1)???2f(x?1,y)?f(x,y?1)??122
在本式中我们采取平方根运算使该处理更为精准，运用此算法实际上是求旋转45两个方向上微分值的和。
通过Roberts边缘算子得到的结果如下图所示：