CN107292844A - 全变差正则化变分随机共振自适应暗图像滤波增强方法 - Google Patents
全变差正则化变分随机共振自适应暗图像滤波增强方法 Download PDFInfo
- Publication number
- CN107292844A CN107292844A CN201710467408.3A CN201710467408A CN107292844A CN 107292844 A CN107292844 A CN 107292844A CN 201710467408 A CN201710467408 A CN 201710467408A CN 107292844 A CN107292844 A CN 107292844A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- msup
- msubsup
- image
- 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 46
- 238000001914 filtration Methods 0.000 title claims abstract description 41
- 230000003044 adaptive effect Effects 0.000 title claims abstract description 26
- 230000002708 enhancing effect Effects 0.000 claims abstract description 47
- 238000005457 optimization Methods 0.000 claims description 15
- 238000004364 calculation method Methods 0.000 claims description 12
- 239000011159 matrix material Substances 0.000 claims description 7
- 238000012935 Averaging Methods 0.000 claims description 2
- 238000007781 pre-processing Methods 0.000 claims description 2
- 230000000007 visual effect Effects 0.000 abstract description 5
- 238000002059 diagnostic imaging Methods 0.000 abstract description 3
- 238000012544 monitoring process Methods 0.000 abstract description 3
- 230000004297 night vision Effects 0.000 abstract description 3
- 230000006978 adaptation Effects 0.000 abstract 1
- 238000000205 computational method Methods 0.000 abstract 1
- 238000012545 processing Methods 0.000 description 9
- 238000002474 experimental method Methods 0.000 description 8
- 238000012360 testing method Methods 0.000 description 6
- 239000003623 enhancer Substances 0.000 description 4
- 230000000694 effects Effects 0.000 description 3
- 230000015572 biosynthetic process Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000010801 machine learning Methods 0.000 description 2
- 238000000926 separation method Methods 0.000 description 2
- 241000208340 Araliaceae Species 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 238000003709 image segmentation Methods 0.000 description 1
- 230000002045 lasting effect Effects 0.000 description 1
- 230000002787 reinforcement Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000005728 strengthening Methods 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
- 238000005303 weighing 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/70—Denoising; Smoothing
-
- 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
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
- Image Processing (AREA)
Abstract
本发明公开了一种全变差正则化变分随机共振自适应暗图像滤波增强方法。该方法在变分框架下将全变差正则化滤波与随机共振增强进行耦合,通过全变差正则化图像去噪和随机共振图像增强两个步骤的交替迭代实现含噪声的低对比度暗图像的滤波去噪和对比度增强。本发明提出的模型自适应计算方法,一方面使本发明可以适用于不同的图像,提高算法适用性,另一方面也可以自适应保持图像亮度阈值在正常范围内。本发明可以在增强图像对比度的同时,有效地抑制噪声并保持图像的细节成分,增强之后的图像具有良好视觉效果,在医学成像、夜视安全监控等领域具有广泛的应用前景。
Description
技术领域
本发明属于图像增强的技术,具体涉及一种全变差正则化变分随机共振自适应暗图像滤波增强方法,
背景技术
低对比度暗图像的增强在医学成像、夜视监测、遥感等领域具有广泛的应用前景。目前,国际上关于暗图像增强问题已有许多研究,如直方图均衡化(HE)(Q.Wang,L.Chenand D.Shen,Fast Histogram Equalization for Medical Image Enhancement,30thAnnual International IEEE EMBS Conference Vancouver,2008),单尺度Retinex(SSR)(Choi D H,Jang I H,Mi H K,et al.Color image enhancement using single-scaleretinex based on an improved image formation model.Signal ProcessingConference,2008,European.IEEE,2008:1-5),多尺度Retinex(MSR)(Rahman Z U,JobsonD J,Woodell G A.Multi-scale retinex for color image enhancement.InternationalConference on Image Processing,1996.Proceedings.IEEE,1996:1003-1006vol.3),动态随机共振(Dynamic Stochastic Resonance,DSR)等,其中DSR方法近年来吸引了越来越多的关注,出现了多种基于DSR的低对比度暗图像增强方法,例如空间域暗图像增强(Chouhan R,Jha R K,Biswas P K.Enhancement of dark and low-contrast imagesusing dynamic stochastic resonance[J].Iet Image Processing,2013,7(7):174-184),小波变换域(Chouhan R,Kumar C P,Kumar R,et al.Contrast Enhancement ofDark Images using Stochastic Resonance in Wavelet Domain.InternationalJournal of Machine Learning and Computing,vol.2,pp.711-715,2012.)等等。传统上,在图像处理中噪声是不受欢迎的。而在基于DSR的随机共振图像增强中,非线性系统中噪声的存在是取得最佳性能必不可少的因素。然而,这些噪声在图像对比度增强的同时也得到增强,导致基于传统DSR图像增强的方法得到的增强图像中有大量的噪声,严重影响图像的质量,并给后续的图像分割等带来很大的困难。
发明内容
本发明目的是提供一种全变差正则化变分随机共振自适应暗图像滤波增强方法。该方法利用全变差正则化与随机共振图像增强变分模型耦合,以达到增强暗图像的同时抑制噪声的目的。
实现本发明目的的技术解决方案为:一种全变差正则化变分随机共振自适应暗图像滤波增强方法,通过下面两个阶段共四个步骤来实现低对比度暗图像的亮度增强和噪声抑制:
阶段一:图像预处理阶段
加噪声步骤:依据当前给定的待处理的低对比度暗图像uo的方差,随机生成高斯白噪声,并加入到待处理的大小为M×N的低对比度暗图像uo中,生成实际处理的大小为M×N的含噪声低对比度暗图像f;
阶段二:图像自适应增强滤波阶段
通过下面三个步骤的交替迭代实现图像的自适应增强滤波:
步骤一:对第n步迭代计算得到的大小为M×N的增强滤波后的图像u(n),通过求解下面的优化问题实现图像的梯度滤波:
公式(1)中各符号含义为:
u(n):第n步迭代计算得到的大小为M×N的增强滤波后的图像;
图像u(n)在像素点(i,j)处的沿垂直方向的一阶差分;
图像u(n)在像素点(i,j)处的沿水平方向的一阶差分;
dx:对于图像u(n)沿垂直方向的一阶差分的近似;
dy:对于图像u(n)沿水平方向的一阶差分的近似;
通过对公式(1)求解,得到的对图像u(n)沿垂直方向的一阶差分的最优近似;
通过对公式(1)求解,得到的对图像u(n)沿水平方向的一阶差分的最优近似;
λ2:用于控制dx和dy对于图像u(n)沿垂直方向和水平方向的一阶差分的近似程度,取为充分大的正数;
步骤二:在步骤一求解得到和后,通过求解下面的随机共振优化问题,更新得到新的增强滤波图像u(n+1)
公式(2)中各符号含义为:
f:在预处理阶段中得到的大小为N×M的实际处理的含噪声低亮度图像;
u:大小为N×M的待求近似真实图像;
fij:图像f在像素点(i,j)处的灰度值;
uij:图像u在像素点(i,j)处的灰度值;
图像u在像素点(i,j)处的沿垂直方向的一阶差分;
图像u在像素点(i,j)处的沿水平方向的一阶差分;
λ1:用于调节控制正则化滤波项的滤波强度的参数;
a:用于调节控制随机共振增强项的增强强度的参数;
像素点(i,j)处的凸优化控制参数;
步骤三:判断增强滤波后的图像u(n+1)是否满足给定的迭代终止条件,如果满足,则迭代终止,输出最终的增强滤波后的图像u(n+1);否则,则依据u(n+1)来自适应更新模型参数,然后转回到阶段二的步骤一继续迭代计算,直到满足迭代终止条件为止。
阶段一加噪声步骤中所述的实际处理的含噪声低对比度暗图像f按照下面公式得到
f=uo+n (3)
其中uo是实际给定的大小为N×M的待处理低对比度暗图像,n是均值为零,方差为的高斯白噪声,为待处理的低对比度暗图像uo的方差。
优化模型(1)的解按照下面的计算公式计算:
其中i=1,2,…,N;j=1,2,…,M,max{x,y}表示取x,y中的最大值,ε>0为一个很小的正数,主要保证(4)式中的分母不等于零。
优化模型(2)的解按照下面的迭代方式进行计算:
对k=0,1,…,K,按照下面的格式迭代计算
其中Δt为迭代步长,K为一个正整数,div(·)为散度算子,为两个二维向量值矩阵,即矩阵的每个元素为一个二维向量。
优化模型(2)中的凸优化控制参数i=1,2,…,N;j=1,2,…,M满足下面的约束条件:
优化模型(2)中的凸优化控制参数的计算方式为
阶段二的步骤三中所采用的迭代终止条件为:
|mean(u(n+1))-TB|<δ (8)
其中mean(·)为平均值算子,TB∈(0,255)为设定的增强后图像的灰度平均值,δ>0为一个很小的数。
本发明与现有技术相比,其显著优点为:本发明方法能够自适应地实现低对比度暗图像增强,在有效增强图像亮度和对比度的同时,也可以有效地抑制图像噪声,得到的增强图像具有良好的视觉效果,本发明方法在医学成像、夜视监测、遥感等领域都有广泛的应用前景。
附图说明
图1为本发明方法整体算法结构及数据流程图。
图2为随机噪声加噪单元算法及数据流程图。
图3为梯度滤波计算单元算法及数据流程图。
图4为随机共振增强计算单元算法及数据流程图。
图5为参数自适应更新计算单元算法及数据流程图。
图6为本发明的实验测试图像,(a)原始Barbara图像,(b)待增强Barbara图像,(c)原始MR图像,(d)待增强MR图像。
图7为本发明方法与现有算法对低亮度含噪声Barbara增强图像的比较:(a)HE,(b)SSR,(c)MSR,(d)DSR,(e)DSR+DWT,(f)VSR+TV。
图8为本发明方法与现有算法对MR图像处理后得到图像的比较:(a)HE,(b)SSR,(c)MSR,(d)DSR,(e)DSR+DWT,(f)VSR+TV。
具体实施方式
下面结合附图对本发明作进一步说明。
附图1至图5给出了本发明方法的整体算法结构以及具体计算单元的算法及数据流程图。附图6为实验的测试图像。
本发明一种全变差正则化变分随机共振自适应暗图像滤波增强方法,记f是大小为N×M的待处理的低亮度含噪声图像,u是大小为N×M的待求近似真实图像,则耦合之后的全变差正则化变分随机共振自适应暗图像滤波增强变分模型可表示为:
其中称为全变差滤波正则项,称为随机共振增强正则项,称为数据保真项。在模型(1)中,和分别表示图像u在像素点(i,j)处的沿垂直方向和水平方向的一阶差分;λ1>0是用于调节控制正则化滤波项的滤波强度的参数,a>0是用于调节控制随机共振增强项的增强强度的参数,ρij则是图像像素点(i,j)处的凸优化控制参数,为保证优化问题(1)为凸优化问题,给出ρij的取值范围为:
在模型(1)中,待处理的低对比度含噪声图像f并不是实际中真正的观测图像uo,待处理图像f是在观测图像g上再加上一个加性的高斯白噪声之后得到的,在得到待处理的低对比度含噪声图像f后,用模型(1)进行同步滤波增强得到最后的图像。
最优化问题(1)的求解等价于下面的约束优化问题:
用拉格朗日乘子法,可将约束优化问题(3)转化为下面的非约束优化问题:
其中λ2>0为拉格朗日乘子。显然,当λ2→+∞时,公式(4)和公式(1)是完全等价的,在实际中通常取充分大的参数λ2,然后通过求解无约束优化问题(4)来进行求解。而利用交替迭代法,极小化问题(4)的解可通过下面的迭代格式来近似求解:
给定初值u(0)=f,对n=1,2,…交替计算下面两式:
其中公式(5)实际上是在全变差正则化基础上利用变量替换得到的,这里起到的作用是通过对图像梯度场的平滑来实现图像去噪滤波,公式(5)则是由随机共振正则化导出的,这里起到的作用就是图像亮度和对比度的增强。通过对(5)-(6)式的交替迭代,就可以实现同步实现图像去噪和滤波,从而得到高质量的增强图像。
本发明所提供的全变差正则化变分随机共振自适应暗图像滤波增强方法由四个处理单元构成,分别是:随机噪声加噪单元、梯度滤波计算单元、随机共振增强计算单元和参数自适应更新计算单元。
1.1随机噪声加噪单元
由于噪声的存在是随机共振图像增强的必要条件,而在实际中给出的低对比度观测图像u(0)中不含噪声或者含噪声非常少,无法起到利用噪声来实现共振增强的作用。因此为保证足够噪声的存在,在给定实际给定的大小为N×M的观测图像u(0)之后,在其中加入一个随机高斯白噪声之后得到真正的待处理低对比度含噪声图像f,因此这个计算单元称为随机噪声加噪单元。
所述随机噪声加噪单元,由下列计算步骤构成:
步骤1.输入一幅大小为N×M的低对比度观测图像uo;
步骤2.生成均值为零,方差为高斯白噪声(记为n),其中为图像uo的像素灰度方差;
步骤3.生成待处理的低对比度含噪声图像f=uo+n;
1.2梯度滤波计算单元
该计算单元的主要任务实现公式(5)所表示的极小化问题的求解,该极小化问题实际上实现了梯度场的滤波,因此成为梯度滤波计算单元。
事实上,该极小化问题是有解析解的,其解析解可表示为
其中i=1,2,…,N;j=1,2,…,M,max{x,y}表示取x,y中的最大值,ε>0为一个很小的正数,主要保证(6)式中的分母不等于零。
所述梯度滤波计算单元,由下列计算步骤构成:
步骤1.输入当前增强图像u(n),计算每个像素点(i,j)处沿垂直方向的一阶差分和沿水平方向的一阶差分i=1,2,…,N;j=1,2,…,M;
步骤2.计算u(n)的每个像素点(i,j)处的梯度模
步骤3.计算大小为N×M的软阈值收缩系数矩阵S,其中
这里max{x,y}表示取x,y中的最大值。
步骤4.计算并输出
1.3随机共振增强计算单元
该计算单元的主要任务实现公式(6)所表示的极小化问题的求解,该极小化问题实际上实现了图像对比度增强,因此成为随机共振增强计算单元。在本计算单元中,利用梯度下降算法求解极小化问题(6)。
所述随机共振增强计算单元,由下列计算步骤构成:
步骤1.输入已得到的增强图像u(n)、由梯度滤波计算单元计算得到的导数近似图像和内循环迭代次数K以及参数矩阵ρ(n);
步骤2.设定内迭代循环的初值u(n+1,0)=u(n),对k=0,1,2,…,K,循环计算步骤3~步骤4
步骤3.计算
步骤4.计算u(n+1,k+1)=u(n+1,k)+Δt[a(u(n+1,k)-ρ(n)·(u(n+1,k))3)+f-λ1λ2·Div],在本发明中,迭代步长默认为Δt=0.1;
步骤5.输出更新增强图像u(n+1)=u(n+1,K+1)
1.4参数自适应更新计算单元
在本发明提出的由公式(1)所表述的图像滤波增强模型中,参数ρij的取值非常关键。一方面,该参数影响随机共振增强的效果,另一方面,该参数也控制公式(1)所表示的最优化问题是否为凸优化问题,从而影响问题解的存在唯一性。在本发明中,从保持图像亮度持续增强和保持问题为凸优化问题这两方面,得到参数ρij的取值范围(i=1,2,…,N;j=1,2,…,M)为(2)式,显然ρij的选择和uij有关。在迭代格式(5)—(6)中,将对该参数进行自适应更新。
所述参数自适应更新计算单元,由下列计算步骤构成:
步骤1.输入当前得到的增强图像u(n+1)、待处理的低对比度含噪声图像f,设定的参数a,δ以及目标增强图像的平均亮度TB;
步骤2.若|mean(u(n+1))-TB|<δ,迭代终止并输出最终的增强图像u(n+1),否则转到步骤3进行参数的更新;
步骤3.输出新增强参数矩阵ρ(n+1),其在每个像素点(i,j)处的计算公式为
下面结合实例及附图对本发明作进一步说明。
1)实验条件:
实验所用的计算环境为Intel Core i5-5200U@2.2GHz双核CPU,内存为4GB,编程平台为Matlab R2012a。实验所用的测试图像为大小为256×256的国际标准测试图像Barbara图像和大小为220×222的MR图像。
2)实验内容:
首先将原始Barbara图像和MR图像的最大灰度值调节至0.2和0.4得到低对比度按图像,然后利用本发明提出的方法来进行处理和现有的一些常用方法进行比较,这些方法主要包括:
(1)HE:Q.Wang,L.Chen and D.Shen,Fast Histogram Equalization forMedical Image Enhancement,30th Annual International IEEE EMBS ConferenceVancouver,2008
(2)SSR:Choi D H,Jang I H,Mi H K,et al.Color image enhancement usingsingle-scale retinex based on an improved image formation model.SignalProcessing Conference,2008,European.IEEE,2008:1-5
(3)MSR:Rahman Z U,Jobson D J,Woodell G A.Multi-scale retinex forcolor image enhancemen.International Conference on Image Processing,1996.Proceedings.IEEE,1996,3:1003-1006
(4)DSR:Chouhan R,Jha R K,Biswas P K.Enhancement of dark and low-contrast images using dynamic stochastic resonance.Iet Image Processing,2013,7(7):174-184
(5)DSR+DWT:Chouhan R,Kumar C P,Kumar R,et al.Contrast Enhancement ofDark Images using Stochastic Resonance in Wavelet Domain.InternationalJournal of Machine Learning and Computing,2012,vol.2:711-715,
本实验中从四个指标对算法进行比较:
(1)PSNR,即峰值信噪比,该指标主要用于衡量算法的去噪性能;
(2)RCEF,即相对对比度增强因子,该指标主要用于衡量算法的对比度增强性能;
(3)IF-SNR,即信噪比改善因子,该指标主要用于衡量算法的复原性能;
(4)DSM,即分部分离测度,该指标主要用于衡量输出图像与原始图像增强的差异。
在本实验中,调节参数TB使得该算法满足迭代终止条件时,能够得到最高的相对对比度增强因子(RCEF)且在正常范围(0~10)。
表1中给出了利用几种方法的性能指标比较。实验结果表明:当输入图像为有原图像的Barbara时,在正常的增强情况下,相较于其它方法,本发明所提出的方法峰值信噪比和信噪比恢复因子都远优于它们,虽然直方图均衡化和MSR的相对对比度增强因子和分部分离测度的值较高,但却使得图像产生了失真的视觉效果。当输入图像为MR图像时,本发明的方法得到的相对对比度增强因子优于其它方法。由于原始MR图像是带有噪声的,因此不计算其增强后的峰值信噪比和信噪比改善因子。
表1本发明方法与其它算法在PSNR、IF-SNR、RCEF和DSM的比较
图7-图8给出了六种方法对于实验图像的处理后图像的比较。实验结果表明,其它方法在增强图像之后噪声也较为明显,本发明的方法在增强图像的同时可以有效抑制噪声保持图像较好的视觉效果。
本发明提出将现有的全变差图像去噪和随机共振图像增强方法进行耦合,将原有的基于偏微分方程形式的随机共振方程改写为变分形式,并将现有最常用的针对图像去噪的全变差正则化方法引入到变分形式的随机共振增强中,提出一种新的全变差正则化变分随机共振自适应暗图像滤波增强模型,通过对模型的求解,在增强图像对比度的同时,也抑制噪声的增强甚至起到对噪声抑制作用。同时,通过对模型参数的控制,可以保证所提出的模型在数学上属于凸优化问题,这也保证模型解的存在性和唯一性,同时也可以有效防止图像亮度的“过增强”现象,从而大大提高增强后的图像的视觉质量。
Claims (7)
1.一种全变差正则化变分随机共振自适应暗图像滤波增强方法,其特征在于:通过下面两个阶段共四个步骤来实现低对比度暗图像的亮度增强和噪声抑制:
阶段一:图像预处理阶段
加噪声步骤:依据当前给定的待处理的低对比度暗图像uo的方差,随机生成高斯白噪声,并加入到待处理的大小为M×N的低对比度暗图像uo中,生成实际处理的大小为M×N的含噪声低对比度暗图像f;
阶段二:图像自适应增强滤波阶段
通过下面三个步骤的交替迭代实现图像的自适应增强滤波:
步骤一:对第n步迭代计算得到的大小为M×N的增强滤波后的图像u(n),通过求解下面的优化问题实现图像的梯度滤波:
<mrow>
<mo>(</mo>
<msubsup>
<mi>d</mi>
<mi>x</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>,</mo>
<msubsup>
<mi>d</mi>
<mi>y</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>)</mo>
<mo>=</mo>
<mi>arg</mi>
<munder>
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
<mrow>
<msub>
<mi>d</mi>
<mi>x</mi>
</msub>
<mo>,</mo>
<msub>
<mi>d</mi>
<mi>y</mi>
</msub>
</mrow>
</munder>
<mfenced open = "{" close = "}">
<mtable>
<mtr>
<mtd>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>M</mi>
</munderover>
<msqrt>
<mrow>
<msubsup>
<mrow>
<mo>(</mo>
<msub>
<mi>d</mi>
<mi>x</mi>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msubsup>
<mrow>
<mo>(</mo>
<msub>
<mi>d</mi>
<mi>y</mi>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
<mn>2</mn>
</msubsup>
</mrow>
</msqrt>
<mo>+</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mfrac>
<msub>
<mi>&lambda;</mi>
<mn>2</mn>
</msub>
<mn>2</mn>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>M</mi>
</munderover>
<mo>&lsqb;</mo>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mo>&dtri;</mo>
<mi>x</mi>
</msub>
<msup>
<mi>u</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
</msup>
</mrow>
<mo>)</mo>
</mrow>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mrow>
<mo>(</mo>
<msub>
<mi>d</mi>
<mi>x</mi>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mo>&dtri;</mo>
<mi>y</mi>
</msub>
<msup>
<mi>u</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
</msup>
</mrow>
<mo>)</mo>
</mrow>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mrow>
<mo>(</mo>
<msub>
<mi>d</mi>
<mi>y</mi>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
</mrow>
公式(1)中各符号含义为:
u(n):第n步迭代计算得到的大小为M×N的增强滤波后的图像;
图像u(n)在像素点(i,j)处的沿垂直方向的一阶差分;
图像u(n)在像素点(i,j)处的沿水平方向的一阶差分;
dx:对于图像u(n)沿垂直方向的一阶差分的近似;
dy:对于图像u(n)沿水平方向的一阶差分的近似;
通过对公式(1)求解,得到的对图像u(n)沿垂直方向的一阶差分的最优近似;
通过对公式(1)求解,得到的对图像u(n)沿水平方向的一阶差分的最优近似;
λ2:用于控制dx和dy对于图像u(n)沿垂直方向和水平方向的一阶差分的近似程度,取为充分大的正数;
步骤二:在步骤一求解得到和后,通过求解下面的随机共振优化问题,更新得到新的增强滤波图像u(n+1)
<mrow>
<msup>
<mi>u</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msup>
<mo>=</mo>
<mi>arg</mi>
<munder>
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
<mi>u</mi>
</munder>
<mfenced open = "{" close = "}">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>a</mi>
<munderover>
<mo>&Sigma;</mo>
<mi>i</mi>
<mi>N</mi>
</munderover>
<munderover>
<mo>&Sigma;</mo>
<mi>j</mi>
<mi>M</mi>
</munderover>
<mrow>
<mo>(</mo>
<mfrac>
<msubsup>
<mi>&rho;</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mn>4</mn>
</mfrac>
<msubsup>
<mi>u</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
<mn>4</mn>
</msubsup>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<msubsup>
<mi>u</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mo>)</mo>
</mrow>
<mo>-</mo>
<munderover>
<mo>&Sigma;</mo>
<mi>i</mi>
<mi>N</mi>
</munderover>
<munderover>
<mo>&Sigma;</mo>
<mi>j</mi>
<mi>M</mi>
</munderover>
<msub>
<mi>f</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<msub>
<mi>u</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>+</mo>
<mfrac>
<mrow>
<msub>
<mi>&lambda;</mi>
<mn>1</mn>
</msub>
<msub>
<mi>&lambda;</mi>
<mn>2</mn>
</msub>
</mrow>
<mn>2</mn>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>M</mi>
</munderover>
<mo>&lsqb;</mo>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mo>&dtri;</mo>
<mi>x</mi>
</msub>
<mi>u</mi>
</mrow>
<mo>)</mo>
</mrow>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mrow>
<mo>(</mo>
<msubsup>
<mi>d</mi>
<mi>x</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
<mrow>
<mi>i</mi>
<mo>.</mo>
<mi>j</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mo>&dtri;</mo>
<mi>y</mi>
</msub>
<mi>u</mi>
</mrow>
<mo>)</mo>
</mrow>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mrow>
<mo>(</mo>
<msubsup>
<mi>d</mi>
<mi>y</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
<mrow>
<mi>i</mi>
<mo>.</mo>
<mi>j</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
公式(2)中各符号含义为:
f:在预处理阶段中得到的大小为N×M的实际处理的含噪声低亮度图像;
u:大小为N×M的待求近似真实图像;
fij:图像f在像素点(i,j)处的灰度值;
uij:图像u在像素点(i,j)处的灰度值;
图像u在像素点(i,j)处的沿垂直方向的一阶差分;
图像u在像素点(i,j)处的沿水平方向的一阶差分;
λ1:用于调节控制正则化滤波项的滤波强度的参数;
a:用于调节控制随机共振增强项的增强强度的参数;
像素点(i,j)处的凸优化控制参数;
步骤三:判断增强滤波后的图像u(n+1)是否满足给定的迭代终止条件,如果满足,则迭代终止,输出最终的增强滤波后的图像u(n+1);否则,则依据u(n+1)来自适应更新模型参数,然后转回到阶段二的步骤一继续迭代计算,直到满足迭代终止条件为止。
2.根据权利要求1所述的全变差正则化变分随机共振自适应暗图像滤波增强方法,其特征在于:阶段一加噪声步骤中所述的实际处理的含噪声低对比度暗图像f按照下面公式得到
f=uo+n (3)
其中uo是实际给定的大小为N×M的待处理低对比度暗图像,n是均值为零,方差为的高斯白噪声,为待处理的低对比度暗图像uo的方差。
3.根据权利要求1所述的全变差正则化变分随机共振自适应暗图像滤波增强方法,其特征在于:优化模型(1)的解按照下面的计算公式计算:
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mrow>
<mo>(</mo>
<msubsup>
<mi>d</mi>
<mi>x</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msub>
<mo>=</mo>
<mi>max</mi>
<mrow>
<mo>{</mo>
<mrow>
<msqrt>
<mrow>
<msubsup>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mo>&dtri;</mo>
<mi>x</mi>
</msub>
<msup>
<mi>u</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
</msup>
</mrow>
<mo>)</mo>
</mrow>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msubsup>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mo>&dtri;</mo>
<mi>y</mi>
</msub>
<msup>
<mi>u</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
</msup>
</mrow>
<mo>)</mo>
</mrow>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
<mn>2</mn>
</msubsup>
</mrow>
</msqrt>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<msub>
<mi>&lambda;</mi>
<mn>2</mn>
</msub>
</mfrac>
<mo>,</mo>
<mn>0</mn>
</mrow>
<mo>}</mo>
</mrow>
<mo>&CenterDot;</mo>
<mfrac>
<msub>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mo>&dtri;</mo>
<mi>x</mi>
</msub>
<msup>
<mi>u</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
</msup>
</mrow>
<mo>)</mo>
</mrow>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mrow>
<mi>&epsiv;</mi>
<mo>+</mo>
<msqrt>
<mrow>
<msubsup>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mo>&dtri;</mo>
<mi>x</mi>
</msub>
<msup>
<mi>u</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
</msup>
</mrow>
<mo>)</mo>
</mrow>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msubsup>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mo>&dtri;</mo>
<mi>y</mi>
</msub>
<msup>
<mi>u</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
</msup>
</mrow>
<mo>)</mo>
</mrow>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
<mn>2</mn>
</msubsup>
</mrow>
</msqrt>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mrow>
<mo>(</mo>
<msubsup>
<mi>d</mi>
<mi>y</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msub>
<mo>=</mo>
<mi>max</mi>
<mrow>
<mo>{</mo>
<mrow>
<msqrt>
<mrow>
<msubsup>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mo>&dtri;</mo>
<mi>x</mi>
</msub>
<msup>
<mi>u</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
</msup>
</mrow>
<mo>)</mo>
</mrow>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msubsup>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mo>&dtri;</mo>
<mi>y</mi>
</msub>
<msup>
<mi>u</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
</msup>
</mrow>
<mo>)</mo>
</mrow>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
<mn>2</mn>
</msubsup>
</mrow>
</msqrt>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<msub>
<mi>&lambda;</mi>
<mn>2</mn>
</msub>
</mfrac>
<mo>,</mo>
<mn>0</mn>
</mrow>
<mo>}</mo>
</mrow>
<mo>&CenterDot;</mo>
<mfrac>
<msub>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mo>&dtri;</mo>
<mi>y</mi>
</msub>
<msup>
<mi>u</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
</msup>
</mrow>
<mo>)</mo>
</mrow>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mrow>
<mi>&epsiv;</mi>
<mo>+</mo>
<msqrt>
<mrow>
<msubsup>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mo>&dtri;</mo>
<mi>x</mi>
</msub>
<msup>
<mi>u</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
</msup>
</mrow>
<mo>)</mo>
</mrow>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msubsup>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mo>&dtri;</mo>
<mi>y</mi>
</msub>
<msup>
<mi>u</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
</msup>
</mrow>
<mo>)</mo>
</mrow>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
<mn>2</mn>
</msubsup>
</mrow>
</msqrt>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
其中i=1,2,…,N;j=1,2,…,M,max{x,y}表示取x,y中的最大值,ε>0为一个很小的正数,保证(4)式中的分母不等于零。
4.根据权利要求1所述的全变差正则化变分随机共振自适应暗图像滤波增强方法,其特征在于:优化模型(2)的解按照下面的迭代方式进行计算:
对k=0,1,…,K,按照下面的格式迭代计算
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msup>
<mi>u</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
<mo>,</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
</msup>
<mo>=</mo>
<msup>
<mi>u</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msup>
<mi>u</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
<mo>,</mo>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msup>
<mo>=</mo>
<msup>
<mi>u</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
<mo>,</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</msup>
<mo>+</mo>
<mi>&Delta;</mi>
<mi>t</mi>
<mo>&lsqb;</mo>
<mi>a</mi>
<mrow>
<mo>(</mo>
<msup>
<mi>u</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
<mo>,</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</msup>
<mo>-</mo>
<mi>&rho;</mi>
<mo>&CenterDot;</mo>
<msup>
<mrow>
<mo>(</mo>
<msup>
<mi>u</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
<mo>,</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</msup>
<mo>)</mo>
</mrow>
<mn>3</mn>
</msup>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>f</mi>
<mo>-</mo>
<msub>
<mi>&lambda;</mi>
<mn>1</mn>
</msub>
<msub>
<mi>&lambda;</mi>
<mn>2</mn>
</msub>
<mo>&CenterDot;</mo>
<mi>d</mi>
<mi>i</mi>
<mi>v</mi>
<mrow>
<mo>(</mo>
<msup>
<mi>d</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
</msup>
<mo>-</mo>
<mo>&dtri;</mo>
<msup>
<mi>u</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
<mo>,</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</msup>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msup>
<mi>u</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msup>
<mo>=</mo>
<msup>
<mi>u</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
<mo>,</mo>
<mi>K</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msup>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
</mrow>
其中Δt为迭代步长,K为一个正整数,div(·)为散度算子,为两个二维向量值矩阵,即矩阵的每个元素为一个二维向量。
5.根据权利要求1所述的全变差正则化变分随机共振自适应暗图像滤波增强方法,其特征在于:优化模型(2)中的凸优化控制参数i=1,2,…,N;j=1,2,…,M满足下面的约束条件:
<mrow>
<mfrac>
<mn>1</mn>
<mrow>
<mn>3</mn>
<msup>
<mrow>
<mo>(</mo>
<msubsup>
<mi>u</mi>
<mi>ij</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo><</mo>
<msubsup>
<mi>&rho;</mi>
<mi>ij</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo><</mo>
<msup>
<mrow>
<mo>(</mo>
<msubsup>
<mi>u</mi>
<mi>ij</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<mfrac>
<msub>
<mi>f</mi>
<mi>ij</mi>
</msub>
<mrow>
<mi>a</mi>
<mo>&CenterDot;</mo>
<msubsup>
<mi>u</mi>
<mi>ij</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
</msubsup>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
<mo>.</mo>
</mrow>
6.根据权利要求1所述的全变差正则化变分随机共振自适应暗图像滤波增强方法,其特征在于:优化模型(2)中的凸优化控制参数的计算方式为
<mrow>
<msubsup>
<mi>&rho;</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<mo>&lsqb;</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mn>3</mn>
<msup>
<mrow>
<mo>(</mo>
<msubsup>
<mi>u</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<msubsup>
<mi>u</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<mfrac>
<msub>
<mi>f</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mrow>
<mi>a</mi>
<mo>&CenterDot;</mo>
<msubsup>
<mi>u</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
</msubsup>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<mo>,</mo>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<mi>N</mi>
<mo>;</mo>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<mi>M</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
</mrow>
7.根据权利要求1所述的全变差正则化变分随机共振自适应暗图像滤波增强方法,其特征在于:阶段二的步骤三中所采用的迭代终止条件为:
|mean(u(n+1))-TB|<δ (8)
其中mean(·)为平均值算子,TB∈(0,255)为设定的增强后图像的灰度平均值,δ>0为一个很小的数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710467408.3A CN107292844B (zh) | 2017-06-20 | 2017-06-20 | 全变差正则化变分随机共振自适应暗图像滤波增强方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710467408.3A CN107292844B (zh) | 2017-06-20 | 2017-06-20 | 全变差正则化变分随机共振自适应暗图像滤波增强方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107292844A true CN107292844A (zh) | 2017-10-24 |
CN107292844B CN107292844B (zh) | 2020-11-13 |
Family
ID=60097904
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710467408.3A Active CN107292844B (zh) | 2017-06-20 | 2017-06-20 | 全变差正则化变分随机共振自适应暗图像滤波增强方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107292844B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110363796A (zh) * | 2018-05-25 | 2019-10-22 | 哈尔滨工程大学 | 基于轻量卷积神经网络结构的前视声纳运动目标跟踪方法 |
CN110570381A (zh) * | 2019-09-17 | 2019-12-13 | 合肥工业大学 | 一种基于高斯总变差的半解耦图像分解暗光图像增强方法 |
CN114897833A (zh) * | 2022-05-13 | 2022-08-12 | 中国科学院光电技术研究所 | 一种多帧激励下的正则化随机共振强背景弱目标增强方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103236046A (zh) * | 2013-04-28 | 2013-08-07 | 南京理工大学 | 基于图像形态模糊隶属度的分数阶自适应相干斑滤波方法 |
CN103729847A (zh) * | 2013-12-26 | 2014-04-16 | 杭州电子科技大学 | 一种基于多重随机共振机制的图像边缘检测方法 |
CN106251315A (zh) * | 2016-08-23 | 2016-12-21 | 南京邮电大学 | 一种基于全变分的图像去噪方法 |
CN106408519A (zh) * | 2015-11-10 | 2017-02-15 | 青岛大学 | 一种基于全变分的非局部图像复原方法 |
-
2017
- 2017-06-20 CN CN201710467408.3A patent/CN107292844B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103236046A (zh) * | 2013-04-28 | 2013-08-07 | 南京理工大学 | 基于图像形态模糊隶属度的分数阶自适应相干斑滤波方法 |
CN103729847A (zh) * | 2013-12-26 | 2014-04-16 | 杭州电子科技大学 | 一种基于多重随机共振机制的图像边缘检测方法 |
CN106408519A (zh) * | 2015-11-10 | 2017-02-15 | 青岛大学 | 一种基于全变分的非局部图像复原方法 |
CN106251315A (zh) * | 2016-08-23 | 2016-12-21 | 南京邮电大学 | 一种基于全变分的图像去噪方法 |
Non-Patent Citations (4)
Title |
---|
NIDHI GUPTA: "Enhancement of Dark Images using Dynamic Stochastic Resonance in Combined DWT and DCT", 《2014 9TH INTERNATIONAL CONFERENCE ON INDUSTRIAL AND INFORMATION SYSTEMS (ICIIS)》 * |
RAJIB KUMAR JHA等: "Enhancement of Digital Images using Stochastic Resonance", 《 TENCON 2005 - 2005 IEEE REGION 10 CONFERENCE》 * |
RAJLAXMI CHOUHAN等: "Enhancement of dark and low-contrast images using dynamic stochastic resonance", 《 IET IMAGE PROCESSING 》 * |
余丽红: "基于全变分的磁共振图像去噪算法的研究", 《中国优秀硕士学位论文全文数据库(电子期刊)信息科技辑》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110363796A (zh) * | 2018-05-25 | 2019-10-22 | 哈尔滨工程大学 | 基于轻量卷积神经网络结构的前视声纳运动目标跟踪方法 |
CN110570381A (zh) * | 2019-09-17 | 2019-12-13 | 合肥工业大学 | 一种基于高斯总变差的半解耦图像分解暗光图像增强方法 |
CN110570381B (zh) * | 2019-09-17 | 2022-04-29 | 合肥工业大学 | 一种基于高斯总变差的半解耦图像分解暗光图像增强方法 |
CN114897833A (zh) * | 2022-05-13 | 2022-08-12 | 中国科学院光电技术研究所 | 一种多帧激励下的正则化随机共振强背景弱目标增强方法 |
CN114897833B (zh) * | 2022-05-13 | 2023-08-11 | 中国科学院光电技术研究所 | 一种多帧激励下的正则化随机共振强背景弱目标增强方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107292844B (zh) | 2020-11-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US9002131B2 (en) | Method and system for enhancing image quality | |
CN109636745B (zh) | 基于分数阶微分图像增强算法的最佳阶图像增强方法 | |
CN107292844A (zh) | 全变差正则化变分随机共振自适应暗图像滤波增强方法 | |
CN103440630A (zh) | 基于引导滤波器的大动态范围红外图像显示与细节增强方法 | |
CN102298774B (zh) | 基于联合相似性的非局部均值去噪方法 | |
CN103049891B (zh) | 基于自适应窗口选择的视频图像去模糊方法 | |
CN110473152B (zh) | 基于改进Retinex算法的图像增强方法 | |
Wang et al. | Low-light image joint enhancement optimization algorithm based on frame accumulation and multi-scale Retinex | |
CN110232670A (zh) | 一种基于高低频分离的图像视觉效果增强的方法 | |
CN102063704A (zh) | 一种机载视景增强方法及其装置 | |
CN103489168A (zh) | 一种红外图像自适应转伪彩增强方法及系统 | |
Ghosh et al. | A novel approach of retinal image enhancement using PSO system and measure of fuzziness | |
CN104616252B (zh) | 基于nsct和pcnn的数字图像增强方法 | |
TWI692970B (zh) | 影像處理電路及相關的影像處理方法 | |
CN109859138A (zh) | 一种基于人眼视觉特性的红外图像增强方法 | |
Liu et al. | Structure-aware texture filtering based on local histogram operator | |
Feng et al. | Low-light color image enhancement based on Retinex | |
CN103218792B (zh) | 一种基于配准的红外图像噪声时域滤波方法 | |
CN107292836B (zh) | 基于外部图像块先验信息与稀疏表示的图像盲去模糊方法 | |
CN106127700B (zh) | 一种凝视红外退化图像自适应复原方法 | |
Huang et al. | Low illumination color image enhancement based on dynamic bistable stochastic resonance | |
CN103810677A (zh) | 基于超经验模式分解的遥感图像增强方法 | |
Dong et al. | Multi-scale total variation with automated regularization parameter selection for color image restoration | |
Zhang et al. | SSIM-based optimal non-local means image denoising with improved weighted Kernel function | |
Han et al. | Image enhancement based on spatial multi-scale homomorphic filtering and local entropy guided image filtering |
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 |