CN114399449B - 一种基于均值滤波分解的形态选通偏振图像融合方法 - Google Patents

一种基于均值滤波分解的形态选通偏振图像融合方法 Download PDF

Info

Publication number
CN114399449B
CN114399449B CN202111388502.2A CN202111388502A CN114399449B CN 114399449 B CN114399449 B CN 114399449B CN 202111388502 A CN202111388502 A CN 202111388502A CN 114399449 B CN114399449 B CN 114399449B
Authority
CN
China
Prior art keywords
image
polarization
source image
detail
polarization angle
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.)
Active
Application number
CN202111388502.2A
Other languages
English (en)
Other versions
CN114399449A (zh
Inventor
刘杰
张耿
杨凡超
蒋鑫
李思远
胡炳樑
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
XiAn Institute of Optics and Precision Mechanics of CAS
Original Assignee
XiAn Institute of Optics and Precision Mechanics of CAS
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by XiAn Institute of Optics and Precision Mechanics of CAS filed Critical XiAn Institute of Optics and Precision Mechanics of CAS
Priority to CN202111388502.2A priority Critical patent/CN114399449B/zh
Publication of CN114399449A publication Critical patent/CN114399449A/zh
Application granted granted Critical
Publication of CN114399449B publication Critical patent/CN114399449B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/50Image enhancement or restoration using two or more images, e.g. averaging or subtraction
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/20Image enhancement or restoration using local operators
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20212Image combination
    • G06T2207/20221Image fusion; Image merging

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Image Processing (AREA)
  • Image Analysis (AREA)

Abstract

本发明涉及偏振图像融合技术,具体涉及一种基于均值滤波分解的形态选通偏振图像融合方法。解决了现有偏振角图像参与度不高,未能有效利用偏振角图像信息以及偏振角图像较强的噪声未得到有效处理,造成融合图像信噪比较低的技术问题。本发明方法的步骤为:S1)利用均值滤波将可见光强度源图像、偏振度源图像、偏振角源图像分解为轮廓图像及细节图像;S2)将偏振度源图像、偏振角源图像进行形态开运算;S3)对各源图像轮廓部分采用加权均值法得到融合轮廓;S4)先对偏振度与偏振角的细节图像信息进行选通,再对选通后的偏振信息与可见光细节图像信息求和得到融合细节;S5)将融合细节部分与融合轮廓部分相加得到最终的融合图像。

Description

一种基于均值滤波分解的形态选通偏振图像融合方法
技术领域
本发明涉及偏振图像融合技术,具体涉及一种基于均值滤波分解的形态选通偏振图像融合方法。
背景技术
物体表面的反射光和散射光会包含其自身特性的偏振信息,通过这些偏振信息,可以反应被测物体的轮廓、纹理、粗糙度等物理信息。偏振度图像及偏振角图像分别根据目标表面偏振度及偏振角进行成像,可实现对目标边缘特征的增强,并对耀光、过曝等影响普通可见光强度图像的因素有较强的抑制作用,有着独特的目标探测能力。
为实现偏振角图像、偏振度图像、可见光强度图像信息的互补以及后续数据利用的便利,需要将上述三种图像进行融合为一。但偏振角图像、偏振度图像中噪声影响较大,尤其偏振角图像,噪声影响较为显著,使得融合过程中想要同时实现滤除噪声和保留源图像信息比较困难。
现有偏振图像融合算法可分为两类:第一类是只利用可见光强度图像和偏振度图像进行融合,噪声影响较大的偏振角图像不参与融合。这类融合方法得到的融合图像信噪比较高,但也同样丢失了偏振角图像中包含的有效信息。第二类是将偏振角图像、偏振度图像、可见光强度图像三者进行融合。这类方法可实现三种源图像信息的融合,但在噪声抑制及信息保留两种目标之间未实现有效的均衡。要么未对偏振角图像中较强的噪声进行有效处理,造成融合图像信噪比较低,要么抑制了噪声,同时保留的有效信息也较少。
发明内容
本发明目的是解决现有偏振角图像、偏振度图像、可见光强度图像三者进行融合时偏振角图像参与度不高,未能有效利用偏振角图像信息以及偏振角图像参与度高的方法中其偏振角图像较强的噪声未得到有效处理,造成融合图像信噪比较低的技术问题,而提供一种基于均值滤波分解的形态选通偏振图像融合方法。本发明利用均值滤波将图像分解为轮廓图像及细节图像,通过形态开运算进行选通,确定区域是否参与细节图像融合,可有效滤除偏振信息图像中的噪声,实现偏振角图像、偏振度图像、可见光强度图像的高质量融合,本发明计算简单,实时性高,有着较高的实用价值。
本发明技术方案为:
一种基于均值滤波分解的形态选通偏振图像融合方法,其特殊之处在于,包括以下步骤:
S1)对可见光强度源图像S、偏振度源图像D、偏振角源图像A进行窗口大小为M×N的均值滤波,得到模糊化的轮廓图像Sf、Df、Af,而后利用源图像灰度值减去对应的轮廓图像灰度值得到可见光强度细节图像Sd、偏振度细节图像Dd、偏振角细节图像Ad
Sd=S-Sf
Dd=D-Df
Ad=A-Af
S2)对偏振度源图像、偏振角源图像进行形态开运算;而后将形态开运算后的偏振度源图像各像素点灰度值、偏振角源图像各像素点灰度值分别与形态开运算前的偏振度源图像像素点灰度均值、偏振角源图像像素点灰度均值进行比较,小于灰度均值置0,大于灰度均值置1,得到偏振度细节图像选通区域模板dopened、偏振角细节图像选通区域模板aopened
S3)对步骤S1)得到的可见光强度源图像、偏振度源图像和偏振角源图像的轮廓图像灰度值采用加权均值法进行加权求和得到融合轮廓Rf
Rf=λ1Sf2Df3Af
式中,λ1,λ2,λ3分别为可见光强度源图像、偏振度源图像和偏振角源图像的轮廓图像在融合图像中的权重;
所述λ1,λ2,λ3均不为0,且λ123=1,λ2=λ3,λ1为0.5~0.6。
S4)先对步骤S1)得到的偏振度细节图像Dd、偏振角细节图像Ad通过步骤S2)得到的选通区域模板进行选通;再对选通后的偏振度细节图像灰度值、偏振角细节图像灰度值与步骤S1)得到的可见光强度细节图像灰度值进行求和得到融合细节Rd
Rd=Sd+dopened×Dd+aopened×Ad
S5)将步骤S3)得到的融合轮廓Rf灰度值与步骤S4)得到的融合细节Rd灰度值相加得到最终的融合图像R:
R=Rf+Rd
进一步地,步骤S1)中,所述窗口大小M×N为30像素*30像素;步骤S2)中,所述形态开运算的窗口为圆形,半径为7个像素;步骤S3)中,所述λ1,λ2,λ3分别为0.5,0.25,0.25。
本发明的有益效果:
1、本发明利用均值滤波将图像分解为轮廓图像及细节图像,通过形态开运算进行选通,确定区域是否参与细节图像融合,可有效滤除偏振信息图像中的噪声,实现偏振角图像、偏振度图像、可见光强度图像的高质量融合。
2、本发明能够较好的保留三种源图像的重要特征信息并对图像噪声有较好的抑制作用,实现了融合图像的高质量表达。
3、本发明计算简单,实时性高,有着较高的实用价值。
附图说明
图1为实施例的三种源图像及采用本发明方法处理后的融合图像(a:偏振度源图像,b:偏振角源图像,c:可见光强度源图像,d:融合图像);
图2为实施例的偏振角选通模板;
图3为实施例的偏振度选通模板。
具体实施方式
下面通过实施例和附图对本发明进行详细的描述。
一种基于均值滤波分解的形态选通偏振图像融合方法,包括以下步骤:
S1)对同一场景的可见光强度源图像S、偏振度源图像D、偏振角源图像A进行窗口大小为M×N的均值滤波,得到模糊化的轮廓图像Sf、Df、Af。本例中窗口为30像素*30像素。而后利用各源图像灰度值减去各自的轮廓图像灰度值得到可见光强度细节图像Sd、偏振度细节图像Dd、偏振角细节图像Ad,实现将源图像分解为轮廓图像及细节图像的目的。
Sd=S-Sf
Dd=D-Df
Ad=A-Af
S2)对偏振度源图像、偏振角源图像进行形态开运算,为较好的区分源图像中目标部分与噪声部分,开运算窗口一般取圆形,半径为7~10个像素。本例中取开运算窗口为圆形,半径为7个像素。而后将形态开运算后的偏振度源图像各像素点灰度值、偏振角源图像各像素点灰度值分别与形态开运算前的偏振度源图像像素点灰度均值、偏振角源图像像素点灰度均值进行比较,小于灰度均值置0,大于灰度均值置1,得到偏振度细节图像选通区域模板dopened、偏振角细节图像选通区域模板aopened
S3)对步骤S1)得到的可见光强度源图像、偏振度源图像和偏振角源图像的轮廓图像灰度值采用加权均值法进行加权求和得到融合轮廓Rf
Rf=λ1Sf2Df3Af
式中,λ1,λ2,λ3分别为可见光强度源图像、偏振度源图像和偏振角源图像的轮廓图像在融合图像中的权重;λ1,λ2,λ3均不为0,且λ123=1,λ2=λ3,λ1为0.5~0.6。可见光强度源图像的轮廓一般较为丰富,故融合图像轮廓以可见光轮廓图像为主要框架,偏振度轮廓图像及偏振角轮廓图像进行补充。本实施例中λ1,λ2,λ3分别取0.5,0.25,0.25。
S4)先对步骤S1)得到的偏振度细节图像Dd、偏振角细节图像Ad通过步骤S2)得到的选通区域模板进行选通;再对选通后的偏振度细节图像灰度值、偏振角细节图像灰度值与步骤S1)得到的可见光强度细节图像灰度值进行求和得到融合细节Rd
Rd=Sd+dopened×Dd+aopened×Ad
S5)将步骤S3)得到的融合轮廓Rf灰度值与步骤S4)得到的融合细节Rd灰度值相加得到最终的融合图像R:
R=Rf+Rd
如图1所示,为室外同一远景的三种源图像及采用本发明方法的融合图像(a:偏振度源图像,b:偏振角源图像,c:可见光强度源图像,d:融合图像)。偏振度源图像a,由于玻璃镜面的强偏振反射,偏振程度显著,造成图像局部亮度高,整体亮度较低,近处树木细节已无法分辨,但对楼房、窗口等目标的边缘反应较好;偏振角源图像b的噪声较大,树木处噪声尤其显著,但对远处电线及楼宇窗口可清晰感知;可见光强度源图像c,远处天空出现过曝情况,电线处已经无法分辨,但近处树木边缘较清晰。上述三种源图像经过本发明融合算法处理后得到融合图像d,融合图像对偏振角源图像中的噪声进行了有效选通抑制,而且树木、楼房、窗口、远处电线等目标的有效信息经过融合均得到有效保留。证明本方法可有效实现对融合图像噪声滤除及信息保留的目标。融合处理过程中只包含均值滤波、形态开运算以及加减乘除等基础运算,计算简单,实时性较好。
如图2所示,根据偏振角选通模板aopened可知,偏振角源图像中信息表达较好的楼房及远处电线处均参与了图像融合,噪声较大的树木部分未参与融合(选通模板亮处参与融合,暗处滤除不参与融合)。
如图3所示,根据偏振度选通模板dopened可知,偏振度源图像中信息表达较好的楼房、窗口等目标边缘均参与了图像融合,源图像中对比度较低部分未参与融合(选通模板亮处参与融合,暗处滤除不参与融合)。
本发明利用均值滤波将源图像分解为轮廓图像及细节图像,通过形态开运算进行选通,确定区域是否参与细节图像融合,可有效滤除偏振信息图像中的噪声,实现偏振角图像、偏振度图像、可见光强度图像的高质量融合,且本发明能够较好的保留三种源图像的重要特征信息实现了融合图像的高质量表达。本发明计算简单,实时性高,有着较高的实用价值。

Claims (2)

1.一种基于均值滤波分解的形态选通偏振图像融合方法,其特征在于,包括以下步骤:
S1)对可见光强度源图像S、偏振度源图像D、偏振角源图像A进行窗口大小为M×N的均值滤波,得到模糊化的轮廓图像Sf、Df、Af,而后利用源图像灰度值减去对应的轮廓图像灰度值得到可见光强度细节图像Sd、偏振度细节图像Dd、偏振角细节图像Ad
Sd=S-Sf
Dd=D-Df
Ad=A-Af
S2)对偏振度源图像、偏振角源图像进行形态开运算;而后将形态开运算后的偏振度源图像各像素点灰度值、偏振角源图像各像素点灰度值分别与形态开运算前的偏振度源图像像素点灰度均值、偏振角源图像像素点灰度均值进行比较,小于灰度均值置0,大于灰度均值置1,得到偏振度细节图像选通区域模板dopened、偏振角细节图像选通区域模板aopened
S3)对步骤S1)得到的可见光强度源图像、偏振度源图像和偏振角源图像的轮廓图像灰度值采用加权均值法进行加权求和得到融合轮廓Rf
Rf=λ1Sf2Df3Af
式中,λ1,λ2,λ3分别为可见光强度源图像、偏振度源图像和偏振角源图像的轮廓图像在融合图像中的权重;
所述λ1,λ2,λ3均不为0,且λ123=1,λ2=λ3,λ1为0.5~0.6;
S4)先对步骤S1)得到的偏振度细节图像Dd、偏振角细节图像Ad通过步骤S2)得到的选通区域模板进行选通;再对选通后的偏振度细节图像灰度值、偏振角细节图像灰度值与步骤S1)得到的可见光细节图像灰度值进行求和得到融合细节Rd
Rd=Sd+dopened×Dd+aopened×Ad
S5)将步骤S3)得到的融合轮廓Rf灰度值与步骤S4)得到的融合细节Rd灰度值相加得到最终的融合图像R:
R=Rf+Rd
2.根据权利要求1所述的一种基于均值滤波分解的形态选通偏振图像融合方法,其特征在于:
步骤S1)中,所述窗口大小M×N为30像素*30像素;
步骤S2)中,所述形态开运算的窗口为圆形,半径为7个像素;
步骤S3)中,所述λ1,λ2,λ3分别为0.5,0.25,0.25。
CN202111388502.2A 2021-11-22 2021-11-22 一种基于均值滤波分解的形态选通偏振图像融合方法 Active CN114399449B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111388502.2A CN114399449B (zh) 2021-11-22 2021-11-22 一种基于均值滤波分解的形态选通偏振图像融合方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111388502.2A CN114399449B (zh) 2021-11-22 2021-11-22 一种基于均值滤波分解的形态选通偏振图像融合方法

Publications (2)

Publication Number Publication Date
CN114399449A CN114399449A (zh) 2022-04-26
CN114399449B true CN114399449B (zh) 2023-04-11

Family

ID=81225682

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111388502.2A Active CN114399449B (zh) 2021-11-22 2021-11-22 一种基于均值滤波分解的形态选通偏振图像融合方法

Country Status (1)

Country Link
CN (1) CN114399449B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115631123B (zh) * 2022-11-22 2023-03-03 北京航空航天大学 仿生视觉融合的恶劣环境成像装置及方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106204510A (zh) * 2016-07-08 2016-12-07 中北大学 一种基于结构相似度约束的红外偏振与光强图像融合方法
WO2017020595A1 (zh) * 2015-08-05 2017-02-09 武汉高德红外股份有限公司 一种可见光图像与红外图像融合处理系统及融合方法
CN110440926A (zh) * 2019-08-13 2019-11-12 北京环境特性研究所 用于动态目标测量的分时型红外偏振成像装置及方法
CN111369484A (zh) * 2020-03-05 2020-07-03 中国铁道科学研究院集团有限公司基础设施检测研究所 钢轨廓形检测方法及装置

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2017020595A1 (zh) * 2015-08-05 2017-02-09 武汉高德红外股份有限公司 一种可见光图像与红外图像融合处理系统及融合方法
CN106204510A (zh) * 2016-07-08 2016-12-07 中北大学 一种基于结构相似度约束的红外偏振与光强图像融合方法
CN110440926A (zh) * 2019-08-13 2019-11-12 北京环境特性研究所 用于动态目标测量的分时型红外偏振成像装置及方法
CN111369484A (zh) * 2020-03-05 2020-07-03 中国铁道科学研究院集团有限公司基础设施检测研究所 钢轨廓形检测方法及装置

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
Active lighting applied to three-dimensional econstruction of specular metallic serfaces by polarization imaging;Morel O 等;《Applied Optics》;20061231;全文 *
Sensitivity study of a polarization imaging-based target detection system;Gu Langlang 等;《IEEE International Conference on Electronic Measurement & Instruments (ICEMI)》;20180122;全文 *
基于偏振信息融合的海洋背景目标检测;梁远安等;《大气与环境光学学报》;20160115(第01期);全文 *
运用小波图像融合技术增强痕迹偏振图像;冯清枝等;《激光与红外》;20160220(第02期);全文 *
雾天偏振成像影响分析及复原方法研究;汪杰君等;《激光技术》;20160725(第04期);全文 *

Also Published As

Publication number Publication date
CN114399449A (zh) 2022-04-26

Similar Documents

Publication Publication Date Title
Zhang et al. Fast haze removal for nighttime image using maximum reflectance prior
Nadimi et al. Physical models for moving shadow and object detection in video
Tripathi et al. Single image fog removal using bilateral filter
CN108876856B (zh) 一种大型建筑火灾火源识别定位方法及系统
US7103227B2 (en) Enhancing low quality images of naturally illuminated scenes
US7359562B2 (en) Enhancing low quality videos of illuminated scenes
US7102638B2 (en) Reducing texture details in images
JP4708343B2 (ja) 背景領域および前景領域をモデリングする方法
US20040183925A1 (en) Stylized imaging using variable controlled illumination
CN111062293B (zh) 一种基于深度学习的无人机森林火焰识别方法
US20040184677A1 (en) Detecting silhouette edges in images
US20040183940A1 (en) Non-photorealistic camera
CN103578084A (zh) 基于亮通道滤波的彩色图像增强方法
CN109523480A (zh) 一种海雾图像的去雾方法、装置、计算机存储介质及终端
Riaz et al. Single image dehazing via reliability guided fusion
CN114399449B (zh) 一种基于均值滤波分解的形态选通偏振图像融合方法
TW201032180A (en) Method and device for keeping image background by multiple gauss models
Mei et al. Single image dehazing using dark channel fusion and haze density weight
Othman et al. Enhanced single image dehazing technique based on HSV color space
Wang et al. Single-image dehazing using color attenuation prior based on haze-lines
CN111368785B (zh) 一种摄像头遮挡判断方法、装置、设备及存储介质
CN113034379A (zh) 天候天时自适应的快速图像清晰化处理方法
CN114399448B (zh) 一种基于非下采样剪切波变换的多偏振信息选通融合方法
Jia et al. A two-step approach to see-through bad weather for surveillance video quality enhancement
CN113052833A (zh) 一种基于红外热辐射的非视域成像方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant