-
文中针对水雷目标图像中的非均匀强噪声提出了一种新的去噪方法。首先,对非均匀强噪声进行统计,并以尺度衡量其强度;其次,提出LEP-EPC滤波算法,使其比LEP具有更优秀的保边缘平滑能力;最后,依据强噪声尺度迭代地将LEP-EPC滤波算法运用到每次的去噪结果上去除非均匀强噪声。图1为去噪算法流程图。
-
水下环境的复杂性决定了强噪声的特殊性。一方面,强噪声空间分布非均匀。水雷成像一般是在30 m以下的水下环境中,由于水体对光照的吸收,水深大于30 m的水下几乎是漆黑一片,所以深水成像必须采用人工照明设备提供光源。人工光源属于汇聚光源[8],虽然一般会在光源前添加准直透镜,但是水体的散射过程会让准直光逐步向漫射光过度,使得光场呈喇叭状分布,出现中间光子密度大,往四周逐渐减弱的现象。受成像距离和视场角限制,在光场的覆盖区域,强噪声会呈现同样的喇叭状分布,而在光场覆盖范围以外的区域,图像会存在暗区,强噪声强度较小甚至没有。另一方面,强噪声大小非均匀。水下里含有的浮游物及有机质相当于一个个微小的多面体,当相干光照射到多面体上发生衍射,而多面体的每个面接收到入射光后会改变该部分光波的相位,使衍射光在空间上产生干涉现象。强噪声在成像上具体表现为大量无规则的颗粒状亮斑。
图2(a)是实地采集的自然水下图像,水下目标是模拟水雷,图2(b)是利用暗通道算法[5]去散射的结果,图中框选了模拟水雷目标区域和两处非均匀强噪声区域,其细节放大图见图2(c)和图2(d)。由图2可以看出,利用去散射方法的确缓解了水下后向散射对图像的干扰,但是同时也放大了强噪声对图像的影响,不利于探测水雷目标。
为了定量描述水下强噪声的非均匀特性,我们对模拟水雷目标图像进行强噪声统计。由于强噪声在图像上表现为无规则的亮斑,因此可以对这些亮斑进行连通域检测,并统计连通域外接矩形的长边长度
$ length $ ,以此来衡量强噪声的强度。为了定量描述强噪声,文中给非均匀强噪声定义不同的尺度$ l\left(l\geqslant 1\right) $ ,最大的强噪声尺度${l}_{{\rm{max}}}$ 取决于统计的连通域外接矩形的长边长度的最大值$ {length}_{{\rm{max}}} $ :$$\mathop {\arg \min }\limits_{{l_{\max }}} \left| {{2^{{l_{\max }} + 1}} - lengt{h_{\max }}} \right|$$ (1) 图3是模拟水雷目标图像的强噪声统计结果。由图3可以看出,
$ {length}_{{\rm {max}}}=58 $ ,因此模拟水雷目标图像中最大的强噪声尺度$ {l}_{{\rm {max}}}=5 $ ,即强噪声尺度$ 1\leqslant l\leqslant 5 $ 。统计得到最大的强噪声尺度后可以对强噪声进行尺度划分,若强噪声的$ length $ 满足$ {2}^{l}< length< $ $ {2}^{l+1} $ ,则该强噪声属于第$ l $ 尺度。 -
强噪声的非均匀特性决定了强噪声具有局部特征,因此文中采用局部平滑的方式去除水下非均匀强噪声。局部保边缘滤波(Local Edge-Preserving, LEP)是依据图像局部梯度特征来区分图像内容中的边缘与非边缘,建立局部线性模型,构造出具有边缘保留特性的权值,从而达到保边缘平滑的目的。运用局部保边缘滤波算法去噪,既能有效地去除多种类型的噪声,又能保留边缘信息[9]。算法假设输入图像为
$ I $ ,输出图像为$ S $ ,在局部窗口$ {\rm{\omega }} $ 中最小化公式(2):$$ \begin{split} &{\rm {min}}\sum\limits_{i \in \omega } {{{(I(i) - {a_\omega }I(i) - {b_\omega })}^2}} + \\ & \alpha {\left| {\nabla I(i)} \right|^{(2 - \beta )}}a_\omega ^2 \end{split}$$ (2) 式中:
$ i=\left(x,y\right) $ ,代表图像的空间坐标;$ {\rm{\alpha }} $ 、$ {\rm{\;\beta }} $ 为已知参数;$ {a}_{\omega } $ 、$ {b}_{\omega } $ 为待求参数。算法假设$ S $ 和$ I $ 在局部窗口$ {\rm{\omega }} $ 里存在公式(3)的线性关系,利用公式(2)求出$ {a}_{\omega } $ 、$ {b}_{\omega } $ 即可得到去噪结果:$$S\left( i \right) = {a_\omega }I\left( i \right) + {b_\omega },{\rm{ }}i \in \omega $$ (3) $${a_\omega } = \frac{{\sigma _\omega ^2\left( i \right)}}{{\sigma _\omega ^2\left( i \right) + \dfrac{1}{N}\alpha {{ \displaystyle\sum\limits_{i \in \omega } {\left| {\nabla I\left( i \right)} \right|} }^{2 - \beta }}}}$$ (4) $${b_\omega } = {\overline I _\omega }\left( i \right) - {a_\omega }{\overline I _\omega }\left( i \right)$$ (5) 式中:
$ {\sigma }_{\omega }^{2}\left(i\right) $ 是$ I $ 在局部窗口内的方差;$ {\bar{I}}_{\omega }\left(i\right) $ 是$ I $ 在局部窗口内的均值。公式(2)的惩罚项中仅仅利用梯度信息去控制图像各部分内容的平滑程度,然而梯度信息容易受噪声干扰,无法很好地表示图像边缘处的像素;由公式(4)易知待求参数$ {a}_{\omega } $ 恒小于1,结合$ \nabla S\left(i\right)={a}_{\omega }\nabla I\left(i\right) $ 可知去噪的同时必然平滑图像的边缘信息。因此,文中根据这两点,对LEP算法模型提出优化。 -
文中针对LEP算法的保边缘平滑特性进行优化,提出了一种基于边缘感知约束的局部保边缘(Local Edge Preserving based on Edge Perception Constraint, LEP-EPC)滤波。为了减小去噪对边缘信息的影响,首先需要更好地区分边缘与噪声。文中采用文献[10]中提出的差分特征值作为边缘感知测度。差分特征值的计算是基于图像的Hessian矩阵,Hessian矩阵是通过计算图像的二阶导数得到。研究表明,相比于一阶导数(梯度),二阶导数能够更好地表征图像的边缘及细节[11]。Hessian矩阵定义为:
$$H{\rm{ = }}\left[ {\begin{array}{*{20}{c}} {{I_{xx}}}&{{I_{xy}}} \\ {{I_{xy}}}&{{I_{yy}}} \end{array}} \right]$$ (6) Hessian矩阵的特征值
$ {\lambda }_{1} $ 和$ {\lambda }_{2} $ 由公式(7)和(8)可以得到:$${\lambda _1} = \frac{1}{2}\left[({I_{x{\rm{x}}}} + {I_{yy}}) + \sqrt {{{({I_{xx}} - {I_{yy}})}^2} + 4{I_{xy}}^2} \right]$$ (7) $${\lambda _2} = \frac{1}{2}\left[({I_{x{\rm{x}}}} + {I_{yy}}) - \sqrt {{{({I_{xx}} - {I_{yy}})}^2} + 4{I_{xy}}^2} \right]$$ (8) 由此,差分特征值
$ D\left(i\right) $ 定义为:$$D\left( i \right) = \left( {{\lambda _1} - {\lambda _2}} \right){\lambda _1}w\left( i \right)$$ (9) 其中
$ w\left(i\right) $ 是用来平衡边缘增强以及噪声抑制的自适应权重,定义为:$$w(i) = \dfrac{{\sigma _\omega ^2(i) - {\rm {min}}(\sigma _\omega ^2)}}{{{\rm {max}}(\sigma _\omega ^2) -{\rm {min}}(\sigma _\omega ^2)}}$$ (10) 式中:
$ {\rm {max}}\left({\sigma }_{\omega }^{2}\right) $ 和$ {\rm {min}}\left({\sigma }_{\omega }^{2}\right) $ 是局部方差$ {\sigma }_{\omega }^{2}\left(i\right) $ 的最大值与最小值。利用
$ D\left(i\right) $ 能够较好地区分目标边缘与噪声:当像素点位于目标边缘处,$ {\lambda }_{1} $ 较大,$ {\lambda }_{2} $ 较小,$ \left({\lambda }_{1}-{\lambda }_{2}\right){\lambda }_{1} $ 较大,因此$ D\left(i\right) $ 也较大;当像素点位于平坦区域,$ \left({\lambda }_{1}-{\lambda }_{2}\right) $ 较小,$ \left({\lambda }_{1}-{\lambda }_{2}\right){\lambda }_{1} $ 较小,因此$ D\left(i\right) $ 也较小。图4(b)是图4(a)与标准差为30的高斯白噪声仿真而成的含噪图像,图4(c)和图4(d)分别是图4(b)的梯度图与差分特征值。分析图4可知,梯度图并没有将噪声和边缘区分开,而差分特征值则能够更好地表示图像边缘信息,具有更强的抗噪声干扰能力。利用
$ D\left(i\right) $ 可以构建出一个空间自适应的正则化系数$ 1/\left(1+\beta D\left(i\right)\right) $ ,那么新的目标函数为:$${\rm {min}}\sum\limits_{i \in \omega } {{{(I(i) - {a_\omega }I(i) - {b_\omega })}^2}} + \frac{{R(i)}}{{1 + \beta D(i)}}$$ (11) 显然,
$ 1/\left(1+\beta D\left(i\right)\right) $ 的取值依赖于$ D\left(i\right) $ ,根据像素点的不同位置给正则化项$ R\left(i\right) $ 施加不同的惩罚权重,从而控制去噪图像各部分内容的平滑程度。当像素点位于图像平坦区域时,$ D\left(i\right) $ 的值较小,$ 1/\left(1+\beta D\left(i\right)\right) $ 会给平坦区域的像素点施加一个较大的惩罚,从而平滑噪声;而当像素点位于边缘时,$ D\left(i\right) $ 的值较大,会给边缘处的像素值施加较小的惩罚,从而保留边缘和细节。为了防止最终的去噪结果中出现块状伪影以及增加算法的鲁棒性,文中对Hessian矩阵进行高斯滤波,即公式(6)修正为:$$H{\rm{ = }}\left[ {\begin{array}{*{20}{c}} {{I_{xx}} * G}&{{I_{xy}} * G} \\ {{I_{xy}} * G}&{{I_{yy}} * G} \end{array}} \right]$$ (12) 其中,
$ G $ 代表高斯核。其次,由公式(3)易知
$ \nabla S\left(i\right)={a}_{\omega }\nabla I\left(i\right) $ ,算法保边缘特性的好坏取决于$ {a}_{\omega } $ 的值。理论上,当某个像素点位于输入图像的边缘处时,$ {a}_{\omega } $ 的值为1,则边缘可以完全保留;当某个像素点位于输入图像的平坦区域时,$ {a}_{\omega } $ 的值为0,则该区域可以得到很好的平滑。若$ R\left(i\right)= $ $ {a}_{\omega }^{2} $ ,则最终求得的$ {a}_{\omega } $ 依然恒小于1。为了让求解出的$ {a}_{\omega } $ 的值更接近于理论值,$ R\left(i\right)={\left({a}_{\omega }-{\gamma }_{\omega }\right)}^{2} $ ,其中$ {\gamma }_{\omega } $ 是边缘感知约束项,定义为:$${\gamma _\omega }{\rm{ = }}1 - \frac{1}{{1 + {\rm{e}^{{\eta _\omega }\left( {{\varphi _\omega } - {\mu _\varphi }} \right)}}}}$$ (13) 式中:
$ {\eta }_{\omega }=\dfrac{n}{\left({\mu _\varphi }-{\rm {min}}\left({\varphi _\omega }\right)\right)} $ ,$ n $ 可以是大于1的任意整数,文中取$ n=2 $ ;${\varphi _\omega }$ 定义为两个局部窗口$ \omega $ 和$ \omega ' $ 内的标准差的乘积$ {\sigma }_{\omega }{\sigma }_{\omega '} $ ;${\mu _\varphi }$ 是${\varphi _\omega }$ 的均值。图5是根据图4(b)的有噪图像得到的边缘感知约束项
$ {\gamma }_{\omega } $ 的灰度分布图,分析图5易知,当图像的像素点位于目标图像边缘处时,$ {\gamma }_{\omega } $ 接近于1;当输入图像的像素点位于平坦区域时,$ {\gamma }_{\omega } $ 接近于0。图 5 边缘感知约束项
$ {\gamma }_{\omega } $ 的灰度分布图Figure 5. Gray distribution of edge perception constraint term
$ {\gamma }_{\omega } $ 综合以上分析,优化后的去噪目标函数可以写成公式(14)的形式:
$$ \begin{split} &{\rm {min}}\sum\limits_{i \in \omega } {{{(I(i) - {a_\omega }I(i) - {b_\omega })}^2}} +\\ & \dfrac{1}{{1 + \beta D(i)}}{\left( {{a_\omega } - {\gamma _\omega }} \right)^2} \\ \end{split} $$ (14) 根据模型(14),可以求得参数如下:
$${a_\omega } = \dfrac{{\sigma _\omega ^2 + \dfrac{1}{{1 + \beta D(i)}}{\gamma _\omega }}}{{\sigma _\omega ^2 + \dfrac{1}{{1 + \beta D(i)}}}}$$ (15) $${b_\omega } = {\overline I _\omega }\left( i \right) - {a_\omega }{\overline I _\omega }\left( i \right)$$ (16) 相比于原始LEP算法模型,LEP-EPC的算法模型求出的参数更接近于理论值。因此,在保边缘平滑效果上LEP-EPC算法比LEP算法更好。
-
给定一幅含有强噪声的图像
$ I $ ,通过单次滤波去噪可以将它分解成两层,一层是去噪结果$ S $ ,另一层是分离出的强噪声$ N $ 。假设把单次滤波去噪抽象为函数$ f $ ,则$$f\left( I \right) = S + N$$ (17) 但是由于噪声是非均匀的,单次滤波只可以分离得到某一个尺度的噪声,去噪不够彻底,且滤除强噪声的尺度取决于局部滤波窗口的大小,因此固定的滤波窗口无法很好地解决强噪声非均匀的问题。为了去除非均匀强噪声,需要迭代地将
$ f $ 应用到每一次的去噪结果中,在迭代过程中依据强噪声尺度逐步增大去噪尺度$ l $ ,从而生成多尺度分解。迭代过程可以用公式(18)表示:$$\left\{ {\begin{array}{*{20}{c}} {{S_l} = {f_l}\left( {{S_{l - 1}}} \right)} \\ {{N_l} = {S_{l - 1}} - {S_l}} \end{array}{\rm{ }}l = 1,2,3 \cdots ,{{m }}} \right.$$ (18) 式中:
$ {S}_{0} $ 初始化为输入图像I,多次迭代后,可以得到每个尺度的去噪图像$ {S}_{1}\sim{S}_{m} $ 和每个尺度的噪声图像$ {N}_{1}\sim{N}_{m} $ ,$ {S}_{m} $ 是最终的去噪结果。 -
文中算法采用MATLAB编程实现,运行环境为Windows10,处理器为Intel(R)Core(TM)i7。文中在这一部分安排如下:首先利用仿真实验来比较LEP滤波算法和LEP-EPC滤波算法在保边缘平滑性能上的优劣;其次对含有真实强噪声的模拟水雷目标图像利用LEP-EPC滤波算法结合多尺度分解策略(下文简称LEP-EPCM算法)进行去噪,并与经典去噪算法在主观评价和客观评价两方面对比去噪结果。主观上,对比的经典算法包括BM3D算法[12]和梯度L0范数最小化算法[13](下文简称L0算法);客观上,文中利用等效视数(Equivalent Number of Looks, ENL)和边缘保持指数(Edge Preserved Index, EPI)分别衡量算法在均匀区域对强噪声的去除能力和在目标区域对边缘信息的保留能力,并且还统计了算法的耗时。最后通过对比去噪前后的图像分割测试结果,验证算法在目标探测方面的有效性。
-
首先利用一组仿真图像来比较LEP滤波算法与LEP-EPC滤波算法在保边缘平滑性能上的优劣,实验结果如图6所示。从结果可以看出,两种算法模型对噪声都有着一定的平滑效果。但是LEP算法的去噪结果中模拟水雷目标的边缘及图像细节信息均发生了丢失,而LEP-EPC算法对模拟水雷和绳索的边缘保持得更好。此外,文中利用峰值信噪比(PSNR)和结构相似度(SSIM)对算法进行定量评估。在实验中高斯白噪声的标准差分别设置为15、20、25、30、35和40,分别利用两种滤波算法对图像进行去噪。去噪结果的PSNR对比如图7所示,结果表明LEP-EPC算法对噪声有着更好的平滑效果,表1中的SSIM结果进一步说明,LEP-EPC算法在边缘保护方面表现良好。
图 7 两种滤波算法的去噪结果在不同噪声水平下的PSNR对比
Figure 7. PSNR comparison of the denoising results of the two filtering algorithms under different noise levels
表 1 不同噪声水平下图像去噪结果的SSIM对比
Table 1. SSIM comparison of image denoising results under different noise levels
-
为了证明文中算法相比于经典算法在去除水下强噪声方面的优势,图8给出了4组含有真实非均匀强噪声的模拟水雷目标图像经过不同算法处理后的结果。主观上看,BM3D算法对目标信息保留较好,但是去除大尺度的噪声效果不理想,且处理结果相较于原图存在着一定的模糊;L0算法有着较明显的去噪效果,但是在目标边界处及强噪声区域有着明显的阶梯效应,而且容易发生目标丢失;LEP-EPCM算法的处理结果中无明显伪影,目标保留较好且去噪较为彻底。说明相较于另外两种经典算法,LEP-EPCM在保边缘去噪的同时有着更好的视觉效果。
为了定量地评价LEP-EPCM与另外两种经典算法的性能,文中采用ENL、EPI以及耗时作为衡量算法性能的指标。ENL是用来衡量图像均匀区域的平滑度,其值越高,说明对均匀区域强噪声的平滑效率越高。定义为:
$$ENL = \frac{{{\mu ^2}}}{{{\sigma ^2}}}$$ (19) 式中:
$ \mu $ 和$ \sigma $ 分别代表均匀区域的像素均值和标准差。EPI反映的是去噪后保留图像边缘细节的程度,定义为:$$EPI = \dfrac{{\displaystyle\sum\nolimits_i {\displaystyle\sum\nolimits_j {\left| {S\left( {i + 1,j} \right) - S\left( {i,j} \right)} \right|} } }}{{\displaystyle\sum\nolimits_i {\displaystyle\sum\nolimits_j {\left| {I\left( {i + 1,j} \right) - I\left( {i,j} \right)} \right|} } }}$$ (20) 式中:
$ S $ 是去噪后的图像;$ I $ 是含噪声图像。EPI越接近1,说明去噪后的图像边缘细节保留越好。文中分别计算了图8中框选的目标区域的EPI值以及强噪声区域的ENL值,对比结果见表2。通过对比表中数据可知,LEP-EPCM算法去噪结果的ENL和EPI值均是最高,说明保边缘去噪效果更好。在耗时方面,由表3数据中不难看出BM3D算法耗时最长,LEP-EPCM算法耗时低于L0算法的耗时,去噪更加高效。表 2 不同算法去噪结果的质量评估
Table 2. Quality assessments of denoising results of different algorithms
Algorithms Mine1 Mine2 Mine3 Mine4 ENL EPI ENL EPI ENL EPI ENL EPI Original image 33.67 1.0 39.10 1.0 31.65 1.0 8.33 1.0 BM3D 48.34 0.84 48.36 0.55 37.43 0.80 6.99 0.52 L0 76.22 0.82 73.75 0.62 61.16 0.62 3.09 0.32 LEP-EPCM 422.75 0.98 307.88 0.64 230.75 0.86 232.16 0.60 表 3 不同算法去噪结果的耗时对比
Table 3. Time cost comparison of denoising results of different algorithms
Algorithms Mine1 Mine2 Mine3 Mine4 BM3D 12 s 11.9 s 11.7 s 11.1 s L0 3.25 s 3.27 s 3.05 s 3.07 s LEP-EPCM 3.06 s 3.08 s 2.95 s 3.01 s -
文中还对模拟水雷目标图像进行了图像分割测试,并对比了不同算法去噪前后的图像分割结果。从图9可以直观地看出,图像去噪之前,目标被淹没在强噪声中,不便于后续的目标探测;通过BM3D算法和L0算法处理后,依然没有将目标从强噪声背景中分割出来,而LEP-EPCM算法将水雷目标从强噪声中分离了出来,强噪声被大量去除,大大减少了强噪声对水雷目标探测的干扰,为水雷目标探测提供了较为清晰的几何特征和候选区域,提高了目标探测的精确度。
Non-uniform strong noise removal method for non-cooperative mine target image
-
摘要: 水雷目标探测会受到水下非均匀强噪声(有机质、悬浮颗粒等)的干扰,为了解决这一问题,提出了一种新的去噪方法。首先优化局部保边缘滤波算法,提出了基于边缘感知约束的局部保边缘滤波,在模型中引入了一个空间自适应的边缘感知约束正则化项,用来更好地表征图像的边缘及细节,使得算法的保边缘平滑特性更好。其次针对强噪声的非均匀特性,采用多尺度策略,迭代地将优化后的模型运用到每个尺度的去噪结果上生成多尺度分解,并在多尺度分解的过程中,逐步增加去噪尺度,将不同尺度的噪声逐步从上一尺度的去噪结果中分离出来。实验结果表明,相较于其他经典去噪方法,提出的算法能够在更好地去除水下非均匀强噪声的同时保留水雷目标信息,对实时水雷作业有着一定的指导意义。Abstract: The detection of mine targets will be interfered by the underwater non-uniform strong noise (organic matter, suspended particles, etc). To solve this problem, a novel denoising method was proposed. Firstly, the local edge preserving filtering algorithm was optimized and the local edge preserving filtering based on edge perception constraint was proposed. A spatially adaptive edge perception constraint regularization term was introduced into the model to better represent the edges and details of the image, so that the edge-preserving and smoothing property could be better. Secondly, the multi-scale strategy was used to solve the heterogeneity of strong noise, the optimized model was iteratively applied to the noise removal results of each scale to generate multi-scale decomposition, and the denoising scale was gradually increaseed in the process of multi-scale decomposition. The noise of different scales was gradually separated from the denoising results of the previous scale. The experimental results show that, compared with other classical denoising methods, the proposed algorithm can better remove the underwater non-uniform strong noise while retaining the mine target information, which also has a certain guiding significance for real-time mine operation.
-
表 1 不同噪声水平下图像去噪结果的SSIM对比
Table 1. SSIM comparison of image denoising results under different noise levels
表 2 不同算法去噪结果的质量评估
Table 2. Quality assessments of denoising results of different algorithms
Algorithms Mine1 Mine2 Mine3 Mine4 ENL EPI ENL EPI ENL EPI ENL EPI Original image 33.67 1.0 39.10 1.0 31.65 1.0 8.33 1.0 BM3D 48.34 0.84 48.36 0.55 37.43 0.80 6.99 0.52 L0 76.22 0.82 73.75 0.62 61.16 0.62 3.09 0.32 LEP-EPCM 422.75 0.98 307.88 0.64 230.75 0.86 232.16 0.60 表 3 不同算法去噪结果的耗时对比
Table 3. Time cost comparison of denoising results of different algorithms
Algorithms Mine1 Mine2 Mine3 Mine4 BM3D 12 s 11.9 s 11.7 s 11.1 s L0 3.25 s 3.27 s 3.05 s 3.07 s LEP-EPCM 3.06 s 3.08 s 2.95 s 3.01 s -
[1] Liang Tianquan, Zhang Xiaoyun, Duan Peng, et al. Underwater target detection under strong scattering medium using improved dark channel method [J]. Infrared and Laser Engineering, 2020, 49(2): 0203012. (in Chinese) doi: 10.3788/IRLA202049.0203012 [2] Hu Haofeng, Li Xiaobo, Liu Tiegen. Recent advances in underwater image restoration technique based on polarimetric imaging [J]. Infrared and Laser Engineering, 2019, 48(6): 0603006. (in Chinese) doi: 10.3788/IRLA201948.0603006 [3] Wang Yibin, Yin Shibai, Lv Zhuowen. Adaptive background light estimation and nonlocal prior underwater image restoration [J]. Optics and Precision Engineering, 2019, 27(2): 499-510. (in Chinese) doi: 10.3788/OPE.20192702.0499 [4] Fan Zhiguo, Song Qiang, Dai Qingqing, et al. Polarization restoration of underwater target based on global parameter estimation [J]. Optics and Precision Engineering, 2018, 26(7): 1621-1632. (in Chinese) doi: 10.3788/OPE.20182607.1621 [5] Kaiming He, Jian Sun, Xiaoou Tang, et al. Single image haze removal using dark channel prior [J]. IEEE Transactions on Pattern Analysis & Machine Intelligence, 2011, 33(12): 2341-2353. [6] Komuro T, Chen K, Enomoto K, et al. Tracking and removal of suspended matter from underwater video images[C]// Thirteenth International Conference on Quality Control by Artificial Vision, 2017. [7] Gu B, Li W, Zhu M, et al. Local edge-preserving multiscale decomposition for high dynamic range image tone mapping [J]. IEEE Transactions on Image Processing, 2013, 22(1): 70-79. doi: 10.1109/TIP.2012.2214047 [8] Quan Xiangqian, Chen Xiangzi, Quan Yongqian, et al. Analysis and development of deep sea optical lighting and imaging system [J]. Chinese Optics, 2018, 11(2): 153-165. (in Chinese) doi: 10.3788/co.20181102.0153 [9] Wang Li. Research on edge preserving image filtering and segmentation under high intensity noise[D]. Nanjing: Nanjing Normal University, 2016 .(in Chinese) [10] Tian Haiying, Cai Hongmin, Lai Jianhuang, et al. Effective image noise removal based on difference eigenvalue[C]//Image Processing (ICIP), 2011 18th IEEE International Conference on, 2011. [11] Tang C, Han L, Ren H, et al. Second-order oriented partial-differential equations for denoising in electronic-speckle-pattern interferometry fringes [J]. Optics Letters, 2008, 33(19): 2179-2181. doi: 10.1364/OL.33.002179 [12] Dabov K, Foi A, Katkovnik V, et al. Image denoising by sparse 3-D transform-domain collaborative filtering [J]. IEEE Transactions on Image Processing, 2007, 16(8): 2080-2095. doi: 10.1109/TIP.2007.901238 [13] Xu L, Lu C, Xu Y, et al. Image smoothing via L0 gradient minimization [J]. ACM Transactions on Graphics, 2012, 30(6): 2024208.