CN103871060A - 基于平稳方向波域概率图模型的图像分割方法 - Google Patents
基于平稳方向波域概率图模型的图像分割方法 Download PDFInfo
- Publication number
- CN103871060A CN103871060A CN201410092303.0A CN201410092303A CN103871060A CN 103871060 A CN103871060 A CN 103871060A CN 201410092303 A CN201410092303 A CN 201410092303A CN 103871060 A CN103871060 A CN 103871060A
- Authority
- CN
- China
- Prior art keywords
- image
- split
- frequency sub
- pixel
- fritter
- 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
本发明公开了一种基于平稳方向波域概率图模型的图像分割方法,克服了现有技术中不能充分挖掘待分割图像中的方向性信息、无法提取包含更多图像信息的训练图像块的问题。本发明的实现步骤为:(1)输入图像;(2)计算特征向量;(3)提取训练图像块;(4)求解先验概率;(5)构建隐马尔科夫链模型参数集;(6)更新隐马尔科夫链模型参数集;(7)求解最大似然值;(8)得到最终的分割图像;(9)输出分割图像。本发明具有边缘准确性高和区域一致性的优点,可应用于在目标识别中对感兴趣区域的获取。
Description
技术领域
本发明属于图像处理技术领域,更进一步涉及图像分割技术领域中的一种基于平稳方向波域概率图模型的图像分割方法。本发明可应用于在目标识别中对感兴趣区域的获取。
背景技术
图像分割是图像处理与计算机视觉领域中的一项关键技术。它是由图像处理到图像分析的关键步骤,分割结果的优劣直接影响到随后的图像分析,理解和求解等问题。图像分割的目的是把图像分成各具特性的不同区域并将其中感兴趣的目标提取出来。这里特性可以是像素的灰度、颜色、纹理,对应的目标可以是单个区域或多个区域。
近年来,作为一种有向无环的概率图模型,隐马尔科夫树模型(HMT)被广泛的应用于图像分割中,该模型可以充分挖掘多尺度变换域系数之间的统计特性,从而将系数分布的求解问题转化为隐马尔科夫树模型中对应的隐状态确定问题。
西安电子科技大学在其专利申请“基于小波域隐马尔科夫树模型的图像分割方法”(专利申请号:200810150924.4,公开号:CN101350099A)中提出了一种基于小波域隐马尔科夫树模型的图像分割方法。该方法从待分割图像中选取训练图像块,对训练图像块进行小波变换得到训练数据,采用最大期望(EM)算法建立隐马尔科夫树模型参数,并利用最大似然算法求解似然值;然后,采用上下文背景进行多尺度融合,得到最终的分割图像。该方法虽然具有区域一致性好和边缘准确的优点,但是仍然存在的不足是,该方法采用小波变换只能描述点奇异性,不能有效地刻画待分割图像中的二维奇异性,因此待分割图像在分割过程中会丢失部分边缘和细节信息。
西安电子科技大学在其专利申请“基于非下采样Contourlet变换的HMT图像分割方法”(专利申请号:200810232336.5,公开号:CN101447080A)中提出了一种基于非下采样轮廓波(Contourlet)变换的隐马尔科夫书模型图像分割方法。该方法从待分割图像中提取训练图像,并对训练图像进行非下采样轮廓波(Coutourlet)变换,采用最大期望(EM)算法建立隐马尔科夫树模型;然后,使用标记树融合得到最大后验概率分类,结合隐马尔科夫随机场模型和像素空间相关性信息,得到最终分割结果。该方法虽然具有边缘准确和区域一致性好的优点,但是仍然存在的不足是,该方法在提取训练图像时没有充分挖掘待分割图像中的纹理信息,导致分割后的图像中误分类效应比较明显。
发明内容
本发明针对上述现有技术的不足,提出一种基于平稳方向波域概率图模型的图像分割方法。本发明能够较好的保持区域的一致性和边缘的准确性,提高了图像分割结果的正确率。
为实现上述目的,本发明的具体步骤如下:
(1)输入任选的一幅待分割图像;
(2)计算特征向量:
(2a)对待分割图像进行多尺度的平稳方向波变换,得到不同尺度的低频子带系数和高频子带系数;
(2b)按照下式,计算待分割图像中像素的特征向量:
其中,T表示待分割图像中像素的特征向量,d表示待分割图像中16×16大小的窗口内像素的平稳方向波子带系数,|·|表示取绝对值操作,∑表示求和操作;
(3)提取训练图像块:
(3a)采用标记分水岭方法,对待分割图像进行分割,得到一副划分为若干图像小块的粗分割图像;
(3b)将待分割图像中像素的特征向量作为粗分割图像中相同位置像素的特征向量,对粗分割图像中图像小块内所有像素的特征向量求和,得到粗分割图像中图像小块的特征向量;
(3c)采用模糊C均值聚类方法,对粗分割图像中图像小块的特征向量进行聚类,分别得到粗分割图像中图像小块的隶属度和类别;
(3d)从粗分割图像中类别相同的图像小块中选取隶属度大于0.8的图像小块,作为待分割图像的训练图像块;
(4)求解先验概率:
(4a)将粗分割图像中每个图像小块的类别值,作为该图像小块内所有像素的灰度值,得到初始分割图像;
(4b)按照下式,求解初始分割图像中像素的先验概率:
P(φ)=exp(-U(φ))
其中,P(φ)表示初始分割图像中像素的先验概率,φ表示初始分割图像中像素的类别,exp(·)表示指数操作,U(φ)表示初始分割图像中3×3大小的邻域窗口内像素点之间的能量函数;
(5)构建隐马尔科夫链模型参数集:
(5a)对待分割图像的高频子带系数,建立隐马尔科夫链模型;
(5b)构建隐马尔科夫链模型参数集;
(6)更新隐马尔科夫链模型参数集:
采用最大期望方法,更新隐马尔科夫链参数集中的期望、方差、概率分布函数和状态转移概率四个参数;
(7)采用最大似然方法,求解待分割图像中像素的最大似然值;
(8)利用最大后验估计公式,得到最终的分割图像;
(9)输出最终的分割图像。
本发明与现有的技术相比具有以下优点:
第一,由于本发明采用了平稳方向波变换对待分割图像进行多尺度变换的方法,克服了现有技术中不能充分挖掘待分割图像中的方向性信息而导致边缘分割不准确的缺点,有效的捕捉图像中的轮廓和边缘信息,使得本发明提高了边缘分割的准确性。
第二,由于本发明采用了基于标记分水岭方法的训练图像块提取的方法,克服了现有技术中不能提取包含更多图像信息的训练图像块而导致错分现象非常严重的缺点,使得本发明保持了区域的一致性,提高了图像分割结果的正确率。
附图说明
图1是本发明的流程图;
图2为本发明与现有技术对两类纹理图像分割效果对比图;
图3为本发明与现有技术对三类纹理图像分割效果对比图。
具体实施方式
下面结合附图对本发明做进一步的描述。
参照附图1,本发明的具体步骤如下:
步骤1:输入图像。
输入任选的一幅待分割图像。本发明实施例中使用的待分割图像分别如附图2(a)和附图3(a)所示。其中,附图2(a)是从Brodatz纹理图像库中选取的二类纹理图像,大小为256×256,附图3(a)是从Brodatz纹理图像库中选取的三类纹理图像,大小为256×256。
步骤2:计算特征向量。
对待分割图像进行多尺度的平稳方向波变换,得到不同尺度的低频子带系数和高频子带系数,其中,平稳方向波变换的实现步骤如下:
第一步,按照下式,构造变换矩阵:
其中,M表示变换矩阵,a1,a2,b1,b2分别表示变换矩阵M的四个整数元素值;
第二步,用待分割图像中像素所在的行数和列数构成1×2维矩阵,将该矩阵与变换矩阵相乘,得到待分割图像的陪集;
第三步,将陪集沿水平和垂直方向分别进行两次一维平稳小波变换和一次一维平稳小波变换,得到待分割图像的平稳方向波低频子带系数和高频子带系数。
按照下式,计算待分割图像中像素的特征向量:
其中,T表示待分割图像中像素的特征向量,d表示待分割图像中16×16大小的窗口内像素的平稳方向波子带系数,|·|表示取绝对值操作,∑表示求和操作。
步骤3:提取训练图像块。
采用标记分水岭方法,对待分割图像进行分割,得到一副划分为若干图像小块的粗分割图像,其中,标记分水岭方法的具体步骤如下:
第一步,对待分割图像进行中值滤波,得到滤波后图像;
第二步,按照下式,得到滤波后图像的梯度图像:
第三步,选取滤波后图像中像素灰度值的绝对值范围在0到5之间的连续区域作为内部标记;
第四步,对内部标记进行分水岭变换,得到内部标记的分水岭线,将这些分水岭线作为外部标记;
第五步,判断梯度图像中内部标记位置上像素的灰度值的绝对值是否大于5,若是,则执行第六步,否则,执行第七步;
第六步,设置像素的邻域窗口大小为3×3,对邻域窗口内其他像素的灰度值的绝对值求和,将求和的结果作为该像素修正后的灰度值;
第七步,对梯度图像进行分水岭变换,得到一副划分为若干图像小块的粗分割图像。
将待分割图像中像素的特征向量作为粗分割图像中相同位置像素的特征向量,对粗分割图像中图像小块内所有像素的特征向量求和,得到粗分割图像中图像小块的特征向量。
采用模糊C均值聚类方法,对粗分割图像中图像小块的特征向量进行聚类,分别得到粗分割图像中图像小块的隶属度和类别,其中,模糊C均值聚类方法的具体步骤如下:
第一步,从粗分割图像中任选A个图像小块的特征向量作为A类图像小块的聚类中心,其中A表示对粗分割图像设定的分割类别数;
第二步,按照下式,更新粗分割图像中图像小块的特征向量的隶属度:
其中,uik表示粗分割图像中序号为k的图像小块特征向量隶属于第i类的隶属度,i表示当前的聚类类别,k表示粗分割图像中图像小块的特征向量的序号,∑表示求和操作,xk表示粗分割图像中序号为k的图像小块特征向量,vi表示第i类的聚类中心,vj表示第j类的聚类中心,j表示聚类类别,j=1,2,...,A,A表示对粗分割图像设定的分割类别数,||·||表示求欧式距离操作,n表示粗分割图像中图像小块的特征向量的个数;
第三步,按照下式,更新聚类中心:
其中,vi表示第i类的聚类中心,i表示当前的聚类类别,∑表示求和操作,uik表示粗分割图像中序号为k的图像小块特征向量隶属于第i类的隶属度,k表示粗分割图像中图像小块的特征向量的序号,k=1,2,...,n,n表示粗分割图像中图像小块的特征向量的个数,xk表示粗分割图像中序号为k的图像小块特征向量;
第四步,重复执行第二步和第三步20次,得到粗分割图像中图像小块的特征向量的隶属矩阵;
第五步,将隶属矩阵中每列隶属度的最大值作为图像小块的隶属度,得到粗分割图像中图像小块的隶属度;
第六步,将隶属矩阵中每列隶属度的最大值的行数作为图像小块的类别,得到粗分割图像中图像小块的类别。
从粗分割图像中类别相同的图像小块中选取隶属度大于0.8的图像小块,作为待分割图像的训练图像块。
步骤4:求解先验概率。
将粗分割图像中每个图像小块的类别值,作为该图像小块内所有像素的灰度值,得到初始分割图像。
按照下式,求解初始分割图像中像素的先验概率:
P(φ)=exp(-U(φ))
其中,P(φ)表示初始分割图像中像素的先验概率,φ表示初始分割图像中像素的类别,exp(·)表示指数操作,U(φ)表示初始分割图像中3×3大小的邻域窗口内像素点之间的能量函数。
步骤5:构建隐马尔科夫链模型参数集。
对待分割图像的高频子带系数,建立隐马尔科夫链模型的步骤如下:
第一步,任意选取待分割图像的两个相邻尺度的高频子带系数;
第二步,将较大尺度的高频子带系数中的系数作为父系数,较小尺度的高频子带系数中的系数作为子系数;
第三步,连结父系数和子系数,得到隐马尔科夫链模型。
按照如下公式,构建隐马尔科夫链模型参数集:
θ={μm,σm,f(D),ε}
其中,θ表示隐马尔科夫链模型参数集,μm表示高频子带系数的隐状态变量为状态值m的期望,m表示高频子带系数的状态值,σm表示高频子带系数的隐状态变量为状态值m方差,f(D)表示高频子带系数的概率分布函数,D表示待分割图像的高频子带系数,ε表示高频子带系数的状态转移概率。
步骤6:更新隐马尔科夫链模型参数集。
采用最大期望方法,更新隐马尔科夫链参数集中的期望、方差、概率分布函数和状态转移概率四个参数,其中,最大期望方法的具体步骤如下:
第一步,将隐马尔科夫链模型参数集中的四个参数全部初始化为0;
第二步,按照下式,更新隐马尔科夫链模型参数集:
其中,θS+1表示循环次数为S+1的隐马尔科夫链模型参数集,S表示循环次数,argmin(·)表示使函数取最小值的操作,θ表示隐马尔科夫链模型参数集,E[·]表示求期望值操作,ln(·)表示取对数操作,p(H,w|θS)表示待分割图像的训练图像块的高频子带系数的后验概率,H表示待分割图像的训练图像块的高频子带系数,w表示待分割图像的训练图像块的高频子带系数的状态值,θS表示循环次数为S的隐马尔科夫链模型参数集;
第三步,将循环次数为S+1的隐马尔科夫链模型参数集与循环次数为S的隐马尔科夫链模型参数集相减,若差值小于1×10-4,则停止循环,得到最终的隐马尔科夫链模型参数集;否则,返回第二步继续循环。
步骤7:采用最大似然方法,求解待分割图像中像素的最大似然值。
第一步,按照下式,求解高频子带系数的似然值:
其中,F(D|θ)表示高频子带系数的似然值,D表示待分割图像的高频子带系数,θ表示隐马尔科夫链模型参数集,∑表示求和操作,m表示高频子带系数的状态值,β(m)表示高频子带系数的隐状态变量为状态值m时的条件似然概率,f(D)表示高频子带系数的概率分布函数;
第二步,将同一尺度下的所有高频子带系数的似然值相乘,将相乘结果作为每个尺度的高频子带系数的似然值;
第三步,求解所有尺度高频子带系数的似然值的最大值,得到待分割图像中像素的最大似然值。
步骤8:利用最大后验估计公式,得到最终的分割图像。
按照下面的公式,得到最终的分割图像:
其中,Y表示最终的分割图像,argmax(·)表示使函数取最大值的操作,φ表示初始分割图像中像素的类别,P(φ)表示初始分割图像中像素的先验概率,F表示待分割图像中像素的最大似然值。
步骤9:输出最终的分割图像。
输出最终的分割图像。
下面结合附图2和附图3的仿真图对本发明的效果做进一步说明。
1.仿真实验条件:
本发明的硬件测试平台是:处理器为Inter Core2Duo CPU E8200,主频为2.67GHz,内存2GB,软件平台为:Windows7旗舰版32位操作系统和MatlabR2012b。本发明的输入图像分别为二类纹理图像和三类纹理图像,大小都为256×256,格式都为PNG。
2.仿真内容:
本发明用到的现有技术进行对比的两个方法分别如下:
Hyeokho Choi等人在文献“Multiscale image segmentation usingwavelet-domain hidden markov models.Image Processing,IEEE Transactions on,10(9):1309-132,2001”中提到的基于小波域隐马尔科夫模型的多尺度图像分割方法,简称HMTseg方法。
Yuheng Sha等人在文献“Unsupervised image segmentation using contourletdomain hidden markov trees model.In Image Analysis and Recognition,pages32-39.Springer,2005.”中提出的基于轮廓波域隐马尔科夫树模型的无监督图像分割方法,简称CHMTseg方法。
3.仿真结果分析:
图2是本发明方法与现有技术对二类纹理图像分割效果对比图。其中,图2(a)为输入的二类纹理图像,其大小为256×256。图2(b)为采用HMTseg方法对二类纹理图像进行分割的效果图,图2(c)为采用CHMTseg方法对二类纹理图像进行分割的效果图,图2(d)为本发明方法对二类纹理图像进行分割的效果图。
图3是本发明方法与现有技术对三类纹理图像分割效果对比图。其中,图3(a)为输入的三类纹理图像,其大小为256×256。图3(b)为采用HMTseg方法对三类纹理图像进行分割的效果图,图3(c)为采用CHMTseg方法对三类纹理图像进行分割的效果图,图3(d)为本发明方法对三类纹理图像进行分割的效果图。
从图2(b)和图3(b)中可以看到,采用HMTseg方法得到的分割图像具有明显的锯齿状边缘,且存在少量大小不一的错分割杂块。从图2(c)和图3(c)中可以看出,采用CHMTseg方法得到的分割图像背景非常凌乱,误分类效应很明显。从图2(d)和图3(d)中可以看出,本发明的方法能够很好地保持区域的一致性和边缘的准确性,提高了图像分割结果的正确率。
Claims (9)
1.一种基于平稳方向波域概率图模型的图像分割方法,包括如下步骤:
(1)输入任选的一幅待分割图像;
(2)计算特征向量:
(2a)对待分割图像进行多尺度的平稳方向波变换,得到不同尺度的低频子带系数和高频子带系数;
(2b)按照下式,计算待分割图像中像素的特征向量:
其中,T表示待分割图像中像素的特征向量,d表示待分割图像中16×16大小的窗口内像素的平稳方向波子带系数,|·|表示取绝对值操作,∑表示求和操作;
(3)提取训练图像块:
(3a)采用标记分水岭方法,对待分割图像进行分割,得到一副划分为若干图像小块的粗分割图像;
(3b)将待分割图像中像素的特征向量作为粗分割图像中相同位置像素的特征向量,对粗分割图像中图像小块内所有像素的特征向量求和,得到粗分割图像中图像小块的特征向量;
(3c)采用模糊C均值聚类方法,对粗分割图像中图像小块的特征向量进行聚类,分别得到粗分割图像中图像小块的隶属度和类别;
(3d)从粗分割图像中类别相同的图像小块中选取隶属度大于0.8的图像小块,作为待分割图像的训练图像块;
(4)求解先验概率:
(4a)将粗分割图像中每个图像小块的类别值,作为该图像小块内所有像素的灰度值,得到初始分割图像;
(4b)按照下式,求解初始分割图像中像素的先验概率:
P(φ)=exp(-U(φ))
其中,P(φ)表示初始分割图像中像素的先验概率,φ表示初始分割图像中像素的类别,exp(·)表示指数操作,U(φ)表示初始分割图像中3×3大小的邻域窗口内像素点之间的能量函数;
(5)构建隐马尔科夫链模型参数集:
(5a)对待分割图像的高频子带系数,建立隐马尔科夫链模型;
(5b)构建隐马尔科夫链模型参数集;
(6)更新隐马尔科夫链模型参数集:
采用最大期望方法,更新隐马尔科夫链参数集中的期望、方差、概率分布函数和状态转移概率四个参数;
(7)采用最大似然方法,求解待分割图像中像素的最大似然值;
(8)利用最大后验估计公式,得到最终的分割图像;
(9)输出最终的分割图像。
2.根据权利要求1所述基于平稳方向波域概率图模型的图像分割方法,其特征在于,步骤(2a)所述的平稳方向波变换,按照如下步骤进行:
第一步,按照下式,构造变换矩阵:
其中,M表示变换矩阵,a1,a2,b1,b2分别表示变换矩阵M的四个整数元素值;
第二步,用待分割图像中像素所在的行数和列数构成1×2维矩阵,将该矩阵与变换矩阵相乘,得到待分割图像的陪集;
第三步,将陪集沿水平和垂直方向分别进行两次一维平稳小波变换和一次一维平稳小波变换,得到待分割图像的平稳方向波低频子带系数和高频子带系数。
3.根据权利要求1所述基于平稳方向波域概率图模型的图像分割方法,其特征在于,步骤(3a)所述的标记分水岭方法,按照如下步骤进行:
第一步,对待分割图像进行中值滤波,得到滤波后图像;
第二步,按照下式,得到滤波后图像的梯度图像:
第三步,选取滤波后图像中像素灰度值的绝对值范围在0到5之间的连续区域作为内部标记;
第四步,对内部标记进行分水岭变换,得到内部标记的分水岭线,将这些分水岭线作为外部标记;
第五步,判断梯度图像中内部标记位置上像素的灰度值的绝对值是否大于5,若是,则执行第六步,否则,执行第七步;
第六步,设置像素的邻域窗口大小为3×3,对邻域窗口内其他像素的灰度值的绝对值求和,将求和的结果作为该像素修正后的灰度值;
第七步,对梯度图像进行分水岭变换,得到一副划分为若干图像小块的粗分割图像。
4.根据权利要求1所述基于平稳方向波域概率图模型的图像分割方法,其特征在于,步骤(3c)所述的模糊C均值聚类方法,按照如下步骤进行:
第一步,从粗分割图像中任选A个图像小块的特征向量作为A类图像小块的聚类中心,其中A表示对粗分割图像设定的分割类别数;
第二步,按照下式,更新粗分割图像中图像小块的特征向量的隶属度:
其中,uik表示粗分割图像中序号为k的图像小块特征向量隶属于第i类的隶属度,i表示当前的聚类类别,k表示粗分割图像中图像小块的特征向量的序号,∑表示求和操作,xk表示粗分割图像中序号为k的图像小块特征向量,vi表示第i类的聚类中心,vj表示第j类的聚类中心,j表示聚类类别,j=1,2,...,A,A表示对粗分割图像设定的分割类别数,||·||表示求欧式距离操作,n表示粗分割图像中图像小块的特征向量的个数;
第三步,按照下式,更新聚类中心:
其中,vi表示第i类的聚类中心,i表示当前的聚类类别,∑表示求和操作,uik表示粗分割图像中序号为k的图像小块特征向量隶属于第i类的隶属度,k表示粗分割图像中图像小块的特征向量的序号,k=1,2,...,n,n表示粗分割图像中图像小块的特征向量的个数,xk表示粗分割图像中序号为k的图像小块特征向量;
第四步,重复执行第二步和第三步20次,得到粗分割图像中图像小块的特征向量的隶属矩阵;
第五步,将隶属矩阵中每列隶属度的最大值作为图像小块的隶属度,得到粗分割图像中图像小块的隶属度;
第六步,将隶属矩阵中每列隶属度的最大值的行数作为图像小块的类别,得到粗分割图像中图像小块的类别。
5.根据权利要求1所述基于平稳方向波域概率图模型的图像分割方法,其特征在于,步骤(5a)所述建立隐马尔科夫链模型的步骤如下:
第一步,任意选取待分割图像的两个相邻尺度的高频子带系数;
第二步,将较大尺度的高频子带系数中的系数作为父系数,较小尺度的高频子带系数中的系数作为子系数;
第三步,连结父系数和子系数,得到隐马尔科夫链模型。
6.根据权利要求1所述基于平稳方向波域概率图模型的图像分割方法,其特征在于,步骤(5b)所述构建隐马尔科夫链模型参数集是指按照如下公式进行:
θ={μm,σm,f(D),ε}
其中,θ表示隐马尔科夫链模型参数集,μm表示高频子带系数的隐状态变量为状态值m的期望,m表示高频子带系数的状态值,σm表示高频子带系数的隐状态变量为状态值m方差,f(D)表示高频子带系数的概率分布函数,D表示待分割图像的高频子带系数,ε表示高频子带系数的状态转移概率。
7.根据权利要求1所述基于平稳方向波域概率图模型的图像分割方法,其特征在于,步骤(6)所述的最大期望方法,按照如下步骤进行:
第一步,将隐马尔科夫链模型参数集中的四个参数全部初始化为0;
第二步,按照下式,更新隐马尔科夫链模型参数集:
其中,θS+1表示循环次数为S+1的隐马尔科夫链模型参数集,S表示循环次数,argmin(·)表示使函数取最小值的操作,θ表示隐马尔科夫链模型参数集,E[·]表示求期望值操作,ln(·)表示取对数操作,p(H,w|θS)表示待分割图像的训练图像块的高频子带系数的后验概率,H表示待分割图像的训练图像块的高频子带系数,w表示待分割图像的训练图像块的高频子带系数的状态值,θS表示循环次数为S的隐马尔科夫链模型参数集;
第三步,将循环次数为S+1的隐马尔科夫链模型参数集与循环次数为S的隐马尔科夫链模型参数集相减,若差值小于1×10-4,则停止循环,得到最终的隐马尔科夫链模型参数集;否则,返回第二步继续循环。
8.根据权利要求1所述基于平稳方向波域概率图模型的图像分割方法,其特征在于,步骤(7)所述最大似然方法,按照如下步骤进行:
第一步,按照下式,求解高频子带系数的似然值:
其中,F(D|θ)表示高频子带系数的似然值,D表示待分割图像的高频子带系数,θ表示隐马尔科夫链模型参数集,∑表示求和操作,m表示高频子带系数的状态值,β(m)表示高频子带系数的隐状态变量为状态值m时的条件似然概率,f(D)表示高频子带系数的概率分布函数;
第二步,将同一尺度下的所有高频子带系数的似然值相乘,将相乘结果作为每个尺度的高频子带系数的似然值;
第三步,求解所有尺度高频子带系数的似然值的最大值,得到待分割图像中像素的最大似然值。
9.根据权利要求1所述基于平稳方向波域概率图模型的图像分割方法,其特征在于,步骤(8)所述的最大后验估计公式如下:
其中,Y表示最终的分割图像,argmax(·)表示使函数取最大值的操作,φ表示初始分割图像中像素的类别,P(φ)表示初始分割图像中像素的先验概率,F表示待分割图像中像素的最大似然值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410092303.0A CN103871060B (zh) | 2014-03-13 | 2014-03-13 | 基于平稳方向波域概率图模型的图像分割方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410092303.0A CN103871060B (zh) | 2014-03-13 | 2014-03-13 | 基于平稳方向波域概率图模型的图像分割方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103871060A true CN103871060A (zh) | 2014-06-18 |
CN103871060B CN103871060B (zh) | 2016-08-17 |
Family
ID=50909565
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410092303.0A Active CN103871060B (zh) | 2014-03-13 | 2014-03-13 | 基于平稳方向波域概率图模型的图像分割方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103871060B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106372663A (zh) * | 2016-08-30 | 2017-02-01 | 北京小米移动软件有限公司 | 构建分类模型的方法及装置 |
US9665945B2 (en) | 2014-07-23 | 2017-05-30 | Xiaomi Inc. | Techniques for image segmentation |
CN106780440A (zh) * | 2016-11-29 | 2017-05-31 | 北京邮电大学 | 爆炸电路板残片图像自动比对识别方法 |
CN109102511A (zh) * | 2018-07-06 | 2018-12-28 | 深圳先进技术研究院 | 一种脑血管分割方法、系统及电子设备 |
CN110782407A (zh) * | 2019-10-15 | 2020-02-11 | 北京理工大学 | 一种基于天空区域概率分割的单幅图像去雾方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101329736A (zh) * | 2008-06-20 | 2008-12-24 | 西安电子科技大学 | 基于特征选择和隐马尔科夫模型的图像分割方法 |
CN101447080A (zh) * | 2008-11-19 | 2009-06-03 | 西安电子科技大学 | 基于非下采样Contourlet变换的HMT图像分割方法 |
US20100310172A1 (en) * | 2009-06-03 | 2010-12-09 | Bbn Technologies Corp. | Segmental rescoring in text recognition |
-
2014
- 2014-03-13 CN CN201410092303.0A patent/CN103871060B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101329736A (zh) * | 2008-06-20 | 2008-12-24 | 西安电子科技大学 | 基于特征选择和隐马尔科夫模型的图像分割方法 |
CN101447080A (zh) * | 2008-11-19 | 2009-06-03 | 西安电子科技大学 | 基于非下采样Contourlet变换的HMT图像分割方法 |
US20100310172A1 (en) * | 2009-06-03 | 2010-12-09 | Bbn Technologies Corp. | Segmental rescoring in text recognition |
Non-Patent Citations (2)
Title |
---|
JING BAI等: "Image Segmentation Using Directionlet-domain Hidden Markov Tree Models", 《RADAR 2011》, 31 December 2011 (2011-12-31), pages 1615 - 1618 * |
侯彪 等: "基于小波域隐马尔科夫树模型的图像分割", 《红外与毫米波学报》, vol. 28, no. 2, 30 April 2009 (2009-04-30), pages 156 - 160 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9665945B2 (en) | 2014-07-23 | 2017-05-30 | Xiaomi Inc. | Techniques for image segmentation |
CN106372663A (zh) * | 2016-08-30 | 2017-02-01 | 北京小米移动软件有限公司 | 构建分类模型的方法及装置 |
CN106372663B (zh) * | 2016-08-30 | 2019-09-10 | 北京小米移动软件有限公司 | 构建分类模型的方法及装置 |
CN106780440A (zh) * | 2016-11-29 | 2017-05-31 | 北京邮电大学 | 爆炸电路板残片图像自动比对识别方法 |
CN106780440B (zh) * | 2016-11-29 | 2019-05-31 | 北京邮电大学 | 爆炸电路板残片图像自动比对识别方法 |
CN109102511A (zh) * | 2018-07-06 | 2018-12-28 | 深圳先进技术研究院 | 一种脑血管分割方法、系统及电子设备 |
CN110782407A (zh) * | 2019-10-15 | 2020-02-11 | 北京理工大学 | 一种基于天空区域概率分割的单幅图像去雾方法 |
CN110782407B (zh) * | 2019-10-15 | 2021-10-19 | 北京理工大学 | 一种基于天空区域概率分割的单幅图像去雾方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103871060B (zh) | 2016-08-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108053417B (zh) | 一种基于混合粗分割特征的3D U-Net网络的肺分割装置 | |
Kim et al. | Improved simple linear iterative clustering superpixels | |
CN110781775A (zh) | 一种多尺度特征支持的遥感影像水体信息精确分割方法 | |
Zeng et al. | Image segmentation using spectral clustering of Gaussian mixture models | |
CN103871060A (zh) | 基于平稳方向波域概率图模型的图像分割方法 | |
CN101493935B (zh) | 基于剪切波隐马尔可夫模型的合成孔径雷达图像分割方法 | |
CN103455991A (zh) | 一种多聚焦图像融合方法 | |
CN104217406B (zh) | 一种基于剪切波系数处理的sar图像降噪方法 | |
CN102402685A (zh) | 基于Gabor特征的三马尔可夫场SAR图像分割方法 | |
CN106296675A (zh) | 一种高噪声灰度不均匀图像的分割方法 | |
CN100461217C (zh) | 一种复杂性测度的图像纹理分割方法 | |
CN106570183B (zh) | 一种彩色图像检索和分类方法 | |
CN103810704A (zh) | 基于支持向量机和判别随机场的sar图像变化检测方法 | |
Vard et al. | Fast texture energy based image segmentation using directional Walsh–Hadamard transform and parametric active contour models | |
CN104732552B (zh) | 基于非平稳条件场的sar图像分割方法 | |
CN104392241A (zh) | 一种基于混合回归的头部姿态估计方法 | |
Doulgeris | A simple and extendable segmentation method for multi-polarisation SAR images | |
CN102496142B (zh) | 基于模糊的三马尔可夫场sar图像分割方法 | |
CN103700089A (zh) | 一种三维医学图像多尺度异构特征的提取与分类方法 | |
CN103761727A (zh) | 基于自适应高斯混合模型的鲁棒图像分割方法 | |
CN108154176B (zh) | 一种针对单一深度图像的3d人体姿态估计算法 | |
CN103578107A (zh) | 一种交互式图像分割方法 | |
CN104392458A (zh) | 基于空间限制邻域混合模型的图像分割方法 | |
Hosseini-Asl et al. | Lung segmentation based on nonnegative matrix factorization | |
CN103413332B (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 |