CN104574296B - 一种用于偏振去雾霾的多小波融合处理图片的方法 - Google Patents
一种用于偏振去雾霾的多小波融合处理图片的方法 Download PDFInfo
- Publication number
- CN104574296B CN104574296B CN201410815711.4A CN201410815711A CN104574296B CN 104574296 B CN104574296 B CN 104574296B CN 201410815711 A CN201410815711 A CN 201410815711A CN 104574296 B CN104574296 B CN 104574296B
- Authority
- CN
- China
- Prior art keywords
- image
- sigma
- polarization
- frequency
- fusion
- 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.)
- Expired - Fee Related
Links
Landscapes
- Image Processing (AREA)
Abstract
本发明涉及一种用于偏振去雾霾的多小波融合处理图片的方法,属于图像处理技术领域。本发明中的去雾方法将偏振探测与多小波融合算法相结合,以目标偏振探测为基础,应用多小波所具有的多个尺度函数和多个小波函数,使其同时满足对称性、正交性和紧支集性等性能,根据分解后高、低频系数所具有的特点分别采用不同的规则对目标偏振度、偏振角等偏振信息及强度信息进行融合。使用本方法得到的偏振融合图像,能够同时综合目标识别的细节与亮度信息,在雾霾天气环境下得到更高识别率、更高质量、具有更高细节信息的目标图像,提高目标探测的去雾能力。
Description
技术领域
本发明涉及一种用于偏振去雾霾的多小波融合处理图片的方法,属于图像信息处理技术领域。
背景技术
偏振成像技术是探测成像领域中的又一全新学科,它将偏振测量与图像处理方法相结合,通过测量目标辐射和反射的偏振度和偏振角等信息,可以弥补传统光度学测量中所存在的缺陷,获取测量所需的高质量、高对比度的图像,在军民应用中都有着广泛的应用前景。
但是由于目前城市污染所导致雾霾天气日益严重的现象,使得环境中某些介质对光的吸收和散射产生衰减作用,则最终探测到的目标图像不清晰、对比度低,严重影响目标的可识别率。为了改善雾霾天气下图像的偏振成像质量,基于分别能反应图像边缘、图像细节及图像强度信息的偏振度、偏振角及强度图像的特点,采用对偏振图像融合的方式得到具有更多图像信息的探测结果,提高图像在复杂环境下的可识别率。
在常用的像素级、特征级和决策级三种融合方法中,像素级融合方法在对图像清晰度的改善方面有着很大的贡献,但对图像的融合多采用单分辨率方法,且以统一的规则进行融合,易产生边缘噪声和丢失信息的现象,应用小波变换的多分辨率融合方法可根据小波多尺度的特性,对图像的高频及低频信息分别应用不同的规则进行融合,能很好地改善单分辨率融合的不足,但这种传统单小波方法仅具有一个尺度和一个小波函数,对于具有图像边缘和细节特点的偏振度和偏振角图像,在分解时很难充分表达其细节信息。基于此许多研究者又将小波的分解算法与其他融合方法相结合,以获得更多的目标信息,但这种方法对图像限制较多,且完成步骤较繁琐。在单小波算法的基础上发展起来的多小波方法,不仅具有与小波方法相同的多分辨率分析特性,而且多小波所能同时具有的紧支撑、对称和正交性等特点,为后续的融合算法提供更多的分解信息。
发明内容
为了在雾霾天气环境下得到更高质量、更高识别率的探测结果,本发明提出了一种用于偏振去雾霾的多小波融合处理图片的方法。
本发明是这样设计的:一种用于偏振去雾霾的多小波融合处理图片的方法,其特征在于:包括下述步骤,
步骤1、对激光器进行定标,测出激光器发射光线的初始偏振态,用偏振片Ⅰ根据初始偏振态起偏,获得的偏振光线照射在雾霾模拟环境中的目标上;在接收端调节偏振片Ⅱ的角度分别为0°、60°和120°,用探测器进行探测,通过图像处理模块获得这三个角度的光强探测结果分别为I′(0°)、I′(60°)和I′(120°);
步骤2、将步骤1中探测到的结果通过公式(1)获得目标的三个斯托克斯线偏振分量
其中,I为总光强,Q为x分量和y分量的强度差,U为在Q的基础上+45°分量和-45°分量的强度差,
将目标的三个斯托克斯线偏振分量通过公式(2)获得偏振度Dolp和偏振角PA
步骤3、将步骤2中得到的偏振度Dolp作为第一图像和总光强作为第二图像分别用多小波融合处理法进行一阶分解,分解顺序为第一顺序,其包含行预处理、列预处理、行多小波变换以及列多小波变换,最终获得含有16个子块的分解结果,其中为低频子块,和为高频子块,i取1或2,表示第一图像或第二图像,其中L1为第一个尺度函数对应的低通滤波器、L2为第二个尺度函数对应的低通滤波器、H1为第一个尺度函数高通滤波器和H2为第二个尺度函数对应的高通滤波器;
步骤4、对第一图像和第二图像中的低频信息进行融合,用LpLq i表示低频的各个子块,p和q分别取1或2,首先通过公式(3)确定两图像的局部区域能量Ei(x,y)
其中,M和N都为3,窗口其次通过公式(4)获得图像间的局部区域能量匹配度M12
其中阈值T1=0.75,用来确定图像间的相关性,获得低频融合系数LpLq F,
当M12<T1时,采用局部区域能量取大方法,通过公式(5)获得融合后的低频系数
当M12≥T1时,通过公式(6)获得融合后的低频系数
步骤5、对步骤3中高频的三个部分12个子块进行融合,通过公式(7)获得图像空间频率SFi
其中i=1时表示第一图像,i=2时表示第二图像,M×N(3×3)窗口下HF为高频系数水平方向上的频率、VF为高频系数垂直方向上的频率和DF为高频系数对角线方向上的频率,通过公式(8)获得
然后在这一窗口下通过公式(9)获得图像间的相关系数R
其中u1为第一图像在此窗口下的像素均值,u2为第二图像在此窗口下的像素均值,
最后根据R值获得融合后高频系数H1L1 F,
当R≥0.8时,表示高度相关,融合后的高频系数为
当0.3<R<0.8时,表示中度相关,有
其中T1为第一图像的空间频率在整个频率中所占权重,T2为第二图像的空间频率在整个空间频率中所占权重,即且T2(x,y)+T1(x,y)=1,
当R≤0.3时,表示低度相关,
重复步骤5操作对12个子块进行融合;
步骤6、对以上步骤得到的融合结果进行一致性检测;
步骤7、对所获得的高频和低频系数按照第二顺序进行重构,获得融合后图像;
步骤8、将第一图像与第二图像的融合结果再按照步骤3~7与偏振角图像进行融合。
所述的接收端包括偏振片Ⅱ、滤光片Ⅱ及探测器,其中探测器为CCD照相机。
所述步骤3中的第一顺序与步骤7中的第二顺序为相反的顺序。
所述步骤6中所述的一致性检测:
a、将融合结果中某一像素点作为中心,在其3×3邻域中统计融合后系数来源于第一图像和第二图像的个数;
b、根据统计结果选择个数占优势的那个图像,并将此区域的中心点的像素值,以该占优势的那个图像相应的像素点的值进行定义。
附图说明
下面结合附图说明和具体实施方式对本发明作进一步说明:
图1是在雾霾模拟环境下目标偏振探测的实验原理图。
图2是多小波分解过程流程图。
如图所示,1为激光器、2为准直扩束***装置、3为滤光片Ⅰ、4为偏振片Ⅰ、5雾霾模拟环境中的目标、6为偏振片Ⅱ、7为滤光片Ⅱ、8为探测器、9为图像处理模块。
具体实施方式
一种采用偏振去雾霾多小波融合处理图片的方法,其所需条件和实现步骤如下:
所需条件:如图1所示,所需条件由4部分组成,
1、起偏部分:包括激光器、准直扩束***、滤光片Ⅰ及偏振片Ⅰ。
2、雾霾模拟环境。
3、检偏部分:包括偏振片Ⅱ、滤光片Ⅱ、CCD相机。
4、计算机图像处理部分为MATLAB软件。
具体实现步骤:
步骤1、对激光器(1)进行定标,测出激光器(1)发射光线的初始偏振态,用偏振片Ⅰ(4)根据初始偏振态起偏,获得的偏振光线照射在雾霾模拟环境中的目标(5)上;在接收端调节偏振片Ⅱ(6)的角度分别为0°、60°和120°,用探测器(8)进行探测,通过图像处理模块(9)获得这三个角度的光强探测结果分别为I′(0°)、I′(60°)和I′(120°);
步骤2、将步骤1中探测到的结果通过公式(1)获得目标的三个斯托克斯线偏振分量
其中,I为总光强,Q为x分量和y分量的强度差,U为在Q的基础上+45°分量和-45°分量的强度差,
将目标的三个斯托克斯线偏振分量通过公式(2)获得偏振度Dolp和偏振角PA
步骤3、将步骤2中得到的偏振度Dolp作为第一图像和总光强作为第二图像分别用多小波融合处理法进行一阶分解,分解顺序为第一顺序,其包含行预处理、列预处理、行多小波变换以及列多小波变换,最终获得含有16个子块的分解结果,其中为低频子块,和为高频子块,i取1或2,表示第一图像或第二图像,其中L1为第一个尺度函数对应的低通滤波器、L2为第二个尺度函数对应的低通滤波器、H1为第一个尺度函数高通滤波器和H2为第二个尺度函数对应的高通滤波器;
步骤4、对第一图像和第二图像中的低频信息进行融合,用LpLq i表示低频的各个子块,p和q分别取1或2,首先通过公式(3)确定两图像的局部区域能量Ei(x,y)
其中,M和N都为3,窗口其次通过公式(4)获得图像间的局部区域能量匹配度M12
其中阈值T1=0.75,用来确定图像间的相关性,获得低频融合系数LpLq F,
当M12<T1时,采用局部区域能量取大方法,通过公式(5)获得融合后的低频系数
当M12≥T1时,通过公式(6)获得融合后的低频系数
步骤5、对步骤3中高频的三个部分12个子块进行融合,通过公式(7)获得图像空间频率SFi
其中i=1时表示第一图像,i=2时表示第二图像,M×N(3×3)窗口下HF为高频系数水平方向上的频率、VF为高频系数垂直方向上的频率和DF为高频系数对角线方向上的频率,通过公式(8)获得
然后在这一窗口下通过公式(9)获得图像间的相关系数R
其中u1为第一图像在此窗口下的像素均值,u2为第二图像在此窗口下的像素均值,
最后根据R值获得融合后高频系数H1L1 F,
当R≥0.8时,表示高度相关,融合后的高频系数为
当0.3<R<0.8时,表示中度相关,有
其中T1为第一图像的空间频率在整个频率中所占权重,T2为第二图像的空间频率在整个空间频率中所占权重,即且T2(x,y)+T1(x,y)=1,
当R≤0.3时,表示低度相关,
重复步骤5操作对12个子块进行融合;
步骤6、对以上步骤得到的融合结果进行一致性检测;
步骤7、对所获得的高频和低频系数按照第二顺序进行重构,获得融合后图像;
步骤8、将第一图像与第二图像的融合结果再按照步骤3~7与偏振角图像进行融合。
所述的接收端包括偏振片Ⅱ6、滤光片Ⅱ7及探测器8,其中探测器8为CCD照相机。
所述步骤3中的第一顺序与步骤7中的第二顺序为相反的顺序。
所述步骤6中所述的一致性检测:
a、将融合结果中某一像素点作为中心,在其3×3邻域中统计融合后系数来源于第一图像和第二图像的个数;
b、根据统计结果选择个数占优势的那个图像,并将此区域的中心点的像素值,以该占优势的那个图像相应的像素点的值进行定义。
利用以上步骤实现了用于偏振去雾霾的多小波融合处理图片的方法。
Claims (3)
1.一种用于偏振去雾霾的多小波融合处理图片的方法,其特征在于:包括下述步骤,
步骤1、对激光器(1)进行定标,测出激光器(1)发射光线的初始偏振态,用偏振片Ⅰ(4)根据初始偏振态起偏,获得的偏振光线照射在雾霾模拟环境中的目标(5)上;在接收端调节偏振片Ⅱ(6)的角度分别为0°、60°和120°,用探测器(8)进行探测,通过图像处理模块(9)获得这三个角度的光强探测结果分别为I′(0°)、I′(60°)和I′(120°);
步骤2、将步骤1中探测到的结果通过公式(1)获得目标的三个斯托克斯线偏振分量
其中,I为总光强,Q为x分量和y分量的强度差,U为在Q的基础上+45°分量和-45°分量的强度差,
将目标的三个斯托克斯线偏振分量通过公式(2)获得偏振度Dolp和偏振角PA
步骤3、将步骤2中得到的偏振度Dolp作为第一图像和总光强作为第二图像分别用多小波融合处理法进行一阶分解,分解顺序为第一顺序,其包含行预处理、列预处理、行多小波变换以及列多小波变换,最终获得含有16个子块的分解结果,其中为低频子块,和为高频子块,i取1或2,表示第一图像或第二图像,其中L1为第一个尺度函数对应的低通滤波器、L2为第二个尺度函数对应的低通滤波器、H1为第一个尺度函数对应高通滤波器和H2为第二个尺度函数对应的高通滤波器;
步骤4、对第一图像和第二图像中的低频信息进行融合,用LpLq i表示低频的各个子块,p和q分别取1或2,首先通过公式(3)确定两图像的局部区域能量Ei(x,y)
其中,M和N都为3,窗口其次通过公式(4)获得图像间的局部区域能量匹配度M12
其中阈值T1=0.75,用来确定图像间的相关性,获得低频融合系数LpLq F,
当M12<T1时,采用局部区域能量取大方法,通过公式(5)获得融合后的低频系数
当M12≥T1时,通过公式(6)获得融合后的低频系数
步骤5、对步骤3中高频的三个部分12个子块进行融合,通过公式(7)获得图像空间频率SFi
其中i=1时表示第一图像,i=2时表示第二图像,M×N(3×3)窗口下HF为高频系数水平方向上的频率、VF为高频系数垂直方向上的频率和DF为高频系数对角线方向上的频率,通过公式(8)获得
然后在这一窗口下通过公式(9)获得图像间的相关系数R
其中u1为第一图像在此窗口下的像素均值,u2为第二图像在此窗口下的像素均值,
最后根据R值获得融合后高频系数H1L1 F,
当R≥0.8时,表示高度相关,融合后的高频系数为
当0.3<R<0.8时,表示中度相关,
其中T1为第一图像的空间频率在整个频率中所占权重,T2为第二图像的空间频率在整个空间频率中所占权重,即且T2(x,y)+T1(x,y)=1,
当R≤0.3时,表示低度相关,
重复步骤5操作对12个子块进行融合;
步骤6、对以上步骤得到的融合结果进行一致性检测;
步骤7、对所获得的高频和低频系数按照第二顺序进行重构,获得融合后图像;
步骤8、将第一图像与第二图像的融合结果再按照步骤3~7与偏振角图像进行融合。
2.根据权利要求1所述的用于偏振去雾霾的多小波融合处理图片的方法,其特征在于:所述步骤3中的第一顺序与步骤7中的第二顺序为相反的顺序。
3.根据权利要求1所述的用于偏振去雾霾的多小波融合处理图片的方法,其特征在于:所述步骤6中所述的一致性检测:
a、将融合结果中某一像素点作为中心,在其3×3邻域中统计融合后系数来源于第一图像和第二图像的个数;
b、根据统计结果选择个数占优势的那个图像,并将此区域的中心点的像素值,以该占优势的那个图像相应的像素点的值进行定义。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410815711.4A CN104574296B (zh) | 2014-12-24 | 2014-12-24 | 一种用于偏振去雾霾的多小波融合处理图片的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410815711.4A CN104574296B (zh) | 2014-12-24 | 2014-12-24 | 一种用于偏振去雾霾的多小波融合处理图片的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104574296A CN104574296A (zh) | 2015-04-29 |
CN104574296B true CN104574296B (zh) | 2017-07-04 |
Family
ID=53090277
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410815711.4A Expired - Fee Related CN104574296B (zh) | 2014-12-24 | 2014-12-24 | 一种用于偏振去雾霾的多小波融合处理图片的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104574296B (zh) |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104899890A (zh) * | 2015-06-23 | 2015-09-09 | 上海一雅实业有限公司 | 基于偏振光成像技术的检测方法 |
CN107942518B (zh) * | 2018-01-05 | 2020-05-19 | 京东方科技集团股份有限公司 | 增强现实装置、控制方法和计算机可读存储介质 |
CN108492274B (zh) * | 2018-04-03 | 2020-08-07 | 中国人民解放军国防科技大学 | 一种长波红外偏振特征提取与融合的图像增强方法 |
CN109300098B (zh) * | 2018-08-17 | 2022-04-05 | 华东师范大学 | 一种基于小波变换的多聚焦显微图像融合方法 |
CN113159051B (zh) * | 2021-04-27 | 2022-11-25 | 长春理工大学 | 一种基于边缘解耦的遥感图像轻量化语义分割方法 |
CN114841885B (zh) * | 2022-05-10 | 2023-03-10 | 中国矿业大学(北京) | 一种基于偏振图像数据的去雾融合处理方法 |
CN116380256B (zh) * | 2023-06-07 | 2023-08-01 | 长春理工大学 | 基于雾霾衰减系数的短波红外全偏振成像装置与方法 |
CN116503704B (zh) * | 2023-06-27 | 2023-09-05 | 长春理工大学 | 一种强背景下目标偏振探测***及其探测方法 |
CN116659520B (zh) * | 2023-08-02 | 2023-10-20 | 中国人民解放军国防科技大学 | 基于仿生偏振视觉增强的匹配定位方法、装置和设备 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102567977A (zh) * | 2011-12-31 | 2012-07-11 | 南京理工大学 | 一种基于小波的红外偏振图像自适应融合方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8649865B2 (en) * | 2009-12-10 | 2014-02-11 | Pacesetter, Inc. | Method and system for hemodynamic optimization using plethysmography |
-
2014
- 2014-12-24 CN CN201410815711.4A patent/CN104574296B/zh not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102567977A (zh) * | 2011-12-31 | 2012-07-11 | 南京理工大学 | 一种基于小波的红外偏振图像自适应融合方法 |
Non-Patent Citations (4)
Title |
---|
光学相关运动目标识别技术的研究;张肃;《中国博士学位论文全文数据库 信息科技辑》;20140715;I138-17 * |
基于区域特性量测的图像融合方法;王江安 等;《光学与光电技术》;20030415;第1卷(第2期);57-64 * |
基于小波变换的图像融合技术研究;吕艳琼;《中国优秀硕士学位论文全文数据库 信息科技辑》;20091115;I138-1308 * |
基于小波的偏振图像融合算法及性能评价;张晶晶 等;《合肥工业大学学报(自然科学版)》;20090728;第32卷(第7期);1101-1105 * |
Also Published As
Publication number | Publication date |
---|---|
CN104574296A (zh) | 2015-04-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104574296B (zh) | 一种用于偏振去雾霾的多小波融合处理图片的方法 | |
CN106504222B (zh) | 一种基于仿生视觉机理的水下偏振图像融合*** | |
CN105959684B (zh) | 基于双目融合的立体图像质量评价方法 | |
CN101609549B (zh) | 视频模糊图像的多尺度几何分析超分辨处理方法 | |
CN101630405B (zh) | 一种利用核Fisher分类与冗余小波变换的多聚焦图像融合方法 | |
CN105139367A (zh) | 一种基于非下采样剪切波的可见光偏振图像融合方法 | |
CN109242888A (zh) | 一种结合图像显著性和非下采样轮廓波变换的红外与可见光图像融合方法 | |
CN106556940A (zh) | 一种tft‑lcd屏自动光学检测中的背景抑制方法 | |
CN107341786A (zh) | 小波变换与联合稀疏表示的红外与可见光图像融合方法 | |
CN103281559A (zh) | 视频质量检测的方法及*** | |
CN102842120B (zh) | 基于超复数小波相位测量的图像模糊程度检测方法 | |
CN104299196A (zh) | 一种图像处理装置及方法、显示设备 | |
CN105654142A (zh) | 基于自然场景统计的无参考立体图像质量评价方法 | |
CN102800070A (zh) | 基于区域和人眼对比敏感特性的异源图像融合方法 | |
Arivazhagan et al. | A modified statistical approach for image fusion using wavelet transform | |
CN101216936A (zh) | 基于成像机理与非采样Contourlet变换多聚焦图像融合方法 | |
CN103971354A (zh) | 低分辨率红外图像重建高分辨率红外图像的方法 | |
CN103065291A (zh) | 基于提升小波变换和像素区域相关性的图像融合方法 | |
CN107590785A (zh) | 一种基于sobel算子的布里渊散射谱图像识别方法 | |
Junwu et al. | An infrared and visible image fusion algorithm based on LSWT-NSST | |
CN104951800A (zh) | 一种面向资源开采型地区的遥感影像融合方法 | |
CN104899890A (zh) | 基于偏振光成像技术的检测方法 | |
CN107341449A (zh) | 一种基于云块特征变化的静止气象卫星降水估算方法 | |
CN101286236A (zh) | 基于多特征图像和均值漂移的红外目标跟踪方法 | |
CN110223273A (zh) | 一种结合离散余弦变换与神经网络的图像修复取证方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170704 Termination date: 20211224 |
|
CF01 | Termination of patent right due to non-payment of annual fee |