CN107203980B - 自适应多尺度暗通道先验的水下目标探测图像增强方法 - Google Patents

自适应多尺度暗通道先验的水下目标探测图像增强方法 Download PDF

Info

Publication number
CN107203980B
CN107203980B CN201710405123.7A CN201710405123A CN107203980B CN 107203980 B CN107203980 B CN 107203980B CN 201710405123 A CN201710405123 A CN 201710405123A CN 107203980 B CN107203980 B CN 107203980B
Authority
CN
China
Prior art keywords
image
dark channel
scale
calculating
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
Application number
CN201710405123.7A
Other languages
English (en)
Other versions
CN107203980A (zh
Inventor
马金祥
朱锡芳
袁洪春
熊超
吴峰
范新南
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Changzhou Institute of Technology
Original Assignee
Changzhou Institute of Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Changzhou Institute of Technology filed Critical Changzhou Institute of Technology
Priority to CN201710405123.7A priority Critical patent/CN107203980B/zh
Publication of CN107203980A publication Critical patent/CN107203980A/zh
Application granted granted Critical
Publication of CN107203980B publication Critical patent/CN107203980B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/70Denoising; Smoothing
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/90Dynamic range modification of images or parts thereof
    • G06T5/92Dynamic range modification of images or parts thereof based on global image properties
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20024Filtering details

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Quality & Reliability (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Image Processing (AREA)
  • Color Image Communication Systems (AREA)
  • Image Analysis (AREA)

Abstract

本发明公开了一种自适应多尺度暗通道先验的水下目标探测图像增强方法,其步骤包括:获取水下目标原始彩色图像;计算原始图像的光照向量;计算原始图像的暗通道图像;计算暗通道图像不同尺度窗口对应的透射图;应用导向滤波方法计算原始彩色图像在不同尺度窗口透射图对应的导向滤波输出图像;导向滤波输出图像与透射图之间误差最小条件,初选出几个最合适的最小滤波尺度窗口;对筛选的导向滤波输出图像加权平均;应用暗通道先验理论,计算原始图像的恢复增强图像;对恢复图像J从均值、方差、对比度和信息熵等方面进行定量评价。本发明中用到的导向滤波函数,能在图像的纹理与平滑之间实现均衡处理,使得处理后的图像视觉质量提高、纹理信息丰富。

Description

自适应多尺度暗通道先验的水下目标探测图像增强方法
技术领域
本发明属于图像信息处理领域,具体涉及一种自适应多尺度暗通道先验的水下目标探测图像增强方法。
背景技术
水下目标探测图像存在非均匀亮度、低信噪比、低对比度等特殊情况,常用的水下目标探测图像增强算法主要分为修改水下图像的光照和抑制图像对比度以保留图像边缘两大类,但不可避免会降低探测图像的视觉质量。基于光照散射物理模型利用先验条件进行图像逆向还原,最经典的为何恺明博士的暗通道先验去雾方法。暗通道先验方法选用一个尺度参数计算彩色图像中对应尺度窗口中像素各分量的最小值,并将此值作为暗通道先验信息,然后通过暗通道先验信息估算场景的深度信息,计算逆向还原恢复图像。暗通道先验理论在单幅图像的去雾清晰化处理中取得了比较好的应用效果。由于受到水的光学特性以及水中各种微粒、浮游生物和水体流动的影响,直接平移和转嫁暗通道先验理论的研究成果进行水下探测图像增强存在一定的困难:一、由于只考虑了一个尺度参数,对于存在较大偏亮区域的探测图像,处理结果误差较大;二、图像增强处理过程是一个开环系统,不具备根据增强效果评价指标参数自适应调整尺度参数的可能性。
发明内容
本发明所要解决的技术问题,面对存在非均匀亮度、低信噪比、低对比度环境下的水下大坝裂缝检测图像的精确定位和准确描述的客观实际需求,研究自适应多尺度暗通道先验的水下目标探测图像增强方法,实现水下目标探测图像去噪处理,提高水下目标探测图像的视觉质量。
本发明采用如下方案实现:
自适应多尺度暗通道先验的水下目标探测图像增强方法,包括以下步骤:
第一步:获取水下目标原始彩色图像I;
第二步:计算图像I的光照向量A,A的具体计算过程为:I对应的灰度图像Gray(I)中最大像素值,即为光照向量A;
第三步:按He方法计算图像I的暗通道图像DC0(I),DC0(I)的具体计算过程为:
Figure BDA0001309506190000011
式中,Ic是I的一个颜色通道;
第四步:计算暗通道图像DC0(I)的m个不同尺度窗口Ωk对应的透射图Tk(I),Tk(I)的具体计算过程如下:
a)设定不同尺度窗口Ωk:Nk×Nk(k=1,2,…,m);
b)计算DC0(I)进行不同尺度窗口Ωk的最小值滤波输出
Figure BDA0001309506190000021
c)对
Figure BDA0001309506190000022
进行反色、归一化处理,即可得到对应的透射图Tk(I)。
第五步:按He方法应用导向滤波方法计算原始彩色图像I在透射图Tk(I)对应的导向滤波投射图Qk(I),Qk(I)的具体计算过程如下:
a)设置导向滤波尺度窗口Wk:4Nk×4Nk(k=1,2,…,m),并设定调整参数ε;
b)设置导向图像为原始图像I的灰度图像Gray(I),设置滤波输入图像为透射图Tk(I);
c)计算导向滤波输出图像Qk(I)。
第六步:根据导向滤波输出图像Qk(I)与透射图Tk(I)之间误差最小条件,初选出n个最合适的最小滤波尺度窗口Ωl′(l=1,2,…n;(n≤(m/2)),以及对应的导向滤波输出图像Ql′(I)。具体计算过程如下:
min∑wh|Qk(I)-Tk(I)|,式中w和h分别是I的宽度和高度像素值;
第七步:对初选的Ql′(I)(l=1,2,…n),按照权值向量d=[d1,d2,…,dn](其中
Figure BDA0001309506190000023
),进行加权平均处理得到Q(I)。具体计算过程如下:
Figure BDA0001309506190000024
第八步:应用暗通道先验理论,计算图像I的恢复图像J,具体计算过程如下:
Figure BDA0001309506190000025
式中q0取常数值,用于降低图像的噪声;
第九步:对恢复图像J从均值、方差、对比度和信息熵等方面进行定量评价,相关定量评价指标函数表示为:
均值:
Figure BDA0001309506190000026
标准差:
Figure BDA0001309506190000031
对比度:
Figure BDA0001309506190000032
式中,n=|i-j|,
Figure BDA0001309506190000033
为灰度共生矩阵
信息熵为:
Figure BDA0001309506190000034
式中,pn为像素值n出现的概率
本发明所达到的有益之处在于:
本发明的方法,可以只利用单幅非均匀亮度、低信噪比、低对比度水下目标探测图像自身的信息,对图像进行自适应多尺度暗通道先验的去噪增强处理。先通过投射图误差最小缩小尺度参数数量范围,自适应确定尺度空间,再进一步通过均值、方差、对比度和信息熵等综合定量评价指标评估暗通道先验的增强图像。本发明中用到的导向滤波函数,能在图像的纹理与平滑之间实现均衡处理,使得处理后的图像视觉质量提高、纹理信息丰富。
附图说明
图1是本发明的方法的控制流程图。
图2是自适应筛选尺度窗口控制流程图。
具体实施方式
下面结合附图对本发明作进一步描述。以下实施例仅用于更加清楚地说明本发明的技术方案,而不能以此来限制本发明的保护范围。
参照图1所示,本发明是一种自适应多尺度暗通道先验的水下目标探测图像增强方法,整体流程图如图1所示,具体实现步骤如下:
第一步:获取水下目标原始彩色图像I;
第二步:计算图像I的光照向量A,A的具体计算过程为:I对应的灰度图像Gray(I)中最大像素值,即为光照向量A;
第三步:按He方法计算图像I的暗通道图像DC0(I),DC0(I)的具体计算过程为:
Figure BDA0001309506190000035
式中,Ic是I的一个颜色通道;
第四步:计算暗通道图像DC0(I)的m个不同尺度窗口Ωk对应的透射图Tk(I),Tk(I)的具体计算过程如下:
d)设定不同尺度窗口Ωk:Nk×Nk(k=1,2,…,m);
e)计算DC0(I)进行不同尺度窗口Ωk的最小值滤波输出
Figure BDA0001309506190000041
f)对
Figure BDA0001309506190000042
进行反色、归一化处理,即可得到对应的透射图Tk(I)。
第五步:按He方法应用导向滤波方法计算原始彩色图像I在透射图Tk(I)对应的导向滤波投射图Qk(I),Qk(I)的具体计算过程如下:
d)设置导向滤波尺度窗口Wk:4Nk×4Nk(k=1,2,…,m),并设定调整参数ε;
e)设置导向图像为原始图像I的灰度图像Gray(I),设置滤波输入图像为透射图Tk(I);
f)计算导向滤波输出图像Qk(I)。
第六步:根据导向滤波输出图像Qk(I)与透射图Tk(I)之间误差最小条件,初选出n个最合适的最小滤波尺度窗口Ωl′(l=1,2,…n;(n≤(m/2)),以及对应的导向滤波输出图像Ql′(I)。具体计算过程如下:
min∑wh|Qk(I)-Tk(I)|,式中w和h分别是I的宽度和高度像素值;
第七步:对初选的Ql′(I)(l=1,2,…n),按照权值向量d=[d1,d2,…,dn](其中
Figure BDA0001309506190000043
),进行加权平均处理得到Q(I)。具体计算过程如下:
Figure BDA0001309506190000044
第八步:应用暗通道先验理论,计算图像I的恢复图像J,具体计算过程如下:
Figure BDA0001309506190000045
式中q0取常数值,用于降低图像的噪声;
第九步:对恢复图像J从均值、方差、对比度和信息熵等方面进行定量评价,相关定量评价指标函数表示为:
均值:
Figure BDA0001309506190000046
标准差:
Figure BDA0001309506190000047
对比度:
Figure BDA0001309506190000048
式中,n=|i-j|,
Figure BDA0001309506190000049
为灰度共生矩阵
信息熵为:
Figure BDA0001309506190000051
式中,pn为像素值n出现的概率
对自适应多尺度参数(或尺度窗口)图像增强有关问题的说明:
(1)暗通道原理存在的前提是图像中存在大片阴影和鲜艳的颜色,显然水下目标探测图像并不总是满足这样的特性。对颜色失真图像的数据分析后发现,当输入图像中某个通道的像素值小于空气光向量中对应值时会出现失真现象,两者差距越大,图像失真情况越严重。为减小失真,应适当增加透射率的取值;
(2)尺度参数的数量应适中,尺度参数太多会增加暗通道信息计算的复杂度,尺度参数太少则会导致自适应选择空间范围太小。经过尺度参数选择后,尺度参数数量范围至少应缩小为初始数量的一半以内;
(3)尺度参数的疏密程度和参数数量,可以根据输入图像的尺寸大小做适当的调整。根据导向滤波输入图像与导向滤波输出图像两者之间的误差最小,缩小尺度参数范围。对恢复图像从均值、方差、对比度和信息熵等方面进行定量评价,可以作为图像进一步增强对比度处理的依据;
(4)对基于暗通道原理的恢复图像进行进一步分析,如果图像灰度值局限在一个较小的区间,灰度层次少、反差较低,可以扩展图像灰度范围,提高图像对比度。可以根据实际情况,采取线性拉伸或者非线性拉伸方法。
综上所述,本发明的自适应多尺度暗通道先验的水下目标探测图像增强方法,包括如下步骤:步骤一:获取水下目标原始彩色图像;步骤二:计算原始图像的光照向量;步骤三:计算原始图像的暗通道图像;步骤四:计算暗通道图像不同尺度窗口对应的透射图;步骤五:应用导向滤波方法计算原始彩色图像在不同尺度窗口透射图对应的导向滤波输出图像;步骤六:导向滤波输出图像与透射图之间误差最小条件,初选出几个最合适的最小滤波尺度窗口;步骤七:对筛选的导向滤波输出图像加权平均;步骤八:应用暗通道先验理论,计算原始图像的恢复增强图像;步骤九:对恢复图像J从均值、方差、对比度和信息熵等方面进行定量评价。本发明中用到的导向滤波函数,能在图像的纹理与平滑之间实现均衡处理,使得处理后的图像视觉质量提高、纹理信息丰富。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明技术原理的前提下,还可以做出若干改进和变形,这些改进和变形也应视为本发明的保护范围。

Claims (5)

1.自适应多尺度暗通道先验的水下目标探测图像增强方法,其特征在于,包括以下步骤:
步骤1:获取水下目标原始彩色图像;
步骤2:计算原始图像的光照向量;
步骤3:计算原始图像的暗通道图像;
步骤4:计算暗通道图像不同尺度窗口对应的透射图;
步骤4具体包括如下步骤:
步骤4-1:设定不同尺度窗口Ωk:Nk×Nk,k=1,2,…,m;
步骤4-2:计算暗通道图像DC0(I)不同尺度窗口Ωk的最小值滤波输出
Figure FDA0002568040420000011
步骤4-3:对
Figure FDA0002568040420000012
进行反色、归一化处理,即可得到对应的透射图Tk(I);
步骤5:应用导向滤波方法计算原始彩色图像在不同尺度窗口透射图对应的导向滤波输出图像;
步骤5中,原始彩色图像I在透射图Tk(I)对应的导向滤波投射图Qk(I)的具体计算过程如下:
步骤5-1:设置导向滤波尺度窗口Wk:4Nk×4Nk,k=1,2,…,m,并设定调整参数ε;
步骤5-2:设置导向图像为原始图像I的灰度图像Gray(I),设置滤波输入图像为透射图Tk(I);
步骤5-3:计算导向滤波输出图像Qk(I);
步骤6:根据导向滤波输出图像与透射图之间误差最小条件,初选出几个最合适的最小滤波尺度窗口以及对应的导向滤波输出图像;
步骤6中,根据导向滤波输出图像Qk(I)与透射图Tk(I)之间误差最小条件,初选出n个最合适的最小滤波尺度窗口Ω′l,l=1,2,…,n;
Figure FDA0002568040420000013
以及对应的导向滤波输出图像Q′l(I);具体计算过程如下:
min∑wh|Qk(I)-Tk(I)|,
式中,w和h分别是I的宽度和高度;
步骤7:对筛选的导向滤波输出图像加权平均;
步骤8:应用暗通道先验理论,计算原始图像的恢复增强图像;
步骤9:对恢复图像从包括均值、方差、对比度和信息熵的方面进行定量评价。
2.根据权利要求1所述的自适应多尺度暗通道先验的水下目标探测图像增强方法,其特征在于,步骤3中暗通道图像的计算方法为:
Figure FDA0002568040420000021
式中,I是水下目标原始彩色图像;Ic是I的一个颜色通道;DC0(I)是暗通道图像。
3.根据权利要求1所述的自适应多尺度暗通道先验的水下目标探测图像增强方法,其特征在于,步骤7具体为对初选的Q′l(I),l=1,2,…,n,按照权值向量d=[d1,d2,…,dn],其中
Figure FDA0002568040420000022
进行加权平均处理得到Q(I),计算方法具体为:
Figure FDA0002568040420000023
式中,n为初选出最合适的最小滤波尺度窗口Ω′l的数量,l=1,2,…,n;
Figure FDA0002568040420000024
Q′l(I)为对应的导向滤波输出图像。
4.根据权利要求3所述的自适应多尺度暗通道先验的水下目标探测图像增强方法,其特征在于,步骤8中恢复增强图像J的具体计算过程如下:
Figure FDA0002568040420000025
式中,A为光照向量;q0取常数值,用于降低图像的噪声。
5.根据权利要求4所述的自适应多尺度暗通道先验的水下目标探测图像增强方法,其特征在于,步骤9中,均值、方差、对比度和信息熵的定量评价指标函数表示为:
均值:
Figure FDA0002568040420000031
标准差:
Figure FDA0002568040420000032
对比度:
Figure FDA0002568040420000033
式中,t=|i-j|,
Figure FDA0002568040420000034
为灰度共生矩阵;
信息熵为:
Figure FDA0002568040420000035
式中,Pt为像素值t出现的概率。
CN201710405123.7A 2017-05-31 2017-05-31 自适应多尺度暗通道先验的水下目标探测图像增强方法 Active CN107203980B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710405123.7A CN107203980B (zh) 2017-05-31 2017-05-31 自适应多尺度暗通道先验的水下目标探测图像增强方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710405123.7A CN107203980B (zh) 2017-05-31 2017-05-31 自适应多尺度暗通道先验的水下目标探测图像增强方法

Publications (2)

Publication Number Publication Date
CN107203980A CN107203980A (zh) 2017-09-26
CN107203980B true CN107203980B (zh) 2020-10-16

Family

ID=59907193

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710405123.7A Active CN107203980B (zh) 2017-05-31 2017-05-31 自适应多尺度暗通道先验的水下目标探测图像增强方法

Country Status (1)

Country Link
CN (1) CN107203980B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108596249B (zh) * 2018-04-24 2022-04-22 苏州晓创光电科技有限公司 图像特征提取及分类的方法和装置
CN110060243B (zh) * 2019-04-11 2022-06-28 国家电网有限公司 一种大坝表面裂缝骨架提取方法
CN110796611B (zh) * 2019-09-29 2022-04-15 深圳市华汉伟业科技有限公司 一种阴影和自由曲面的图像补正方法及其图像补正装置
CN111968062B (zh) * 2020-09-07 2022-12-09 新疆大学 基于暗通道先验镜面高光图像增强方法、装置及存储介质
CN115063404B (zh) * 2022-07-27 2022-11-08 建首(山东)钢材加工有限公司 基于x射线探伤的耐候钢焊缝质量检测方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104008528A (zh) * 2014-05-21 2014-08-27 河海大学常州校区 基于阈值分割的非均匀光场水下目标探测图像增强方法
CN106056559A (zh) * 2016-06-30 2016-10-26 河海大学常州校区 基于暗通道先验的非均匀光场水下目标探测图像增强方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101448164B1 (ko) * 2013-04-22 2014-10-14 금오공과대학교 산학협력단 파라미터 최적화를 이용한 영상의 안개 제거 방법
US9349170B1 (en) * 2014-09-04 2016-05-24 The United States Of America As Represented By The Secretary Of The Navy Single image contrast enhancement method using the adaptive wiener filter

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104008528A (zh) * 2014-05-21 2014-08-27 河海大学常州校区 基于阈值分割的非均匀光场水下目标探测图像增强方法
CN106056559A (zh) * 2016-06-30 2016-10-26 河海大学常州校区 基于暗通道先验的非均匀光场水下目标探测图像增强方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
雾天降质图像增强处理技术研究;曾宇骏 等;《中国优秀硕士学位论文全文数据库 信息科技辑》;20160315;参见9-25页 *

Also Published As

Publication number Publication date
CN107203980A (zh) 2017-09-26

Similar Documents

Publication Publication Date Title
CN107203980B (zh) 自适应多尺度暗通道先验的水下目标探测图像增强方法
CN105654436B (zh) 一种基于前景背景分离的背光图像增强去噪方法
CN107767354B (zh) 一种基于暗原色先验的图像去雾算法
CN107527332B (zh) 基于改进Retinex的低照度图像色彩保持增强方法
CN108765336B (zh) 基于暗亮原色先验与自适应参数优化的图像去雾方法
CN111292257B (zh) 一种基于Retinex的暗视觉环境下图像增强方法
CN102750674B (zh) 基于自适应容差的视频图像去雾方法
Zhou et al. Underwater image enhancement method with light scattering characteristics
CN103218778B (zh) 一种图像和视频的处理方法及装置
CN108876743A (zh) 一种图像快速去雾方法、系统、终端及存储介质
CN109118446B (zh) 一种水下图像复原及去噪方法
CN107798661B (zh) 一种自适应的图像增强方法
CN109949247A (zh) 一种基于yiq空间光学成像模型的梯度域自适应增益水下图像增强方法
CN110675340A (zh) 一种基于改进的非局部先验的单幅图像去雾方法及介质
CN108182671B (zh) 一种基于天空区域识别的单幅图像去雾方法
CN105447825B (zh) 图像去雾方法及其系统
CN110807742A (zh) 一种基于一体式网络的微光图像增强方法
CN112862721B (zh) 基于暗通道与Retinex的地下管道图像去雾方法
CN109961415A (zh) 一种基于hsi空间光学成像模型的自适应增益水下图像增强方法
CN108648160B (zh) 一种水下海参图像去雾增强方法及系统
CN114693548B (zh) 一种基于亮区域检测的暗通道去雾方法
CN107977941B (zh) 一种明亮区域色彩保真及对比度增强的图像去雾方法
Liang et al. Underwater image quality improvement via color, detail, and contrast restoration
CN115937019A (zh) 一种lsd二次分割和深度学习相结合的不均匀去雾方法
CN107292844B (zh) 全变差正则化变分随机共振自适应暗图像滤波增强方法

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