发明名称 一种基于模糊增强和曲面拟合的图像边缘特征提取方法
摘要 本发明公开了一种基于模糊增强和曲面拟合的图像边缘特征提取方法,主要为了克服现有方法泛化性能差及参数难以确定等不足。其步骤为:首先,对原始待处理图像进行模糊增强处理;然后,采用最小二乘支持向量回归函数对模糊增强后的图像进行曲面拟合;最后,通过求拟合图像的一阶和二阶导数提取图像边缘,得到最终的边缘特征图像。采用本发明方法能够检测清晰、细致的图像边缘信息,并且方法适应性强,适用于对不同类型的图像进行边缘特征提取。
申请公布号 CN101976336A 申请公布日期 2011.02.16
申请号 CN201010517199.7 申请日期 2010.10.21
申请人 西北工业大学 发明人 郭雷;杜亚勤
分类号 G06K9/46(2006.01)I;G06T5/00(2006.01)I 主分类号 G06K9/46(2006.01)I
代理机构 西北工业大学专利中心 61204 代理人 王鲜凯
主权项 1.一种基于模糊增强和曲面拟合的图像边缘特征提取方法,其特征在于步骤如下:步骤一:利用模糊隶属函数<img file="FSA00000315334200011.GIF" wi="581" he="266" />将待处理图像转换到模糊特征平面,得到模糊隶属度图像;其中,x<sub>mn</sub>是待处理图像m行n列像素点的灰度值,u<sub>mn</sub>为模糊隶属度图像中m行n列像素点的值,L为待处理图像的最大灰度级,x<sub>T</sub>为灰度级阈值;所述的灰度级阈值x<sub>T</sub>利用最大类间方差法计算得到;所述的m的取值范围为[1,M],所述的n的取值范围为[1,N],M为待处理图像的总行数,N为待处理图像的总列数;步骤二:利用增强函数<img file="FSA00000315334200012.GIF" wi="735" he="145" />对模糊隶属度图像进行增强处理,得到增强后的模糊隶属度图像;其中,u′<sub>mn</sub>为增强后的模糊隶属度图像中m行n列像素点的值;步骤三:利用<img file="FSA00000315334200013.GIF" wi="805" he="139" />计算得到模糊增强后的灰度图像;其中,x′<sub>mn</sub>是模糊增强后灰度图像中m行n列像素点的灰度值;步骤四:利用最小二乘支持向量回归函数<img file="FSA00000315334200014.GIF" wi="992" he="123" />对模糊增强后的灰度图像进行曲面拟合,得到拟合图像;其中,m表示图像的第m行,n表示图像的第n列;(m<sub>i</sub>,n<sub>i</sub>)为观测像素坐标,f(m,n)为拟合图像m行n列像素点的灰度值;a为按a=Ω<sup>-1</sup>(Y-bI)计算的高维空间加权矢量;b为按<img file="FSA00000315334200015.GIF" wi="241" he="119" />计算的偏移量;σ<sup>2</sup>是高斯核函数的参数,取值范围为[1,10];Y=[x′<sub>11</sub>,x′<sub>12</sub>,…,x′<sub>MN</sub>]<sup>T</sup>为由模糊增强后图像中各行像素首尾相接组成的向量;Ω=exp{-(|m-m<sub>i</sub>|<sup>2</sup>+|n-n<sub>i</sub>|<sup>2</sup>)/σ<sup>2</sup>}+γ<sup>-1</sup>I,γ为惩罚因子,取值范围为[0,1];I=[1,1,…,1]<sup>T</sup>为单位向量;exp{·}表示指数运算;步骤五:分别按<img file="FSA00000315334200016.GIF" wi="1071" he="123" />和<img file="FSA00000315334200017.GIF" wi="1080" he="122" />计算拟合图像水平方向和垂直方向的梯度,如果拟合图像中某像素点的水平方向梯度值<img file="FSA00000315334200018.GIF" wi="71" he="110" />或垂直方向梯度值<img file="FSA00000315334200019.GIF" wi="59" he="112" />大于阈值一T<sub>1</sub>,则令边缘图像一中对应位置像素点的值为255,否则为0,得到边缘图像一;所述的阈值一T<sub>1</sub>的取值范围为[0.6,0.9];按<img file="FSA00000315334200021.GIF" wi="494" he="118" />计算拟合图像的二阶导数,如果拟合图像中某像素点的二阶导数值<img file="FSA00000315334200022.GIF" wi="210" he="57" />为负且与其直接相邻的8个像素点的二阶导数值中最大值为正且该最大值与该像素点的二阶导数值之差大于阈值二T<sub>2</sub>,则令边缘图像二中对应位置像素点的值为255,否则为0,得到边缘图像二;所述的阈值二T<sub>2</sub>的取值范围为[0.2,0.6];其中,<maths num="0001"><![CDATA[<math><mrow><mfrac><mrow><msup><mo>&PartialD;</mo><mn>2</mn></msup><mi>f</mi></mrow><mrow><mo>&PartialD;</mo><msup><mi>m</mi><mn>2</mn></msup></mrow></mfrac><mo>=</mo><mo>-</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>L</mi></munderover><mfrac><msub><mrow><mn>2</mn><mi>a</mi></mrow><mi>i</mi></msub><msup><mi>&sigma;</mi><mn>2</mn></msup></mfrac><mo>[</mo><mn>1</mn><mo>-</mo><mfrac><mn>2</mn><msup><mi>&sigma;</mi><mn>2</mn></msup></mfrac><msup><mrow><mo>(</mo><mi>m</mi><mo>-</mo><msub><mi>m</mi><mi>i</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>]</mo><mi>exp</mi><mo>{</mo><mo>-</mo><mrow><mo>(</mo><msup><mrow><mo>|</mo><mi>m</mi><mo>-</mo><msub><mi>m</mi><mi>i</mi></msub><mo>|</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>|</mo><mi>n</mi><mo>-</mo><msub><mi>n</mi><mi>i</mi></msub><mo>|</mo></mrow><mn>2</mn></msup><mo>)</mo></mrow><mo>/</mo><msup><mi>&sigma;</mi><mn>2</mn></msup><mo>}</mo><mo>,</mo></mrow></math>]]></maths><maths num="0002"><![CDATA[<math><mrow><mfrac><mrow><msup><mo>&PartialD;</mo><mn>2</mn></msup><mi>f</mi></mrow><msup><mrow><mo>&PartialD;</mo><mi>n</mi></mrow><mn>2</mn></msup></mfrac><mo>=</mo><mo>-</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>L</mi></munderover><mfrac><msub><mrow><mn>2</mn><mi>a</mi></mrow><mi>i</mi></msub><msup><mi>&sigma;</mi><mn>2</mn></msup></mfrac><mo>[</mo><mn>1</mn><mo>-</mo><mfrac><mn>2</mn><msup><mi>&sigma;</mi><mn>2</mn></msup></mfrac><msup><mrow><mo>(</mo><mi>n</mi><mo>-</mo><msub><mi>n</mi><mi>i</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>]</mo><mi>exp</mi><mo>{</mo><mo>-</mo><mrow><mo>(</mo><msup><mrow><mo>|</mo><mi>m</mi><mo>-</mo><msub><mi>m</mi><mi>i</mi></msub><mo>|</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>|</mo><mi>n</mi><mo>-</mo><msub><mi>n</mi><mi>i</mi></msub><mo>|</mo></mrow><mn>2</mn></msup><mo>)</mo></mrow><mo>/</mo><msup><mi>&sigma;</mi><mn>2</mn></msup><mo>}</mo><mo>;</mo></mrow></math>]]></maths>步骤六:对边缘图像一和边缘图像二进行与运算,得到最终的边缘特征图像。
地址 710072 陕西省西安市友谊西路127号
您可能感兴趣的专利