本发明涉及偏振图像融合技术，具体涉及一种基于非下采样剪切波变换的多偏振信息选通融合方法。解决了现有图像融合时偏振角图像参与度不高无法保留有效信息以及偏振角图像噪声滤除较差，造成融合图像信噪比较低的技术问题。本发明方法包括以下步骤：S1)利用非下采样剪切波变换对可见光强度源图像、偏振度源图像、偏振角源图像进行分解得到低频分量A0，B0，C0及高频方向分量S2)利用加权平均法对源图像的低频分量进行融合得到低频融合分量；S3)对源图像的高频方向分量进行选通，取邻域能量最大的像素点灰度值得到融合图像的高频方向分量；S4)对低频融合分量和融合图像高频方向分量进行非下采样剪切波变换的逆变换，得到融合图像。
1.一种基于非下采样剪切波变换的多偏振信息选通融合方法，其特征在于，包括以下步骤：S1)利用非下采样剪切波变换对可见光强度源图像、偏振度源图像、偏振角源图像进行多级多方向分解，得到可见光强度源图像、偏振度源图像、偏振角源图像的低频分量A0，B0，C0及高频方向分量 k为非下采样剪切波变换的分解级数，l为各分解级下的方向级数；
S2)利用加权求和方法对步骤S1)得到的可见光强度源图像、偏振度源图像、偏振角源图像的低频分量灰度值进行融合得到低频融合分量R0；
λ1，λ2，λ3分别为可见光强度源图像、偏振度源图像、偏振角源图像的低频分量在融合图像中的权重；
所述λ1，λ2，λ3均不为0，且λ1+λ2+λ3＝1，λ2＝λ3，λ1为0.5～0.6；
S3)对步骤S1)得到的可见光强度源图像、偏振度源图像、偏振角源图像的高频方向分量中各像素点计算邻域窗口ω1的邻域方差大小，并根据邻域方差大小进行选通；若邻域方差大于阈值，则对应像素点不参与后续融合；若邻域方差小于阈值，则对相应的像素点再取邻域窗口ω2，并计算该邻域窗口ω2的能量，选取可见光强度源图像、偏振度源图像、偏振角源图像中邻域窗口ω2能量最大的像素点灰度值作为融合图像的相应高频方向分量像素点灰度值；各像素点灰度值确定后，便得到融合图像高频方向分量分别为可见光强度源图像高频方向分量、偏振度源图像高频方向分量、偏振角源图像高频方向分量；
为经过高频融合 算法选通后得到的融合图像的高频方向分量；
所述阈值为图像整体邻域方差平均值的1.2～1.5倍；
S4)对步骤S2)得到的低频融合分量R0及步骤S3)得到的融合图像的高频方向分量 进行非下采样剪切波变换的逆变换，得到最终的融合图像R。
2.根据权利要求1所述的一种基于非下采样剪切波变换的多偏振信息选通融合方法，其特征在于：步骤S3)中，方差Var计算公式为：
g(i+p,j+q)为图像(i+p,j+q)处的灰度值；
3.根据权利要求2所述的一种基于非下采样剪切波变换的多偏振信息选通融合方法，其特征在于：步骤S1)中，k＝4；l在各分解级上分别为8，8，16，16；
步骤S2)中，λ1，λ2，λ3分别取0.5，0.25，0.25。
4.根据权利要求3所述的一种基于非下采样剪切波变换的多偏振信息选通融合方法，其特征在于：步骤S3)中，所述邻域窗口ω1为21×21；
所述邻域窗口ω2为3×3，邻域窗口各像素权重(ω(p,q))为1/16*[1 2 1；2 4 2；1 2
一种基于非下采样剪切波变换的多偏振信息选通融合方法
技术领域
[0001] 本发明涉及偏振图像融合技术，具体涉及一种基于非下采样剪切波变换的多偏振信息选通融合方法。
背景技术
[0002] 物体表面的反射光和散射光会包含其自身特性的偏振信息，通过这些偏振信息，可以反应被测物体的轮廓、纹理、粗糙度等物理信息。偏振度图像及偏振角图像分别根据目标表面偏振度及偏振角进行成像，可实现对目标边缘特征的增强，并对耀光、过曝等影响普通可见光强度图像的因素有较强的抑制作用，有着独特的目标探测能力。
[0003] 为实现偏振角图像、偏振度图像、可见光强度图像信息的互补以及后续数据利用的便利，需要将上述三种图像进行融合为一。但偏振角图像、偏振度图像中噪声影响较大，尤其偏振角图像，噪声影响较为显著，使得融合过程中想要同时实现滤除噪声和保留源图
像信息比较困难。
[0004] 现有偏振图像融合算法可分为两类：第一类是只利用可见光强度图像和偏振度图像进行融合，噪声影响较大的偏振角图像不参与融合。这类融合方法得到的融合图像信噪
比较高，但也同样丢失了偏振角图像中包含的有效信息。第二类是将偏振角图像、偏振度图像、可见光强度图像三者进行融合。这类方法可实现三种源图像信息的融合，但在噪声抑制及信息保留两种目标之间未实现有效的均衡。要么未对偏振角图像中较强的噪声进行有效
处理，造成融合图像信噪比较低，要么抑制了噪声，同时保留的有效信息也较少。
发明内容
[0005] 本发明目的是解决现有偏振角图像、偏振度图像、可见光强度图像三者进行融合时，在噪声抑制及信息保留两种目标之间未实现有效的均衡，造成融合图像质量较差的技
术问题，而提供一种基于非下采样剪切波变换的多偏振信息选通融合方法。该方法利用非
下采样剪切波变换进行滤波分解，通过邻域方差大小对源图像窗口邻域进行选通，判定是
否噪声较大，是否参与融合，可有效抑制源图像中的噪声并保留其有效信息。该方法可对图像各分量进行精细融合调控，实现了多偏振信息的高质量融合。
[0006] 本发明技术方案为：
[0007] 一种基于非下采样剪切波变换的多偏振信息选通融合方法，其特殊之处在于，包括以下步骤：
[0008] S1)利用非下采样剪切波变换对可见光强度源图像、偏振度源图像、偏振角源图像进行多级多方向分解，得到可见光强度源图像、偏振度源图像、偏振角源图像的低频分量
A0，B0，C0及高频方向分量 k为非下采样剪切波变换的分解级数，l为各分解
级下的方向级数；
[0009] 所述k为3～4；
[0010] 所述l为4的倍数；
[0011] S2)利用加权求和方法对步骤S1)得到的可见光强度源图像、偏振度源图像、偏振角源图像的低频分量灰度值进行融合得到低频融合分量R0；
[0012] R0＝λ1A0+λ2B0+λ3C0
[0013] λ1，λ2，λ3分别为可见光强度源图像、偏振度源图像、偏振角源图像的低频分量在融合图像中的权重；
[0014] 所述λ1，λ2，λ3均不为0，且λ1+λ2+λ3＝1，λ2＝λ3，λ1为0.5～0.6；
[0015] S3)对步骤S1)得到的可见光强度源图像、偏振度源图像、偏振角源图像的高频方向分量中各像素点计算邻域窗口ω1的邻域方差大小，并根据邻域方差大小进行选通；若邻域方差大于阈值，则对应像素点不参与后续融合；若邻域方差小于阈值，则对相应的像素点再取邻域窗口ω2，并计算该邻域窗口ω2的能量，选取可见光强度源图像、偏振度源图像、偏振角源图像中邻域窗口ω2能量最大的像素点灰度值作为融合图像的相应高频方向分量像
素点灰度值；各像素点灰度值确定后，便得到融合图像高频方向分量
[0016]
[0017] 分别为可见光强度源图像高频方向分量、偏振度源图像高频方向分量、偏振角源图像高频方向分量；
[0018] F{}为高频融合算法；
[0019] 为经过高频算法选通后得到的融合图像的高频方向分量；
[0020] 所述邻域窗口ω1为21～31个像素；
[0021] 所述邻域窗口ω2为3～5个像素；
[0022] 所述阈值为图像整体邻域方差平均值的1.2～1.5倍；
[0023] S4)对步骤S2)得到的低频融合分量R0及步骤S3)得到的融合图像的高频方向分量进行非下采样剪切波变换的逆变换，得到最终的融合图像R。
[0024] 进一步地，步骤S3)中，方差Var计算公式为：
[0025]
[0026] MV为邻域灰度均值；
[0027] (i,j)为邻域坐标；
[0028] g(i,j)为图像(i,j)处的灰度值；
[0029] M×N为窗口大小；
[0030] 邻域能量Ar计算公式为：
[0031]
[0032] ω(p,q)为邻域窗口各像素权重；
[0033] g(i+p,j+q)为图像(i+p,j+q)处的灰度值；
[0034] 窗口横向宽度为2m+1，纵向宽度为2n+1。
[0035] 进一步地，步骤S1)中，k＝4；l在各分解级上分别为8，8，16，16；
[0036] 步骤S2)中，λ1，λ2，λ3分别取0.5，0.25，0.25。
[0037] 进一步地，步骤S3)中，所述邻域窗口ω1为21×21；
[0038] 所述阈值为图像整体邻域方差平均值的1.2倍；
[0039] 所述邻域窗口ω2为3×3，邻域窗口各像素权重(ω(p,q))为1/16*[1 2 1；2 4 2；1 2 1]。
[0040] 本发明的有益效果：
[0041] 1、该方法利用非下采样剪切波变换进行滤波分解，通过邻域方差大小对源图像高频方向分量窗口邻域进行选通，判定是否噪声较大，是否参与融合，可有效抑制源图像中的噪声并保留其有效信息。
[0042] 2、该方法可对各源图像高频方向分量进行精细融合调控，实现了多偏振信息的高质量融合。
[0043] 3、该方法能够较好的保留三种源图像的重要特征信息实现了融合图像的高质量表达。
附图说明
[0044] 图1为本发明实施例的偏振度源图像；
[0045] 图2为本发明实施例的偏振角源图像；
[0046] 图3为本发明实施例的可见光强度源图像；
[0047] 图4为本发明实施例的融合图像。
具体实施方式
[0048] 一种基于非下采样剪切波变换的多偏振信息选通融合方法，包括以下步骤：
[0049] S1)利用非下采样剪切波变换对可见光强度源图像、偏振度源图像、偏振角源图像进行多级多方向分解，得到可见光强度源图像、偏振度源图像、偏振角源图像的低频分量
A0，B0，C0及高频方向分量 k为非下采样剪切波变换的分解级数，l为各分解
级下的方向级数。分解级数太少，融合不够精细；分解级数太多，计算量太大，因此k一般为3～4。同样分解方向太少，不利用融合分析；分解方向太多，计算量太大，l为4的倍数，一般取
8或16。本实施例中分解级数k＝4，在各分解级上的方向级数l分别为8，8，16，16。可见光强度源图像、偏振度源图像、偏振角源图像均被分解为1幅低频分量图像及48幅高频方向分量图像。
[0050] S2)利用加权求和方法对步骤S1)得到的可见光强度源图像、偏振度源图像、偏振角源图像的低频分量灰度值进行融合得到低频融合分量R0；
[0051] R0＝λ1A0+λ2B0+λ3C0
[0052] λ1，λ2，λ3分别为可见光强度源图像、偏振度源图像、偏振角源图像的低频分量在融合图像中的权重；λ1，λ2，λ3均不为0，且λ1+λ2+λ3＝1，λ2＝λ3，λ1为0.5～0.6。本实施例λ1，λ2，λ3分别取0.5，0.25，0.25。
[0053] S3)先对步骤S1)得到的可见光强度源图像、偏振度源图像、偏振角源图像的高频方向分量中各像素点计算邻域窗口ω1的邻域方差大小，并根据邻域方差大小进行选通；若邻域方差大于阈值，则对应像素点不参与后续融合；若邻域方差小于阈值，则对相应的像素点再取邻域窗口ω2，并计算该邻域窗口ω2的能量，选取可见光强度源图像、偏振度源图像、偏振角源图像中邻域窗口ω2能量最大的像素点灰度值作为融合图像的相应高频方向分量
像素点灰度值；各像素点灰度值确定后，便得到融合图像高频方向分量
[0054]
[0055] 为各源图像高频方向分量；
[0056] F{}为高频融合算法；
[0057] 为经过高频算法选通后得到的融合图像高频方向分量。
[0058] 方差Var计算公式为：
[0059]
[0060] MV为邻域灰度均值；
[0061] (i,j)为邻域坐标；
[0062] g(i,j)为图像(i,j)处的灰度值；
[0063] M×N为窗口大小。
[0064] 邻域能量Ar计算公式为：
[0065]
[0066] ω(p,q)为邻域窗口各像素权重，权重大小一般以窗口中间像素为主，兼顾窗口边缘；
[0067] g(i+p,j+q)为图像(i+p,j+q)处的灰度值；
[0068] 窗口横向宽度为2m+1，纵向宽度为2n+1。
[0069] 阈值为图像整体邻域方差平均值的1.2～1.5倍，阈值较大，各源图像噪声部分易引入到融合图像中，阈值较小，各源图像的有效信息可能被滤除，本实施例的选通阈值设为图像整体邻域方差平均值的1.2倍。为区分像素点区域是否噪声较大，窗口ω1的取值一般
在21～31个像素之间，而窗口ω2一般取3～5个像素，其是在窗口ω1内且噪声较小的像素点区域部分，选取较小领域内质量最好的像素值进行融合。本实施例中邻域窗口ω1为21×
21，计算能量的邻域窗口ω2为3×3，m为1，n为1。邻域窗口ω2的各像素权重为1/16*[1 2 1；
2 4 2；1 2 1]。
[0070] S4)对步骤S2)得到的低频融合分量R0及步骤S3)得到的融合图像的高频方向分量进行非下采样剪切波变换的逆变换，得到最终的融合图像R。
[0071] 图1、图2和图3分别为室外同一远景的偏振度源图像、偏振角源图像和可见光强度源图像。如图1所示，偏振度源图像由于玻璃镜面的强偏振反射，偏振程度显著，造成图像局部亮度高，整体亮度较低，近处树木细节已无法分辨，但对楼房、窗口等目标的边缘反应较好。如图2所示，偏振角源图像噪声较大，树木处噪声尤其显著，但对远处电线及楼宇窗口可清晰感知。如图3所示，可见光强度源图像，远处天空出现过曝情况，电线处已经无法分辨，但近处树木边缘较清晰。
[0072] 如图4所示，经过本发明融合算法处理后得到融合图像，融合图像对偏振角源图像中的噪声进行了有效选通抑制，而且树木、楼房、窗口、远处电线等目标的有效信息经过融合均得到有效保留，证明本方法可有效实现对源图像的噪声滤除及信息保留。
[0073] 此外，本发明方法将可见光强度源图像、偏振度源图像、偏振角源图像进行了多层多方向分解，最终各源图像均被分解为1幅低频分量图像及48幅高频方向分量图像，实现了精细融合操作。
