基于剪切波系数概率密度分布的图像边缘检测方法
【专利摘要】本发明公开一种基于剪切波系数概率密度分布的图像边缘检测方法，首先抑制噪声、增强图像的边缘，然后统计图像剪切波变换后的概率密度分布规律，进而利用高斯混合模型对剪切波系数所对应的大、小状态建模，分离出图像的候选边缘像素，最后采用大津法和细化方法从候选像素中提取图像的边缘，实验结果表明本发明方法提取出的边缘连续性好，人眼主观质量较高。
【专利说明】
基于剪切波系数概率密度分布的图像边缘检测方法
技术领域
[0001] 本发明设及图像处理领域，尤其是一种边缘连续性好、人眼视觉质量高的基于剪 切波系数概率密度分布的图像边缘检测方法。
【背景技术】
[0002] 边缘是图像奇异点、线的集合，反映了图像灰度的不连续性，是图像的重要特性之 一。边缘检测结果的好坏将影响目标分割、特征提取和图像识别等高层处理的质量，是图像 分析与识别的关键预处理步骤之一。
[0003] 虽然早期的化nny算子W及Robe;rts算子、Sobel算子、Prewitt算子、Lap lace算子 等均能实现图像的边缘检测，但是检测出的边界往往存在不连续的地方，且易受到噪声的 干扰。小波分析理论的不断完善，使得它在图像的边缘检测中得W成功应用，研究表明小波 变换的模极大值点可W定位图像的奇异信号，若某一点的模值比其同幅角方向上相邻像素 大时，可确定此点为局部极大值点并保留，从而检测出图像的边缘信息。同时，众多研究者 还对图像边缘多尺度分析的数学理论框架进行了深入探讨，主要是根据LipscMtz正则性 判断小波模极大值的衰减速度，进而得知边缘的奇异性，运为基于小波分析的边缘检测奠 定了严密的数学基础。但是，图像的小波系数模极大值点还有可能是噪声或纹理信息等，由 于小波变换只能捕获水平、竖直和对角线方向的线奇异，不具有较为完备的方向分析能力， 无法进一步有效辨别噪声和纹理信息对应的模极大值点，运将不可避免地造成边缘丢失或 受到噪声干扰。而理论表明，剪切波(Shearlet)变换能够更加准确地检测图像边缘的几何 信息和方向特征，是目前最优的图像稀疏表示方法之一。Kan曲ui Guo等人建立了剪切波变 换用于边缘检测的严密的数学理论；Glenn Easley等人设计了利用剪切波变换系数的 Lipschitz正则性从有噪声的图像中区分出真正边缘点的方法；侯彪等人利用改进的 化ear let变换并结合模糊C均值、Snake模型进行边缘检测的新方法，克服了边界追踪法、 Markovian分割法、Snake算法等传统检测算法存在的计算速度慢、易受噪声干扰、难处理轮 廓曲线分离或合并的缺点；Pengfei Xu等人则结合剪切波和小波变换，借助不同频率的信 号特征更好地从噪声图像中提取出了边缘信息。然而，现有的基于剪切波变换的边缘检测 算法均没有很好地考虑变换后的系数分布特征。实际上，噪声与图像边缘的变换系数概率 密度分布特征迴异，现有方法却未在边缘检测过程中引入概率密度分布的约束，在一定程 度上影响了边缘检测效率。
【发明内容】
[0004] 本发明是为了解决现有技术所存在的上述技术问题，提供一种可改善边缘连续 性、提高视觉效果的基于剪切波系数概率密度分布的图像边缘检测方法。
[0005] 本发明的技术解决方案是:一种基于剪切波系数概率密度分布的图像边缘检测方 法，其特征在于按照如下步骤进行： 步骤1.利用蒙特卡罗方法估算输入图像X的噪声方差巧2,再采用方差为。的高斯平 滑滤波器对图像i进行平滑处理； 步骤2.采用Kirsch算子对输入图像1进行锐化预处理，设处理后的图像为?'，并且
,其中減基哲(6,???載黨巧:獲^ 模板，表示卷积运算;对f进行剪切波变换； 步骤3.采用坎尼算子提取最低频子带的边缘系数，并保持运些边缘系数不变，将最低 频子带中的非边缘系数全部置为0; 步骤4.用高斯混合模型对变换后的剪切波高频子带系数集合采=:麵鴻,搞|进 行建模，根据大、小状态获取图像的边缘系数，所述N表示高频子带系数总数： 步骤4.1采用高斯混合分布函数对剪切波高频子带系数进行建模：
其中，游二策为剪切波高频子带系数取大状态的状态变量，緻二::;2为剪切波高频子带 系数取小状态的状态变量,3?分别表示高频子带系数取大、小两个状态的先验概率且满足
吸与I;，霉是待估计的高斯混合分布参数，茶滅館:!端，巧3 为高频子带系数取大、小两个状态所对应的高斯分布的概率密度函数，>媒为高斯分布的均 值，巧^为高斯分布的方差； 步骤4.2建立似然函数
并利用期望最大化极大似然估计 算法训练参数撲二谭装进而判定每个剪切波高频子带系数所处的状态：
步骤4.3保持那些处于大状态的概率超过ο. 5的高频子带系数不变，而将其余高频子 带系数全部置成0; 步骤4.4利用邻域分析法对大状态高频子带系数进行细化处理，将不在细化边缘上的 大状态高频子带系数也置成0; 步骤5.对处理后的变换系数进行逆剪切波变换，得到一幅新的图像1% 步骤6.利用大津法计算r的自适应阔值了，采用该阔值将图像r进行二值化； 步骤7.利用邻域分析法对二值化后的图像Γ进行细化处理，输出单像素宽的边缘检 测结果。
[0006] 与现有技术相比，本发明具有W下优点：首先，利用剪切波变换对悼连续的线奇 异信号的优秀捕获能力及其变换系数的概率密度分布，实现了更加准确的边缘系数定位； 其次，通过高斯平滑、在低频子带提取候选边缘和两次细化处理，能够尽可能地排除图像中 的噪声干扰，提取出的边缘具有更好的连续性，整体视觉效果好。
【附图说明】
[0007] 图1是本发明实施例与同类方法的边缘检测结果对比。
【具体实施方式】
[000引一种供基于剪切波系数概率密度分布的图像边缘检测方法，其特征在于按照如下 步骤进行： 步骤1.利用蒙特卡罗(Monte-Carlo)方法估算输入图像I的噪声方差再采用方 差为勺高斯平滑滤波器对图像I进行平滑处理； 步骤2.采用Kirsch算子对输入图像1进行锐化预处理，设处理后的图像为當，并且
，其中蛛；(裏爲蘇卑畢夺，:;制 模板，表示卷积运算;对f进行剪切波变换； 步骤3.采用坎尼(化nny)算子提取最低频子带的边缘系数，并保持运些边缘系数不变， 将最低频子带中的非边缘系数全部置为0; 步骤4.用高斯混合模型对变换后的剪切波高频子带系数集合苯二鶴勒'·',％!进 行建模，根据大、小状态获取图像的边缘系数，所述Ν表示高频子带系数总数： 步骤4.1采用高斯混合分布函数对剪切波高频子带系数进行建模：
其中，嫌二I为剪切波高频子带系数取大状态的状态变量，游二为剪切波高频子带 系数取小状态的状态变量，馬分别表示高频子带系数取大、小两个状态的先验概率且满足
为高频子带系数取大、小两个状态所对应的高斯分布的概率密度函数，&为高斯分布的均 值，巧为高斯分布的方差； 步骤4.2建立似然函蠻
并利用期望最大化化xpectation Maximization, EM)极大似然估计算法训练参数梦二避袋猫1扱去揉I，进而判定每个剪 切波高频子带系数所处的状态：
步骤4.3保持那些处于大状态的概率超过0.5的高频子带系数不变，而将其余高频子 带系数全部置成0; 步骤4.4利用邻域分析法对大状态高频子带系数进行细化处理，将不在细化边缘上的 大状态高频子带系数也置成0; 步骤5.对处理后的变换系数进行逆剪切波变换，得到一幅新的图像Γ ; 步骤6.利用大津法计算f的自适应阔值巫，采用该阔值将图像Γ进行二值化； 步骤7.利用邻域分析法对二值化后的图像r进行细化处理，输出单像素宽的边缘检 测结果，算法结束。
[0009]本发明实施例与同类方法的边缘检测结果对比如图1所示，其中(a)为原图像;（b) 为Canny算子的检测结果；（C)为基于小波的高斯混合模型的检测结果；（d)为基于轮廓波的 高斯混合模型的检测结果；（e)为本发明实施例的检测结果。从图1可W看出，本发明方法提 取出的边缘连续性好，人眼主观质量较高。
【主权项】
1. 一种基于剪切波系数概率密度分布的图像边缘检测方法，其特征在于按照如下步骤 进行： 步骤1.利用蒙特卡罗方法估算输入图像1的噪声方差再采用方差为α2的高斯平 滑滤波器对图像1进行平滑处理； 步骤2.采用Kirsch算子对输入图像丨进行锐化预处理，设处理后的图像为Γ，并且表示Kir sch算子的8个 模板，表示卷积运算;对f进行剪切波变换； 步骤3.采用坎尼算子提取最低频子带的边缘系数，并保持这些边缘系数不变，将最低 频子带中的非边缘系数全部置为〇; 步骤4.用高斯混合模型对变换后的剪切波高频子带系数集合iT …,?)进 行建模，根据大、小状态获取图像的边缘系数，所述N表示高频子带系数总数： 步骤4.1采用高斯混合分布函数对剪切波高频子带系数进行建模：其中，=4为剪切波高频子带系数取大状态的状态变量，·为剪切波高频子带系 数取小状态的状态变量，分别表示高频子带系数取大、小两个状态的先验概率且满足是待估计的高斯混合分布参数，彳 为高频子带系数取大、小两个状态所对应的高斯分布的概率密度函数，.?为高斯分布的均 值，ffi为高斯分布的方差；并利用期望最大化极大似然估计 算法训练参数康二ail:通輯康0Γ1瑪，进而判定每个剪切波高频子带系数所处的状态：步骤4.3保持那些处于大状态的概率超过0.5的高频子带系数不变，而将其余高频子 带系数全部置成0; 步骤4.4利用邻域分析法对大状态高频子带系数进行细化处理，将不在细化边缘上的 大状态高频子带系数也置成〇; 步骤5.对处理后的变换系数进行逆剪切波变换，得到一幅新的图像?"; 步骤6.利用大津法计算γ的自适应阈值，，采用该阈值将图像r进行二值化； 步骤7.利用邻域分析法对二值化后的图像进行细化处理，输出单像素宽的边缘检 测结果。
【文档编号】G06T7/00GK105825514SQ201610164812
【公开日】2016年8月3日
【申请日】2016年3月22日
【发明人】王相海, 宋传鸣, 苏欣, 沈雨桐
【申请人】辽宁师范大学