《哈尔滨工程大学 模式识别实验报告》由会员分享，可在线阅读，更多相关《哈尔滨工程大学 模式识别实验报告（24页珍藏版）》请在金锄头文库上搜索。
1、实 验 报 告实验课程名称：模式识别姓名：班级： 20120811 学号：实验名称规范程度原理叙述实验过程实验结果实验成绩图像的贝叶斯分类K 均值聚类算法神经网络模式识别平均成绩折合成绩注：1、每个实验中各项成绩按照 5 分制评定，实验成绩为各项总和2、平均成绩取各项实验平均成绩3、折合成绩按照教学大纲要求的百分比进行折合2015 年 4 月实验 1 图像的贝叶斯分类1.1 实验目的将模式识别方法与图像处理技术相结合，掌握利用最小错分概率贝叶斯分类器进行图像分类的基本方法，通过实验加深对基本概念的理解。1.2 实验仪器设备及软件HP D538、MATLAB1.3 实验原理1.3.1 基本原理阈
2、值化分割算法是计算机视觉中的常用算法， 对灰度图象的阈值分割就是先确定一个处于图像灰度取值范围内的灰度阈值， 然后将图像中每个像素的灰度值与这个阈值相比较。 并根据比较的结果将对应的像素划分为两类， 灰度值大于阈值的像素划分为一类， 小于阈值的划分为另一类，等于阈值的可任意划分到两类中的任何一类。此过程中，确定阈值是分割的关键。对一般的图像进行分割处理通常对图像的灰度分布有一定的假设， 或者说是基于一定的图像模型。最常用的模型可描述如下：假设图像由具有单峰灰度分布的目标和背景组成，处于目标和背景内部相邻像素间的灰度值是高度相关的， 但处于目标和背景交界处两边的像素灰度值有较大差别， 此时， 图
3、像的灰度直方图基本上可看作是由分别对应于目标和背景的两个单峰直方图混合构成。而且这两个分布应大小接近，且均值足够远，方差足够小，这种情况下直方图呈现较明显的双峰。类似地，如果图像中包含多个单峰灰度目标，则直方图可能呈现较明显的多峰。上述图像模型只是理想情况， 有时图像中目标和背景的灰度值有部分交错。 这时如用全局阈值进行分割必然会产生一定的误差。 分割误差包括将目标分为背景和将背景分为目标两大类。实际应用中应尽量减小错误分割的概率，常用的一种方法为选取最优阈值。这里所谓的最优阈值， 就是指能使误分割概率最小的分割阈值。 图像的直方图可以看成是对灰度值概率分布密度函数的一种近似。 如一幅图像中只
4、包含目标和背景两类灰度区域， 那么直方图所代表的灰度值概率密度函数可以表示为目标和背景两类灰度值概率密度函数的加权和。 如果概率密度函数形式已知，就有可能计算出使目标和背景两类误分割概率最小的最优阈值。假设目标与背景两类像素值均服从正态分布且混有加性高斯噪声， 上述分类问题可用模式识别中的最小错分概率贝叶斯分类器来解决。以1p与2p分别表示目标与背景的灰度分布概率密度函数，1P与2P分别表示两类的先验概率，则图像的混合概率密度函数用下式表示1122( )( )( )p xPp xP px式中1p和2p分别为2 1 2 1()2 1 11( )2x p xe 2 2 2 2()2 2 21( )
5、2x pxe 121PP1、2是针对背景和目标两类区域灰度均值1与2的标准差。若假定目标的灰度较亮，其灰度均值为2，背景的灰度较暗，其灰度均值为1，因此有12现若规定一门限值T对图像进行分割， 势必会产生将目标划分为背景和将背景划分为目标这两类错误。通过适当选择阈值T，可令这两类错误概率为最小，则该阈值T即为最佳阈值。把目标错分为背景的概率可表示为12( )( )TE Tpx dx 把背景错分为目标的概率可表示为21( )( ) TE Tp x dx总的误差概率为2112( )( )( )E TP E TPE T为求得使误差概率最小的阈值T，可将( )E T对T求导并令导数为零，可得1122(
6、 )( )Pp TP p T代换后，可得22 1212 22 2111()()ln22PTT P 此时，若设12，则有2 122121ln2PTP 若还有12PP的条件，则12 2T这时的最优阈值就是两类区域灰度均值1与2的平均值。上面的推导是针对图像灰度值服从正态分布时的情况，如果灰度值服从其它分布，依理也可求出最优阈值来。一般情况下，在不清楚灰度值分布时，通常可假定灰度值服从正态分布。因此，本课题中亦可使用此方法来求得最优阈值，来对实验图像进行分割。1.3.2 最优阈值的迭代算法在实际使用最优阈值进行分割的过程中， 需要利用迭代算法来求得最优阈值。 设有一幅数字图像( , )f x y，混
7、有加性高斯噪声，可表示为( , )( , )( , )g x yf x yn x y此处假设图像上各点的噪声相互独立，且具有零均值，如果通过阈值分割将图像分为目标与背景两部分，则每一部分仍然有噪声点随机作用于其上，于是，目标1( , )g x y和2( , )gx y可表示为11( , )( , )( , )g x yf x yn x y22( , )( , )( , )gx yfx yn x y迭代过程中，会多次地对1( , )g x y和2( , )gx y求均值，则111( , ) ( , )( , ) ( , )E g x yE f x yn x yE f x y222( , )( ,
8、 )( , )( , )E gx yE fx yn x yE fx y可见，随着迭代次数的增加，目标和背景的平均灰度都趋向于真实值。因此，用迭代算法求得的最佳阈值不受噪声干扰的影响。利用最优阈值对实验图像进行分割的迭代步骤为：（1）确定一个初始阈值0T，0T可取为minmax 02SST式中，minS和maxS为图像灰度的最小值和最大值。（2）利用第 k 次迭代得到的阈值将图像分为目标1R和背景2R两大区域，其中1 ( , )|( , )kRf x yf x yT2 ( , )|0( , )kRf x yf x yT（3）计算区域1R和2R的灰度均值1S和2S。（4）计算新的阈值1kT，其中1
9、2 12kSST（5）如果1|kkTT小于允许的误差，则结束，否则1kk，转步骤（2） 。利用迭代法求得最优阈值后，仍需进行一些人工调整才能将此阈值用于实验图像的分割，这是因为，这种最优阈值仍然属于全局阈值，它利用了图像中所有像素点的信息，但当光照不均匀时， 图像中部分区域的灰度值可能差距较大， 造成计算出的最优阈值分割效果不理想，此时，可设一人工经验因子进行校正。1.4、实验步骤及程序实验步骤： 1、读取指定图像，取矩阵的最大值和最小值，并以最大值、最小值的平均值为初始阈值 A。 2、比较所有的矩阵因子和初始阈值的大小，若某矩阵因子较大，则有效区域的像素点数增 加 1，该点灰度值需计入有效区
10、域的灰度总值。反之，背景的像素点增加 1，该点灰度值需 计入背景的灰度值。 3、所有的矩阵因子都比较完以后，计算有效区域的像素平均灰度值和背景的平均灰度值。 取这两个平均值的平均，记为 B，若 A=B,则循环结束，该值为最优阈值。否则，令 A=B， 重复步骤 2、3。读取指定图像，取矩阵的最 大值和最小值， 并以最大值、 最小值的平均值作为初始的 阈值 A矩阵因子A?有效区域的像素点数增 加 1，该点灰度值需计入 有效区域的灰度总值背景的像素点增加 1，该 点灰度值需计入背景的 灰度值计算有效区域的像素平均灰度值和背景的平 均灰度值。取这两个平均值的平均，记为 B比较完毕？A-B=TK)iFo
11、reground=iForeground+1;ForegroundSum=ForegroundSum+double(tmp);elseiBackground=iBackground+1;BackgroundSum=BackgroundSum+double(tmp);endendendZO=ForegroundSum/iForeground;ZB=BackgroundSum/iBackground;TKTmp=double(ZO+ZB)/2);if(TKTmp=TK)bCal=0;elseTK=TKTmp;endenddisp(strcat(diedaihoudeyuzhi,num2str(TK
12、);newI=im2bw(I,double(TK)/255);imshow(I)figureimshow(newI)1.5、实验结果与分析实验得到的迭代后的分割阈值：94.8064分割效果图如下所示。图 1.5.1原始图像图 1.5.2分割后的图像020040060080010001200050100150200250图 1.5.3原始图像的灰度直方图实验分析：对灰度图象的阈值分割就是先确定一个处于图像灰度取值范围内的灰度阈值， 然后将图像中每个像素的灰度值与这个阈值相比较。 并根据比较的结果将对应的像素划分为两类， 灰度值大于阈值的像素划分为一类， 小于阈值的划分为另一类， 等于阈值的可任意
13、划分到两类中的任何一类。其中确定阈值是分割的关键。最优阈值的求得需要使用迭代算法。它将会影响到迭代的次数和结果精度。实验 2 K 均值聚类算法2.1 实验目的将模式识别方法与图像处理技术相结合，掌握利用 K 均值聚类算法进行图像分类的基本方法，通过实验加深对基本概念的理解。2.2 实验仪器设备及软件HP D538、MATLAB、WIT2.3 实验原理K 均值聚类法分为如下几个步骤：一、 初始化聚类中心1、 根据具体问题，凭经验从样本集中选出 C 个比较合适的样本作为初始聚类中心。2、 用前 C 个样本作为初始聚类中心。3、 将全部样本随机地分成 C 类，计算每类的样本均值，将样本均值作为初始聚
14、类中心。二、 初始聚类1、 按就近原则将样本归入各聚类中心所代表的类中。2、 取一样本，将其归入与其最近的聚类中心的那一类中，重新计算样本均值，更新聚类中心。然后取下一样本，重复操作，直至所有样本归入相应类中。三、 判断聚类是否合理采用误差平方和准则函数判断聚类是否合理， 不合理则修改分类。 循环进行判断、修改直至达到算法终止条件。2.4 实验步骤及程序实验步骤及程序实验步骤： 1、读取原始图像，确定四个初始聚类中心。2、计算各点与聚类中心的距离，以及各点到不同聚类中心的距离之差，选取距离最近的聚 类中心作为该点的聚类中心，依据此原理将属于不同聚类中心的元素聚类。 3、求各类的平均值作为新的聚
15、类中心，检验是否满足精度条件。 4、输出的四个聚类中心值，将图像分成四类输出。开始读取原始图像， 确定四 个初始聚类中心计算各点与聚类中心的距离，以及各点到 不同聚类中心的距离之差，选取距离最近 的聚类中心作为该点的聚类中心，据此原 理将属于不同聚类中心的元素聚类将各类的平均值作 为新的聚类中心满足误差条件？新的聚类中心即 为最终结果NY结束图 2.4.1实验程序流程图实验程序： clccleartic% A=imread(peppers.jpg);A=imread(N: 模式识别实验资料实验图片peppers.bmp); figure,imshow(A)figure,imhist(A)A=d
16、ouble(A);for i=1:200c1(1)=25;c2(1)=75;c3(1)=120;c4(1)=200;r=abs(A-c1(i);g=abs(A-c2(i);b=abs(A-c3(i);y=abs(A-c4(i);r_g=r-g;g_b=g-b;r_b=r-b;b_y=b-y;r_y=r-y;g_y=g-y;n_r=find(r_g0i=i+1;c1(i)=sum(A(n_r)/length(n_r);c2(i)=sum(A(n_g)/length(n_g);c3(i)=sum(A(n_b)/length(n_b);c4(i)=sum(A(n_y)/length(n_y);d1(i)=sqrt(abs(c1(i)-c1(i-1);d2(i)=sqrt(abs(c2(i)-c2(i-
