CN105869129B - 针对非均匀校正后的热红外图像剩余非均匀噪声去除方法 - Google Patents

针对非均匀校正后的热红外图像剩余非均匀噪声去除方法 Download PDF

Info

Publication number
CN105869129B
CN105869129B CN201610235714.XA CN201610235714A CN105869129B CN 105869129 B CN105869129 B CN 105869129B CN 201610235714 A CN201610235714 A CN 201610235714A CN 105869129 B CN105869129 B CN 105869129B
Authority
CN
China
Prior art keywords
image
value
black matrix
correction
pixel
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
CN201610235714.XA
Other languages
English (en)
Other versions
CN105869129A (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.)
Shanghai Institute of Technical Physics of CAS
Original Assignee
Shanghai Institute of Technical Physics 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 Shanghai Institute of Technical Physics of CAS filed Critical Shanghai Institute of Technical Physics of CAS
Publication of CN105869129A publication Critical patent/CN105869129A/zh
Application granted granted Critical
Publication of CN105869129B publication Critical patent/CN105869129B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/70Denoising; Smoothing
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10048Infrared image

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

本发明公开了一种针对非均匀校正后的热红外图像剩余非均匀噪声去除方法。其实施步骤如下:1)利用红外成像***在不同温度下对黑体成像得到红外图像;2)将得到的黑体图像进行非均匀校正;3)求得校正后黑体图像的像素平均值;4)将校正后的黑体图像中每个像素值与其平均值相减存储结果;5)拍摄外景图像并进行非均匀校正;6)将校正后图像像素值与不同温度下图像平均值做比较找到最接近值对应的相减数组,将像素值与数组中对应位置的值相减得到新的图像。本发明弥补了传统红外图像噪声去除方法不足,提供了一种更有效、更精确去除剩余非均匀空间噪声的方法。本发明能够有效去除红外图像中的剩余非均匀空间噪声,并不损失图像原信息,获得高质量的复原图像。

Description

针对非均匀校正后的热红外图像剩余非均匀噪声去除方法
技术领域:
本发明属于高灵敏度红外成像及应用技术领域,基于一种制冷型红外成像硬件***,涉及红外图像非均匀校正后的剩余空间噪声的去除。
背景技术:
随着红外成像技术的发展,人们对热红外成像***灵敏度水平的要求越来越高,获取更高灵敏度红外图像的需求也更加迫切,经过非均匀性校正的红外图像已不能满足高灵敏度场合的需要,因为即使经过非均匀性校正,红外图像仍然存在剩余非均匀性空间噪声,在高灵敏度场合,这会极大的影响对信号的分辨能力,必须采取新的措施以进一步提高红外图像的灵敏度水平。高灵敏度红外成像***一般采用制冷型红外焦平面器件,由于焦平面阵列中各探测元响应特性不完全一致,不同的探测光敏元件对相同的红外辐射会产生不同的输出信号,使得红外图像出现非均匀性噪声。常规的非均匀性校正算法无法完全去除该噪声,但在一般应用场合该剩余噪声不影响图像质量,所以可见文献中尚无针对非均匀处理后的红外图像进一步去除非均匀性噪声的方法。本方法基于此问题而提出,基于红外焦平面硬件特性,针对性的进一步削弱非均匀性噪声,可以在不损失图像信息的前提下有效抑制剩余非均匀性噪声。
发明内容:
发明目的:本发明提供了一种有效、精确去除剩余非均匀空间噪声的全新方法。这种方法根据像元的响应特性以及在定标图像校正后像元与理想值偏差的固定性出发,创造性的提出了针对剩余非均匀性空间噪声的去除方法:
本发明提出的红外成像***剩余非均匀空间噪声去除方法,包括以下步骤:
1)利用基于M×N规模的焦平面阵列的红外成像***,分别在A,A+e,A+2e,……B的温度下对均匀的黑体进行成像,得到H幅图像Y1、Y2……YH,A为设定的起始温度值,B为设定的终止温度值,e为温度间隔值;
2)选取C1、C2、C3、C4四个温度下的四幅黑体图像F1、F2、F3、F4作为进行非均匀校正的标准图像,用这四幅图像分别对步骤1)中所得到的H幅黑体图像进行非均匀校正,得到H幅校正后的黑体图像X1、X2……XH;非均匀性校正所用系数由如下公式求得:
其中,φk为取定的定标点,Vi,jk)为测得的第i行,第j列个光敏元对应的响应输出;对红外焦平面阵列所有光敏元的输出信号求平均得;
3)对步骤2)中非均匀校正后的每幅图像求像素点平均值,计为aver1,aver2,……averH,将计算得到的平均值结果作为每幅黑体图像的标准值:
其中X(i,j)是图像第i行,第j列元素的像素值;
4)将步骤2)中非均匀校正后的每幅图像中的每个像素点值与步骤3)中求得的对应图像的平均值求差,得到H个M×N的二维数组D1,D2……DH,每个数组中存储着对应图像求得的差值;
5)拍摄外景图像得到红外图像f,用步骤2)中选取的四幅黑体图像F1、F2、F3、F4对图像f进行非均匀校正,得到校正后的图像F;
6)将校正后图像F中的每个像素点值F(i,j)分别与步骤3)中得到的平均值进行比较,找到与像素点的值最接近的平均值averX,找到最接近的平均值后,再找出与该平均值所对应的步骤4)中求出的差值二维矩阵DX,其中X为1~H中的值;
7)将图像F中的像素点F(i,j)与步骤6)中找到的DX数组中与该像素位置对应的点DX(i,j)值求差,将求得的结果作为新的像素值,得到最终去除空间噪声后的红外图像。
本发明相对于现有方法具有如下的优点及效果:
1本发明紧接传统非均匀校正方法之后,创造性的进一步降低了剩余非均匀性噪声,为降低红外图像空间噪声提供了新的手段。
2本发明针对传统非均匀校正方法处理后的红外图像的剩余非均匀性噪声,在实验室以不同像元在不同辐射量于固定定标点校正后对均值的偏差作为补偿,降低了在高灵敏度运用场合剩余非均匀性噪声对图像的干扰,降低了空间噪声。,
3本发明针对非均匀性噪声的固定性得到的偏差补偿数据集,可广泛用于各种场景对固定图案噪声的消除,提高图像的非均匀性,无场景使用限制。
附图说明:
图1是本发明的红外成像***剩余非均匀的空间噪声去除方法框图。
图2是本发明中(10,40)位置上的黑体图像像素在各温度下与其均值求得的差值分布曲线。
图3(a,b)是本发明实施例中非均匀性校正后经过直方图2/255拉伸后的黑体图像(a)与用本文方法降低剩余非均匀空间噪声的同样按照2/255同样比例拉伸后的黑体图像(b)的对比。
图4(a,b)是本发明实施例中待去噪外景图像(a)与经过处理后的外景图像(b)的对比。
具体实施方式:
下面结合实施例及附图对本发明作进一步详细的描述,但本发明的实施方式不限于此。
如图1,以320×256焦平面为硬件的红外成像***为例,改进的红外***剩余非均匀空间噪声去除方法,包括以下步骤:
1)利用红外成像***,分别在6,8,10,12,14,15,17,20,22,24,25,27摄氏度下对均匀黑体进行成像,得到12幅黑体图像Y1、Y2、……Y12,温度的间隔不等。
2)选取10、15、20、25摄氏度四个温度下的四幅黑体图像F1、F2、F3、F4作为进行非均匀校正的标准图像。用这四幅图像分别对1)中所得到的12幅黑体图像进行非均匀校正,得到12幅校正后的黑体图像。
3)对2)中非均匀校正后的每幅图像求像素点平均值,计为aver1=36108,aver2=3737,aver3=38853,aver4=40313…aver7=44174…aver10=49857…aver12=52449。将计算得到的平均值结果作为每幅黑体图像的标准值。
4)将2)中非均匀校正后的每幅图像中的每个像素点值与3)中求得的对应图像的平均值求差。得到12个256×320的二维数组D1,D2……D46,每个数组中存储着对应图像求得的差值,图2所示是各温度下(10,40)位置的像素与各平均值求得的差。
5)拍摄外景图像得到红外图像f,用2)中选取的四幅图像对图像f进行非均匀校正,得到校正后的图像F,如图4(a)所示。
6)将校正后图像F中的每个像素点值F(i,j)分别与3)中得到的平均值进行比较,找到与像素点的值最接近的平均值averX,其中i的值为1~256,j的值为1~320,找到最接近的平均值后,再找出与该平均值所对应的4)中求出的差值二维矩阵DX,其中X为1~12中的值。
7)将图像F中的像素点F(i,j)与6)中找到的DX数组中与该像素位置对应的点DX(i,j)值求差,将求得的结果作为新的像素值,得到最终去除空间噪声后的红外图像。
如图3(a)所示,是没有经过本文方法处理的非均匀性校正拉伸后黑体图像,与原始图像拉伸比例为2/255。其中可以明显的看到存在着剩余非均匀性空间噪声,剩余非均匀性为0.044052%。图3(b)为经过本文方法处理后的图像,可以看到剩余非均匀性噪声降低为0.039722%。
用本文的去除剩余非均匀性空间噪声的方法对外景图像进行降噪,如图4(b)所示为本实施例的去剩余非均匀性空间噪声结果,将图4(b)与图4(a)中的原始空间噪声比较,可见,其中的空间固定噪声被大幅降低,大幅提高了图像质量,凸显了弱信号目标。
以上所述仅为本发明的优选实施方式,本发明的保护范围并不仅限于上述实施方式,凡是属于本发明原理的技术方案均属于本发明的保护范围。对于本领域的技术人员而言,在不脱离本发明的原理的前提下进行的若干改进和润饰,这些改进也应视为本发明的保护范围。

Claims (1)

1.一种针对非均匀校正后的热红外图像剩余非均匀噪声去除方法,包括如下步骤:
1)利用基于M×N规模的焦平面阵列的红外成像***,分别在A,A+e,A+2e,......B的温度下对均匀的黑体进行成像,得到H幅图像Y1、Y2......YH,A为设定的起始温度值,B为设定的终止温度值,e为温度间隔值;
2)在步骤1)A,A+e,A+2e,......B的温度下获得的图像中选取C1、C2、C3、C4四个温度下的四幅黑体图像F1、F2、F3、F4作为进行非均匀校正的标准图像,用这四幅图像分别对步骤1)中所得到的H幅黑体图像进行非均匀校正,得到H幅校正后的黑体图像X1、X2......XH;非均匀性校正所用系数由如下公式求得:
其中,φk为取定的定标点,Vi,jk)为测得的第i行,第j列个光敏元对应的响应输出;为红外焦平面阵列所有光敏元的输出信号的平均值;
3)对步骤2)中非均匀校正后的每幅图像求像素点平均值,计为aver1,aver2,……averH,将计算得到的平均值结果作为每幅黑体图像的标准值:
其中X(i,j)是图像第i行,第j列元素的像素值;
4)将步骤2)中非均匀校正后的每幅图像中的每个像素点值与步骤3)中求得的对应图像的平均值求差,得到H个M×N的二维数组D1,D2……DH,每个数组中存储着对应图像求得的差值;
5)拍摄外景图像得到红外图像f,用步骤2)中选取的四幅黑体图像F1、F2、F3、F4对图像f进行非均匀校正,得到校正后的图像F;
6)将校正后图像F中的每个像素点值F(i,j)分别与步骤3)中得到的平均值进行比较,找到与像素点的值最接近的平均值averX,找到最接近的平均值后,再找出与该平均值所对应的步骤4)中求出的差值二维矩阵DX,其中X为1~H中的值;
7)将图像F中的像素点F(i,j)与步骤6)中找到的DX数组中与该像素位置对应的点DX(i,j)值求差,将求得的结果作为新的像素值,得到最终去除空间噪声后的红外图像。
CN201610235714.XA 2015-12-01 2016-04-15 针对非均匀校正后的热红外图像剩余非均匀噪声去除方法 Active CN105869129B (zh)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN2015108689252 2015-12-01
CN201510868925 2015-12-01

Publications (2)

Publication Number Publication Date
CN105869129A CN105869129A (zh) 2016-08-17
CN105869129B true CN105869129B (zh) 2019-04-02

Family

ID=56632183

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610235714.XA Active CN105869129B (zh) 2015-12-01 2016-04-15 针对非均匀校正后的热红外图像剩余非均匀噪声去除方法

Country Status (1)

Country Link
CN (1) CN105869129B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108254083B (zh) * 2018-01-18 2019-06-14 武汉格物优信科技有限公司 一种红外图像固定图案噪声抑制和评估方法
CN111524057B (zh) * 2020-04-14 2023-06-02 烟台艾睿光电科技有限公司 一种红外图像生成方法、装置、设备及红外热成像***
CN113436088B (zh) * 2021-06-09 2022-07-26 浙江兆晟科技股份有限公司 一种红外图像热窗效应实时抑制方法及***
CN113489907B (zh) * 2021-07-16 2023-07-14 北京富吉瑞光电科技股份有限公司 红外热像仪的成像控制方法、装置及红外成像设备
CN115375590B (zh) * 2022-10-27 2023-04-07 之江实验室 一种亮度非均匀性在线补偿的图像处理方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1430048A (zh) * 2001-12-29 2003-07-16 中国科学院上海技术物理研究所 用于红外焦平面探测器非均匀性校正的方法及装置
CN1834689A (zh) * 2006-04-05 2006-09-20 中国科学院上海技术物理研究所 基于统计约束模型的线阵扫描图像非均匀性校正方法
CN101226639A (zh) * 2008-01-29 2008-07-23 航天东方红卫星有限公司 一种星载tdiccd相机的相对辐射校正方法
CN102589707A (zh) * 2012-02-29 2012-07-18 华中光电技术研究所中国船舶重工集团公司第七一七研究所 红外焦平面阵列探测器非均匀性校正残差的实时补偿方法
CN103335716A (zh) * 2013-06-21 2013-10-02 中国科学院西安光学精密机械研究所 基于变积分时间的面阵红外相机定标与非均匀性校正方法
CN104346784A (zh) * 2014-10-31 2015-02-11 北京津同利华科技有限公司 基于msr的短波红外成像辐射量校正方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1430048A (zh) * 2001-12-29 2003-07-16 中国科学院上海技术物理研究所 用于红外焦平面探测器非均匀性校正的方法及装置
CN1834689A (zh) * 2006-04-05 2006-09-20 中国科学院上海技术物理研究所 基于统计约束模型的线阵扫描图像非均匀性校正方法
CN101226639A (zh) * 2008-01-29 2008-07-23 航天东方红卫星有限公司 一种星载tdiccd相机的相对辐射校正方法
CN102589707A (zh) * 2012-02-29 2012-07-18 华中光电技术研究所中国船舶重工集团公司第七一七研究所 红外焦平面阵列探测器非均匀性校正残差的实时补偿方法
CN103335716A (zh) * 2013-06-21 2013-10-02 中国科学院西安光学精密机械研究所 基于变积分时间的面阵红外相机定标与非均匀性校正方法
CN104346784A (zh) * 2014-10-31 2015-02-11 北京津同利华科技有限公司 基于msr的短波红外成像辐射量校正方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
Research on spatial noise of high resolution thermal infrared imaging system;Hongzhen Ji等;《Proceedings of SPIE》;20141118;第9263卷;第1-10页 *
焦平面红外成像设备非均匀性校正残差分析;樊宏杰 等;《红外技术》;20130731;第35卷(第7期);第409-413页 *
辐射源定标红外焦平面阵列非均匀性校正算法研究;殷世民 等;《光子学报》;20080531;第37卷(第5期);第992-995页 *

Also Published As

Publication number Publication date
CN105869129A (zh) 2016-08-17

Similar Documents

Publication Publication Date Title
CN105869129B (zh) 针对非均匀校正后的热红外图像剩余非均匀噪声去除方法
US11570423B2 (en) System and methods for calibration of an array camera
EP1080443B1 (en) Improved dark frame subtraction
Lin et al. Determining the radiometric response function from a single grayscale image
JP4885179B2 (ja) 映像歪曲補正方法及びその補正方法を採用した映像処理装置
US8306274B2 (en) Methods for estimating peak location on a sampled surface with improved accuracy and applications to image correlation and registration
WO2020224518A1 (zh) 基于局部中值直方图的自适应红外图像去条纹算法
US8463068B2 (en) Methods, systems and apparatuses for pixel value correction using multiple vertical and/or horizontal correction curves
US10931901B2 (en) Method and apparatus for selectively correcting fixed pattern noise based on pixel difference values of infrared images
CN109272520B (zh) 一种联合运动指导与边缘检测的自适应红外焦平面非均匀校正方法
CN106017695B (zh) 基于运动状态估计的自适应红外非均匀性校正方法
CN107093196B (zh) 视频卫星面阵相机在轨相对辐射定标方法
CN108171735B (zh) 基于深度学习的十亿像素视频对齐方法及***
US8089534B2 (en) Multi illuminant shading correction using singular value decomposition
CN113379636A (zh) 一种红外图像非均匀性校正方法、装置、设备及存储介质
CN116934833A (zh) 基于双目视觉水下结构病害检测方法、设备及介质
CN111932478A (zh) 一种非制冷红外焦平面的自适应非均匀校正方法
CN115002360A (zh) 一种基于稳健估计的红外视频非均匀性校正方法
CN109636740B (zh) 一种红外图像多尺度智能化非均匀性校正方法
Zhou et al. Local spatial correlation-based stripe non-uniformity correction algorithm for single infrared images
JP2018160024A (ja) 画像処理装置、画像処理方法及びプログラム
CN112819738A (zh) 红外图像融合方法、装置、计算机设备和存储介质
KR20190113251A (ko) 열화상 카메라의 불균일 상태 보정 방법
Chen et al. An adaptive regression method for infrared blind-pixel compensation
CN111076815A (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