CN110246099A - 一种保持结构边缘的图像去纹理方法 - Google Patents
一种保持结构边缘的图像去纹理方法 Download PDFInfo
- Publication number
- CN110246099A CN110246099A CN201910497118.2A CN201910497118A CN110246099A CN 110246099 A CN110246099 A CN 110246099A CN 201910497118 A CN201910497118 A CN 201910497118A CN 110246099 A CN110246099 A CN 110246099A
- Authority
- CN
- China
- Prior art keywords
- image
- filter
- filtering
- texture
- joint
- 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
- 238000000034 method Methods 0.000 title claims abstract description 71
- 238000001914 filtration Methods 0.000 claims abstract description 83
- 230000002146 bilateral effect Effects 0.000 claims abstract description 51
- 238000002156 mixing Methods 0.000 claims abstract description 10
- 230000002708 enhancing effect Effects 0.000 claims abstract description 8
- 230000000694 effects Effects 0.000 claims abstract description 6
- 238000012545 processing Methods 0.000 claims abstract description 6
- 230000003044 adaptive effect Effects 0.000 claims description 14
- 230000008859 change Effects 0.000 claims description 10
- 230000009466 transformation Effects 0.000 claims description 10
- 230000008569 process Effects 0.000 claims description 5
- 239000004576 sand Substances 0.000 claims description 5
- 238000004364 calculation method Methods 0.000 claims description 3
- 238000009795 derivation Methods 0.000 claims description 3
- 238000013507 mapping Methods 0.000 claims description 3
- 238000005070 sampling Methods 0.000 claims description 3
- 230000001131 transforming effect Effects 0.000 claims description 3
- 238000005728 strengthening Methods 0.000 claims 1
- 230000006870 function Effects 0.000 description 25
- 230000007547 defect Effects 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 125000001475 halogen functional group Chemical group 0.000 description 1
- 238000003709 image segmentation Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 230000008447 perception Effects 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 230000000717 retained effect Effects 0.000 description 1
- 238000000844 transformation Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/20—Image enhancement or restoration using local operators
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/73—Deblurring; Sharpening
-
- 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
- G06T2207/20028—Bilateral filtering
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Image Processing (AREA)
Abstract
本发明提供了一种保持结构边缘的图像去纹理方法,对给定输入图像,计算自适应空间滤波核尺度、增强处理,增强图像中的结构边缘、进行自适应空间滤波核尺度的联合双边滤波,获得引导图像;根据引导图像,局部拉普拉斯图像滤波器,构造联合局部拉普拉斯滤波器,与联合双边滤波方法的线性混合;计算滤波器线性混合的插值系数;结合输入图像和引导图像,计算联合双边滤波的结果;使用插值系数将输入图像、联合双边滤波结果进行线性混合,得到纹理滤波后的图像。本发明充分考虑了输入图像结构边缘形状,保证滤波后的结果图像具有光滑的边缘和跟输入图像尽可能一致的边缘形状,在滤波结果中能够很好地保持输入图像的色彩效果,提升图像编辑处理的质量。
Description
技术领域
本发明涉及数字图像编辑及计算机视觉技术领域,具体涉及一种保持结构边缘的图像去纹理方法。
背景技术
图像通常包含结构信息和纹理信息。图像去纹理是指采用适当的方法分离图像的结构与纹理信息,在去除纹理信息时要尽可能保护图像结构信息。图像去纹理方法在图像分割、物体识别、显著性检测、图像增强和图像风格化等领域都具有广泛的应用。例如图像去纹理方法能够用于去除非机动车道图像中包含的地砖纹理,可以作为预处理操作来提高盲道的自动化检测效果等。因此,图像去纹理在计算机视觉和计算摄影学等领域具有重要的作用。
近年来,图像去纹理方法被广泛研究。现有的方法包括基于边缘感知的局部滤波方法、基于优化的方法和基于图像块的方法。
基于边缘感知的局部滤波方法,这些方法依靠像素梯度来区分纹理和结构边缘,并采用加权平均机制剔除图像纹理。在这类方法中双边滤波器和引导滤波器是两种最著名的边缘感知图像滤波器。然而,这类方法没有采用明确的度量方法来区分边缘和纹理,因此它们通常无法很好地去除纹理。另一方面,基于全局优化的方法通常需要求解一个大型线性系统,导致它们在计算量上劣势明显。基于图像块的方法根据图像块的特征统计,采用加权平均或者联合双边滤波器对图像进行滤波,获得纹理滤波后的图像。此类方法会破坏去纹理后的图像结构边缘形状,一些细小的结构无法保留,导致在图像后续增强处理时在边缘附近产生明显的瑕疵。
现在多采用联合双边滤波进行图像去纹理,联合双边滤波函数的表达式如下(Johannes Kopf,Michael F.Cohen,Dani Lischinski,and Matt Uyttendaele.Jointbilateral upsampling.ACM Transactions on Graphics.2007,26(3),Article 96):
其中:p表示图像内的当前像素,q是p邻域内的像素,是二维空间定义域滤波高斯核函数,其中σs是位置方差;是值域滤波高斯核函数,其中σr是色彩方差;M表示引导图像;I表示输入图像;J表示输出图像;是归一化系数。
为了保证能保留图像的小结构,一些方法在联合双边滤波中使用自适应空间尺度的引导图像,获得纹理滤波结果。然而这些方法只考虑了图像中纹理的空间尺度,可能会使滤波后输入图像中色彩差异不大的结构边缘也被去除,导致结构形状的缺失。此外基于联合双边滤波器的图像去纹理方法可能会导致边缘的梯度反置而产生诸如光晕等瑕疵。
为了取得在后续图像编辑操作中取得满意的结果,输入图像的色彩效果应该被保护,而且在纹理滤波后输出图像的结构边缘应当是光滑的,其边缘形状应该跟输入图像尽可能一致。因此,本发明提出了一种新的结构边缘保护的图像去纹理方法。该方法基于局部拉普拉斯滤波器,可以很好地保持图像滤波后结构边缘的形状以及保持输入图像的色彩变化。
发明内容
本发明的目的:本发明提出了一种新的保持结构边缘的图像去纹理方法。该方法提出了一种新的计算引导图像方法,结合引导图像和局部拉普拉斯滤波器可以很好地利用各自优点,有效解决传统联合双边滤波纹理滤波方法的问题,较好地在去除纹理信息同时能保护输入图像的结构边缘形状。
一种新的保持结构边缘的图像去纹理方法,其步骤如下:
步骤1:给定一幅输入图像I,I中像素p的强度值记为Ip。给定奇数k值,像素p的邻域N(p)是以像素p为中心的大小为k×k的二维矩形区域。计算自适应纹理滤波空间尺度σs(p)。
步骤2:对输入图像进行增处理,增强图像中较弱的结构边缘。
步骤3:根据步骤一求得的自适应纹理滤波空间尺度和步骤二增强后图像,计算纹理滤波的引导图像M。
步骤4:根据步骤3中求得的引导图像M,结局部拉普拉斯图像滤波器,构造联合局部拉普拉斯滤波器,变换该滤波器使其等价于输入图像与联合双边滤波方法的线性混合。
步骤5:根据步骤1求得的自适应纹理滤波空间尺度、用户给定的色彩方差和步骤3求得的引导图像,计算联合局部拉普拉斯滤波器的插值系数。
步骤6:结合输入图像I和步骤3求得的引导图像,计算联合双边滤波的结果。
步骤7:结合输入图像I、步骤5求得的插值系数和步骤6求得的联合双边滤波结果,采用线性混合方法,得到图像去纹理的结果。
所述步骤1中计算自适应纹理滤波空间尺度σs(p)时,首先计算结构方向相对总变化图像,记作dRTV。计算方法见式(1):
其中,gσ(p,q)表示方差为σ2的高斯函数,ε=10-6用来防止公式中分母为零。表示像素q沿着结构方向角度φ的方向偏微分算子,见式(2):
其中,和分别表示水平和垂直方向的微分算子。在[0,2π]中均匀采样12个不同方向的φ值,使用上述方法计算像素p的12个方向对应的值,然后使用式(3)求得像素p的结构方向θp:
根据求得的结构方向θp,获得输入图像的结构方向相对变化图像
根据求得的结构方向相对变化图像,使用式(4)求得自适应纹理滤波空间尺度σs(p):
Round(·)表示四舍五入运算;λ的值默认为0.005;|N(p)|表示像素p的邻域N(p)内的像素数量。
步骤2中所述增强图像较弱的结构边缘,使用引导图像滤波器(Guided imagefilter,详细内容参见Kaiming He,Jian Sun,and Xiaoou Tang.Guided ImageFiltering.European Conference on ComputerVision.2010)对输入图像I进行增强处理,增强图像对比度,强化输入图像I内强度较弱的结构边缘,求得增强后图像D。
步骤3中所述计算纹理滤波的引导图像M,根据步骤1求得的自适应空间尺度σs(p)和步骤2求得的增强后图像D,对输入图像I进行联合双边滤波,得到引导图像M。联合双边滤波方法见式(5):
式中是归一化系数。
步骤4中所述构造联合局部拉普拉斯滤波器时,以边缘感知的局部拉普拉斯金字塔图像滤波方法(详细内容参见Paris S,Hasinoff S W,Kautz J.Local Laplacianfilters[J].Communications of the ACM,2015,58(3):81-91)为基础,考虑图像两层拉普拉斯金字塔,引入步骤3求得的引导图像M,修改局部拉普拉斯滤波器的变换函数空间r为式(6):
r(p)=p-(p-g)f(pm-gm) (6);
p表示输入图像I的像素,g表示该图像高斯金字塔系数,f(*)表示为连续函数,pm表示引导图像M的像素,gm表示引导图像的高斯金字塔系数。对于两层的滤波器,需要计算输出图像J的拉普拉斯金字塔的两层图像L0[J]和L1[J]。假设输出图像的残差L1[J]保持不变,即L1[J]=L1[I]。输出图像J的拉普拉斯金字塔的第0层图像是变换图像r(I)和其对应低通滤波图像之间的差见式(7):
p表示像素,表示构建拉普拉斯金字塔的归一化高斯核函数,*表示卷积操作。输入图像I的金字塔最精细层表示为金字塔系数g=Ip,引导图像M的金字塔最精细层的系数gm=Mp。将重新定义的映射函数r代入式(7),求得式(8):
将残差层图像L1[·]上采样,并加入到式(8),推导后求得联合局部拉普拉斯滤波器的输出图像J,表示为式(9):
q表示邻域Ωp内的像素。结合联合双边滤波和式(9),发现联合局部拉普拉斯滤波器跟联合双边滤波器有相似之处:两层的联合局部拉普拉斯滤波器输出图像J的第二项是使用和连续函数f在空间邻域的加权平均。将f定义为引导图像M像素值范围偏差的高斯核函数,看到输出图像J的定义与联合双边滤波方法是高度相关的。
所述步骤4中变换联合局部拉普拉斯滤波器,其等价于输入图像与联合双边滤波方法的线性插值,令连续函数f为值域高斯核函数则变换过程表示为式(10)~(14):
Jp=Ip+μp(JBFp-Ip) (13);
Jp=(1-μp)Ip+μpJBFp (14);
其中表示归一化高斯核函数,是二维空间定义域滤波高斯核函数,其中σs是位置方差,其值为步骤1中求得的自适应纹理滤波空间尺度σs(p);是值域滤波高斯核函数,其中σr是色彩方差;μ表示插值系数,值为JBF为输入图像I与引导图像M的联合双边滤波结果。
步骤5中所述计算插值系数μ,根据步骤4中推导结果,使用式(15)计算滤波器插值系数μ:
其中q是p邻域N(p)内的像素,是归一化系数。
步骤6中所述计算联合双边滤波结果,根据自适应空间滤波尺度σs和用户指定的色彩方差σr,使用引导图像M,对输入图像进行联合双边滤波,得到相应的结果JBF,见式(16):
步骤7中所述计算图像去纹理结果,结合输入图像I、步骤5求得的插值系数μ和步骤6求得的联合双边滤波结果JBF,采用步骤4中推导出的线性混合方法,得到图像去纹理的结果J见式(17):
J=(1-μ)I+μ·JBF (17);
如果一次滤波效果不满意,可以将步骤7中求得的J赋值为新的输入图像,重复上述步骤,得到新的图像去纹理结果。一般重复迭代3-5次就可以获得满意的图像去纹理结果。
至此就完成了输入图像的去纹理操作。
本发明涉及的基于联合局部拉普拉斯滤波图像去纹理的优势在于,采用新技术手段,保证在纹理滤波后保护输入图像结构边缘的形状。结合引导图像和局部拉普拉斯滤波方法能够很好地利用各自的优点,有效解决传统联合双边纹理滤波方法的问题,实现高质量的纹理滤波。
附图说明
图1为本发明的基本流程示意图。
具体实施方式
下面结合附图对本发明进行详细说明;本实施例是以本发明技术方案为前提,进行实施的,并结合了详细的实施方式和过程。
如图1所示,本实施例所描述的一种新的结构边缘保护的图像去纹理方法包括如下步骤:
步骤1:给定一幅输入图像I,I中像素p的强度值记为Ip。给定奇数k值,像素p的邻域N(p)是以像素p为中心的大小为k×k的二维矩形区域。计算自适应纹理滤波空间尺度σs(p)。,首先计算结构方向相对总变化图像,记作dRTV。计算方法见式(1):
其中,gσ(p,q)表示方差为σ2的高斯函数,ε=10-6用来防止公式中分母为零。表示像素q沿着结构方向角度φ的方向偏微分算子,见式(2):
其中,和分别表示水平和垂直方向的微分算子。在[0,2π]中均匀采样12个不同方向的φ值,使用上述方法计算像素p的12个方向对应的值,然后使用式(3)求得像素p的结构方向θp:
根据求得的结构方向θp,获得输入图像的结构方向相对变化图像
根据求得的结构方向相对变化图像,使用式(4)求得自适应纹理滤波空间尺度σs(p):
Round(·)表示四舍五入运算;λ的值默认为0.005;|N(p)|表示像素p的邻域N(p)内的像素数量。
步骤2:对输入图像进行增处理,增强图像中较弱的结构边缘。增强图像较弱的结构边缘,使用引导图像滤波器(Guided image filter,详细内容参见Kaiming He,JianSun,and Xiaoou Tang.Guided Image Filtering.European Conference on ComputerVision.2010)对输入图像I进行增强处理,增强图像对比度,强化输入图像I内强度较弱的结构边缘,求得增强后图像D。
步骤3:根据步骤一求得的自适应纹理滤波空间尺度和步骤二增强后图像,计算纹理滤波的引导图像M。根据步骤1求得的自适应空间尺度σs(p)和步骤2求得的增强后图像D,对输入图像I进行联合双边滤波,得到引导图像M。联合双边滤波方法见式(5):
式中是归一化系数。
步骤4:根据步骤3中求得的引导图像M,结局部拉普拉斯图像滤波器,构造联合局部拉普拉斯滤波器,变换该滤波器使其等价于输入图像与联合双边滤波方法的线性混合。
构造联合局部拉普拉斯滤波器时,以边缘感知的局部拉普拉斯金字塔图像滤波方法(详细内容参见Paris S,Hasinoff S W,Kautz J.Local Laplacian filters[J].Communications of the ACM,2015,58(3):81-91)为基础,考虑图像两层拉普拉斯金字塔,引入步骤3求得的引导图像M,修改局部拉普拉斯滤波器的变换函数空间r为式(6):
r(p)=p-(p-g)f(pm-gm) (6);
p表示输入图像I的像素,g表示该图像高斯金字塔系数,f(*)表示为连续函数,pm表示引导图像M的像素,gm表示引导图像的高斯金字塔系数。对于两层的滤波器,需要计算输出图像J的拉普拉斯金字塔的两层图像L0[J]和L1[J]。假设输出图像的残差L1[J]保持不变,即L1[J]=L1[I]。输出图像J的拉普拉斯金字塔的第0层图像是变换图像r(I)和其对应低通滤波图像之间的差见式(7):
p表示像素,表示构建拉普拉斯金字塔的归一化高斯核函数,*表示卷积操作。输入图像I的金字塔最精细层表示为金字塔系数g=Ip,引导图像M的金字塔最精细层的系数gm=Mp。将重新定义的映射函数r代入式(7),求得式(8):
将残差层图像L1[·]上采样,并加入到式(8),推导后求得联合局部拉普拉斯滤波器的输出图像J,表示为式(9):
q表示邻域Ωp内的像素。结合联合双边滤波和式(9),发现联合局部拉普拉斯滤波器跟联合双边滤波器有相似之处:两层的联合局部拉普拉斯滤波器输出图像J的第二项是使用和连续函数f在空间邻域的加权平均。将f定义为引导图像M像素值范围偏差的高斯核函数,看到输出图像J的定义与联合双边滤波方法是高度相关的。
变换联合局部拉普拉斯滤波器,其等价于输入图像与联合双边滤波方法的线性插值,令连续函数f为值域高斯核函数则变换过程表示为式(10)~(14):
Jp=Ip+μp(JBFp-Ip) (13);
Jp=(1-μp)Ip+μpJBFp (14);
其中表示归一化高斯核函数,是二维空间定义域滤波高斯核函数,其中σs是位置方差,其值为步骤1中求得的自适应纹理滤波空间尺度σs(p);是值域滤波高斯核函数,其中σr是色彩方差;μ表示插值系数,值为JBF为输入图像I与引导图像M的联合双边滤波结果。
步骤5:根据步骤1求得的自适应纹理滤波空间尺度、用户给定的色彩方差和步骤3求得的引导图像,计算联合局部拉普拉斯滤波器的插值系数。计算插值系数μ,根据步骤4中推导结果,使用式(15)计算滤波器插值系数μ:
其中q是p邻域N(p)内的像素,是归一化系数。
步骤6:结合输入图像I和步骤3求得的引导图像,计算联合双边滤波的结果。根据自适应空间滤波尺度σs和用户指定的色彩方差σr,使用引导图像M,对输入图像进行联合双边滤波,得到相应的结果JBF,见式(16):
步骤7:结合输入图像I、步骤5求得的插值系数和步骤6求得的联合双边滤波结果,采用线性混合方法,得到图像去纹理的结果。计算图像去纹理结果,结合输入图像I、步骤5求得的插值系数μ和步骤6求得的联合双边滤波结果JBF,采用步骤4中推导出的线性混合方法,得到图像去纹理的结果J见式(17):
J=(1-μ)I+μ·JBF (17);
如果一次滤波效果不满意,可以将步骤7中求得的J赋值为新的输入图像,重复上述步骤,得到新的图像去纹理结果。一般重复迭代3-5次就可以获得满意的图像去纹理结果。
至此就完成了输入图像的去纹理操作。
Claims (8)
1.一种新的保持结构边缘的图像去纹理方法,其特征在于:包括如下步骤:
步骤1:给定一幅输入图像I,I中像素p的强度值记为Ip;给定奇数k值,像素p的邻域N(p)是以像素p为中心的大小为k×k的二维矩形区域;计算自适应纹理滤波空间尺度σs(p);
步骤2:对输入图像进行增处理,增强图像中较弱的结构边缘;
步骤3:根据步骤一求得的自适应纹理滤波空间尺度和步骤二增强后图像,计算纹理滤波的引导图像M;
步骤4:根据步骤3中求得的引导图像M,结局部拉普拉斯图像滤波器,构造联合局部拉普拉斯滤波器,变换该滤波器使其等价于输入图像与联合双边滤波方法的线性混合;
步骤5:根据步骤1求得的自适应纹理滤波空间尺度、用户给定的色彩方差和步骤3求得的引导图像,计算联合局部拉普拉斯滤波器的插值系数;
步骤6:结合输入图像I和步骤3求得的引导图像,计算联合双边滤波的结果;
步骤7:结合输入图像I、步骤5求得的插值系数和步骤6求得的联合双边滤波结果,采用线性混合方法,得到图像去纹理的结果。
2.如权利要求1所述的一种新的保持结构边缘的图像去纹理方法,其特征在于:步骤1中所述计算自适应纹理滤波空间尺度σs(p),首先计算结构方向相对总变化图像,记作dRTV;计算方法见式⑴:
其中,gσ(p,q)表示方差为σ2的高斯函数,ε=10-6用来防止公式中分母为零;表示像素q沿着结构方向角度φ的方向偏微分算子,见式⑵:
其中,和分别表示水平和垂直方向的微分算子;在[0,2π]中均匀采样12个不同方向的φ值,使用上述方法计算像素p的12个方向对应的值,然后使用式⑶求得像素p的结构方向θp:
根据求得的结构方向θp,获得输入图像的结构方向相对变化图像
根据求得的结构方向相对变化图像,使用式⑷求得自适应纹理滤波空间尺度σs(p):
Round(·)表示四舍五入运算;λ的值默认为0.005;|N(p)|表示像素p的邻域N(p)内的像素数量。
3.如权利要求1所述的一种新的保持结构边缘的图像去纹理方法,其特征在于:步骤2中所述增强图像较弱的结构边缘,使用引导图像滤波器对输入图像I进行增强处理,增强图像对比度,强化输入图像I内强度较弱的结构边缘,求得增强后图像D。
4.如权利要求1所述的一种新的保持结构边缘的图像去纹理方法,其特征在于:步骤3中所述计算纹理滤波的引导图像M,根据步骤1求得的自适应空间尺度σs(p)和步骤2求得的增强后图像D,对输入图像I进行联合双边滤波,得到引导图像M;联合双边滤波方法见式⑸:
式中是归一化系数。
5.如权利要求1所述的一种新的保持结构边缘的图像去纹理方法,其特征在于:步骤4中所述构造联合局部拉普拉斯滤波器,以边缘感知的局部拉普拉斯金字塔图像滤波方法为基础,考虑图像两层拉普拉斯金字塔,引入步骤3求得的引导图像M,修改局部拉普拉斯滤波器的变换函数空间r为式⑹:
r(p)=p-(p-g)f(pm-gm) ⑹;
p表示输入图像I的像素,g表示该图像高斯金字塔系数,f(*)表示为连续函数,pm表示引导图像M的像素,gm表示引导图像的高斯金字塔系数;对于两层的滤波器,需要计算输出图像J的拉普拉斯金字塔的两层图像L0[J]和L1[J];假设输出图像的残差L1[J]保持不变,即L1[J]=L1[I];输出图像J的拉普拉斯金字塔的第0层图像是变换图像r(I)和其对应低通滤波图像之间的差见式⑺:
p表示像素,表示构建拉普拉斯金字塔的归一化高斯核函数,*表示卷积操作;输入图像I的金字塔最精细层表示为金字塔系数g=Ip,引导图像M的金字塔最精细层的系数gm=Mp;将重新定义的映射函数r代入式⑺,求得式⑻:
将残差层图像L1[·]上采样,并加入到式⑻,推导后求得联合局部拉普拉斯滤波器的输出图像J,表示为式⑼:
q表示邻域Ωp内的像素;结合联合双边滤波和式⑼,发现联合局部拉普拉斯滤波器跟联合双边滤波器有相似之处:两层的联合局部拉普拉斯滤波器输出图像J的第二项是使用和连续函数f在空间邻域的加权平均;将f定义为引导图像M像素值范围偏差的高斯核函数,看到输出图像J的定义与联合双边滤波方法是高度相关的;
变换联合局部拉普拉斯滤波器,其等价于输入图像与联合双边滤波方法的线性插值,令连续函数f为值域高斯核函数则变换过程表示为式⑽~⒁:
Jp=Ip+μp(JBFp-Ip) ⒀;
Jp=(1-μp)Ip+μpJBFp ⒁;
其中表示归一化高斯核函数,是二维空间定义域滤波高斯核函数,其中σs是位置方差,其值为步骤1中求得的自适应纹理滤波空间尺度σs(p);是值域滤波高斯核函数,其中σr是色彩方差;μ表示插值系数,值为JBF为输入图像I与引导图像M的联合双边滤波结果。
6.如权利要求1所述的一种新的保持结构边缘的图像去纹理方法,其特征在于:步骤5中所述计算插值系数μ,根据步骤4中推导结果,使用式⒂计算滤波器插值系数μ:
其中q是p邻域N(p)内的像素,是归一化系数。
7.如权利要求1所述的一种新的保持结构边缘的图像去纹理方法,其特征在于:步骤6中所述计算联合双边滤波结果,根据自适应空间滤波尺度σs和用户指定的色彩方差σr,使用引导图像M,对输入图像进行联合双边滤波,得到相应的结果JBF,见式⒃:
8.如权利要求1所述的一种新的保持结构边缘的图像去纹理方法,其特征在于:步骤7中所述计算图像去纹理结果,结合输入图像I、步骤5求得的插值系数μ和步骤6求得的联合双边滤波结果JBF,采用步骤4中推导出的线性混合方法,得到图像去纹理的结果J见式⒄:
J=(1-μ)I+μ·JBF ⒄;
如果一次滤波效果不满意,将步骤7中求得的J赋值为新的输入图像,重复上述步骤一至步骤七,得到新的图像去纹理结果;一般重复迭代3-5次就可以获得满意的图像去纹理结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910497118.2A CN110246099B (zh) | 2019-06-10 | 2019-06-10 | 一种保持结构边缘的图像去纹理方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910497118.2A CN110246099B (zh) | 2019-06-10 | 2019-06-10 | 一种保持结构边缘的图像去纹理方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110246099A true CN110246099A (zh) | 2019-09-17 |
CN110246099B CN110246099B (zh) | 2021-09-07 |
Family
ID=67886459
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910497118.2A Active CN110246099B (zh) | 2019-06-10 | 2019-06-10 | 一种保持结构边缘的图像去纹理方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110246099B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112819733A (zh) * | 2021-01-29 | 2021-05-18 | 成都国科微电子有限公司 | 一种定向双边图像滤波方法及装置 |
CN112862715A (zh) * | 2021-02-08 | 2021-05-28 | 天津大学 | 一种实时且可控的尺度空间滤波方法 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050100241A1 (en) * | 2003-11-07 | 2005-05-12 | Hao-Song Kong | System and method for reducing ringing artifacts in images |
US20120308153A1 (en) * | 2011-06-03 | 2012-12-06 | Sunghyun Hwang | Device and method of removing noise in edge area |
CN104458766A (zh) * | 2014-12-31 | 2015-03-25 | 江南大学 | 一种基于结构纹理法的布匹表面瑕疵检测方法 |
CN104991269A (zh) * | 2015-06-04 | 2015-10-21 | 中国科学技术大学 | 一种边缘引导和结构约束的全波形反演快速方法 |
CN105913396A (zh) * | 2016-04-11 | 2016-08-31 | 湖南源信光电科技有限公司 | 一种噪声估计的图像边缘保持混合去噪方法 |
CN107133924A (zh) * | 2017-03-31 | 2017-09-05 | 长安大学 | 一种利用颜色二阶变化信息的保结构特征图像滤波方法 |
CN107274363A (zh) * | 2017-06-02 | 2017-10-20 | 北京理工大学 | 一种具有尺度敏感特性的边缘保持图像滤波方法 |
CN109859145A (zh) * | 2019-02-27 | 2019-06-07 | 长安大学 | 一种基于多级权重相对总变差的图像去纹理方法 |
-
2019
- 2019-06-10 CN CN201910497118.2A patent/CN110246099B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050100241A1 (en) * | 2003-11-07 | 2005-05-12 | Hao-Song Kong | System and method for reducing ringing artifacts in images |
US20120308153A1 (en) * | 2011-06-03 | 2012-12-06 | Sunghyun Hwang | Device and method of removing noise in edge area |
CN104458766A (zh) * | 2014-12-31 | 2015-03-25 | 江南大学 | 一种基于结构纹理法的布匹表面瑕疵检测方法 |
CN104991269A (zh) * | 2015-06-04 | 2015-10-21 | 中国科学技术大学 | 一种边缘引导和结构约束的全波形反演快速方法 |
CN105913396A (zh) * | 2016-04-11 | 2016-08-31 | 湖南源信光电科技有限公司 | 一种噪声估计的图像边缘保持混合去噪方法 |
CN107133924A (zh) * | 2017-03-31 | 2017-09-05 | 长安大学 | 一种利用颜色二阶变化信息的保结构特征图像滤波方法 |
CN107274363A (zh) * | 2017-06-02 | 2017-10-20 | 北京理工大学 | 一种具有尺度敏感特性的边缘保持图像滤波方法 |
CN109859145A (zh) * | 2019-02-27 | 2019-06-07 | 长安大学 | 一种基于多级权重相对总变差的图像去纹理方法 |
Non-Patent Citations (2)
Title |
---|
JUNHO JEON等: "Scale-aware Structure-Preserving Texture Filtering", 《PACIFIC GRAPHICS 2016》 * |
PARIS S等: "Local Laplacian filters:Edge-aware Image Processing with a Laplacian Pyramid", 《COMMUNICATIONS OF THE ACM》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112819733A (zh) * | 2021-01-29 | 2021-05-18 | 成都国科微电子有限公司 | 一种定向双边图像滤波方法及装置 |
CN112819733B (zh) * | 2021-01-29 | 2024-04-16 | 成都国科微电子有限公司 | 一种定向双边图像滤波方法及装置 |
CN112862715A (zh) * | 2021-02-08 | 2021-05-28 | 天津大学 | 一种实时且可控的尺度空间滤波方法 |
Also Published As
Publication number | Publication date |
---|---|
CN110246099B (zh) | 2021-09-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Guo et al. | LIME: Low-light image enhancement via illumination map estimation | |
Min et al. | Fast global image smoothing based on weighted least squares | |
CN103942758B (zh) | 基于多尺度融合的暗通道先验图像去雾方法 | |
Li et al. | Fast guided global interpolation for depth and motion | |
US8509559B2 (en) | Hierarchical motion deblurring method for single image | |
Xie et al. | Single depth image super resolution and denoising via coupled dictionary learning with local constraints and shock filtering | |
CN104933679B (zh) | 一种放大图像的方法及其对应系统 | |
CN107292834B (zh) | 红外图像细节增强方法 | |
Tan et al. | Multipoint filtering with local polynomial approximation and range guidance | |
CN110246099B (zh) | 一种保持结构边缘的图像去纹理方法 | |
Noori et al. | A bilateral image inpainting | |
CN107516302A (zh) | 一种基于OpenCV的混合图像增强的方法 | |
CN112215776A (zh) | 人像磨皮方法、电子设备和计算机可读存储介质 | |
Liu et al. | A generalized relative total variation method for image smoothing | |
Seo | Image denoising and refinement based on an iteratively reweighted least squares filter | |
CN105243647A (zh) | 一种基于线性空间滤波的图像增强方法 | |
CN102509266B (zh) | 一种图像快速保边滤波方法 | |
Chang et al. | A self-adaptive single underwater image restoration algorithm for improving graphic quality | |
CN105741255A (zh) | 图像融合方法及装置 | |
CN108717699B (zh) | 一种基于连续最小割的超声波图像分割方法 | |
Naqash et al. | Edge sharpening in grayscale images using modified Sobel technique | |
KR101715266B1 (ko) | 그래픽 가속기에 기반한 3d 모델의 라인 드로잉 방법 및 이를 이용한 프로그램이 기록된 컴퓨터 판독가능 기록매체 | |
Liu et al. | Infrared ship target image smoothing based on adaptive mean shift | |
CN107798663A (zh) | 一种基于偏微分方程和bm3d的无参数化图像恢复方法 | |
CN110738612B (zh) | X光透视图像的降噪方法及计算机可读存储介质 |
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 |