基于两通道纹理分割主动轮廓模型的图像分割方法.pdf

(19)中华人民共和国国家知识产权局*CN103413332A*(12)发明专利申请(10)申请公布号(10)申请公布号CNCN103413332103413332A(43)申请公布日2013.11.27(21)申请号201310371336.4(22)申请日2013.08.23(71)申请人华北电力大学地址102206北京市昌平区回龙观朱辛庄2号(72)发明人许刚马爽史巍刘坤(74)专利代理机构北京麟保德和知识产权代理事务所(普通合伙)11428代理人周恺丰(51)Int.Cl.G06T7/40(2006.01)权权利要求书3页利要求书3页说明书11页说明书11页附图3页附图3页(54)发明名称基于两通道纹理分割主动轮廓模型的图像分割方法(57)摘要本发明公开了数字图像处理技术领域中的一种基于两通道纹理分割主动轮廓模型的图像分割方法。包括，提取图像中每个像素的灰度值、水平梯度场和垂直梯度场；计算图像中每个像素的灰度值、水平梯度场和垂直梯度场对应的纹理特征；根据所述纹理特征获取灰度特征通道和边缘特征通道；建立两通道纹理分割主动轮廓线模型；通过水平集函数的演化最小化纹理分割模型完成图像分割。本发明提高了算法效率并避免由灰度信息引起的误分割，提高了算法的准确性。CN103413332ACN10342ACN103413332A权利要求书1/3页1.一种基于两通道纹理分割主动轮廓模型的图像分割方法，其特征是所述方法包括：步骤1：提取图像中每个像素的灰度值、水平梯度场和垂直梯度场；步骤2：计算图像中每个像素的灰度值、水平梯度场和垂直梯度场对应的纹理特征；步骤3：根据所述纹理特征获取灰度特征通道和边缘特征通道；步骤4：建立两通道纹理分割主动轮廓线模型；步骤5：通过水平集函数的演化最小化纹理分割模型完成图像分割。2.根据权利要求1所述的方法，其特征是所述提取图像中每个像素的水平梯度场具体为采用公式计算图像中的第i行第j列像素的水平梯度场其中，I(i,j)是图像中的第i行第j列像素的灰度值，I(i+1,j)是图像中的第i+1行第j列像素的灰度值。3.根据权利要求1所述的方法，其特征是所述提取图像中每个像素的垂直梯度场具体为采用公式计算图像中的第i行第j列像素的水平梯度场其中，I(i,j)是图像中的第i行第j列像素的灰度值，I(i,j+1)是图像中的第i行第j+1列像素的灰度值。4.根据权利要求2或3所述的方法，其特征是所述步骤3具体为：步骤301：根据图像中每个像素的水平梯度场和垂直梯度场对应的纹理特征u2(x,y)和u3(x,y)，采用公式提取边缘特征uedge；其中，I是像素的灰度值，是像素的灰度值的梯度；步骤302：根据公式计算灰度特征通道u′1(x,y)和边缘特征通道u′2(x,y)；其中，i=1,2，L1=u1(x,y)，L2=uedge(x,y)。5.根据权利要求4所述的方法，其特征是所述建立两通道纹理分割主动轮廓线模型为：其中，和分别是灰度特征通道和边缘特征通道中的曲线C内部和外部的均值；曲线C满足C={(x,y):φ(x,y)=0}，φ(x,y)为水平函数集；μ、α和β分别为长度项、灰度项和边缘项的参数；Ω为积分区域，即图像区域；δ(·)为Dirac函数；为水平函数集φ(x,y)的梯度；2CN103413332A权利要求书2/3页和分别为灰度特征通道中的曲线C内外区域的灰度均值；和分别为灰度特征通道的参数，且u′1(x,y)为像素(x,y)对应的灰度特征通道的值；H(·)为Heaviside函数；a为用于调整函数形状的常数且a>0；和分别为边缘特征通道中的曲线C内外区域的灰度均值；和分别为边缘特征通道的参数，且u′2(x,y)为像素(x,y)对应的边缘特征通道的值；H(·)为Heaviside函数。6.根据权利要求5所述的方法，其特征是所述步骤5包括：步骤501：随机给定初始闭合分割曲线C0，并计算初始闭合分割曲线C0对应的初始水平集函数φ0(x,y)；步骤502：设定模型参数μ,α,β,步骤503：令k=0，分别计算初始闭合分割曲线C0内部和外部的均值；初始闭合分割曲线C0内部的均值的计算公式为：初始闭合分割曲线C0外部的均值的计算公式为：上述两个公式中，i=1,2；Ω为积分区域，即图像区域；u′1(x,y)为像素(x,y)对应的灰度特征通道的值；u′2(x,y)为像素(x,y)对应的边缘特征通道的值；H(·)为Heaviside函数；步骤504：根据公式φk+1(x,y)-φk(x,y)=Δt×L(φk(x,y))迭代计算φk+1(x,y)；3CN103413332A权利要求书3/3页即Δt为设定时间步长，δε(·)为Dirac函数；步骤505：从水平函数集φk+1(x,y)中提取零水平集，该提取的零水平集即演化曲线；步骤