CN107274363B - 一种具有尺度敏感特性的边缘保持图像滤波方法 - Google Patents
一种具有尺度敏感特性的边缘保持图像滤波方法 Download PDFInfo
- Publication number
- CN107274363B CN107274363B CN201710408799.1A CN201710408799A CN107274363B CN 107274363 B CN107274363 B CN 107274363B CN 201710408799 A CN201710408799 A CN 201710408799A CN 107274363 B CN107274363 B CN 107274363B
- Authority
- CN
- China
- Prior art keywords
- image
- gradient
- pixel
- point
- filtering
- 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.)
- Active
Links
- 238000001914 filtration Methods 0.000 title claims abstract description 68
- 238000000034 method Methods 0.000 title claims abstract description 45
- 230000035945 sensitivity Effects 0.000 title claims abstract description 37
- 238000005457 optimization Methods 0.000 claims abstract description 24
- 239000011159 matrix material Substances 0.000 claims description 21
- 238000009499 grossing Methods 0.000 claims description 19
- 238000004364 calculation method Methods 0.000 claims description 10
- 239000013598 vector Substances 0.000 claims description 10
- 238000012937 correction Methods 0.000 claims description 9
- 238000005259 measurement Methods 0.000 claims description 7
- 239000000126 substance Substances 0.000 claims description 2
- 230000000694 effects Effects 0.000 description 7
- 238000012545 processing Methods 0.000 description 7
- 125000001475 halogen functional group Chemical group 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 4
- 238000000354 decomposition reaction Methods 0.000 description 3
- 230000006835 compression Effects 0.000 description 2
- 238000007906 compression Methods 0.000 description 2
- 238000009792 diffusion process Methods 0.000 description 2
- 238000004321 preservation Methods 0.000 description 2
- 238000000926 separation method Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000002146 bilateral effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 150000001875 compounds Chemical class 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000004927 fusion Effects 0.000 description 1
- 238000012804 iterative process Methods 0.000 description 1
- 238000012417 linear regression Methods 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 230000014759 maintenance of location Effects 0.000 description 1
- 239000002184 metal Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000035772 mutation Effects 0.000 description 1
- 238000003786 synthesis reaction Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20024—Filtering details
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20172—Image enhancement details
- G06T2207/20192—Edge enhancement; Edge preservation
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Image Processing (AREA)
Abstract
本发明提供了一种具有尺度敏感特性的边缘保持图像滤波方法,不仅能够通过全局优化很好地保持边缘,还能够通过在优化方程中引入尺度敏感性度量其具有良好的尺度敏感特性。本发明采用迭代优化的方式计算输入图像I的滤波结果图像S。在每一次优化过程中,通过最小化当前的能量函数求解出该次优化得到的滤波图像St,最后一次优化得到的滤波图像即为最后的滤波结果图像S。
Description
技术领域
本发明属于数字图像非线性滤波领域,具体涉及一种具有尺度敏感特性的边缘保持图像滤波方法。
背景技术
数字图像的边缘保持滤波属于一种非线性滤波,目的是在滤除图像细节信息的同时使边缘信息不被模糊,即保持图像里物体的边缘。由于这种特性,边缘保持滤波在图像和计算机视觉处理中被广泛应用,能够使数字图像和视频的增强、合成、高动态压缩和实现各种风格化等处理中获得良好的效果,并避免出现光晕现象。目前,已经提出和发明了很多边缘保持滤波方法和技术。其中,双边滤波是一种使用最为广泛的图像边缘保持滤波技术。它利用图像空间域和亮度域内的两个高斯滤波核的乘积对图像进行卷积运算,能够平滑和模糊掉各像素邻域内小的细节特征,而保持由大的像素值突变而产生的边缘特征。各向异性扩散也是一种经典的边缘保持滤波器,它模仿热量的物理扩散过程对图像细节进行模糊,并通过在此过程的计算中引入一个边缘停止函数,使得遇边缘而终止,从而起到保护边缘的作用。
为提高滤波器计算速度和进一步改善边缘保持效果,近些年来又出现了一些新的滤波器。例如,引导图像滤波以一幅引导图像和原图像构造局部线性回归模型,通过模型的最优值求解可得到图像滤波结果。由于该模型求解的形式简单、易于计算,这种滤波的计算速度较快,且该方法可以有效避免梯度反转等现象。基于域转换的DTF滤波器通过将图像的二维滤波转换到一维空间上进行处理,能够使滤波的计算速度更快,同时又能取得很好的边缘保持效果。上述滤波方法的计算过程都是基于图像局部信息的操作,这样做的缺点是会在边缘处造成一定的模糊和引入一定程度的光晕。近几年,又出现了基于全局优化的边缘保持滤波技术,通过最小化一个全局的能量函数,使滤波后的结果中模糊效果尽量全局扩散,而不会积聚在边缘附近。
目前,绝大多数的边缘保持滤波算法都是首先滤除低对比度图像信息,而保持较为强烈的边缘信息。而很多情况下希望边缘保持滤波器能够具有尺度敏感特性,也就是图像中尺寸较小(小尺度)的信息首先被滤除,同时保持大尺度信息及其边缘不被模糊。实际上,尺度敏感特性是传统图像滤波器的基本特性之一。例如,最常用的高斯滤波就是根据图像中物体尺寸的大小,优先滤除小尺度信息,
其优点是能对图像中的信息进行尺度分离,从而能很好地在图像信息的多尺度分解、分析及处理中发挥作用。然而,高斯滤波同时也会模糊大尺度信息的边缘,引起光晕,以及在尺度分离时不能很好地处理物体边缘。传统的边缘保持滤波虽然能够保持边缘,但很难具备上述尺度敏感的基本滤波特性;而要使滤波器具有良好的尺度敏感特性,又难以获得理想的边缘保持效果。
发明内容
有鉴于此,本发明提供了一种具有尺度敏感特性的边缘保持图像滤波方法,不仅能够通过全局优化很好地保持边缘,还能够通过在优化方程中引入尺度敏感性度量其具有良好的尺度敏感特性。
一种具有尺度敏感特性的边缘保持图像滤波方法,该方法的具体步骤包括:
通过迭代的方式计算输入图像I的滤波结果图像S;第t次迭代优化中,求解满足能量函数最小的图像,获得此次迭代优化得到的滤波图像St;
所述能量函数的值为所有像素的能量之和;像素p的能量为像素差异的平方与平滑梯度之和;其中,所述像素差异为滤波图像St与输入图像I在像素p处的差异,所述平滑梯度为滤波图像St中像素p的梯度乘以像素p的平滑参数
优选地,所述能量函数为:
优选地,p点处r尺度敏感性度量值Qr,p
第三步:在图像F中以每一像素点p为中心、半径为ri的圆周内,找出过p点且方向为θp的直线段;找到所述直线段上的所有像素点,以从线段一端到另一端的顺序,将所有像素点记为:{p-n,…,p-1,p0,p1,…,pn},n为正整数;
其中,ε为设定的小值;
优选地,所述求解满足能量函数最小的图像的具体过程为:
将所述能量函数写成矩阵表达式的形式:
式中,VS表示由滤波图像St中各元素组成的列向量;VI表示由输入图像I中各元素组成的列向量;Dx和Dy分别为x和y方向上一阶梯度算子和的矩阵形式;上标T表示对应矩阵的转置;W为对角线矩阵,且对角线元素值由所有像素点处的平滑参数ωt-1给定;
对能量函数的矩阵表达式求导,并令导数为零,获得:
(o+λL)VS=VI (2)
通过解算式(2),获得列向量VS,即得到滤波图像St。
优选地,当迭代次数达到设定的总次数N时,停止迭代,得到最优值。
优选地,迭代优化的总次数N=3。
优选地,全局平滑参数λ=0.1。
有益效果:
1)传统的数字图像滤波方法难以同时具有良好的边缘保持和尺度敏感特性。本发明提供了一种具有尺度敏感特性的边缘保持滤波技术,在保持图像边缘特征的同时,又能够使得滤波过程具有尺度敏感特性,使尺度较小的图像信息首先被滤除。基于这种滤波技术,可以实现边缘保持的多尺度图像信息分解,为各种数字图像的多尺度处理提供一种更加精准、不会引入光晕的滤波和分解手段。
2)本发明可以作为一种通用的底层图像滤波和分解工具,取代传统的边缘保持滤波方法,广泛应用于现有的图像和视频增强、高动态压缩、融合、去噪和风格化等处理任务中,得到视觉效果更好或真实感更强的处理结果。
附图说明
图1为滤波计算流程图。
图2为不同大小方块的图像的不同尺度滤波结果。
具体实施方式
下面结合附图并举实施例,对本发明进行详细描述。
本发明提供了一种具有尺度敏感特性的边缘保持图像滤波方法,其采用迭代优化的方式计算输入图像I的滤波结果图像S;在每一次优化过程中,通过最小化当前的能量函数求解出该次迭代优化得到的滤波图像St,最后一次优化得到的滤波图像即为最后的滤波结果图像S,迭代流程如图1所示。
其中,在上述迭代过程中,能量函数采用图像中所有像素的能量之和表示。任意像素p的能量为像素差异的平方与平滑梯度之和;其中,所述像素差异为滤波图像St与输入图像I在像素p处的差异,所述平滑梯度为滤波图像St中像素p的梯度乘以像素p的平滑参数
可见,本发明采用多尺度敏感性度量值构建平滑参数,在保持图像边缘特征的同时,又能够使得滤波过程具有尺度敏感特性,使尺度较小的图像信息首先被滤除。基于这种滤波技术,可以实现边缘保持的多尺度图像信息分解,为各种数字图像的多尺度处理提供一种更加精准、不会引入光晕的滤波和分解手段。
下面对本发明具有尺度敏感特性的边缘保持图像滤波方法的实施过程进行详细描述,其具体步骤如下:
步骤1:设定迭代优化的总次数N、全局平滑参数λ和尺度参数r;第t次优化的能量函数表达式为:
式中,Ω为图像中所有像素点集合,下角标p表示相应图像中的像素点p;St为第t次迭代优化所求解的滤波图像,I为输入图像;分别为St在x和y方向上的梯度图像;λ为预先设定的全局平滑参数;为t-1次迭代后在p点处的平滑参数;t为迭代优化次数,t的初始值为1。
步骤3:判断优化次数t是否达到预设的总次数N;如果t=N,则t=N对应获得的滤波图像SN为最终的滤波结果图像S;如果t<N,则t自加1,重复步骤2-3。
下面对上述步骤2中的技术细节问题进行详细说明:
式中,为t-1次迭代后获得的p点处梯度强度;为t-1次迭代后得到的p点处梯度强度校正量;Qr,p为p点处r尺度敏感性度量值;ε为一个非常小的常值,用以避免上式中分母为零。下面详细说明上述变量的计算方法:
式中,Gp表示输入图像I在p点处梯度强度,其中表示输入图像I在x和y方向上的梯度强度图像;表示第t-1次迭代得到的滤波图像St-1在p点处梯度强度,其中 为第t-1次迭代得到的滤波图像St-1在x和y方向上的梯度图像,且令S0=I。
(3)在图像F中以每一像素点p为中心、半径为ri的圆周内,找出过p点且方向为θp的直线段。然后,找到该线段上的所有像素点,以从线段一端到另一端的顺序,将这些像素点记为:{p-n,…,p-1,p0,p1,…,pn},n为正整数,并以N(p)表示该像素点集合。在该集合中,中心点p0(也就是p)两侧的像素点个数相等。
式中,i∈[1,m],m取正整数;0<r1<r2<…<rm且rm=r;min(·)表示取最小值运算。本实施例中取r1=1,r2=2,…,rm-1=r-1,rm=r。
二、通过最小化能量函数求解滤波图像的具体方法为:
首先,将能量函数写成矩阵的形式:
式中,VS表示由滤波图像St中各元素组成的列向量;VI表示由输入图像I中各元素组成的列向量;Dx和Dy分别为x和y方向上一阶梯度算子和的矩阵形式;上标T表示对应矩阵的转置;W为对角线矩阵,且对角线元素值由所有像素点处的平滑参数ωt-1给定。
然后,对上述矩阵表达式求导,并令导数为零,得到
(o+λL)VS=VI
一般情况下,令全局平滑参数λ=0.1、迭代优化的总次数N=3即可得到效果较好的滤波效果。图2给出了利用上述滤波步骤,对包含不同大小方块的图像进行的在不同尺度参数r条件下滤波得到的结果。从图中可以看到,较小尺寸的方块首先被滤除,且随着尺度r的增大,较大尺寸的方块也相继被滤除,表现出了很好尺度敏感特性。同时,每一滤波结果中剩下来的方块及其边缘很好地保持了下来,没有被模糊,说明该滤波技术具有很好的边缘保持特性。图2给出了两幅自然图像的不同尺度滤波结果,同样可以看到,小尺度的细节信息逐渐被抑制,而剩下来的图像信息的边缘被很好地保持了。
综上所述,以上仅为本发明的较佳实施例而已,并非用于限定本发明的保护范围。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (8)
1.一种具有尺度敏感特性的边缘保持图像滤波方法,其特征在于,该方法的具体步骤包括:
通过迭代的方式计算输入图像I的滤波结果图像S;第t次迭代优化中,求解满足能量函数最小的图像,获得此次迭代优化得到的滤波图像St;
所述能量函数的值为所有像素的能量之和;像素p的能量为像素差异的平方与平滑梯度之和;其中,所述像素差异为滤波图像St与输入图像I在像素p处的差异,所述平滑梯度为滤波图像St中像素p的梯度乘以像素p的平滑参数
第三步:在图像F中以每一像素点p为中心、半径为ri的圆周内,找出过p点且方向为θp的直线段;找到所述直线段上的所有像素点,以从线段一端到另一端的顺序,将所有像素点记为:{p-n,…,p-1,p0,p1,…,pn},n为正整数;
其中,ε为设定的小值;
6.如权利要求1所述的方法,其特征在于,当迭代次数达到设定的总次数N时,停止迭代,得到最优值。
7.如权利要求6所述的方法,其特征在于,迭代优化的总次数N=3。
8.如权利要求2所述的方法,其特征在于,全局平滑参数λ=0.1。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710408799.1A CN107274363B (zh) | 2017-06-02 | 2017-06-02 | 一种具有尺度敏感特性的边缘保持图像滤波方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710408799.1A CN107274363B (zh) | 2017-06-02 | 2017-06-02 | 一种具有尺度敏感特性的边缘保持图像滤波方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107274363A CN107274363A (zh) | 2017-10-20 |
CN107274363B true CN107274363B (zh) | 2020-09-22 |
Family
ID=60065673
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710408799.1A Active CN107274363B (zh) | 2017-06-02 | 2017-06-02 | 一种具有尺度敏感特性的边缘保持图像滤波方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107274363B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110159868B (zh) * | 2019-05-27 | 2021-07-30 | 北京奥蓝仕技术有限公司 | 智能套筒的数据处理方法、装置及智能套筒 |
CN110246099B (zh) * | 2019-06-10 | 2021-09-07 | 浙江传媒学院 | 一种保持结构边缘的图像去纹理方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101908208A (zh) * | 2010-07-27 | 2010-12-08 | 浙江大学 | 面向图像边缘检测的平滑滤波空间尺度的自适应确定方法 |
CN104240212A (zh) * | 2014-09-03 | 2014-12-24 | 西安电子科技大学 | 基于目标特征的isar图像融合方法 |
CN104915929A (zh) * | 2014-03-14 | 2015-09-16 | 电子科技大学 | 一种图像滤波处理方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8897378B2 (en) * | 2013-03-12 | 2014-11-25 | Tandent Vision Science, Inc. | Selective perceptual masking via scale separation in the spatial and temporal domains using intrinsic images for use in data compression |
-
2017
- 2017-06-02 CN CN201710408799.1A patent/CN107274363B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101908208A (zh) * | 2010-07-27 | 2010-12-08 | 浙江大学 | 面向图像边缘检测的平滑滤波空间尺度的自适应确定方法 |
CN104915929A (zh) * | 2014-03-14 | 2015-09-16 | 电子科技大学 | 一种图像滤波处理方法 |
CN104240212A (zh) * | 2014-09-03 | 2014-12-24 | 西安电子科技大学 | 基于目标特征的isar图像融合方法 |
Non-Patent Citations (4)
Title |
---|
Edge-aware Gradient Domain Optimization Framework for Image Filtering by Local Propagation;Miao Hua 等;《IEEE》;20141231;第1-8页 * |
Rolling Guidance Filter;Qi Zhang 等;《ECCV 2014》;20141231;第815-830页 * |
Smooth and iteratively Restore: A simple and fast edge-preserving smoothing model;Philipp Kniefacz 等;《arXiv》;20150525;第1-9页 * |
基于双边与高斯滤波混合分解的图像融合方法;周志强 等;《系统工程与电子技术》;20160131;第38卷(第1期);第8-13页 * |
Also Published As
Publication number | Publication date |
---|---|
CN107274363A (zh) | 2017-10-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Zhang et al. | Learning fully convolutional networks for iterative non-blind deconvolution | |
Pan et al. | Kernel estimation from salient structure for robust motion deblurring | |
Xu et al. | An improved anisotropic diffusion filter with semi-adaptive threshold for edge preservation | |
CN101902547A (zh) | 图像处理方法以及图像处理装置 | |
Pan et al. | Spatially variant linear representation models for joint filtering | |
Yu et al. | Efficient patch-wise non-uniform deblurring for a single image | |
Zhou et al. | Scale-aware edge-preserving image filtering via iterative global optimization | |
CN105427259A (zh) | 多方向加权tv和非局部自相似性正则化图像去模糊方法 | |
CN107274363B (zh) | 一种具有尺度敏感特性的边缘保持图像滤波方法 | |
CN109410246B (zh) | 基于相关滤波的视觉跟踪的方法及装置 | |
CN111882555B (zh) | 基于深度学习的网衣检测方法、装置、设备及存储介质 | |
Song et al. | Structure-preserving texture filtering for adaptive image smoothing | |
Ji et al. | Image recovery via geometrically structured approximation | |
CN112308872A (zh) | 基于多尺度Gabor一阶导数的图像边缘检测方法 | |
Hu et al. | Image deblurring based on enhanced salient edge selection | |
CN109635809B (zh) | 一种面向视觉退化图像的超像素分割方法 | |
Guo et al. | Structure-texture decomposition via joint structure discovery and texture smoothing | |
US20140363090A1 (en) | Methods for Performing Fast Detail-Preserving Image Filtering | |
Zhao et al. | Constant time texture filtering | |
CN116740501A (zh) | 一种图像模糊区域修复补偿模型的训练方法及应用 | |
CN108615224B (zh) | 图像校正设备和方法以及非易失性计算机可读存储介质 | |
Pang et al. | Improved L0 gradient minimization with L1 fidelity for image smoothing | |
Zhang et al. | Non-blind deblurring of structured images with geometric deformation | |
Umnov et al. | Ringing artifact suppression using sparse representation | |
Barbu | Nonlinear fourth-order diffusion-based model for image denoising |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |