西安交大模式识别实验报告

西安交大模式识别实验报告

模式识别实验报告

实 验 报 告

实验课程名称： 模式识别

姓名： 王宇 班级：号： 2011081325

注：1、每个实验中各项成绩按照5分制评定，实验成绩为各项总和

2、平均成绩取各项实验平均成绩

3、折合成绩按照教学大纲要求的百分比进行折合

2014年 6月

实验一、 图像的贝叶斯分类

一、实验目的

将模式识别方法与图像处理技术相结合，掌握利用最小错分概率贝叶斯分类器进行图像分类的基本方法，通过实验加深对基本概念的理解。 二、实验仪器设备及软件 HP D538、MATLAB 三、实验原理 概念：

阈值化分割算法是计算机视觉中的常用算法，对灰度图象的阈值分割就是先确定一个处于图像灰度取值范围内的灰度阈值，然后将图像中每个像素的灰度值与这个阈值相比较。并根据比较的结果将对应的像素划分为两类，灰度值大于阈值的像素划分为一类，小于阈值的划分为另一类，等于阈值的可任意划分到两类中的任何一类。

最常用的模型可描述如下：假设图像由具有单峰灰度分布的目标和背景组成，处于目标和背景内部相邻像素间的灰度值是高度相关的，但处于目标和背景交界处两边的像素灰度值有较大差别，此时，图像的灰度直方图基本上可看作是由分别对应于目标和背景的两个单峰直方图混合构成。而且这两个分布应大小接近，且均值足够远，方差足够小，这种情况下直方图呈现较明显的双峰。类似地，如果图像中包含多个单峰灰度目标，则直方图可能呈现较明显的多峰。

上述图像模型只是理想情况，有时图像中目标和背景的灰度值有部分交错。这时如用全局阈值进行分割必然会产生一定的误差。分割误差包括将目标分为背景和将背景分为目标两大类。实际应用中应尽量减小错误分割的概率，常用的一种方法为选取最优阈值。这里所谓的最优阈值，就是指能使误分割概率最小的分割阈值。图像的直方图可以看成是对灰度值概率分布密度函数的一种近似。如一幅图像中只包含目标和背景两类灰度区域，那么直方图所代表的灰度值概率密度函数可以表示为目标和背景两类灰度值概率密度函数的加权和。如果概率密度函数形式已知，就有可能计算出使目标和背景两类误分割概率最小的最优阈值。 假设目标与背景两类像素值均服从正态分布且混有加性高斯噪声，上述分类问题可以使用模式识别中的最小错分概率贝叶斯分类器来解决。以p1与p2分别表示目标与背景的灰度分布概率密度函数，P1与P2分别表示两类的先验概率，则图像的混合概率密度函数可用下式表示为

p(x)?P1p1(x)?P2p2(x)

式中p1和p2分别为

p1(x)?

p2(x)?

?

(x??1)22?1

?

(x??2)22?22

P1?P2?1

?1、?2是针对背景和目标两类区域灰度均值?1与?2的标准差。若假定目标的灰度较亮，

其灰度均值为?2，背景的灰度较暗，其灰度均值为?1，因此有

?1??2

现若规定一门限值T对图像进行分割，势必会产生将目标划分为背景和将背景划分为目标这两类错误。通过适当选择阈值T，可令这两类错误概率为最小，则该阈值T即为最佳阈值。 把目标错分为背景的概率可表示为

E1(T)??p2(x)dx

??

T

把背景错分为目标的概率可表示为

E2(T)??

总的误差概率为

??

T

p1(x)dx

E(T)?P2E1(T)?PE12(T)

为求得使误差概率最小的阈值T，可将E(T)对T求导并令导数为零，可得

P1p1(T)?P2p2(T)

代换后，可得

P(T??1)2(T??2)21?2

ln???

P2?12?122?12

此时，若设?1??2??，则有

T?

?1??2

2

?P??2

?ln?2??1??2?P1?

若还有P1?P2的条件，则

T?

?1??2

2

这时的最优阈值就是两类区域灰度均值?1与?2的平均值。

上面的推导是针对图像灰度值服从正态分布时的情况，如果灰度值服从其它分布，依理也可求出最优阈值来。一般情况下，在不清楚灰度值分布时，通常可假定灰度值服从正态分布。 在实际使用最优阈值进行分割的过程中，需要利用迭代算法来求得最优阈值。设有一幅数字图像f(x,y)，混有加性高斯噪声，可表示为

g(x,y)?f(x,y)?n(x,y)

此处假设图像上各点的噪声相互独立，且具有零均值，如果通过阈值分割将图像分为目标与背景两部分，则每一部分仍然有噪声点随机作用于其上，于是，目标g1(x,y)和g2(x,y)可表示为

g1(x,y)?f1(x,y)?n(x,y) g2(x,y)?f2(x,y)?n(x,y)

迭代过程中，会多次地对g1(x,y)和g2(x,y)求均值，则

E{g1(x,y)}?E{f1(x,y)?n(x,y)}?E{f1(x,y)} E{g2(x,y)}?E{f2(x,y)?n(x,y)}?E{f2(x,y)}

可见，随着迭代次数的增加，