CN107203980B - 自适应多尺度暗通道先验的水下目标探测图像增强方法 - Google Patents
自适应多尺度暗通道先验的水下目标探测图像增强方法 Download PDFInfo
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 47
- 238000001514 detection method Methods 0.000 title claims abstract description 23
- 238000001914 filtration Methods 0.000 claims abstract description 34
- 230000005540 biological transmission Effects 0.000 claims abstract description 23
- 238000012545 processing Methods 0.000 claims abstract description 11
- 238000005286 illumination Methods 0.000 claims abstract description 10
- 238000012935 Averaging Methods 0.000 claims abstract description 3
- 238000004364 calculation method Methods 0.000 claims description 19
- 238000010586 diagram Methods 0.000 claims description 9
- 230000003044 adaptive effect Effects 0.000 claims description 8
- 230000002708 enhancing effect Effects 0.000 claims description 8
- 238000011158 quantitative evaluation Methods 0.000 claims description 4
- 239000011159 matrix material Substances 0.000 claims description 3
- 238000010606 normalization Methods 0.000 claims description 3
- 230000000007 visual effect Effects 0.000 abstract description 5
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 3
- ZXQYGBMAQZUVMI-QQDHXZELSA-N [cyano-(3-phenoxyphenyl)methyl] (1r,3r)-3-[(z)-2-chloro-3,3,3-trifluoroprop-1-enyl]-2,2-dimethylcyclopropane-1-carboxylate Chemical compound CC1(C)[C@@H](\C=C(/Cl)C(F)(F)F)[C@H]1C(=O)OC(C#N)C1=CC=CC(OC=2C=CC=CC=2)=C1 ZXQYGBMAQZUVMI-QQDHXZELSA-N 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000011084 recovery Methods 0.000 description 2
- 239000003086 colorant Substances 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000010365 information processing Effects 0.000 description 1
- 230000002401 inhibitory effect Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000013519 translation Methods 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
- 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/90—Dynamic range modification of images or parts thereof
- G06T5/92—Dynamic range modification of images or parts thereof based on global image properties
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
-
- 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
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)的具体计算过程为:
第四步:计算暗通道图像DC0(I)的m个不同尺度窗口Ωk对应的透射图Tk(I),Tk(I)的具体计算过程如下:
a)设定不同尺度窗口Ωk:Nk×Nk(k=1,2,…,m);
第五步:按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∑w∑h|Qk(I)-Tk(I)|,式中w和h分别是I的宽度和高度像素值;
第八步:应用暗通道先验理论,计算图像I的恢复图像J,具体计算过程如下:
第九步:对恢复图像J从均值、方差、对比度和信息熵等方面进行定量评价,相关定量评价指标函数表示为:
本发明所达到的有益之处在于:
本发明的方法,可以只利用单幅非均匀亮度、低信噪比、低对比度水下目标探测图像自身的信息,对图像进行自适应多尺度暗通道先验的去噪增强处理。先通过投射图误差最小缩小尺度参数数量范围,自适应确定尺度空间,再进一步通过均值、方差、对比度和信息熵等综合定量评价指标评估暗通道先验的增强图像。本发明中用到的导向滤波函数,能在图像的纹理与平滑之间实现均衡处理,使得处理后的图像视觉质量提高、纹理信息丰富。
附图说明
图1是本发明的方法的控制流程图。
图2是自适应筛选尺度窗口控制流程图。
具体实施方式
下面结合附图对本发明作进一步描述。以下实施例仅用于更加清楚地说明本发明的技术方案,而不能以此来限制本发明的保护范围。
参照图1所示,本发明是一种自适应多尺度暗通道先验的水下目标探测图像增强方法,整体流程图如图1所示,具体实现步骤如下:
第一步:获取水下目标原始彩色图像I;
第二步:计算图像I的光照向量A,A的具体计算过程为:I对应的灰度图像Gray(I)中最大像素值,即为光照向量A;
第三步:按He方法计算图像I的暗通道图像DC0(I),DC0(I)的具体计算过程为:
第四步:计算暗通道图像DC0(I)的m个不同尺度窗口Ωk对应的透射图Tk(I),Tk(I)的具体计算过程如下:
d)设定不同尺度窗口Ωk:Nk×Nk(k=1,2,…,m);
第五步:按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∑w∑h|Qk(I)-Tk(I)|,式中w和h分别是I的宽度和高度像素值;
第八步:应用暗通道先验理论,计算图像I的恢复图像J,具体计算过程如下:
第九步:对恢复图像J从均值、方差、对比度和信息熵等方面进行定量评价,相关定量评价指标函数表示为:
对自适应多尺度参数(或尺度窗口)图像增强有关问题的说明:
(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的最小值滤波输出
步骤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:根据导向滤波输出图像与透射图之间误差最小条件,初选出几个最合适的最小滤波尺度窗口以及对应的导向滤波输出图像;
min∑w∑h|Qk(I)-Tk(I)|,
式中,w和h分别是I的宽度和高度;
步骤7:对筛选的导向滤波输出图像加权平均;
步骤8:应用暗通道先验理论,计算原始图像的恢复增强图像;
步骤9:对恢复图像从包括均值、方差、对比度和信息熵的方面进行定量评价。
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)
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)
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)
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 |
-
2017
- 2017-05-31 CN CN201710405123.7A patent/CN107203980B/zh active Active
Patent Citations (2)
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)
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 |