CN105184740A - 一种红外焦平面图像非均匀性条带校正方法 - Google Patents
一种红外焦平面图像非均匀性条带校正方法 Download PDFInfo
- Publication number
- CN105184740A CN105184740A CN201510261452.XA CN201510261452A CN105184740A CN 105184740 A CN105184740 A CN 105184740A CN 201510261452 A CN201510261452 A CN 201510261452A CN 105184740 A CN105184740 A CN 105184740A
- Authority
- CN
- China
- Prior art keywords
- dtri
- lambda
- image
- subproblem
- plane 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.)
- Pending
Links
Abstract
本发明提出了一种红外焦平面图像非均匀性条带校正方法,具体为:1)采集含条带噪声的红外焦平面图像的像素点位置;2)对含条带噪声的红外焦平面图像建立单方向总变分去条带能量泛函模型;3)引入辅助变量替换能量泛函模型中的不可微项,进而得到新的能量泛函;4)依据新的能量泛将求解问题拆解为三个子问题,对所述三个子问题进行交替迭代求解,输出得到校正后图像f。本方法利用了条带只会对垂直于其的梯度产生影响,而不影响沿其方向的梯度这一特性知识,能快速有效地去除单方向条带噪声的方法,具有兼顾条带去除和图像细节保存的优点,适应性强,计算复杂度低。
Description
技术领域
本发明属于红外图像处理领域,更具体地,涉及到一种红外焦平面图像非均匀性条带校正方法,适用于快速去除多传感元成像数据中的单方向条带噪声。
背景技术
由于红外传感器材料、工艺水平限制和外界环境的影响,红外焦平面阵列中各探测元间响应不一致,导致图像数据中出现非均匀性条带噪声,严重影响了后续的数据处理。条带噪声校正的困难在于有效地去除各种不同类型的条带噪声的同时较完整地保存图像原有的结构信息。条带校正算法可以大致分为三大类:
第一大类基于统计匹配的方法,以直方图匹配和矩匹配方法为代表。统计匹配的方法简单快速,易于实现。缺点是需要假设图像的行或列的统计特性一致,并且还需事先已知条带的位置,不符合实际情况,去条带效果差。
第二大类基于图像滤波技术,以低通滤波、功率滤波器为代表,优点在于只对条带进行滤波处理,而不处理非条带内容,简单易于实现。但,严重依赖条带检测的准确性,漏检的条带被忽略不被处理,错检则会将图像内容也滤除掉,导致信息丢失。实际情况中很难准确检测条带。
第三大类基于变分正则化方法,将条带校正问题看作估计不含条带图像的逆问题,建立包含数据项和正则化项的能量函数,通过最小化能量函数迭代求解得到不含条带图像。Bouali等人的文章《TowardoptimaldestripingofMODISdatausingaunidirectionalvariationalmodel》提出了一种遥感卫星图像条带校正的变分方法,其能量泛函只包含x方向的梯度域数据保真项TVx(f)和y方向上的数据约束项αTVy(f-g)的总变分能量泛函:
其中TVx和TVy分别表示对图像x与y方向求总变分。第二项表示我们寻找一个和条带图像在竖直方向有几乎一样变化的图像,因为条带噪声几乎只对水平方向的图像梯度产生影响,而第一项表示只对条带图像在水平方向进行平滑,该模型有效的反映了条带的单方向性质,但是该模型结果图像与原始图像在灰度值上整体出现较大偏差。
发明内容
本发明的目的在于提出一种红外焦平面图像非均匀性条带校正方法,将条带校正问题表示为从条带图像中估计真实图像的逆问题,利用了条带只会对垂直于其的梯度产生影响,而不影响沿其方向的梯度这一特性知识,能快速有效地去除单方向条带噪声的方法,具有兼顾条带去除和图像细节保存的优点,适应性强,计算复杂度低。
本发明与现有技术相比,具有如下优点:
第一,兼顾条带去除和图像细节保存。本发明将条带校正问题表示为从条带图像中估计真实图像的逆问题,利用了条带只会对垂直于其的梯度产生影响,而不影响沿其方向的梯度这一特性知识。既能有效地去除条带噪声,又能较好的保存图像细节。
第二,计算复杂度低。本发明利用分裂Bregman方法进行数值优化求解,有效地规避了不可微项。同时采用Guass_Sidel迭代进行求解,具有更快的收敛速度。
第三,适应性强。本发明中提供算法参数调节接口,适应多种条带噪声。
附图说明
图1为本发明方法流程图;
图2为真实红外图像轻度条带去除效果,图2a为原始的噪声图像;图2b为非均匀性校正后的图像;图2c为原始噪声图像与非均匀性校正后图像差值;
图3为真实红外图像中度条带去除结果,图3a为原始的噪声图像;图3b为非均匀性校正后的图像;图3c为原始噪声图像与非均匀性校正后图像差值;
图4为真实红外图像重度条带去除效果,图4a为原始的噪声图像;图4b为非均匀性校正后的图像;图4c为原始噪声图像与非均匀性校正后图像差值。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。
本发明根据单方向条带噪声的特点,在提出包括灰度域数据保真项y方向的梯度域数据保真项βTVy(f-g)和x方向上的数据约束项αTVx(f)的能量泛函,通过分裂Bregman方法求解和离散化实现。
实现本发明的方法是:首先利用单方向条带噪声的性质提出合理的能量泛函模型
g(x,y)为探测器输出的条带图像,f(x,y)为原始无条带图像,E(f)表示关于待估计图像f的能量泛函;能量泛函包括灰度域数据保真项y方向梯度数据保真项βTVy(f-g)和x方向梯度惩罚项αTVx(f);其中TVx和TVy分别表示对图像x与y方向求总变分;α,β是正则化参数,用于调节正则化强度。(x、y)表示图像中像素点的位置,Ω表示图像像素坐标集合。第二行为具体展开表达式。
分裂Bregman方法求解引入辅助变量dx,dy,bx,by,建立新的能量泛函:
dx,dy是替代项,dx替代αTVx(f),dy替代βTVy(f-g),bx,by则是数据找回项,使 更快的接近于0,加快迭代收敛速度。α,β是正则化参数,调节正则化约束强度。λ1,λ2为惩罚参数,用于约束引入的辅助变量。
对单方向总变分去条带能量泛函模型进行迭代求解,输出得到校正后图像f。依据新的能量泛将求解问题拆解为三个子问题,对所述三个子问题进行交替迭代求解,输出得到校正后图像f;其中,
第一子问题为:固定d,b,求解f,得到
第二子问题为:固定f,b,求解d,得到
第三子问题为:固定f,d,求解b,得到
为了方便硬件实现,采用离散化实现。如图1所示,本发明迭代求解的具体流程如下:
1)输入含条带噪声的红外焦平面图像g,图像大小为M×N;初始化迭代次数k=1,图像fk-1=g,辅助变量为M×N大小的零矩阵。
2)将辅助变量代入能量泛函迭代目标函数求解得到校正后图像fk;具体求解过程为:
固定b,d求解f的数值实现
3)将校正后图像fk代入第二子问题求解d和第三子问题求解b对应的迭代目标函数得到辅助变量
固定b,f求解d的数值实现
固定d,f求解b的数值实现
4)如果‖fk+1-fk‖2/‖fk+1‖2<ε或者迭代次数大于要求的迭代次数,ε为预定阈值,则输出fk+1作为最终处理图像,结束;否则,进入步骤5)。
5)更新k=k+1,返回步骤3)。
在本方法中引入了正则化参数α和β,在分裂Bregman优化求解过程中又引入惩罚参数λ1和λ2。本发明针对不同强度的条带示例列出多组与之对应的参数,其仅仅作为示例,不作为对本发明限制。迭代次数固定在40次,在轻条带下的α=20,β=350,λ1=20,λ2=90;在中度条带下的α=80,β=800,λ1=30,λ2=110;在重度条带下的α=200,β=800,λ1=15,λ2=130。
图2是真实红外噪声图像(轻度条带噪声)的条带去除结果。图2a为原始的噪声图像,含有轻度的条带噪声;图2b为非均匀性校正后的图像;图2c为原始噪声图像与非均匀性校正后图像差值(用于描述条带噪声去除情况),α=20,β=350,λ1=20,λ2=90。这组图像清晰地展示了本发明对轻度噪声图像的条带具有良好的去除效果,即保留了图像本身的细节,又最大程度的去除了噪声。
图3是真实红外噪声图像(中度条带噪声)的条带去除结果。图3a为原始的噪声图像,含有中度的条带噪声;图3b为非均匀性校正后的图像;图3c为原始噪声图像与非均匀性校正后图像差值(用于描述条带噪声去除情况),α=80,β=800,lambda1=30,lambda2=110。这组图像清晰地展示了本发明对中度噪声图像的条带具有良好的去除效果,即保留了图像本身的细节(图像中的汉字依然很清晰),又最大程度的去除了噪声。
图4是真实红外噪声图像(重度条带噪声)的条带去除结果。图4a为原始的噪声图像,含有重度的条带噪声;图4b为非均匀性校正后的图像;图4c为原始噪声图像与非均匀性校正后图像差值(用于描述条带噪声去除情况),α=200,β=800,lambda1=15,lambda2=130。这组图像清晰地展示了本发明对重度噪声图像的条带具有良好的去除效果,同时保留了图像本身的细节(树叶和建筑的边缘依然清晰),条带噪声被很好的抑制。
通过图2、图3、图4对比,本专利发明的方法对不同强度的条带噪声红外焦平面图像都有很好的条带去除效果,对细节的保存也是很好的。
Claims (5)
1.一种红外焦平面图像非均匀性条带校正方法,其特征在于,该方法具体为:
1)采集含条带噪声的红外焦平面图像;
2)对含条带噪声的红外焦平面图像建立单方向总变分去条带能量泛函模型:
其中,E(f)表示关于校正后图像f的能量泛函,g为含条带噪声的红外焦平面图像;能量泛函包括灰度域数据保真项TVx(f)和TVy(f)分别表示对校正后图像f的x和y方向的求总变分;α,β是正则化参数,表示求偏导,(x、y)表示图像中像素点的位置,Ω表示图像像素坐标集合;
3)引入辅助变量dx,dy,bx,by替换能量泛函模型中的不可微项,进而得到新的能量泛函:
惩罚参数λ1、λ2为用于约束引入的辅助变量;
4)依据新的能量泛将求解问题拆解为三个子问题,对所述三个子问题进行交替迭代求解,输出得到校正后图像f;其中,
第一子问题为:固定d,b,求解f,得到:
第二子问题为:固定f,b,求解d,得到:
第三子问题为:固定f,d,求解b,得到:
2.根据权利要求1所述的红外焦平面图像非均匀性条带校正方法,其特征在于,所述步骤3)的具体实施方式为:
1)初始化辅助变量为零矩阵和迭代次数k=1;
2)将辅助变量代入第一子问题求解f对应的迭代目标函数,求解得到校正后图像fk;
3)将校正后图像fk代入第二子问题求解d和第三子问题求解b对应的迭代目标函数得到辅助变量
4)若或k达到预定迭代次数最大值,ε为预定阈值,则校正后图像fk即为最终校正后图像f,结束;否则,进入步骤4);
5)更新k=k+1,返回步骤3)。
3.根据权利要求2所述的红外焦平面图像非均匀性条带校正方法,其特征在于,所述第一子问题对应后的迭代目标函数为
fk(i,j)是第k次迭代(i,j)位置处的像素值;代表第k-1次迭代时的辅助变量;分别代表第k-1次迭代时的Bregman变量;g代表原始图像。
4.根据权利要求3所述的红外焦平面图像非均匀性条带校正方法,其特征在于,所述第二子问题对应后的迭代目标函数为:
其中,是从到的线性算子,其中表示N2维的实数维空间,类似地,也是从到的线性算子。shrink是软阈值函数。
5.根据权利要求4所述的红外焦平面图像非均匀性条带校正方法,其特征在于,所述第三子问题对应的迭代目标函数为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510261452.XA CN105184740A (zh) | 2015-05-20 | 2015-05-20 | 一种红外焦平面图像非均匀性条带校正方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510261452.XA CN105184740A (zh) | 2015-05-20 | 2015-05-20 | 一种红外焦平面图像非均匀性条带校正方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN105184740A true CN105184740A (zh) | 2015-12-23 |
Family
ID=54906795
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510261452.XA Pending CN105184740A (zh) | 2015-05-20 | 2015-05-20 | 一种红外焦平面图像非均匀性条带校正方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105184740A (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106462957A (zh) * | 2016-05-19 | 2017-02-22 | 深圳大学 | 一种红外图像中条纹噪声的去除方法及系统 |
CN106952254A (zh) * | 2017-03-17 | 2017-07-14 | 北京思比科微电子技术股份有限公司 | 图像坏行坏列的检测方法 |
CN109767403A (zh) * | 2019-01-23 | 2019-05-17 | 西安科技大学 | 一种基于场景的红外焦平面阵列条状噪声消除方法 |
CN112150377A (zh) * | 2020-08-27 | 2020-12-29 | 北京空间机电研究所 | 一种红外图像非均匀性校正系数解析解交替迭代优化方法 |
CN116402725A (zh) * | 2023-06-08 | 2023-07-07 | 四川工程职业技术学院 | 一种斜条带去除方法、装置、设备及介质 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102183307A (zh) * | 2011-03-17 | 2011-09-14 | 北京航空航天大学 | 红外成像制导中梯度温度气体介质的热辐射噪声获取方法 |
WO2014198625A1 (en) * | 2013-06-13 | 2014-12-18 | Basf Se | Optical detector and method for manufacturing the same |
-
2015
- 2015-05-20 CN CN201510261452.XA patent/CN105184740A/zh active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102183307A (zh) * | 2011-03-17 | 2011-09-14 | 北京航空航天大学 | 红外成像制导中梯度温度气体介质的热辐射噪声获取方法 |
WO2014198625A1 (en) * | 2013-06-13 | 2014-12-18 | Basf Se | Optical detector and method for manufacturing the same |
Non-Patent Citations (2)
Title |
---|
YI CHANG等: "Simultaneous Destriping and Denoising for Remote Sensing Images With Unidirectional Total Variation and Sparse Representation", 《IEEE GEOSCIENCE AND REMOTE SENSING LETTERS》 * |
彭小军等: "基于DSP+FPGA+ASIC 的实时红外图像处理系统", 《微计算机信息》 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106462957A (zh) * | 2016-05-19 | 2017-02-22 | 深圳大学 | 一种红外图像中条纹噪声的去除方法及系统 |
CN106462957B (zh) * | 2016-05-19 | 2019-03-08 | 深圳大学 | 一种红外图像中条纹噪声的去除方法及系统 |
CN106952254A (zh) * | 2017-03-17 | 2017-07-14 | 北京思比科微电子技术股份有限公司 | 图像坏行坏列的检测方法 |
CN109767403A (zh) * | 2019-01-23 | 2019-05-17 | 西安科技大学 | 一种基于场景的红外焦平面阵列条状噪声消除方法 |
CN109767403B (zh) * | 2019-01-23 | 2022-12-16 | 西安科技大学 | 一种基于场景的红外焦平面阵列条状噪声消除方法 |
CN112150377A (zh) * | 2020-08-27 | 2020-12-29 | 北京空间机电研究所 | 一种红外图像非均匀性校正系数解析解交替迭代优化方法 |
CN112150377B (zh) * | 2020-08-27 | 2024-03-29 | 北京空间机电研究所 | 一种红外图像非均匀性校正系数解析解交替迭代优化方法 |
CN116402725A (zh) * | 2023-06-08 | 2023-07-07 | 四川工程职业技术学院 | 一种斜条带去除方法、装置、设备及介质 |
CN116402725B (zh) * | 2023-06-08 | 2023-08-22 | 四川工程职业技术学院 | 一种斜条带去除方法、装置、设备及介质 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Yang et al. | Multi-temporal remote sensing image registration using deep convolutional features | |
CN109102477B (zh) | 一种基于非凸低秩稀疏约束的高光谱遥感图像恢复方法 | |
CN105184740A (zh) | 一种红外焦平面图像非均匀性条带校正方法 | |
CN108230264B (zh) | 一种基于ResNet神经网络的单幅图像去雾方法 | |
CN105825200B (zh) | 基于字典学习和结构稀疏表示的高光谱异常目标检测方法 | |
CN104599292B (zh) | 一种基于低秩矩阵分解的抗噪运动目标检测算法 | |
CN102663693B (zh) | 一种基于最小二乘法的线阵推扫式影像自适应辐射校正方法 | |
CN104112261A (zh) | 基于范数比值正则化的快速图像盲去模糊方法 | |
CN104867150A (zh) | 遥感影像模糊聚类的波段修正变化检测方法及系统 | |
CN109063663B (zh) | 一种由粗到精的时序遥感影像厚云检测与去除方法 | |
CN105654430B (zh) | 一种反差约束的气动热辐射校正方法 | |
CN104240264A (zh) | 一种运动物体的高度检测方法和装置 | |
CN108230249B (zh) | 基于各向异性的l1范数全变分正则化非均匀性校正方法 | |
CN110135344B (zh) | 基于加权固定秩表示的红外弱小目标检测方法 | |
CN103984955A (zh) | 基于显著性特征和迁移增量学习的多摄像机目标识别方法 | |
Peng et al. | Active contours driven by normalized local image fitting energy | |
Huang et al. | Iteratively reweighted unidirectional variational model for stripe non-uniformity correction | |
CN106815820B (zh) | 一种红外图像条状噪声消除方法 | |
CN103824260A (zh) | 一种高分辨率遥感图像薄雾快速去除技术 | |
CN111815620B (zh) | 一种基于卷积特征和低秩表示的织物疵点检测方法 | |
Li et al. | A structure prior weighted hybrid ℓ2–ℓp variational model for single infrared image intensity nonuniformity correction | |
CN113034371A (zh) | 一种基于特征嵌入的红外与可见光图像融合方法 | |
Shit et al. | An encoder‐decoder based CNN architecture using end to end dehaze and detection network for proper image visualization and detection | |
CN104574269A (zh) | 一种月基天文望远镜观测图像行方向拖曳的矫正方法 | |
CN109948571B (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 | ||
WD01 | Invention patent application deemed withdrawn after publication | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20151223 |