【技术实现步骤摘要】
本专利技术涉及光学成像，尤其涉及一种基于穆勒矩阵极分解的水下目标清晰化成像方法。
技术介绍
1、水下成像环境复杂，特别是在浑浊水域中，受到水体吸收和大量微粒的散射作用，目标细节信息损失严重，导致水下成像质量急剧下降，严重影响后续对水下目标的识别能力。因此，如何提高水下光学成像质量，成为浑浊水体下清晰化成像的迫切需求。
2、为了提高水下光学成像的质量，学者们提出了许多研究方法，包括：
3、1、基于图像处理的图像增强方法。该方法采用数字图像处理算法实现，虽然简单高效，但是该类算法以数据统计为基础，往往忽略掉图像质量降低的原因，处理后的图像极易导致感兴趣的细节部分模糊。
4、2、基于深度学习的偏振成像技术利用偏振信息特征提取实现清晰成像，将神经网络与偏振信息相结合，是偏振水下成像技术的新探索，此方法往往需要大量的样本数据以优化模型参数，泛化能力不强，并不适用于陌生的场景。
5、3、基于偏振成像技术的去散射方法。偏振成像技术提升了目标信息维度，是目前实现水下清晰化成像的有效方法之一。采用偏振成像去散射的方法一般有两类，第一类是将偏振度引入到水下光学成像物理模型，其核心思想是求解目标与背景的偏振度，并利用两者之间的偏振度差异实现清晰化成像，但该方法不可避免地存在偏振度的计算难题；第二类是偏振成像结合深度学习实现去散射，利用增加的偏振维度信息实现清晰成像，是偏振水下成像技术的新探索，但是该方法往往需要大量的样本数据以优化模型参数，泛化能力不强，而且很难适用于陌生场景。
6、因此，现
技术实现思路
1、本专利技术的目的在于提供一种基于穆勒矩阵极分解的水下目标清晰化成像方法，通过穆勒矩阵极分解的方式获取退偏矩阵图像，以图像质量为目标优化水下光学成像物理模型中透射率图的估计，利用此模型反演出目标的清晰图像，无需估计目标和场景的偏振度，解决现有计算方法存在偏振度的计算难、样本量要求高、参数计算难等问题。
2、为达成上述目的，本专利技术提供一种基于穆勒矩阵极分解的水下目标清晰化成像方法，包括以下步骤：
3、步骤1、搭建穆勒矩阵水下光学成像物理模型，选取穆勒矩阵的第一个元素为总光强信号，进而根据目标场景选择无穷远处的散射光强；将穆勒矩阵极分解得到退偏图像，使用退偏图像初步估计光学成像模型的透射率图；
4、步骤2、使用一修正系数对初步估计的透射率图进行修正，以优化图像质量；
5、步骤3、选择增强度量eme作为图像质量的评价指标，使用粒子群算法遍寻修正系数，直至eme为最大值，选取eme为最大值时的修正系数对初步估计的透射率图进行修正，获得清晰图像的透射率图；
6、步骤4、根据获得的清晰图像透射率图、总光强信号和无穷远处的散射光强计算得到水下目标清晰图像。
7、进一步，在步骤4中，计算水下目标清晰图像的方法如下：
8、根据现有水下光学成像原理，水体的吸收和微粒的散射作用，探测器接受到的信号包括三个部分：
9、i(x)＝ia(x)+ib(x)+if(x)#(1)
10、式中，i(x)为探测器像素x处收到的总光强信号；ia(x)为衰减的目标信号光；ib(x)为后向散射光，是图像退化的主要原因；if(x)为前向散射光，忽略不计，ia(x)用目标的清晰图像和整个场景的透射率图表达，ib(x)用整个场景的透射率图和无穷远处的散射光强表达，由此，式(1)改写为：
11、i(x)＝iobj(x)t(x)+(1-t(x))#(2)
12、式中，iobj(x)为目标的清晰图像；t(x)为整个场景的透射率图，取值范围为0-1；b∞为无穷远处的散射光强，不同的取值会影响成像效果，b∞选择背景区域的均值，或者选择i(x)最亮的0.1％像素的均值，根据式(2)，反演出目标的清晰图像：
13、
14、进一步，在步骤1中，根据穆勒矩阵极分解的原理，一个归一化的4×4穆勒矩阵表示为三个子矩阵的乘积：
15、m＝mδmrmd#(4)
16、式中，mδ、mr、md分别代表退偏、相位延迟和二向衰减特性，穆勒矩阵中的第一个元素为总光强信号i(x)，b∞选择i(x)最亮的0.1％像素的均值，退偏矩阵mδ的表达式为：
17、
18、式中，为3×1的偏振向量，mδ是表征场景退偏信息的关键参数，表达式为：
19、
20、将mδ的元素用于估计光学成像物理模型中的透射率图t(x)′：
21、
22、式中，
23、
24、式中，将mδ平方处理是以增强对比度且不损失任何信息；将mδ最大值处理以使初步估计的透射率图t(x)′的取值范围为0-1。
25、进一步，在步骤2中，使用一修正系数μ来对t(x)′来进行修正，以优化图像质量，得到t(x)，则式(3)改写为：
26、
27、式中，μ在一个场景下为固定常数，取值范围为0-1/max(t(x))。
28、进一步，在步骤3中，eme的表达式为：
29、
30、式中，iobj(x)图像被划分为m×n块，smax；m,n和smin；m,n表示每块图像中像素灰度值的最大值与最小值，同时为了避免被零划分，参数q取接近于0的正值，设为0.0001；使用粒子群算法寻找μ，将寻找的μ值代入式(9)获得eme值，判断获得的eme值是否为最大值，若为最大值，则代入的μ值为最终选取的修正系数，若获得的eme值不是最大值，则继续寻找μ值代入式(9)中重复计算，直至获得eme的最大值。
31、进一步，在步骤1中，所述穆勒矩阵水下光学成像物理模型采用一套双dofp型穆勒矩阵成像装置搭建，所述双dofp型穆勒矩阵成像装置包括偏振态产生器、偏振态分析器、led点光源、准直透镜、带通滤光片、扩束镜，所述偏振态产生器包括一个固定角度的偏振片p和一个可旋转的1/4波片r1，所述led点光源、准直透镜、带通滤光片、p、r1、扩束镜依次排列设置在一条直线上以照射目标场景，所述偏振态分析器包括两个线偏振相机dofp1、dofp2、一个50:50的非偏振分光棱镜和一个固定角度的1/4波片r2，所述dofp1和dofp2分别固定在分光棱镜垂直的两侧，r2固定在dofp1和分光棱镜之间，以p的透光轴固定为水平方向，作为成像物理模型的基准轴，所述dofp1的相机主轴与基准轴之间具有夹角；
32、led点光源经过准直透镜和带通滤光片后形成准直单色光源，准直单色光源通过偏振态产生器和扩束镜照射目标场景，后经偏振态分析器检测，获得场景的穆勒矩阵，以斯托克斯矢量描述偏振光束，穆勒矩阵描述光学元件和被测场景，两个相机采集的8×4强度矩阵i表示为：
33、i＝apsa·m·apsg#(11)
【技术保护点】
1.一种基于穆勒矩阵极分解的水下目标清晰化成像方法，其特征在于：包括以下步骤：
2.如权利要求1所述的一种基于穆勒矩阵极分解的水下目标清晰化成像方法，其特征在于：在步骤4中，计算水下目标清晰图像的方法如下：
3.如权利要求2所述的一种基于穆勒矩阵极分解的水下目标清晰化成像方法，其特征在于：在步骤1中，根据穆勒矩阵极分解的原理，一个归一化的4×4穆勒矩阵表示为三个子矩阵的乘积：
4.如权利要求3所述的一种基于穆勒矩阵极分解的水下目标清晰化成像方法，其特征在于：在步骤2中，使用一修正系数μ来对t(x)′来进行修正，以优化图像质量，得到t(x)，则式(3)改写为：
5.如权利要求4所述的一种基于穆勒矩阵极分解的水下目标清晰化成像方法，其特征在于：在步骤3中，EME的表达式为：
6.如权利要求1-5任意一项所述的一种基于穆勒矩阵极分解的水下目标清晰化成像方法，其特征在于：在步骤1中，所述穆勒矩阵水下光学成像物理模型采用一套双DoFP型穆勒矩阵成像装置搭建，所述双DoFP型穆勒矩阵成像装置包括偏振态产生器、偏振态分析器、LED点光
7.如权利要求6所述的一种基于穆勒矩阵极分解的水下目标清晰化成像方法，其特征在于：所述DoFP1的相机主轴与基准轴之间的夹角为21°。
...【技术特征摘要】
1.一种基于穆勒矩阵极分解的水下目标清晰化成像方法，其特征在于：包括以下步骤：
2.如权利要求1所述的一种基于穆勒矩阵极分解的水下目标清晰化成像方法，其特征在于：在步骤4中，计算水下目标清晰图像的方法如下：
3.如权利要求2所述的一种基于穆勒矩阵极分解的水下目标清晰化成像方法，其特征在于：在步骤1中，根据穆勒矩阵极分解的原理，一个归一化的4×4穆勒矩阵表示为三个子矩阵的乘积：
4.如权利要求3所述的一种基于穆勒矩阵极分解的水下目标清晰化成像方法，其特征在于：在步骤2中，使用一修正系数μ来对t(x)′来进行修正，以优化图像质量，得到t(x)，则式(3)改写为：
5.如权利要求4所述的一种基于穆勒矩阵极分解的水下目标清晰化成像方法，其特征在于：在步骤3中，eme的表达式为：
6.如权利要求1-5任意一项所述的一种基于穆勒矩阵极分解的水下目标清晰化成像方法，其特征...
【专利技术属性】
技术研发人员：李慧慧，俞万能，杨静雷，廖卫强，陈垚，郭必成，
申请(专利权)人：集美大学，
类型：发明
国别省市：
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。