CN110706165B - 一种基于EME与Mask的水下图像匀光算法 - Google Patents
一种基于EME与Mask的水下图像匀光算法 Download PDFInfo
- Publication number
- CN110706165B CN110706165B CN201910863140.4A CN201910863140A CN110706165B CN 110706165 B CN110706165 B CN 110706165B CN 201910863140 A CN201910863140 A CN 201910863140A CN 110706165 B CN110706165 B CN 110706165B
- Authority
- CN
- China
- Prior art keywords
- eme
- image
- value
- underwater
- values
- 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
- 230000011218 segmentation Effects 0.000 claims abstract description 13
- 230000000694 effects Effects 0.000 claims abstract description 7
- 238000000034 method Methods 0.000 claims description 13
- 238000004364 calculation method Methods 0.000 claims description 7
- 238000003709 image segmentation Methods 0.000 claims description 6
- 238000012545 processing Methods 0.000 claims description 6
- 230000008859 change Effects 0.000 claims description 4
- 238000001914 filtration Methods 0.000 claims description 4
- 230000000903 blocking effect Effects 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 3
- 238000005286 illumination Methods 0.000 abstract description 8
- 230000009286 beneficial effect Effects 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 4
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 4
- 230000008569 process Effects 0.000 description 3
- 238000011161 development Methods 0.000 description 2
- 230000018109 developmental process Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012546 transfer Methods 0.000 description 2
- 230000000007 visual effect Effects 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 239000000126 substance Substances 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/90—Dynamic range modification of images or parts thereof
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/50—Image enhancement or restoration using two or more images, e.g. averaging or subtraction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/11—Region-based segmentation
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Image Analysis (AREA)
- Image Processing (AREA)
Abstract
本发明公开了一种基于EME与Mask的水下图像匀光算法,对原始图像F(x,y)进行图像块分割处理,将以不同步长分割后的图像的EME值进行加权平均,原始图像F(x,y)与背景图像G(x,y)进行相减运算,以得到亮度分布均匀的理想图像I(x,y),将图像I(x,y)进行拉伸处理,进一步改善输出图像的效果。本发明的有益效果是可以用于光照不均匀的水下图像的匀光处理,具有自适应性和准确性的特点。
Description
技术领域
本发明属于水下图像处理技术领域,涉及一种基于EME与Mask的水下图像匀光算法。
背景技术
随着社会经济的发展,人类对于海洋资源的探测与开发的需求在逐步上升。对于水下目标物的识别、分类、分割以及目标识别等任务也有了更高的需求。输入图像的质量严重影响水下任务的准确率,所以获取高质量的水下图像数据对于提升水下任务的准确率非常重要。由于水下图像成像条件复杂,水体中悬浮颗粒对光线的吸收及散射作用以及不同波长的光在水下的衰减程度不同,水下图像通常存在细节模糊、对比度低、颜色失真、光照不均匀等问题。为了改善水下图像质量、提升水下任务的准确率,我们需要对原始的水下图像进行相应的处理。匀光算法能够有效的改善水下光学图像光照不均匀的问题,有重要的实际应用价值。目前,现存的匀光算法有Mask匀光算法、Retinex匀光算法、Wallis匀光算法等。其中,Mask匀光算法由于模型简单,结果视觉效果好而受到广泛应用。但传统的Mask匀光算法存在着一些缺陷,即无法自适应的调整滤波器参数。通过高斯低通滤波器获取背景图像是Mask匀光算法中的重要一步,高斯低通滤波器参数的选择尤为重要。但传统的Mask匀光算法需要人工选择高斯低通滤波器的参数,这个过程有较强的主观性,会对图像的匀光结果造成一定的影响。
发明内容
本发明的目的在于提供一种基于EME与Mask的水下图像匀光算法,本发明的有益效果是可以用于光照不均匀的水下图像的匀光处理,提升水下图像的视觉效果。依据EME值与滤波系数的正比关系,将EME值作为Mask匀光算法中高斯低通滤波器的滤波系数,避免了人工选择参数的主观因素,又引入了图像分块的方法以及加权平均的思想,避免了计算EME过程的偶然性,具有自适应性和准确性的特点。
本发明所采用的技术方案是按照以下步骤进行:
(1)对原始图像F(x,y)进行分块处理:将水下原始RGB图像分别分割为x×x大小相等的图像块,其中x为分割步长,进行步长分别为L、M、N的三次图像块分割,将分割后的三类图像再分别进行后续计算;
(2)对不同步长分割后的图像分别计算其EME值:EME值全称为图像增强测度(Measure of Image Enhancement),其描述了图像灰度的平均动态范围,其定义为:
其中,Imax为图像块最大灰度值,Imin为图像块的最小灰度值,k1、k2为指定步长将原图像分为块k1×k2,c为任意常数以防止分母为0。计算得出三种不同步长的分割方式的EME值,分别为EME1、EME2、EME3;
(3)将不同步长分割后的图像的EME值进行加权平均:加权平均计算的引入能够减小EME值得误差。首先分别计算出以三种步长分割方式的EME值的权重:
再计算出以不同步长分割后的图像EME值的加权平均数:
EME=EME1×w1+EME2×w2+EME3×w3
上式结果即为最终得到的EME值;
(4)加权平均后的EME值作为高斯低通滤波器的截止频率,将原始图像F(x,y)进行滤波以得到背景图像G(x,y):EME值体现了图像灰度差变化情况,其与滤波系数成正比关系。对不均匀光照的水下图像进行数学建模如下:
F(x,y)=I(x,y)+G(x,y)
其中I(x,y)为亮度分布均匀的理想图像,G(x,y)为反映噪声的背景图像。将原始水下RGB图像F(x,y)经过傅里叶变换后,输入到截止频率为EME值的高斯低通滤波器进行滤波,以得到背景图像G(x,y);
(5)原始图像F(x,y)与背景图像G(x,y)进行相减运算:
I(x,y)=F(x,y)-G(x,y)+Imean
其中Imean为代表像素值的偏移量,以调节图像的灰度值分布在0-255之间,且使结果与原始图像整体亮度一致,一般取Imean为原始图像各通道灰度值平均值的均值,即为:
其中IR、IG、IB分别为R、G、B三个通道的灰度平均值;
(6)将图像I(x,y)进行拉伸处理,进一步改善输出图像的效果:经过相减运算的图像灰度值会在动态范围内缩小,使得图像反差变小。将图像的灰度值两端进行拉伸:
其中,β为拉伸系数,取值范围为(-127,127),处理后Z(x,y)的值如果小于0,则Z(x,y)的值为0;如果大于255,则Z(x,y)的值为255。Z(x,y)为EME-Mask匀光算法得到的最终结果;
进一步,步骤(1)采用三种不同步长的图像分割方式,且步长为L、M、N,其中L、M、N之间存在约束关系M≥2L,N≥2M。
进一步,步骤(2)在匀光算法中引入EME值,并且针对三种分割方式分别计算其EME值。
进一步,步骤(3)引入加权平均的思想,将三种分割方式所得到的三个EME值进行加权平均的计算,分别算出三种分割方式EME值的权重再计算出这三种分割方式EME值的加权平均数EME=EME1×w1+EME2×w2+EME3×w3,将此结果作为最终的EME值,避免偶然性,减小了EME值的误差。
进一步,步骤(5)将加权平均后的EME值作为传统Mask匀光算法中高斯低通滤波器的截止频率。
附图说明
图1是传统Mask匀光算法流程示意图;
图2是本发明图像分割方法示意图;
图3是本发明图像分割误差说明示意图;
图4是本发明EME-Mask匀光算法流程示意图;
图5是实例中针对水下光照不均匀图像的EME-Mask匀光算法的实施结果图;
图6是本发明与传统算法处理之后PSNR值的对比结果图。
具体实施方式
下面结合具体实施方式对本发明进行详细说明。
图1为传统Mask匀光算法流程示意图;本发明实施例如下:
第一步:本发明实施例采用Bumblebee XB3立体摄像头进行图像数据采集,该摄像头出厂时已针对镜头扭曲变形和摄像头失准问题进行校准,从而确保所有摄像头校准的一致性,再将该摄像头水密封舱,适用于采集水下高质量图像。将目标物放在暗室的消声水池内,采用水密封舱的摄像头,通过提供外部光源来采集光照不均匀的水下原始RGB图像。本实施例截取摄像头拍摄图像中光照最不均匀的一张为F(x,y),其分辨率为1280×960;
第二步:将水下原始RGB彩色图像F(x,y)分割步长分别为L、M、N的像素块,其中M≥2L,N≥2M,本实施例取L=8、M=16、M=32;
第三步:根据Agaian基于韦伯定理提出的图像增强测度,即EME值,对三种不同步长分割后的图像分别计算其EME值:
其中Imax为图像块最大灰度值,Imin为图像块的最小灰度值,k1、k2为指定步长将原图像分为块k1×k2,c为任意常数防止分母为0。计算得出三种分割方式的EME值分别为EME1=0.0015、EME2=0.0025、EME3=0.0040;
第四步:为了避免偶然性,减小误差,引入加权平均的思想。首先分别计算出三种分割方式EME值的权重:
再根据加权平均的公式计算出这三种分割方式EME值的加权平均数:
EME=EME1×w1+EME2×w2+EME3×w3=0.0031;
第五步:将原始图像F(x,y)经过傅里叶变换,变换到频域:
B(x,y)=F{F(x,y)}
高斯低通滤波器的传递函数为:
其中,u是距离傅里叶变换原点的距离,v是高斯曲线的膨胀程度,EME值则为高斯低通滤波器的截止频率。由于高斯低通滤波器的传递函数经过傅里叶变换是其本身,则将图像输入到截止频率为EME值的高斯低通滤波器进行滤波:
L(x,y)=B(x,y)*H(u,v)
再将L(x,y)进行傅里叶逆变换,变回到空间域:
G(x,y)=F-1{L(x,y)}
滤除高频部分后得到的低频部分即为背景图像G(x,y);
第六步:将原始图像与反映噪声的背景图像作相减运算即可得到亮度分布均匀的理想图像I(x,y):
I(x,y)=F(x,y)-G(x,y)+Imean
即可得到亮度分布均匀的理想图像I(x,y),其中Imean为代表偏移量,调节图像的灰度值分布在0-255之间。为了使结果与原始图像亮度一致,取Imean为原始图像灰度值的均值:
其中IR、IG、IB分别为R、G、B三个通道的灰度平均值;
第七步:为了使经过相减运算的灰度值缩小的图像变回原始平均灰度值,将图像的灰度值两端进行拉伸:
其中β为拉伸系数,取值范围为(-127,127),处理后Z(x,y)的值如果小于0,则Z(x,y)的值为0;如果大于255,则Z(x,y)的值为255。Z(x,y)为EME-Mask匀光算法得到的最终结果;
第八步:比较基于EME与Mask的水下图像匀光算法匀光算法与传统Mask匀光算法的PSNR值。PSNR值为峰值信噪比,是一种评价图像处理前后质量好坏的参数,它的计算公式为:
其中MSE为均方误差,是一种常见的损失函数。
PSNR值越大,则其图像经过处理后的失真越小,图像质量越高。计算基于EME与Mask的水下图像匀光算法的PSNR值:
计算传统Mask匀光算法的PSNR值:
通过实验对比,经过基于EME与Mask的水下图像匀光算法匀光算法的水下RGB图像的PSNR值达到19.2265,经过传统Mask匀光算法的水下图像PSNR值为16.6251,两者差值2.6014。证明了基于EME与Mask的水下图像匀光算法对于提升图像质量的效果,对比了基于EME与Mask的水下图像匀光算法匀光算法的性能相较于传统匀光算法更佳。
图2是本发明图像切割方法示意图;图3是本发明图像分割误差说明示意图;图4是本发明EME-Mask匀光算法流程示意图;图5是实例中针对水下光照不均匀图像的EME-Mask匀光算法的实施结果图;图5(a)为原图像,图5(b)为传统算法匀光后图像,图5(c)为EME-Mask算法匀光后图像。图6是本发明与传统算法处理之后PSNR值的对比结果图。
本发明的优点还在于:
(1)本发明具有较强的自适应性,对于不同的水下图像无需手动选择高斯滤波器的参数,通过引入EME值以及图像切割、加权平均的思想,得到了效果良好、误差较小的高斯滤波器参数,相较于传统算法更加智能。
(2)本发明具有较强的适用性,不仅可以用于水下图像的匀光,通过实验表明也适用于遥感图像的匀光且效果良好。
以上所述仅是对本发明的较佳实施方式而已,并非对本发明作任何形式上的限制,凡是依据本发明的技术实质对以上实施方式所做的任何简单修改,等同变化与修饰,均属于本发明技术方案的范围内。
Claims (5)
1.一种基于EME与Mask的水下图像匀光方法,其特征在于按照以下步骤进行:
(1)对原始图像F(x,y)进行图像分块处理:将水下原始RGB图像分别分割为x×x的大小相等的图像块,其中x为分割步长;进行步长分别为L、M、N的三次图像分割,将分割后的三类图像再分别进行后续计算;
(2)对不同步长分割后的图像分别计算其EME值,EME值描述了图像灰度的平均动态范围,其定义为:
其中Imax为图像块最大灰度值,Imin为图像块的最小灰度值,k1、k2为指定步长将原图像分为块k1×k2,c为任意常数以防止分母为0,计算得出三种不同步长的分割方式的EME值,分别为EME1、EME2、EME3;
(3)将不同步长分割后的图像的EME值进行加权平均:加权平均计算的引入能够减小EME值的误差;首先分别计算出三种步长的分割方式EME值的权重:
再计算出以不同步长分割后的图像EME值的加权平均数:
EME=EME1×w1+EME2×w2+EME3×w3
上式结果即为最终得到的EME值;
(4)加权平均后的EME值作为高斯低通滤波器的截止频率,将原始图像F(x,y)进行滤波以得到背景图像G(x,y):EME值体现了图像的色差变化情况,其与滤波系数成正比关系;对不均匀光照的水下图像进行数学建模如下:
F(x,y)=I(x,y)+G(x,y)
其中I(x,y)为亮度分布均匀的理想图像,G(x,y)为反映噪声的背景图像,将原始水下RGB图像F(x,y)经过傅里叶变换后,输入到截止频率为EME值的高斯低通滤波器进行滤波,以得到背景图像G(x,y);
(5)将原始图像F(x,y)与背景图像G(x,y)进行相减运算:
I(x,y)=F(x,y)-G(x,y)+Imean
其中Imean为代表偏移量,调节图像的灰度值分布在0-255之间,为了使结果与原始图像亮度一致,取Imean为原始图像各通道灰度值平均值的均值,即为:
其中IR、IG、IB分别为R、G、B三个通道的灰度平均值;
(6)将图像I(x,y)进行拉伸处理,进一步改善输出图像的效果:经过相减运算的图像灰度值会在动态范围内缩小,使得图像反差变小,将图像的灰度值两端进行拉伸:
其中β为拉伸系数,取值范围为(-127,127),处理后Z(x,y)的值如果小于0,则Z(x,y)的值为0;如果大于255,则Z(x,y)的值为255;Z(x,y)为EME-Mask匀光算法得到的最终结果。
2.按照权利要求1所述一种基于EME与Mask的水下图像匀光方法,其特征在于:所述步骤(1)采用三种不同步长的图像分割方式,且步长为L、M、N,其中L、M、N之间存在约束关系M≥2L,N≥2M。
3.按照权利要求1所述一种基于EME与Mask的水下图像匀光方法,其特征在于:所述步骤(2)在匀光算法中引入EME值,并且针对三种分割方式分别计算其EME值。
4.按照权利要求1所述一种基于EME与Mask的水下图像匀光方法,其特征在于:所述步骤(3)将三种分割方式所得到的三个EME值进行加权平均计算,分别算出三种分割方式EME值的权重再计算出这三种分割方式EME值的加权平均数EME=EME1×w1+EME2×w2+EME3×w3,将此结果作为最终的EME值,避免参数计算的偶然性,减小了EME值的误差。
5.按照权利要求1所述一种基于EME与Mask的水下图像匀光方法,其特征在于:所述步骤(5)将加权平均后的EME值作为传统Mask匀光算法中高斯低通滤波器的截止频率。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910863140.4A CN110706165B (zh) | 2019-09-12 | 2019-09-12 | 一种基于EME与Mask的水下图像匀光算法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910863140.4A CN110706165B (zh) | 2019-09-12 | 2019-09-12 | 一种基于EME与Mask的水下图像匀光算法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110706165A CN110706165A (zh) | 2020-01-17 |
CN110706165B true CN110706165B (zh) | 2023-05-12 |
Family
ID=69194780
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910863140.4A Active CN110706165B (zh) | 2019-09-12 | 2019-09-12 | 一种基于EME与Mask的水下图像匀光算法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110706165B (zh) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106934775A (zh) * | 2017-03-08 | 2017-07-07 | 中国海洋大学 | 一种基于低秩恢复的非局部图像去噪方法 |
CN108053374A (zh) * | 2017-12-05 | 2018-05-18 | 天津大学 | 一种结合双边滤波与Retinex的水下图像增强方法 |
CN108932700A (zh) * | 2018-05-17 | 2018-12-04 | 常州工学院 | 基于目标成像模型的自适应梯度增益水下图像增强方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1345172A1 (en) * | 2002-02-26 | 2003-09-17 | Sony International (Europe) GmbH | Contrast enhancement for digital images |
CN102968769B (zh) * | 2012-11-27 | 2015-07-22 | 宁波盈芯信息科技有限公司 | 一种图像一致性增强装置 |
-
2019
- 2019-09-12 CN CN201910863140.4A patent/CN110706165B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106934775A (zh) * | 2017-03-08 | 2017-07-07 | 中国海洋大学 | 一种基于低秩恢复的非局部图像去噪方法 |
CN108053374A (zh) * | 2017-12-05 | 2018-05-18 | 天津大学 | 一种结合双边滤波与Retinex的水下图像增强方法 |
CN108932700A (zh) * | 2018-05-17 | 2018-12-04 | 常州工学院 | 基于目标成像模型的自适应梯度增益水下图像增强方法 |
Non-Patent Citations (2)
Title |
---|
shenghui rong.An improved non-uniformity correction algorithm and its GPU parallel implementation.《web of science》.2018,全文. * |
冉红格 ; 何波 ; 刘贵杰.自主式水下机器人中数字电子罗盘的误差分析及校正.中国海洋大学学报(自然科学版).2009,(第0S1期),全文. * |
Also Published As
Publication number | Publication date |
---|---|
CN110706165A (zh) | 2020-01-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111986120A (zh) | 一种基于帧累加和多尺度Retinex的低光照图像增强优化方法 | |
US8319843B2 (en) | Image processing apparatus and method for blur correction | |
CN105046663B (zh) | 一种模拟人类视觉感知的自适应低照度图像增强方法 | |
CN103973990B (zh) | 宽动态融合方法及装置 | |
CN112734650A (zh) | 一种基于虚拟多曝光融合的不均匀光照图像增强方法 | |
Xiong et al. | Application of histogram equalization for image enhancement in corrosion areas | |
CN110163807B (zh) | 一种基于期望亮通道的低照度图像增强方法 | |
Chen et al. | Weighted sparse representation and gradient domain guided filter pyramid image fusion based on low-light-level dual-channel camera | |
CN117274113B (zh) | 基于图像增强的碎硅片清洗效果视觉检测方法 | |
CN115272303A (zh) | 基于高斯模糊的纺织物缺陷程度评估方法、装置及*** | |
Feng et al. | Research on defect detection method for high-reflective-metal surface based on high dynamic range imaging | |
Mu et al. | Low and non-uniform illumination color image enhancement using weighted guided image filtering | |
Cao et al. | Vignetting correction based on a two-dimensional Gaussian filter with harmony for area array sensors | |
CN107478656B (zh) | 基于机器视觉的纸浆搅拌效果检测评价方法、装置、*** | |
CN110706165B (zh) | 一种基于EME与Mask的水下图像匀光算法 | |
CN115249211A (zh) | 一种基于水下非均匀入射光模型的图像复原方法 | |
Xie et al. | Image defogging method combining light field depth estimation and dark channel | |
Sur et al. | Measuring the noise of digital imaging sensors by stacking raw images affected by vibrations and illumination flickering | |
CN112825189A (zh) | 一种图像去雾方法及相关设备 | |
CN108364274B (zh) | 微纳米尺度下光学图像的无损清晰重构方法 | |
CN115829967A (zh) | 一种工业金属表面缺陷图像去噪和增强方法 | |
CN115937024A (zh) | 基于Retinex理论的多帧融合低照度图像增强方法 | |
CN115293989A (zh) | 一种综合MsRcR和自动色阶的图像增强方法 | |
CN111612720B (zh) | 一种广角红外图像优化方法、***及相关组件 | |
CN109544481B (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 |