模式识别实验报告.doc
- 文档编号:163290
- 上传时间:2022-10-04
- 格式:DOC
- 页数:21
- 大小:804.50KB
模式识别实验报告.doc
《模式识别实验报告.doc》由会员分享,可在线阅读,更多相关《模式识别实验报告.doc(21页珍藏版)》请在冰豆网上搜索。
实验1图像的贝叶斯分类
1.1实验目的
将模式识别方法与图像处理技术相结合,掌握利用最小错分概率贝叶斯分类器进行图像分类的基本方法,通过实验加深对基本概念的理解。
1.2实验仪器设备及软件
HPD538、MATLAB
1.3实验原理
1.3.1基本原理
阈值化分割算法是计算机视觉中的常用算法,对灰度图象的阈值分割就是先确定一个处于图像灰度取值范围内的灰度阈值,然后将图像中每个像素的灰度值与这个阈值相比较。
并根据比较的结果将对应的像素划分为两类,灰度值大于阈值的像素划分为一类,小于阈值的划分为另一类,等于阈值的可任意划分到两类中的任何一类。
此过程中,确定阈值是分割的关键。
对一般的图像进行分割处理通常对图像的灰度分布有一定的假设,或者说是基于一定的图像模型。
最常用的模型可描述如下:
假设图像由具有单峰灰度分布的目标和背景组成,处于目标和背景内部相邻像素间的灰度值是高度相关的,但处于目标和背景交界处两边的像素灰度值有较大差别,此时,图像的灰度直方图基本上可看作是由分别对应于目标和背景的两个单峰直方图混合构成。
而且这两个分布应大小接近,且均值足够远,方差足够小,这种情况下直方图呈现较明显的双峰。
类似地,如果图像中包含多个单峰灰度目标,则直方图可能呈现较明显的多峰。
上述图像模型只是理想情况,有时图像中目标和背景的灰度值有部分交错。
这时如用全局阈值进行分割必然会产生一定的误差。
分割误差包括将目标分为背景和将背景分为目标两大类。
实际应用中应尽量减小错误分割的概率,常用的一种方法为选取最优阈值。
这里所谓的最优阈值,就是指能使误分割概率最小的分割阈值。
图像的直方图可以看成是对灰度值概率分布密度函数的一种近似。
如一幅图像中只包含目标和背景两类灰度区域,那么直方图所代表的灰度值概率密度函数可以表示为目标和背景两类灰度值概率密度函数的加权和。
如果概率密度函数形式已知,就有可能计算出使目标和背景两类误分割概率最小的最优阈值。
假设目标与背景两类像素值均服从正态分布且混有加性高斯噪声,上述分类问题可以使用模式识别中的最小错分概率贝叶斯分类器来解决。
以与分别表示目标与背景的灰度分布概率密度函数,与分别表示两类的先验概率,则图像的混合概率密度函数可用下式表示
式中和分别为
、是针对背景和目标两类区域灰度均值与的标准差。
若假定目标的灰度较亮,其灰度均值为,背景的灰度较暗,其灰度均值为,因此有
现若规定一门限值对图像进行分割,势必会产生将目标划分为背景和将背景划分为目标这两类错误。
通过适当选择阈值,可令这两类错误概率为最小,则该阈值即为最佳阈值。
把目标错分为背景的概率可表示为
把背景错分为目标的概率可表示为
总的误差概率为
为求得使误差概率最小的阈值,可将对求导并令导数为零,可得
代换后,可得
此时,若设,则有
若还有的条件,则
这时的最优阈值就是两类区域灰度均值与的平均值。
上面的推导是针对图像灰度值服从正态分布时的情况,如果灰度值服从其它分布,依理也可求出最优阈值来。
一般情况下,在不清楚灰度值分布时,通常可假定灰度值服从正态分布。
因此,本课题中亦可使用此方法来求得最优阈值,来对实验图像进行分割。
1.3.2最优阈值的迭代算法
在实际使用最优阈值进行分割的过程中,需要利用迭代算法来求得最优阈值。
设有一幅数字图像,混有加性高斯噪声,可表示为
此处假设图像上各点的噪声相互独立,且具有零均值,如果通过阈值分割将图像分为目标与背景两部分,则每一部分仍然有噪声点随机作用于其上,于是,目标和可表示为
迭代过程中,会多次地对和求均值,则
可见,随着迭代次数的增加,目标和背景的平均灰度都趋向于真实值。
因此,用迭代算法求得的最佳阈值不受噪声干扰的影响。
利用最优阈值对实验图像进行分割的迭代步骤为:
(1)确定一个初始阈值,可取为
式中,和为图像灰度的最小值和最大值。
(2)利用第k次迭代得到的阈值将图像分为目标和背景两大区域,其中
(3)计算区域和的灰度均值和。
(4)计算新的阈值,其中
(5)如果小于允许的误差,则结束,否则,转步骤
(2)。
利用迭代法求得最优阈值后,仍需进行一些人工调整才能将此阈值用于实验图像的分割,这是因为,这种最优阈值仍然属于全局阈值,它利用了图像中所有像素点的信息,但当光照不均匀时,图像中部分区域的灰度值可能差距较大,造成计算出的最优阈值分割效果不理想,此时,可设一人工经验因子进行校正。
四、实验步骤及程序
1、实验步骤:
(1)利用最优阈值对实验图像进行分割的迭代步骤编写程序流程图;
(2)编写程序,用Matlab语言实现此算法,完成选择图像的分割。
理解最优阈值迭代算法,设计程序实现对自选图像的最优阈值分割。
2、程序流程图:
3、程序如下:
a=imread('e:
/lena.bmp');
figure
(1)
imshow(a)
b=a(:
);
Smax=max(b);
Smin=min(b);
T0=(Smax+Smin)/2;
delta=1;
whiledelta>=0.1
clearR1R2
m=1;
l=1;
forn=1:
65536;
ifb(n)<=T0
R1(m)=b(n);
m=m+1;
else
R2(l)=b(n);
l=l+1;
end
end
Tnext=0.5*(mean(R1)+mean(R2));
delta=abs(Tnext-T0);
T0=Tnext;
end
forn=1:
256
form=1:
256
ifa(n,m)>=T0
c(n,m)=1;
else
c(n,m)=0;
end
end
end
figure
(2)
imshow(c)
五、实验结果与分析
分割域值:
T0=124.3653
分割前原图:
分割后效果图:
结论:
利用迭代法求得最优阈值仍然属于全局阈值。
由图可知分割后的效果图有两个灰度级,即分为目标背景两个部分。
将原图中的每一个像素与Tnest作比较,大于的为目标,小于的为背景。
实验2K均值聚类算法
2.1实验目的
将模式识别方法与图像处理技术相结合,掌握利用均值聚类算法进行图像分类的基本方法,通过实验加深对基本概念的理解。
2.2实验仪器设备及软件
HPD538、MATLAB、WIT
2.3实验原理
均值聚类法分为如下几个步骤:
一、初始化聚类中心
1、凭经验选择。
根据具体问题,凭经验从样本集中选出个比较合适的样本作为初始聚类中心。
2、用前个样本作为初始聚类中心。
3、将全部样本随机地分成类,计算每类的样本均值,将样本均值作为初始聚类中心。
4、密度法。
以每个样本为球心,用某个正数为半径作一个球形邻域,落在邻域内的样本数为该点密度,选密度最大点为第一初始聚类中心。
在离开第一点规定距离范围外确定次大密度点,以避免初始聚类中心聚集。
5、从个聚类划分的解中产生个聚类划分的初始聚类中心。
先把全部样本看作一个聚类,其聚类中心为样本的总均值;然后确定两聚类问题的聚类中心是一聚类问题的总均值和离它最远的点;以此类推。
二、初始聚类
1、按就近原则将样本归入各聚类中心所代表的类中。
2、取一样本,将其归入与其最近的聚类中心的那一类中,重新计算样本均值,更新聚类中心。
然后取下一样本,重复操作,直至所有样本归入相应类中。
三、判断聚类是否合理
采用误差平方和准则函数判断聚类是否合理,不合理则修改分类。
循环进行判断、修改直至达到算法终止条件。
2.4实验步骤及程序
1、实验步骤
(1)利用K-均值聚类法的算法步骤编写程序流程图;
(2)编写程序,用Matlab语言实现此算法,完成选择图像的分割。
(3)利用WIT实现K均值聚类算法的图像分割
2、程序流程图
WTI聚类程序如图:
3、程序如下:
tic
img=imread('e:
/lena.bmp');
[m,n]=size(img);
subplot(2,2,1),imshow(img);title('原始图')
subplot(2,2,2),imhist(img);title('灰度直方图')
holdoff;
img=double(img);
fori=1:
200
c1
(1)=25;
c2
(1)=125;
c3
(1)=200;
r=abs(img-c1(i));
g=abs(img-c2(i));
b=abs(img-c3(i));
r_g=r-g;
g_b=g-b;
r_b=r-b;
n_r=find(r_g<=0&r_b<=0);
n_g=find(r_g>0&g_b<=0);
n_b=find(g_b>0&r_b>0);
i=i+1;
c1(i)=sum(img(n_r))/length(n_r);c2(i)=sum(img(n_g))/length(n_g);c3(i)=sum(img(n_b))/length(n_b);d1(i)=(c1(i)-c1(i-1))^2;
d2(i)=(c2(i)-c2(i-1))^2;
d3(i)=(c3(i)-c3(i-1))^2;
ifd1(i)<=0.001&&d2(i)<=0.001&&d3(i)<=0.001
R=c1(i);
G=c2(i);
B=c3(i);
k=i;
break;
end
end
R
G
B
img=uint8(img);
img(find(img img(find(img>R&img img(find(img>G))=255; toc subplot(2,2,3),imshow(img);title('聚类后图象') subplot(2,2,4),imhist(img);title('聚类后灰度直方图') 2.5实验结果与分析 1、MATLAB实验结果: 聚类类别数: K=3 聚类中心: R=68.3649 G=129.3
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 模式识别 实验 报告