发明名称 利用流形谱聚类进行图像分割的方法
摘要 本发明公开了一种利用流形谱聚类进行图像分割的方法,主要解决现有方法的存储规模大、计算效率和分割精度低的问题。实现步骤为:(1)输入一幅图像,提取输入图像的颜色和纹理特征,并利用分水岭算法获得输入图像的流形集;(2)计算流形特征集,构造距离矩阵,由弗洛伊德算法得流行距离矩阵;(3)计算相似度矩阵,进而构造度矩阵以及归一化拉普拉斯矩阵;(4)特征分解归一化拉普拉斯矩阵,进而构造谱矩阵;(5)对谱矩阵进行归一化,得到归一化谱矩阵,由k均值算法得流形集的标签向量,输出分割结果。本发明存储规模小、计算效率和分割精度高,可应用于医学图像检测病灶区、精密零件表面缺陷检测、卫星拍摄的地形地貌照片的处理。
申请公布号 CN102024262B 申请公布日期 2012.07.04
申请号 CN201110001586.X 申请日期 2011.01.06
申请人 西安电子科技大学 发明人 郑喆坤;刘娟;沈彦波;焦李成;尚荣华;李阳阳;马文萍;王爽;公茂果
分类号 G06T7/00(2006.01)I 主分类号 G06T7/00(2006.01)I
代理机构 陕西电子工业专利中心 61205 代理人 王品华;朱红星
主权项 1.一种利用流形谱聚类进行图像分割的方法,包括如下步骤:1)输入一幅图像,在Luv颜色空间中获得该图像的颜色特征F<sub>C</sub>={f<sub>L</sub>,f<sub>u</sub>,f<sub>v</sub>},并提取图像的十维小波纹理特征F<sub>W</sub>={f<sub>1</sub>,f<sub>2</sub>,...,f<sub>10</sub>},采用特征融合的方法将颜色特征与小波纹理特征进行融合,得到图像的融合特征F=aF<sub>C</sub>+bF<sub>W</sub>,其中F的每一行表示一个像素点的融合特征,f<sub>L</sub>、f<sub>u</sub>、f<sub>v</sub>分别表示Luv颜色空间的亮度分量L、色度坐标分量u和色度坐标分量v的特征,f<sub>i</sub>表示小波纹理特征F<sub>W</sub>第i维特征i=1,2...10,a和b分别表示颜色特征和纹理特征所占的比例,且a+b=1;2)采用分水岭算法将该输入图像预分割为流形集<img file="FDA0000138380210000011.GIF" wi="258" he="61" />将各流形所包含的所有像素点融合特征的平均值分别作为各流形的流形特征,得到流形特征集<img file="FDA0000138380210000012.GIF" wi="325" he="74" />其中流形m<sub>i</sub>为分水岭算法预分割输入图像所得的第i个区域,i=1,2,...,n,n表示流形集M的大小,mf<sup>i</sup>为第i个流形m<sub>i</sub>的流形特征;3)计算每个流形特征mf<sup>i</sup>与其他流形特征之间的欧氏距离,依据每个流形特征mf<sup>i</sup>的所有欧氏距离寻找到流形m<sub>i</sub>的k个近邻流形,构造一个大小为n*n的距离矩阵:d=[d<sub>ij</sub>]<sub>n*n</sub>当第j个流形m<sub>j</sub>是m<sub>i</sub>的k个近邻中的一个时,d<sub>ij</sub>等于mf<sup>i</sup>与mf<sup>j</sup>之间的欧氏距离,否则为-1表示不连通,且d<sub>ii</sub>=0,其中i,j=1,2,...,n,k的取值范围为5到9之间的整数,依据不同的图像人为选择不同的值,d<sub>ij</sub>表示距离矩阵d的第i行第j列元素,d<sub>ii</sub>表示距离矩阵d的第i行第i列元素;4)使用弗洛伊德算法根据距离矩阵d计算各个流形之间的最短路径,得到大小为n*n的流形距离矩阵:D=[D<sub>ij</sub>]<sub>n*n</sub>,其中D<sub>ij</sub>为第i个流形m<sub>i</sub>与第j个流形m<sub>j</sub>之间的最短路径,其中D<sub>ij</sub>表示流形距离矩阵D中第i行第j列元素;5)构造大小为n*n相似度矩阵:W=[W<sub>ij</sub>]<sub>n*n</sub>,其中第i行第j列的元素<img file="FDA0000138380210000021.GIF" wi="249" he="124" />i≠j,对角线元素为0,且当D<sub>ij</sub>=-1时W<sub>ij</sub>等于0;6)计算度矩阵A=[A<sub>ij</sub>]<sub>n*n</sub>,其中第i行第i列的元素<img file="FDA0000138380210000022.GIF" wi="241" he="119" />其它元素为0,并计算归一化拉普拉斯矩阵L=A<sup>-1/2</sup>WA<sup>-1/2</sup>;7)对归一化拉普拉斯矩阵L进行特征分解,并取该L的前K个最大的特征值对应的特征向量构造谱矩阵:V={v<sub>1</sub>,v<sub>2</sub>,...,v<sub>K</sub>},其中K为人为设定的图像分类数,v<sub>i</sub>为L的第i个特征向量i=1,2,...K;8)计算归一化谱矩阵Y=[Y<sub>ij</sub>]<sub>n*n</sub>,其中第i行第j列的元素<img file="FDA0000138380210000023.GIF" wi="392" he="98" />V<sub>ij</sub>为谱矩阵V的第i行第j列的元素;9)采用k均值聚类算法,将归一化谱矩阵Y的各行聚为K类,得到聚类标签向量<img file="FDA0000138380210000024.GIF" wi="225" he="61" />其中c<sub>i</sub>表示第i个流形m<sub>i</sub>被分割为第c<sub>i</sub>类,c<sub>i</sub>∈{1,2,...K};10)将图像中有相同标签的流形分配同一种颜色,输出分割后的图像。
地址 710071 陕西省西安市太白南路2号