CN106709876B - 一种基于暗像元原理的光学遥感图像去雾方法 - Google Patents
一种基于暗像元原理的光学遥感图像去雾方法 Download PDFInfo
- Publication number
- CN106709876B CN106709876B CN201510533132.5A CN201510533132A CN106709876B CN 106709876 B CN106709876 B CN 106709876B CN 201510533132 A CN201510533132 A CN 201510533132A CN 106709876 B CN106709876 B CN 106709876B
- Authority
- CN
- China
- Prior art keywords
- image
- dark
- defogging
- transmittance
- original image
- 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
Links
Images
Landscapes
- Image Processing (AREA)
Abstract
本发明一种基于暗像元原理的光学遥感图像去雾方法,先读入原始图像,规范化原始图像;自适应缩小原始图像,逐点获取缩小图像的暗像素,得到暗像素图;根据暗像素获取图像的透射率数据,并得到估算的透射率图;利用引导滤波器对图像的透射率进行引导滤波处理;对图像的透射率进行内插放大,得到与原始图像一样大小的透射率图像;对图像进行去雾处理,得到去雾后图像;对去雾后图像进行亮度增强处理,得到增强后图像。本发明对航空及航天遥感图像具备较好的去雾效果;效率高,改善了去雾后图像偏暗的现象。
Description
技术领域
本发明涉及数字图像处理领域,尤其涉及一种基于暗像元原理的光学遥感图像去雾方法。
背景技术
在雾、霭等恶劣气象条件下,大气散射大大降低了光学遥感图像的对比度和色彩保真度,成为制约可见光遥感数据获取与应用的重要因素之一。因此,如何滤除云雾的干扰,清晰地还原地物的原貌,提高数据的利用率,便成为相当重要的课题。
图像去雾处理技术是最近数字图像处理的一个热点,国内外学者取得了一系列研究成果,总体来说,图像去雾方法主要包括两大类:图像增强和图像复原。
由于增强的图像去雾技术并不考虑图像质量下降原因,不能针对图像退化的原因进行补偿,所以只能在一定程度上改善视觉效果,并不能有效去除雾的影响,而且不可避免会产生颜色失真。
近几年来国内外学者基于物理成像模型对雾天成像的退化过程进行了深入分析。在大气介质对光传输影响的研究中,将物理学中的传输理论应用于对雾天图像退化的建模当中,取得较好效果。通过多次散射模型,利用传输方程数值来构建退化模型,也能够取得较好的恢复效果。但以上图像恢复算法所需参数获取困难,因此在实际应用中,无法进行广泛推广使用。
近几年,基于先验信息的图像复原算法快速发展,何凯明的暗像元原理去雾技术最具代表性,其去雾技术在普通户外图像的去雾处理中获得了成功,达到非常好的去雾效果,但直接将其应用于大幅面的遥感图像处理仍然存在暗像元先验失效、计算效率低下和结果影像较暗等突出问题。
发明内容
本发明的目的是为了解决现有技术的不足,提供一种基于暗像元原理的光学遥感图像去雾方法。本发明的目的是通过以下技术方法实现的:
至少包括以下步骤:
(1)读入原始图像,规范化原始图像;
(2)自适应缩小原始图像,逐点获取缩小图像的暗像素,得到暗像素图;
(3)根据暗像素获取图像的透射率数据,并得到估算的透射率图;
(4)利用引导滤波器对图像的透射率进行引导滤波处理,消除块效应;
(5)对图像的透射率进行内插放大,得到与原始图像一样大小的透射率图像;
(6)对图像进行去雾处理,得到去雾后图像;
(7)对去雾后图像进行亮度增强处理,得到增强后图像。
上述的一种基于暗像元原理的光学遥感图像去雾方法,其中,所述步骤(1)中规范化原始图像的方法为:统计原始图像中每个波段的像素值,获取每个波段的最大像素值,将每个波段的像素除以对应的最大像素值,得到规范化的图像。
上述的一种基于暗像元原理的光学遥感图像去雾方法,其中,所述步骤(2)还包括:统计原始图像的均值和方差,根据方差,计算图像缩放系数,方差大,缩小系数小,方差小,缩小系数大。
上述的一种基于暗像元原理的光学遥感图像去雾方法,其中,所述步骤(2)中,暗像素的计算公式为:
上述的一种基于暗像元原理的光学遥感图像去雾方法,其中,所述步骤(6)中,对图像进行去雾处理所采用的公式为:
上述的一种基于暗像元原理的光学遥感图像去雾方法,其中,所述步骤(7)中,对去雾后图像进行亮度增强处理采用的公式为:
综上所述,由于采用了上述技术方案,本发明与现有技术相比有以下优点:
(1)本发明具有更好的适应性,不仅对传统的户外图像有较好的去雾效果,而且对航空及航天遥感图像也具备较好的去雾效果,有效的减少天气状况的干扰和影响;
(2)在效率方面,由于暗像元的提取及优化是在缩小图像上进行的,使效率提高,从而使得大幅面遥感图像的高效处理成为现实,以5倍的缩放系数来算,上述两种功能计算效率提高了近25倍;
(3)通过对去雾后图像亮度的增强,使去雾后图像偏暗的现象也得到较大改善。
附图说明
图1为本发明一种基于暗像元原理的光学遥感图像去雾方法的流程图;
图2为本发明与现有技术的去雾处理比较图。
具体实施方式
下面结合附图对本发明的具体实施方式作详细介绍。
请参见图1,本发明至少包括包括以下步骤:
(2)求取图像的均值及方差,根据均值和方差,计算图像缩放系数,方差大,缩小系数小,方差小,缩小系数大;根据缩小系数自适应缩小原始图像,图像缩小得到图像,逐点对图像计算暗像素: ()=,得到暗像素图;式中,()为缩小图像某点的暗像素;是以为中心的一片方形小区域。
(3)根据暗像素获取图像的透射率数据,并得到估算的透射率图;透射率公式为:
(4)利用引导滤波器对图像的透射率进行引导滤波处理,消除块效应。
(5)对图像的透射率进行内插放大,得到与原始图像一样大小的透射率图像。
(6)对图像进行去雾处理,得到去雾后图像。
(7)对去雾后图像进行亮度增强处理,得到增强后图像。
本发明基于暗像元原理,具体为:光学遥感图像的绝大部分区域中总会有一些像素的值很低。在这些区域中,若无雾气干扰,这些像素的光强最小值会很小并趋于零;在被雾气干扰后,这些暗像素的灰度值会变得较高。因此,这些暗像素能够直接用来评估雾光的透射信息。
对一个影像,定义:
有雾影像的结构方程可描述如下:
实际上,是规范化的有雾图像的暗像素。它直接提供了透射率的估计值。利用抠图技术,应用计算Laplacian矩阵对进行平滑修正处理,,修正后的图能够准确估计图像中场景深度的变换情况,得到透射率的最终估计值。
本发明基于暗像元原理,对去雾结构方程进行了改进,去雾结构方程的改进主要是为了解决暗像元先验失效的问题。暗像元先验有一定的物理基础,但作为一种统计规律,在某些特殊的图像区域(如反射率较高的云、白色建筑物等),暗像元先验失效,利用式(6)会出现严重偏色的问题。
根据式(3),实际的大气透射率计算过程如下式所示:
由于大气透射率低估,即使原始图像中通道间相差几个灰度级,但除以一个很小的数,通道间的差异会被放大很多倍,使得最终图像与原始图像的色彩存在较大的落差。为消除暗像元先验失效引起的色彩失真问题,必须调整失效区域的大气透射率方程,使更加符合。
对此,本发明引入一容差机制来调整复原过程中的大气透射率数值,大气光照A基本上落在暗像元先验失效的区域,记光照为A的透射率为,引入一个容差参数K,对于|I-A|小于K的区域(失效区域),重新计算透射率,适当增加其值;对于对于|I-A|大于K的区域(满足先验区域),保持原始透射率不变。相应的公式表示为如下:
当K=0是,式(8)等价于式(6)。
上述步骤(6)根据式(8)对图像进行去雾处理,得到去雾后图像J。
本发明还针对去雾后图像偏暗的问题进行了改进,上述步骤(7)利用此方法增强去雾后图像的亮度。先利用去雾后各个通道的最大、最小值进行线性拉伸处理;然后,用下列式(13)来调整增个影像的色彩,将去雾后图像的各个通道的灰度均值调整到原始图像相应通道的均值,使得去雾后图像处理与原始图像色彩近似。
如图2显示,图(a)为原始图像,图(b)为直接利用暗像元原理处理的图像,图(c)为采用本发明方法处理的图像。原始图像(a)受雾气干扰严重,直接利用暗像元原理对原始图像(a)进行去雾处理后,结果图像(b)存在严重偏色现象,且雾气去除不明显。从中可看出,图像(c)地物色彩正常,不存在偏色现象,且去雾效率及效果较图像(a)都有了较大提高。
综上所述,本发明与现有技术相比有以下优点:
(1)本发明具有更好的适应性,不仅对传统的户外图像有较好的去雾效果,而且对航空及航天遥感图像也具备较好的去雾效果,有效的减少天气状况的干扰和影响;
(2)在效率方面,由于暗像元的提取及优化是在缩小图像上进行的,使效率提高,从而使得大幅面遥感图像的高效处理成为现实,以5倍的缩放系数来算,上述两种功能计算效率提高了近25倍;
(3)通过对去雾后图像亮度的增强,使去雾后图像偏暗的现象也得到较大改善。
以上所述的实施例仅用于说明本专利的技术思想及特点,其目的在于使本领域内的技术人员能够了解本专利的内容并据以实施,不能仅以本实施例来限定本专利的专利范围,即凡依本专利所揭示的精神所作的同等变化或修饰,仍落在本专利的专利范围内。
Claims (3)
1.一种基于暗像元原理的光学遥感图像去雾方法,其特征在于,至少包括以下步骤:
(1)读入原始图像,规范化原始图像;
(2)统计原始图像的均值和方差,根据方差,计算图像缩放系数,方差大,缩小系数小,方差小,缩小系数大,自适应缩小原始图像,逐点获取缩小图像的暗像素,得到暗像素图,暗像素的计算公式为:
(4)利用引导滤波器对图像的透射率进行引导滤波处理,消除块效应;
(5)对图像的透射率进行内插放大,得到与原始图像一样大小的透射率图像;
(6)对图像进行去雾处理,得到去雾后图像,对图像进行去雾处理所采用的公式为:
(7)对去雾后图像进行亮度增强处理,得到增强后图像。
2.根据权利要求1所述的一种基于暗像元原理的光学遥感图像去雾方法,其特征在于,所述步骤(1)中规范化原始图像的方法为:统计原始图像中每个波段的像素值,获取每个波段的最大像素值,将每个波段的像素除以对应的最大像素值,得到规范化的图像。
3.根据权利要求1所述的一种基于暗像元原理的光学遥感图像去雾方法,其特征在于,所述步骤(7)中,对去雾后图像进行亮度增强处理采用的公式为:
f(x,y)=g(x,y)-mg+mf,
式中,g(x,y)和f(x,y)分别为原始图像和去雾后图像的灰度值;mg为原始图像的灰度均值;mf为去雾后图像的灰度均值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510533132.5A CN106709876B (zh) | 2015-08-27 | 2015-08-27 | 一种基于暗像元原理的光学遥感图像去雾方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510533132.5A CN106709876B (zh) | 2015-08-27 | 2015-08-27 | 一种基于暗像元原理的光学遥感图像去雾方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106709876A CN106709876A (zh) | 2017-05-24 |
CN106709876B true CN106709876B (zh) | 2020-04-14 |
Family
ID=58918631
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510533132.5A Active CN106709876B (zh) | 2015-08-27 | 2015-08-27 | 一种基于暗像元原理的光学遥感图像去雾方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106709876B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108309708A (zh) * | 2018-01-23 | 2018-07-24 | 李思霈 | 盲人拐杖 |
CN111539891A (zh) * | 2020-04-27 | 2020-08-14 | 高小翎 | 单张遥感图像的波段自适应除雾优化处理方法 |
CN112750089B (zh) * | 2020-12-27 | 2023-03-28 | 同济大学 | 基于局部块最大和最小像素先验的光学遥感影像去雾方法 |
CN115829887B (zh) * | 2023-02-24 | 2023-07-07 | 执鼎医疗科技(杭州)有限公司 | 血管图像的处理方法、装置和电子设备 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0645736A2 (en) * | 1993-09-27 | 1995-03-29 | Canon Kabushiki Kaisha | Image processing apparatus |
CN1458791A (zh) * | 2002-04-25 | 2003-11-26 | 微软公司 | 分段分层的图像系统 |
CN103761720A (zh) * | 2013-12-13 | 2014-04-30 | 中国科学院深圳先进技术研究院 | 图像去雾方法以图像去雾装置 |
CN104063848A (zh) * | 2014-06-19 | 2014-09-24 | 中安消技术有限公司 | 一种低照度图像增强方法和装置 |
CN104303208A (zh) * | 2012-05-15 | 2015-01-21 | Sk电信有限公司 | 用于去除包含在视频中的雾的图像处理装置及其方法 |
CN104318535A (zh) * | 2014-11-20 | 2015-01-28 | 广东欧珀移动通信有限公司 | 图像去雾的方法、装置及移动终端 |
-
2015
- 2015-08-27 CN CN201510533132.5A patent/CN106709876B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0645736A2 (en) * | 1993-09-27 | 1995-03-29 | Canon Kabushiki Kaisha | Image processing apparatus |
CN1458791A (zh) * | 2002-04-25 | 2003-11-26 | 微软公司 | 分段分层的图像系统 |
CN104303208A (zh) * | 2012-05-15 | 2015-01-21 | Sk电信有限公司 | 用于去除包含在视频中的雾的图像处理装置及其方法 |
CN103761720A (zh) * | 2013-12-13 | 2014-04-30 | 中国科学院深圳先进技术研究院 | 图像去雾方法以图像去雾装置 |
CN104063848A (zh) * | 2014-06-19 | 2014-09-24 | 中安消技术有限公司 | 一种低照度图像增强方法和装置 |
CN104318535A (zh) * | 2014-11-20 | 2015-01-28 | 广东欧珀移动通信有限公司 | 图像去雾的方法、装置及移动终端 |
Also Published As
Publication number | Publication date |
---|---|
CN106709876A (zh) | 2017-05-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106157267B (zh) | 一种基于暗通道先验的图像去雾透射率优化方法 | |
CN110148095B (zh) | 一种水下图像增强方法及增强装置 | |
CN108596849B (zh) | 一种基于天空区域分割的单幅图像去雾方法 | |
CN102750674B (zh) | 基于自适应容差的视频图像去雾方法 | |
CN108876743B (zh) | 一种图像快速去雾方法、系统、终端及存储介质 | |
Gao et al. | Sand-dust image restoration based on reversing the blue channel prior | |
CN104240194B (zh) | 一种基于抛物线函数的低照度图像增强算法 | |
Fang et al. | Improved single image dehazing using segmentation | |
WO2016206087A1 (zh) | 一种低照度图像处理方法和装置 | |
CN112837233B (zh) | 一种基于差分偏振获取透射率的偏振图像去雾方法 | |
CN107451966B (zh) | 一种采用灰度图引导滤波实现的实时视频去雾方法 | |
TWI808406B (zh) | 圖像去霧方法和使用圖像去霧方法的圖像去霧設備 | |
CN111598791B (zh) | 一种基于改进动态大气散射系数函数的图像去雾方法 | |
CN106709876B (zh) | 一种基于暗像元原理的光学遥感图像去雾方法 | |
WO2020124873A1 (zh) | 图像处理方法 | |
TWI489416B (zh) | 影像還原方法 | |
CN108133462B (zh) | 一种基于梯度场区域分割的单幅图像的复原方法 | |
CN105023246B (zh) | 一种基于对比度和结构相似度的图像增强方法 | |
CN111598800B (zh) | 基于空间域同态滤波和暗通道先验的单幅图像去雾方法 | |
CN115456905A (zh) | 一种基于明暗区域分割的单幅图像去雾方法 | |
CN115393216A (zh) | 基于偏振特性与大气传输模型的图像去雾方法和装置 | |
CN110349113B (zh) | 一种基于暗原色先验改进的自适应图像去雾方法 | |
CN111598814A (zh) | 基于极端散射通道的单图像去雾方法 | |
CN106504216B (zh) | 基于变分模型的单幅图像去雾方法 | |
CN111598811A (zh) | 一种单幅全色遥感雾霾图像清晰化方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |