CN106934771B - 一种基于局部相关性的红外图像条纹噪声去除方法 - Google Patents
一种基于局部相关性的红外图像条纹噪声去除方法 Download PDFInfo
- Publication number
- CN106934771B CN106934771B CN201710084056.3A CN201710084056A CN106934771B CN 106934771 B CN106934771 B CN 106934771B CN 201710084056 A CN201710084056 A CN 201710084056A CN 106934771 B CN106934771 B CN 106934771B
- Authority
- CN
- China
- Prior art keywords
- correction
- expected output
- pixel
- image
- pixel point
- 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
- 238000000034 method Methods 0.000 title claims abstract description 25
- 238000001914 filtration Methods 0.000 claims abstract description 8
- 230000003044 adaptive effect Effects 0.000 claims description 8
- 238000004364 calculation method Methods 0.000 claims description 6
- 238000003331 infrared imaging Methods 0.000 description 5
- 230000000694 effects Effects 0.000 description 2
- 230000003595 spectral effect Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 238000003702 image correction Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 238000002945 steepest descent method Methods 0.000 description 1
- 230000001131 transforming effect Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/70—Denoising; Smoothing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10048—Infrared image
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20024—Filtering details
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Transforming Light Signals Into Electric Signals (AREA)
- Image Processing (AREA)
Abstract
本发明属于红外图像处理领域,具体公开了一种基于局部相关性的红外图像条纹噪声去除方法,其包括:通过一维水平均值滤波计算预校正像素点的期望输出值;计算预校正像素点的预期输出,使用预校正像素点上面一个像素点对应的校正参数对预校正像素点进行校正,得到预期输出;满足期望输出和预期输出的差平方和最小的条件下,通过列间迭代的方式实现校正参数的自适应更新;采用新的校正参数对预校正像素点进行校正;对同一行下一个像素点进行校正,直至一帧图像校正完成。本发明的方法不需要存储一帧图像来计算图像的统计特性,只需要缓存当前像素点的相邻元素,结合条纹噪声的列间局部相关性,就可以实现单帧图像的条纹噪声去除,更适合红外图像数据流水式的特点。
Description
技术领域
本发明属于红外图像处理领域,更具体地,涉及一种基于局部相关性的红外图像条纹噪声去除方法。
背景技术
在红外成像系统中,红外焦平面阵列通常是同一列像元共享同一个输出电路,由于行输出电路偏置电压不完全一致,同时输出电路存在非均匀性,使得红外成像系统产生的图像中会包含以条纹为主要特征的非均匀噪声,称为条纹噪声。传统的基于场景的非均匀校正算法或者基于定标的校正算法,无法有效地去除这种条纹噪声。
目前,解决红外图像条纹噪声的方法主要分为两类:一类是基于统计特征的去除算法,该类算法往往对同一列像素应用相同的校正参数,由于外界随机因素和一些预校正步骤的影响,条纹的统计特性受到了影响,也就是说同一列不同像素对应的校正系数相关,但并不相同(称之为条纹的局部相关性)。因此,该类算法不仅不能有效的去除条纹噪声,而且容易产生“伪影”。另一类是与统计特征无关的去除算法,但这些方法适应性和广泛性较差,往往要求图像中地物类型单一,且该类算法容易造成图象退化和图像“伪影”。还有一类是基于变换的条纹去除算法,比如基于傅里叶变换和小波变化的条纹噪声校正算法,例如专利文献CN104580937A中公开了一种红外成像系统条纹噪声去除方法,其实施步骤如下:1)利用红外成像系统对黑体成像并将得到的红外图像转换为一维数列;2)将得到的一维数列进行FFT傅里叶变换;3)对变换后频谱图中代表条纹噪声的谱线进行标定;4)以标定的谱线位置对同样转换为一维数列的待降噪图像FFT转换后的结果进行处理;5)将处理后的结果逆变换为一维数列,然后拼接回二维图像。该方案可以去除红外图像中的条纹噪声,并不损失图像原信息,但这些算法计算量较大,不利于实时地硬件实现。因此,本领域亟需寻找一种有效、实用的条纹噪声去除技术,既保证红外成像质量,又具有实时性。
发明内容
针对现有技术的以上缺陷或不足,本发明提供一种基于局部相关性的红外图像条纹噪声去除方法,该方案充分利用条纹噪声的局部相关特性,在上方像素校正参数的基础上迭代计算出预校正像素点的实时校正参数,从而实现单帧红外图像的条纹噪声去除。本发明的方法无需计算图像的统计特性,无需储存一整幅图像,只需缓存与预校正像素点相邻的同一行的像素点,对流水式的特点的红外图像数据具有良好的噪声去除效果。
为实现上述目的,按照本发明,提出了一种基于局部相关性的红外图像条纹噪声去除方法,其特征在于,该方法包括:
S1:通过一维水平均值滤波计算图像任一行中的预校正像素点的期望输出;
S2:计算预校正像素点的预期输出,具体为:使用预校正像素点的上一行对应像素点的校正参数对预校正像素点进行校正,得到预期输出;
S3:对校正参数进行自适应更新,具体为:满足期望输出和预期输出的差平方和最小的条件下,通过列间迭代的方式实现校正参数的自适应更新,其中所述实现校正参数自适应更新的公式为:
其中,d(i,j)是像素点(i,j)的期望输出,z(i,j)是预校正像素点的预期输出,g(i,j)和o(i,j)是像素点(i,j)的增益校正系数和偏置校正参数,u是可调节的迭代步长;
S4:采用更新后的校正参数对预校正像素点进行校正;
S5:对同一行下一个像素点进行校正,并依次逐行校正,直至一帧图像校正完成,即可实现该帧图像的条纹噪声去除。
作为本发明的进一步优选,步骤S1中所述利用一维水平均值滤波计算预校正像素点的期望输出具体为:
其中,(i,j)表示任一像素点坐标,x(i,j)是预校正像素点,d(i,j)是像素点(i,j)的期望输出,S是一维水平均值滤波器的窗口大小,N是窗口S中像素点的个数。
作为本发明的进一步优选,步骤S2中利用同一列上方相邻像素的校正参数计算预校正像素的预期输出具体为:
z(i,j)=g(i-1,j)*x(i,j)+o(i-1,j)
其中,g(i-1,j)和o(i-1,j)是像素点(i-1,j)的增益校正系数和偏置校正参数,z(i,j)是预校正像素点的预期输出。
作为本发明的进一步优选,步骤S4中利用更新的校正参数对预校正像素点进行校正具体为:
y(i,j)=g(i,j)*x(i,j)+o(i,j)
其中,x(i,j)是预校正像素点,y(i,j)校正输出灰度值。
总体而言,通过本发明所构思的以上技术方案与现有技术相比,具有以下有益效果:
1)本发明的方法充分考虑并利用条纹噪声的局部相关性,在上方像素校正参数的基础上迭代计算出预校正像素点的实时校正参数,从而实现单帧红外图像的条纹噪声去除;
2)本发明的方法在保证红外图像校正效果的基础上,并未存在过多的计算复杂度,不需要存储一帧图像来计算图像的统计特性,只需缓存与预校正像素点相邻的同一行的像素点(灰度值),更适合红外图像数据流水式的特点;
3)本发明的方法通过一维水平滤波计算期望输出,同时利用列间自适应迭代的方式对参数进行更新,计算量小,实时性高。
附图说明
图1是按照本发明实施例的基于局部相关性的红外图像条纹噪声去除方法流程图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
图1是按照本发明实施例的基于局部相关性的红外图像条纹噪声去除方法流程图。
如图1所示,本发明实施例的基于局部相关性的红外图像条纹噪声去除方法,可以实现对任意帧图像的条纹去除。该方法通过对一帧图像的第一行的预校正像素点开始,逐行进行校正,且与校正像素点相邻的同一列的上方相邻的像素点作为校正参数的依据,实现单帧图像的条纹去除。
具体地,该方法的具体实施过程如下:
步骤S1计算预校正像素点的期望输出。
首先利用一维水平均值滤波计算预校正像素点的期望输出值,具体为:
其中,(i,j)表示像素坐标位置,x(i,j)是预校正像素,d(i,j)是像元(i,j)的期望输出,S是一维水平均值滤波器的窗口大小,N是窗口S中像素点的个数。在一个优选实施例中,窗口大小S可以为1×7,当然本发明中并不限于此。
步骤S2计算预校正像素点的预期输出值。
利用同一列上方相邻像素的校正参数来计算当前预校正像素的预期输出,具体为:
z(i,j)=g(i-1,j)*x(i,j)+o(i-1,j)
其中,g(i-1,j)和o(i-1,j)是像元(i-1,j)的增益校正系数和偏置校正参数,z(i,j)是预校正像素点的预期输出值。
第一行的校正参数即增益校正系数和偏置校正参数通过初始化预设得到,本实施例中,优选地,第一行的增益校正系数和偏置校正参数分别初始化为1和0。
步骤S3校正参数的自适应更新
在使期望输出值和预期输出值的差平方和最小的条件下,通过列间迭代的方式实现校正参数的自适应更新。
本实施例中通过最陡下降法来更新校正参数,以达到期望输出值和预期输出值的差平方和最小的条件。校正参数自适应更新的公式为:
其中g(i,j)和o(i,j)是像元(i,j)的增益校正系数和偏置校正参数,u是可调节的迭代步长。
步骤S4条纹非均匀性校正
利用更新的校正参数对预校正像素点进行校正:
y(i,j)=g(i,j)*x(i,j)+o(i,j)
其中x(i,j)是预校正像素点的输入灰度值,y(i,j)校正输出灰度值。
步骤S5以此类推,对同一行下一个像素点进行校正,并依次逐行校正,直至一帧图像校正完成,即可实现该帧图像的条纹噪声去除。
本发明的方法与常规的现有技术相比,不需要存储一帧图像来计算图像的统计特性,只需要缓存当前像素点的相邻元素,结合条纹噪声的列间局部相关性,就可以实现单帧图像的条纹噪声去除,更适合红外图像数据流水式的特点。另外,本发明方法通过一维水平滤波计算期望输出,利用列间自适应迭代的方式对参数进行更新,计算量小,实时性高。
本领域的技术人员容易理解,以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (2)
1.一种基于局部相关性的红外图像条纹噪声去除方法,其特征在于,该方法包括:
S1:通过一维水平均值滤波计算图像中预校正像素点的期望输出值;
S2:计算预校正像素点的预期输出值,具体为:使用所述预校正像素点上一行对应的像素点的校正参数对预校正像素点进行校正,得到预期输出值;
S3:对校正参数进行自适应更新,具体为:满足期望输出值和预期输出值的差平方和最小的条件下,通过列间迭代的方式实现校正参数的自适应更新,其中所述实现校正参数自适应更新的公式为:
其中,d(i,j)是像素点(i,j)的期望输出,z(i,j)是预校正像素点的预期输出,g(i,j)和o(i,j)是像素点(i,j)的增益校正系数和偏置校正参数,u是可调节的迭代步长;
S4:采用更新后的校正参数对预校正像素点再进行校正;
S5:对同一行下一个像素点进行校正,并依次逐行校正,直至一帧图像校正完成,即可实现该帧图像的条纹噪声去除;
其中,步骤S1中所述利用一维水平均值滤波计算预校正像素点的期望输出值具体为:
其中,(i,j)表示任一像素点坐标,x(i,j)是预校正像素点灰度值,d(i,j)是像素点(i,j)的期望输出值,S是一维水平均值滤波器的窗口大小,N是窗口S中像素点的个数;
步骤S2中利用同一列上方相邻像素的校正参数计算预校正像素的预期输出值具体为:
z(i,j)=g(i-1,j)*x(i,j)+o(i-1,j)
其中,g(i-1,j)和o(i-1,j)是像素点(i-1,j)的增益校正系数和偏置校正参数,z(i,j)是预校正像素点的预期输出值。
2.根据权利要求1所述的一种基于局部相关性的红外图像条纹噪声去除方法,其中,步骤S4中利用更新的校正参数对预校正像素点进行校正具体为:
y(i,j)=g(i,j)*x(i,j)+o(i,j)
其中,x(i,j)是预校正像素点灰度值,y(i,j)校正输出灰度值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710084056.3A CN106934771B (zh) | 2017-02-16 | 2017-02-16 | 一种基于局部相关性的红外图像条纹噪声去除方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710084056.3A CN106934771B (zh) | 2017-02-16 | 2017-02-16 | 一种基于局部相关性的红外图像条纹噪声去除方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106934771A CN106934771A (zh) | 2017-07-07 |
CN106934771B true CN106934771B (zh) | 2020-01-21 |
Family
ID=59424217
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710084056.3A Active CN106934771B (zh) | 2017-02-16 | 2017-02-16 | 一种基于局部相关性的红外图像条纹噪声去除方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106934771B (zh) |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108665422B (zh) * | 2017-08-30 | 2022-05-10 | 西安电子科技大学 | 傅里叶域内逆向感知的单帧红外非均匀性检测方法 |
WO2019183843A1 (zh) * | 2018-03-28 | 2019-10-03 | 西安电子科技大学 | 基于帧间配准和自适应步长的红外图像非均匀性校正方法 |
CN109636755B (zh) * | 2018-12-12 | 2023-03-28 | 电子科技大学 | 一种通过加权估计实现红外热图像增强的方法 |
CN109903235A (zh) * | 2019-01-21 | 2019-06-18 | 天津大学 | 一种红外图像条纹噪声的消除方法 |
CN109767403B (zh) * | 2019-01-23 | 2022-12-16 | 西安科技大学 | 一种基于场景的红外焦平面阵列条状噪声消除方法 |
CN110910324B (zh) * | 2019-11-19 | 2023-04-14 | 山东神戎电子股份有限公司 | 红外视频的去竖条纹方法 |
CN111383196B (zh) * | 2020-03-13 | 2023-07-28 | 浙江大华技术股份有限公司 | 红外图像条纹消除方法、红外探测器及存储装置 |
CN112037301A (zh) * | 2020-08-27 | 2020-12-04 | 上海联影医疗科技股份有限公司 | 重建图像的校正方法、装置和存储介质 |
CN114841900B (zh) * | 2022-07-01 | 2022-10-11 | 南京智谱科技有限公司 | 红外图像横纹去除方法、装置及固定式红外成像设备 |
CN115375545B (zh) * | 2022-08-22 | 2023-08-18 | 杭州微影软件有限公司 | 一种图像校正方法及装置 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102968776A (zh) * | 2012-11-13 | 2013-03-13 | 华中科技大学 | 一种线性与非线性滤波器相结合的非均匀性校正方法 |
CN104296876A (zh) * | 2013-07-15 | 2015-01-21 | 南京理工大学 | 基于fpga的场景非均匀性校正方法及其装置 |
CN104599261A (zh) * | 2014-11-19 | 2015-05-06 | 北京泰和云信科技有限公司 | 短波红外焦平面的均值直方图神经网络非均匀性校正算法 |
CN105092043A (zh) * | 2014-05-20 | 2015-11-25 | 南京理工大学 | 一种基于场景的变积分时间的非均匀性校正方法 |
-
2017
- 2017-02-16 CN CN201710084056.3A patent/CN106934771B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102968776A (zh) * | 2012-11-13 | 2013-03-13 | 华中科技大学 | 一种线性与非线性滤波器相结合的非均匀性校正方法 |
CN104296876A (zh) * | 2013-07-15 | 2015-01-21 | 南京理工大学 | 基于fpga的场景非均匀性校正方法及其装置 |
CN105092043A (zh) * | 2014-05-20 | 2015-11-25 | 南京理工大学 | 一种基于场景的变积分时间的非均匀性校正方法 |
CN104599261A (zh) * | 2014-11-19 | 2015-05-06 | 北京泰和云信科技有限公司 | 短波红外焦平面的均值直方图神经网络非均匀性校正算法 |
Non-Patent Citations (1)
Title |
---|
Minimum mean square error method for stripe nonuniformity correction;Weixian Qian et al.;《CHINESE OPTICS LETTERS》;20110510;正文第1-2页 * |
Also Published As
Publication number | Publication date |
---|---|
CN106934771A (zh) | 2017-07-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106934771B (zh) | 一种基于局部相关性的红外图像条纹噪声去除方法 | |
JP4571670B2 (ja) | 画像モデルにおける色成分の修復のための方法、システム、プログラムモジュールおよびコンピュータプログラム製品 | |
CN108846805B (zh) | 一种基于场景自适应的红外热图像两点非均匀校正方法 | |
US8331714B2 (en) | Methods and systems for image processing | |
US9307208B2 (en) | System and method for solving inverse imaging problems | |
JP2007188493A (ja) | 動きボケ画像中の動きボケを減らす方法、それぞれが各自のボケ・パラメータを持つ複数の動きボケ画像を用いて動きボケの減らされた画像を生成する方法、動きボケ画像中の動きボケを減らす装置、およびそれぞれが各自のボケ・パラメータを持つ複数の動きボケ画像を用いて動きボケの減らされた画像を生成する装置 | |
US9965832B2 (en) | Method for performing super-resolution on single images and apparatus for performing super-resolution on single images | |
CN108230249B (zh) | 基于各向异性的l1范数全变分正则化非均匀性校正方法 | |
US8611691B2 (en) | Automated video data fusion method | |
JP2008146643A (ja) | 動きでぶれた画像における動きのぶれを低減する方法、動きでぶれた画像における動きのぶれを低減するための装置、および動きでぶれた画像における動きのぶれを低減するコンピュータ・プログラムを具現するコンピュータ読み取り可能な媒体 | |
CN109978774B (zh) | 多帧连续等曝光图像的去噪融合方法及装置 | |
KR20100064369A (ko) | 영상 처리 방법 및 장치 | |
JP2017010092A (ja) | 画像処理装置、撮像装置、画像処理方法、画像処理プログラム、および、記憶媒体 | |
CN110033414B (zh) | 一种基于均值化处理的红外图像非均匀性校正方法及系统 | |
CN109813442B (zh) | 一种基于多帧处理的内部杂散辐射非均匀性校正方法 | |
CN107727231B (zh) | 一种对光谱成像仪条带噪声进行修正的方法 | |
CN108961163A (zh) | 一种高分辨率卫星影像超分辨重建方法 | |
Foi | Pointwise shape-adaptive DCT image filtering and signal-dependent noise estimation | |
CN113870110A (zh) | 遥感影像的图像融合方法、装置、电子设备和存储介质 | |
US20140320496A1 (en) | Methods for Reducing Row and Column Patterns in a Digital Image | |
JP6344934B2 (ja) | 画像処理方法、画像処理装置、撮像装置、画像処理プログラムおよび記録媒体 | |
JP6541454B2 (ja) | 画像処理装置、撮像装置、画像処理方法、画像処理プログラム、および、記憶媒体 | |
WO2008086037A2 (en) | Color filter array interpolation | |
CN103997611A (zh) | 一种基于噪声模板的图像固定模式噪声的抑制方法 | |
US11145033B2 (en) | Method and device for image correction |
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 |