CN102332163A - 基于三层fcm聚类的小波域多尺度非监督纹理分割方法 - Google Patents
基于三层fcm聚类的小波域多尺度非监督纹理分割方法 Download PDFInfo
- Publication number
- CN102332163A CN102332163A CN201110324778A CN201110324778A CN102332163A CN 102332163 A CN102332163 A CN 102332163A CN 201110324778 A CN201110324778 A CN 201110324778A CN 201110324778 A CN201110324778 A CN 201110324778A CN 102332163 A CN102332163 A CN 102332163A
- Authority
- CN
- China
- Prior art keywords
- scale
- yardstick
- segmentation result
- fcm
- layer
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Landscapes
- Image Analysis (AREA)
Abstract
本发明涉及图像处理领域,公开了一种基于三层FCM聚类的小波域多尺度非监督纹理分割方法。采用第一层FCM分别获取每一尺度上低频部分和高频部分的分割结果;计算每一分割结果的局部标记直方图,并将其作为特征向量进行第二层FCM聚类,从而获取每一尺度上的分割结果;从大尺度到小尺度逐尺度进行分割结果融合,即计算相邻尺度分割结果的局部标记直方图,根据小波变换的采样关系,将对应位置的标记直方图组织在一起形成新的特征向量进行第三层FCM聚类,从而获取每一尺度最终的分割结果。最细尺度的分割结果作为本算法最终的分割结果。本发明简单实用,可用于纹理影像和航拍图像的分割。
Description
技术领域
本发明属于图像处理技术领域,涉及一种图像分割方法,可用于纹理图像和航拍图像的分割。
背景技术
计算机视频主要任务就是通过对采集的图片或视频进行处理以获得相应场景的三维信息,涉及图像处理、模式识别和图像理解。图像分割就是指把影像分成各具特性的区域并提取感兴趣目标的过程。图像分割是由图像处理到图像分析的关键步骤,也是一种基本的计算机视觉技术。这是因为图像的分割、目标的分离、特征的提取和参数的测量将原始图像转化为更为紧凑的形式,使得更高层次的分析和理解成为可能。然而,图像分割本身是一个病态问题,分割的目的是为了理解,但是理想的分割结果往往需要理解后的结果作为先验知识,这种病态性给问题的解决带来极大的困难,成为阻碍计算机视觉发展的一个主要瓶颈。但正因为图像分割是计算机视觉的一个基本问题,分割结果对***的性能影响很大,因此图像分割始终是图像处理领域的一个研究热点。
图像分割在诸如医学、安全、遥感和工业自动化等应用许多应用中,它都是具有关键作用的一种低层图像处理技术。但纹理图像分割依然存在很大的困难。目前来看,纹理分割算法主要分为基于特征的算法、基于模型的算法和基于结构的算法三种。
基于特征的方法首先提取纹理特征,然后再进行纹理的分类与识别,最终结果的好坏主要取决于纹理特征的优劣。常见的纹理特征有灰度共生矩阵、纹理能量、Gabor特征、多尺度多方向特征MPMS等。分割结果的好坏完全取决于纹理特征是否能准确描述纹理特性。
基于模型的方法采用统计模型对纹理进行建模,并在贝叶斯框架下实现纹理分割,所采用模型是否能准确描述纹理的分布特性直接决定了分割结果的好坏,而且越复杂的模型实现起来越耗时。常见的统计模型有:MRF、MRMRF、HMT、MSMRF等。
基于结构的方法采用语法规则将纹理基元组织在一起描述纹理,分割结果的好坏受到基元定义和语法规则准确性的双重约束,因该方法不够灵活而在应用受到很大的限制。
发明内容
为了避免了复杂纹理模型或纹理特征的设计,本发明提出了一种基于三层FCM聚类的小波域多尺度非监督纹理分割方法,以得到较好的分割结果,并将其应用于纹理图像和航拍图像。
FCM聚类算法,是Bezkek于1981年提出的一种模糊聚类算法,是一种逐步迭代的算法,本发明是基于三层FCM聚类的小波域多尺度大非监督纹理分割方法。
实现本发明目的的技术方案是:
第一层FCM用于计算每一尺度上不同子带的分割结果;第二层FCM用以计算子带分割结果的融合结果;第三层FCM用以计算相邻尺度分割结果的融合结果。具体过程如下:
(1)输入待分割图像、给定图像的分类数Nc和小波变换的最大尺度J;
(2)对图像进行J层小波变换。将三个高频子带系数矩阵对应位置的小波系数构成向量,作为该尺度高频部分的特征fj,H;低频部分的特征fj,L直接采用低频子带小波系数矩阵wj,LL;
(3)采用第一层FCM分别对低频和高频部分的特征进行聚类,获取每一尺度上低频与高频部分的分割结果,分别记为segj,L和segj,H;
(4)每一尺度上,采用第二层FCM对segj,L和segj,H进行融合,获取单尺度分割结果xj;
本发明与现有技术相比具有以下优点:
1、本发明设计了一个全新的纹理分割框架,它不需要设计复杂的纹理模型和纹理特征。基本的思路是融合不同尺度、不同频带上的纹理信息分割结构。
2、本发明设计的纹理分割框架能够集成所有的聚类算法以及它们的组合,也可以方便地扩展到其它的多尺度空间中进行分析,具有较好的扩展性。
3、仿真方法表明,本发明设计的方法对于纹理图像和遥感图像具有较好的分割结果。
附图说明
图1是本发明的流程示意图
图2是本发明应用于合成纹理影像1的分割结果
图3是本发明应用于合成纹理影像2的分割结果
图4是本发明应用于航拍遥感图像1的分割结果
图5是本发明应用于航拍遥感图像2的分割结果
具体实施方式
下面通过具体实施例对本发明作进一步详细描述。
参照图1,本发明的具体实现过程如下:
步骤一输入图像:输入待分割图像、给定图像的分类数Nc和小波变换的最大尺度J。
步骤二 离散小波变换(DWT,discrete wavelet transform):对图像进行J层小波变换,每一尺度上包含一个低频子带小波系数矩阵wj,LL和三个高频子带小波系数矩阵wj,LL、wj,LH和wj,HH。将三个高频子带系数矩阵对应位置的小波系数构成向量,作为该尺度高频部分的特征fj,H;低频部分的特征fj,L直接采用低频子带小波系数矩阵wj,LL。具体为:
(2.1)对输入图像进行J层小波变换;
(2.2)对每一尺度的低频和高频部分的特征分别按照如下形式构造低频特征和高频特征:
其中(r,c)为存储位置。
步骤三 第一层FCM:采用第一层FCM分别对低频和高频部分的特征进行聚类,获取每一尺度上低频与高频部分的分割结果,分别记为segj,L和segj,H。
设定低频部分和高频部分的聚类数K1L和K1H,在每一尺度j上执行如下操作:
(3.1)对特征矩阵fj,L和fj,H分别进行第二层FCM聚类。为描述方便,采用如下符号表示有关信息:特征矩阵大小M×N、聚类数K1j,b、特征矩阵fj,b(b=L,H)。聚类过程描述如下:
(a)初始化:初始化大小为M×N×K1j,b的隶属度矩阵Uj,b,每一元素均在[0,1]上取值;
(b)迭代过程:
更新聚类中心:
更新隶属度矩阵:
(c)判断是否满足终止条件:迭代次数是否达到最大次数或者两次迭代的隶属度之间的最大差值是否小于给定阈值epsilon。如果满足终止条件则算法结束;否则继续转步骤(3b2)继续执行。
步骤四 第二层FCM聚类:每一尺度上,采用第二层FCM对segj,L和segj,H进行融合,获取单尺度分割结果xj。
在每一尺度j上执行如下操作:
(a)在给定的窗口大小内,计算segj,b(b=L,H)的局部标记直方图fj,b:
其中,k1=1,2,…,K1j,b, u为以(r,c)为中心的窗口Nr,c内像素的个数。
(b)将fj,L和fj,H对应位置的值串在一起构成尺度j上的特征图像fj:
(c)给定聚类数K,对fj中每一位置的特征采用第二层FCM进行聚类,获取尺度j上的单尺度分割结果xj。具体过程与步骤三类似。
(a)令j=J-1,
(f)设定j=j-1,转(b)继续执行。
本发明的内容可以通过以下的仿真结果进一步进行说明。
1、仿真内容:应用本发明方法与MRMRF方法分别对两幅幅合成灰度纹理影像和两幅航拍影像进行影像分割。两种方法均采用3层“haar”小波变换。MRMRF算法的特征场采用GMRF进行建模,标记场采用只考虑双点势团的二阶MLL模型进行建模,每一尺度上的势函数取值为1.0。本发明方法设定K1L=2×K,K1H=6×K,在计算局部标记直方图时采用的窗口大小随尺度而改变,尺度0(原始分辨率)、1、2、3上采用的窗口大小依次为3、5、7、3。合成纹理影像采用错分率作为客观评价标准。错分率的定义为:
2、仿真结果
图2为本发明方法应用于合成灰度纹理影像的分割结果。其中图2(a)为第一幅合成灰度纹理影像,含有两类纹理;图3(a)为第二幅纹理影像,含有4类纹理。图2(b)和图3(b)分别为图2(a)和图3(a)的真实分割结果;图2(c)和图3(c)为MRMRF算法分割结果;图2(d)和图3(d)为本文算法分割结果。
从图2可以看出,本发明方法的分割结果与MRMRF算法具有相近的能力。这表明本发明的所设计的纹理分割框架能够达到基于模型的分割算法的分割精度。而从图3的分割结果还可以看出,本发明方法的结果要明显优于MRMRF算法。这也表明本发明方法具有不错的分类能力。从对应的错分率上也可明显看出,本发明方法能够能够达到甚至超越MRMRF算法的分割效果
图4和图5为本发明方法应用于航拍影像的分割结果。其中图4(a)为第一幅航拍影像,含有两类地物覆盖。图4(b)和图4(c)分别为图4(a)的MRMRF算法分割结果和本发明方法的分割结果。其中图5(a)为第二幅航拍影像,含有三类地物覆盖。图5(b)和图5(c)分别为图5(a)的MRMRF算法分割结果和本发明方法的分割结果。
从图4和图5可以看出,MRMRF分割方法分割结果的区域一致性有所欠缺。而本发明方法则在区域性和边界性两方面均有较好的结果。这进一步表明本发明方法具有一定的实用性。
综上,本发明设计了一种基于三层FCM聚类的小波域多尺度非监督纹理分割算法,作为一种全新的纹理分割框架,能够避免设计复杂的纹理模型或纹理特征。这个结论从对合成纹理影像和航拍图像的分割结果可以明显看出。
Claims (4)
1.一种基于三层FCM聚类的小波域多尺度非监督纹理分割方法,主要包括如下过程:
(1)输入待分割图像、给定图像的分类数Nc和小波变换的最大尺度J;
(2)对图像进行J层小波变换,每一尺度上包含一个低频子带小波系数矩阵wj,LL和三个高频子带小波系数矩阵wj,LL、wj,LH和wj,HH,将三个高频子带系数矩阵对应位置的小波系数构成向量,作为该尺度高频部分的特征fj,H;低频部分的特征fj,L直接采用低频子带小波系数矩阵wj,LL;
(3)采用第一层FCM分别对低频和高频部分的特征进行聚类,获取每一尺度上低频与高频部分的分割结果,分别记为segj,L和segj,H;
(4)每一尺度上,采用第二层FCM对segj,L和segj,H进行融合,获取单尺度分割结果xj;
3.根据权利要求1所述的分割方法,其中步骤(3)按如下过程进行:
(3a)设fj,L和fj,H的聚类数分别为K1L和K1H;
(3b)对特征矩阵fj,L和fj,H分别进行第二层FCM聚类,采用如下符号表示有关信息:特征矩阵大小M×N、聚类数K1j,b、特征矩阵fj,b(b=L,H),聚类过程描述如下:
(3b1)初始化:初始化大小为M×N×K1j,b的隶属度矩阵Uj,b,每一元素均在[0,1]上取值;
(3b2)迭代过程:
更新聚类中心:
更新隶属度矩阵:
(3b3)判断是否满足终止条件:迭代次数是否达到最大次数或者两次迭代的隶属度之间的最大差值是否小于给定阈值epsilon,如果满足终止条件则算法结束;否则继续转步骤(3b2)继续执行;
4.根据权利要求1所述的分割方法,其中步骤(4)按如下过程进行:
(4a)在给定的窗口大小内,计算segj,b(b=L,H)的局部标记直方图fj,b:
(4b)将fj,L和fj,H对应位置的值串在一起构成尺度j上的特征图像fj:
(4c)给定聚类数K,对fj中每一位置的特征采用第二层FCM进行聚类,获取尺度j上的单尺度分割结果xj。
5.根据权利要求1所述的分割方法,其中步骤(5)按如下过程执行:
(5f)设定j=j-1,转(5b)继续执行。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110324778 CN102332163B (zh) | 2011-10-24 | 2011-10-24 | 基于三层fcm聚类的小波域多尺度非监督纹理分割方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110324778 CN102332163B (zh) | 2011-10-24 | 2011-10-24 | 基于三层fcm聚类的小波域多尺度非监督纹理分割方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102332163A true CN102332163A (zh) | 2012-01-25 |
CN102332163B CN102332163B (zh) | 2013-07-31 |
Family
ID=45483926
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201110324778 Expired - Fee Related CN102332163B (zh) | 2011-10-24 | 2011-10-24 | 基于三层fcm聚类的小波域多尺度非监督纹理分割方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102332163B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108021920A (zh) * | 2017-11-09 | 2018-05-11 | 华南理工大学 | 一种图像对象协同发现的方法 |
CN108694720A (zh) * | 2018-04-25 | 2018-10-23 | 南京邮电大学 | 一种基于多尺度纹理空间聚类的高分影像建筑区提取方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2001041071A1 (en) * | 1999-12-03 | 2001-06-07 | Samsung Electronics Co., Ltd. | Texture description method and texture-based image retrieval method using gabor filter in frequency domain |
CN1822046A (zh) * | 2006-03-30 | 2006-08-23 | 上海电力学院 | 基于区域特征模糊的红外与可见光图像融合方法 |
CN101551905A (zh) * | 2009-05-08 | 2009-10-07 | 西安电子科技大学 | 基于模糊c均值聚类和空间信息的多尺度纹理图像分割方法 |
CN102081791A (zh) * | 2010-11-25 | 2011-06-01 | 西北工业大学 | 一种基于多尺度特征融合的sar图像分割方法 |
-
2011
- 2011-10-24 CN CN 201110324778 patent/CN102332163B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2001041071A1 (en) * | 1999-12-03 | 2001-06-07 | Samsung Electronics Co., Ltd. | Texture description method and texture-based image retrieval method using gabor filter in frequency domain |
CN1822046A (zh) * | 2006-03-30 | 2006-08-23 | 上海电力学院 | 基于区域特征模糊的红外与可见光图像融合方法 |
CN101551905A (zh) * | 2009-05-08 | 2009-10-07 | 西安电子科技大学 | 基于模糊c均值聚类和空间信息的多尺度纹理图像分割方法 |
CN102081791A (zh) * | 2010-11-25 | 2011-06-01 | 西北工业大学 | 一种基于多尺度特征融合的sar图像分割方法 |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108021920A (zh) * | 2017-11-09 | 2018-05-11 | 华南理工大学 | 一种图像对象协同发现的方法 |
CN108694720A (zh) * | 2018-04-25 | 2018-10-23 | 南京邮电大学 | 一种基于多尺度纹理空间聚类的高分影像建筑区提取方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102332163B (zh) | 2013-07-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Dong et al. | High quality multi-spectral and panchromatic image fusion technologies based on curvelet transform | |
CN101551863B (zh) | 基于非下采样轮廓波变换的遥感图像道路提取方法 | |
Tzotsos et al. | Object-based image analysis through nonlinear scale-space filtering | |
CN108491849A (zh) | 基于三维稠密连接卷积神经网络的高光谱图像分类方法 | |
CN102169584B (zh) | 基于分水岭和treelet的遥感图像变化检测方法 | |
CN100573559C (zh) | 基于小波和均值漂移的自适应多尺度纹理图像分割方法 | |
Wang et al. | Multifocus image fusion using convolutional neural networks in the discrete wavelet transform domain | |
CN107067405B (zh) | 基于尺度优选的遥感影像分割方法 | |
Alidoost et al. | A CNN-based approach for automatic building detection and recognition of roof types using a single aerial image | |
Abdollahi et al. | Extraction of road features from UAV images using a novel level set segmentation approach | |
Wu et al. | Classification of defects with ensemble methods in the automated visual inspection of sewer pipes | |
CN109766858A (zh) | 结合双边滤波的三维卷积神经网络高光谱影像分类方法 | |
CN105069796B (zh) | 基于小波散射网络的sar图像分割方法 | |
CN104794477B (zh) | 基于3-d小波变换和稀疏张量的高光谱图像特征抽取方法 | |
CN103700089B (zh) | 一种三维医学图像多尺度异构特征的提取与分类方法 | |
CN101587189B (zh) | 用于合成孔径雷达图像的纹理基元特征提取方法 | |
CN106408030A (zh) | 基于中层语义属性和卷积神经网络的sar图像分类方法 | |
CN103324940A (zh) | 基于多示例多标签学习的皮肤病理图片特征识别方法 | |
Dey et al. | Voronoi-based feature curves extraction for sampled singular surfaces | |
CN102122353A (zh) | 利用增量字典学习与稀疏表示进行图像分割的方法 | |
CN104182985A (zh) | 遥感图像变化检测方法 | |
CN104318563A (zh) | 一种基于医学图像的器官骨架提取方法 | |
CN112115871B (zh) | 适用于行人目标检测的高低频交织边缘特征增强方法 | |
CN104156723B (zh) | 一种具有尺度不变性的最稳定极值区域的提取方法 | |
CN104732551A (zh) | 基于超像素和图割优化的水平集图像分割方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20130731 Termination date: 20171024 |