1.School of Physics and Optoelectronic Engineering, Xidian University, Xi’an 710071, China 2.Xi’an Key Laboratory of Computational Imaging, Xidian University, Xi’an 710071, China 3.Key Laboratory of Optical Engineering, Institute of Optics and Electronics, Chinese Academic of Science, Chengdu 610209, China 4.Academic of Advanced Interdisciplinary Research, Xidian University, Xi’an 710071, China
Fund Project:Project supported by the National Natural Science Foundation of China (Grant Nos. 62075175, 62005203), and the Foundation of the Key Laboratory of Optical Engineering, Chinese Academic of Sciences (Grant No. QC20191097)
Received Date:10 February 2021
Accepted Date:17 April 2021
Available Online:07 June 2021
Published Online:20 August 2021
Abstract:Underwater imaging plays a critical role in marine rescue, seabed resource exploration, underwater archaeology, etc. by providing human-vision-system-friendly information. A variety of approaches have been exploited to realize clear underwater imaging. Noticeably, underwater polarization imaging has attracted attention due to its simple imaging system and clear vision. It can remove the backscattered light from degraded image and recover abundant high-fidelity information of target. Descattering is conducted by using the difference in polarization characteristics between the target and background. A classical underwater polarization imaging method is presented by Schechner [Tali T, Schechner Y Y 2009 IEEE Trans. Pattern Anal. Mach. Intell.31 385], in which the differential polarization characteristics of backscattered light and target light are used to recover clear image. More researches were conducted including Huang et al.’s research [Huang B J, Liu T G, Hu H F, et al. 2016 Optics Express24 9826], Liu et al.’s study [Liu F, Han P L, Wei Y, et al. 2018 Opt. Lett.43 4903], etc.However, in the polarization imaging methods, the uniform underwater backscattered light and polarization parameters over the whole image are usually assumed. In most practical applications, these assumptions cannot hold true. Therefore, the inaccurate estimation of backscattered light makes it difficult to completely descatter an image, leading many methods to fail to detect the target in non-uniform turbid water.In this study, we propose a low-rank-and-sparse-decomposition-based polarization imaging combined with common mode rejection feature of polarization information in scattered light field to eliminate non-uniformity and scattering caused by severe scattering during active polarization imaging of turbid water. The backscattered light is highly reduced and the information contained in background is single and highly correlated. It conforms to the low-rank characteristics of the image. What is more, the target in underwater scene occupies a relatively small proportion, which conforms to the sparsity characteristics of the image. Therefore, combining the low-rank characteristics of backscattered light with the sparse characteristics of target information light, we separate them through low-rank and sparse matrix decomposition to recover clear underwater image. Both experimental and objective image quality evaluation results demonstrate the validity of the proposed method.The proposed method works well in improving polarization vision in non-uniform turbid water, which is due to its ability to make the underwater scene uniform and the target and background information separated through their distribution difference of polarization characteristics. It possesses potential applications in turbid water imaging. Keywords:underwater imaging/ polarization characteristic/ common-mode rejection/ scattering
全文HTML
--> --> -->
2.浑浊水体散射光场强度分布特性研究光在水下传播过程中的散射情况如图1所示, 随着传播距离增加, 散射程度不断增强. 通常可根据水中散射程度的差异性将成像区域分为弹道光区域、散射增强区域和随机游走区域; 弹道光区域距离短, 直接成像即可获得较好观测效果; 随机游走区域内光波散射作用强, 故无法直接获取任何目标信息[26]. 在如图1(a)所示的散射增强区内, 随着传输距离或水体浑浊度的增加, 光波的散射程度急剧增强, 散射光场中背景散射光强度分布${I_{{\rm{scat}}}}(x, y)$如下式所示[22]: 图 1 浑浊水体散射光场特性分析 Figure1. Scattering conditions of light at different turbidity.
$\begin{split}{I_{{\rm{scat}}}}\left( {x,y} \right) =\;& \int_{{R_{{\rm{cam}} = 0}}}^{{R_{{\rm{cam}}}}({x_{{\rm{obj}}}})} b\left[ {\theta \left( z \right)} \right]{I_{{\rm{source}}}}\left( z \right)\\&\times\exp \left[ { - c{R_{{\rm{cam}}}}\left( z \right)} \right]d{R_{{\rm{cam}}}}, \end{split}$