CN103868601B - Irfpa探测器非均匀响应的双边全变分正则化校正方法 - Google Patents
Irfpa探测器非均匀响应的双边全变分正则化校正方法 Download PDFInfo
- Publication number
- CN103868601B CN103868601B CN201410116715.3A CN201410116715A CN103868601B CN 103868601 B CN103868601 B CN 103868601B CN 201410116715 A CN201410116715 A CN 201410116715A CN 103868601 B CN103868601 B CN 103868601B
- Authority
- CN
- China
- Prior art keywords
- value
- parameter
- bilateral
- scene
- correction parameter
- 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.)
- Expired - Fee Related
Links
Landscapes
- Transforming Light Signals Into Electric Signals (AREA)
Abstract
本发明公开了一种IRFPA探测器非均匀响应的双边全变分正则化校正方法,其特征在于,针对IRFPA探测器非均匀响应的校正输出对伪影抑制和细节保持的要求,在能量泛函中引入双边全变分正则项,并与校正参数的回归估计方法相结合,充分利用了双边全变分各向异性的扩散特性,有效地降低了校正后响应输出的伪影效应,更好地保持了原始场景中的细节信息,从而使得IRFPA探测器的温度分辨率得到有效的提升。
Description
技术领域
本发明涉及一种IRFPA探测器非均匀性校正方法,具体涉及一种将双边全变分正则化与校正参数回归估计相结合的IRFPA探测器非均匀性校正方法,属于数字信号处理技术领域。
背景技术
红外焦平面阵列探测器(简称IRFPA探测器)响应不均匀所引起的固定图案噪声已严重影响到IRFPA成像***的信噪比、温度分辨率和空间分辨率。而消除固定模式噪声的最根本途径是在探测器制造过程中采用新的材料和工艺以提升各探测元响应的一致性,然而,新材料和工艺的开发需要大量的经费投入和较长的研究周期,难以满足工程应用对成像***输出信噪比提升的迫切需求。鉴于此,将数字信号处理技术应用于修正IRFPA探测器的非均匀响应,确是一条见效快且效费比高的技术实现途径。
现存的IRFPA探测器非均匀响应校正方法主要有:参考辐射源定标校正法和场景自适应滤波校正法。
参考辐射源定标校正法利用参考辐射源为IRFPA探测器提供均匀辐照,通过对各探测单元的响应进行测量以计算出对应单元的校正参数,进而用于对探测器的实际场景辐照响应进行校正,此类方法由于在实际场景的校正过程中只需将已通过定标获取的校正参数直接用于计算校正输出,因此实时性高,但参考辐射源定标校正法的局限性在于无法自适应地更新校正参数。
场景自适应滤波校正法则利用当前场景辐照响应数据自适应地估计校正参数,并用于后续场景辐照响应的校正,此类方法不断地对校正参数进行迭代估计,不受器件工作状态变化引起参数漂移的影响,因此校正结果的时间稳定性好、精度高。鉴于此,在校正精度要求较高的应用中多采用基于场景的自适应滤波校正方法。
现有的各种基于场景的自适应校正方法中,采用参数回归估计的神经网络校正方法(Adaptive nonuniformity correction for IR FocalPlane Arrays Using Neural Networks,Infrared Sensors:Detectors,Electronics,and Signal Processing,SPIE,1991,Vol.1541,pp.100-109)较之其它方法具有校正效果好、稳定性高以及计算量小的优点。然而,在非均匀响应较为严重的情况下,该方法未能很好地解决伪影现象且对真实场景的结构信息保护欠佳。
综上前述,针对如何在有效保持场景细节信息的前提下消除探测器的非均匀响应,寻求一种可靠的方法对IRFPA非均匀响应进行自适应校正就显得尤为重要。
发明内容
本发明的目的在于将双边全变分正则化加入校正参数的回归估计过程,提出一种新颖的IRFPA探测器非均匀响应的双边全变分正则化校正方法,以减小校正后响应值与实际场景辐照值之间的偏差,进而提高IRFPA探测器的温度分辨率。
为了实现上述目标,本发明采用如下的技术方案:
一种IRFPA探测器非均匀响应的双边全变分正则化校正方法,其特征在于,将双边全变分正则化引入校正参数回归估计中,具体包括如下步骤:
(1)、IRFPA探测器采集得到t时刻的场景辐射观测值Yt;
(2)、利用当前增益校正参数的估计值和偏置校正参数的估计值计算出当前场景辐射观测值Yt的校正值
(3)、利用邻域平均法计算出当前场景的校正值的期望值Tt;
(4)、利用当前场景的校正值和期望值Tt,获得由保真项和双边全变分正则项构成的能量泛函Et;
(5)、采用最陡下降优化方法求解能量泛函Et最小化问题,获得t+1时刻的增益校正参数和偏置校正参数
(6)、根据t+1时刻的增益校正参数和偏置校正参数计算出t+1时刻IRFPA探测器场景辐射观测值Yt+1的校正值
(7)、对IRFPA探测器采集到的后续场景辐射观测值Yn依次执行步骤(3)到步骤(6),便可得到相应的校正值前述n≥t+2。
前述的方法,其特征在于,在步骤(2)和步骤(6)中,基于该时刻增益校正参数和偏置校正参数的场景辐射观测值Y的校正值的计算公式为:
前述的方法,其特征在于,在步骤(3)中,场景校正值的期望值Tt的计算公式为:
式中,符号代表卷积运算。
前述的方法,其特征在于,在步骤(4)中,能量泛函Et表示为:
式中,和分别代表将图像沿x和y方向平移l和m个像素的平移算子。
前述的方法,其特征在于,在步骤(5)中,获取t+1时刻的增益校正参数和偏置校正参数的迭代计算公式分别为:
式中,λ为迭代步长,sign(·)代表符号函数,δ为尺度参数,正则项中的
本发明的有益之处在于:本发明的方法在能量泛函中引入双边全变分正则项,并与校正参数的回归估计方法相结合,充分利用了双边全变分各向异性的扩散特性,有效地降低了校正后响应输出的伪影效应,更好地保持了原始场景中的细节信息,从而使得IRFPA探测器的温度分辨率得到有效的提升。
附图说明
图1是本发明的实现框图;
图2是一幅室外场景非均匀响应图像;
图3是图2中的图像经神经网络方法校正后输出的图像;
图4是图2中的图像经本发明的方法校正后输出的图像。
具体实施方式
以下结合附图和具体实施例对本发明作具体的介绍。
参照图1,本发明的IRFPA探测器非均匀响应的双边全变分正则化校正方法,将双边全变分正则化引入校正参数回归估计中,具体包括如下步骤:
(1)、IRFPA探测器采集得到t时刻的场景辐射观测值Yt。
(2)、利用当前增益校正参数的估计值和偏置校正参数的估计值计算出当前场景辐射观测值Yt的校正值
(3)、利用邻域平均法计算出当前场景校正值的期望值Tt:
式中,符号代表卷积运算。
(4)、利用当前场景的校正值和期望值Tt,获得由保真项和双边全变分正则项构成的能量泛函Et:
式中,为保真项,用于控制估计图像与期望值Tt之间的逼近度;
为双边全变分正则项,用于控制平滑度;
和分别代表将图像沿x和y方向平移l和m个像素的平移算子。
(5)、采用最陡下降优化方法求解能量泛函Et最小化问题,获得t+1时刻的增益校正参数和偏置校正参数计算公式分别为:
式中,sign(·)代表符号函数,λ为迭代步长,δ为尺度参数,在保真项和正则项间起平衡作。而正则项中:
(6)、根据步骤(5)所获得的t+1时刻的增益校正参数和偏置校正参数计算出t+1时刻IRFPA探测器场景辐射观测值Yt+1的校正值
(7)、对IRFPA探测器采集到的后续场景辐射观测值Yn(n≥t+2)依次执行步骤(3)到步骤(6),便可得到相应的校正值
采用现在使用的神经网络方法校正图2中的图像后,输出的图像见图3。
采用本发明的方法校正图2中的图像后,输出的图像见图4。
通过图3和图4的对比可以看出,经本发明方法校正后图像较之神经网络校正方法的校正结果,细节更清晰,锐度更高,残余非均匀性更小,基本观察不到伪影现象。
下面,分别采用粗糙度参数(ρ)、峰值信噪比(PSNR)和结构相似度(SSIM),来量化评估本发明提出的双边全变分正则化非均匀性校正方法与现有神经网络校正方法的性能。
(1)、粗糙度参数
粗糙度参数ρ的数学表达式为:
式中,f表示校正后输出图像,h1表示水平差分模板[1,-1],h2=h1 T表示垂直差分模板,符号代表卷积运算,||·||1则表示L1范数。ρ取值越小则表明图像更平滑,非均匀性越小。
(2)、峰值信噪比
峰值信噪比PSNR的数学表达式为:
式中,I和I0分别表示尺寸为M×N的校正后图像和未受噪声污染图像,PSNR取值越大表明对噪声抑制越好。
(3)、结构相似度
结构相似度SSIM的数学表达式为:
其中,亮度比较项μX和μY分别表示参考图像X和估计图像Y的均值;对比度比较项σX和σY分别表示参考图像X和估计图像Y的方差;结构比较项用参考图像X和估计图像Y的归一化信号(X-μX)/σX和(X-μY)/σY的相关系数来表示,SSIM取值越大说明校正后图像与原始场景结构相似度越高,细节保持能力越好。
对比测试的实验结果见表1。
表1对比测试结果的量化参数对比表
由表1可见:
(1)、经本发明提出的方法校正后的图像的粗糙度参数(ρ)明显低于神经网络校正方法校正结果,说明校正后图像均匀性更高。
(2)、经本发明提出的方法校正后的图像的峰值信噪比(PSNR)明显高于神经网络校正方法校正结果,说明校正后图像残留噪声更小。
(3)、经本发明提出的方法校正后的图像的结构相似性系数(SSIM)明显高于神经网络校正方法校正结果,说明校正后图像保留了原始场景更多的结构信息。
上述结果充分说明了,本发明的方法可有效提升IRFPA探测器响应的均匀性,使校正输出值更加接近实际的响应值。
综上,本发明所提出的双边全变分正则化非均匀性校正方法可有效降低校正后的响应值与实际场景辐照值之间的偏差,对提升IRFPA探测器的温度分辨率起到了重要作用。
需要说明的是,上述实施例不以任何形式限制本发明,凡采用等同替换或等效变换的方式所获得的技术方案,均落在本发明的保护范围内。
Claims (5)
1.IRFPA探测器非均匀响应的双边全变分正则化校正方法,其特征在于,将双边全变分正则化引入校正参数回归估计中,具体包括如下步骤:
(1)、IRFPA探测器采集得到t时刻的场景辐射观测值Yt;
(2)、利用当前增益校正参数的估计值和偏置校正参数的估计值计算出当前场景辐射观测值Yt的校正值
(3)、利用邻域平均法计算出当前场景的校正值的期望值Tt;
(4)、利用当前场景的校正值和期望值Tt,获得由保真项和双边全变分正则项构成的能量泛函Et;
(5)、采用最陡下降优化方法求解能量泛函Et最小化问题,获得t+1时刻的增益校正参数和偏置校正参数
(6)、根据t+1时刻的增益校正参数和偏置校正参数计算出t+1时刻IRFPA探测器场景辐射观测值Yt+1的校正值
(7)、对IRFPA探测器采集到的后续场景辐射观测值Yn依次执行步骤(3)到步骤(6),便可得到相应的校正值所述n≥t+2。
2.根据权利要求1所述的方法,其特征在于,在步骤(2)和步骤(6)中,基于该时刻增益校正参数和偏置校正参数的场景辐射观测值Y的校正值的计算公式为:
3.根据权利要求1所述的方法,其特征在于,在步骤(3)中,场景校正值的期望值Tt的计算公式为:
式中,符号代表卷积运算。
4.根据权利要求1所述的方法,其特征在于,在步骤(4)中,能量泛函Et表示为:
式中,和分别代表将图像沿x和y方向平移l和m个像素的平移算子,δ为尺度参数,α为双边距离权调节参数,0<α<1。
5.根据权利要求1所述的方法,其特征在于,在步骤(5)中,获取t+1时刻的增益校正参数和偏置校正参数的迭代计算公式分别为:
式中,λ为迭代步长,sign(·)代表符号函数,δ为尺度参数,α为双边距离权调节参数,0<α<1,正则项中的
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410116715.3A CN103868601B (zh) | 2014-03-26 | 2014-03-26 | Irfpa探测器非均匀响应的双边全变分正则化校正方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410116715.3A CN103868601B (zh) | 2014-03-26 | 2014-03-26 | Irfpa探测器非均匀响应的双边全变分正则化校正方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103868601A CN103868601A (zh) | 2014-06-18 |
CN103868601B true CN103868601B (zh) | 2016-08-17 |
Family
ID=50907377
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410116715.3A Expired - Fee Related CN103868601B (zh) | 2014-03-26 | 2014-03-26 | Irfpa探测器非均匀响应的双边全变分正则化校正方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103868601B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106780358B (zh) * | 2016-11-15 | 2019-09-17 | 西安电子科技大学 | 基于全变分模型的irfpa非均匀性神经网络校正方法 |
CN107194889B (zh) * | 2017-05-11 | 2021-02-02 | 西安电子科技大学 | 分块双边全变分正则化图像噪声消除方法 |
CN108550122B (zh) * | 2017-12-29 | 2021-10-29 | 西安电子科技大学 | 基于自主路径分块传导滤波的图像去噪方法 |
CN110033414B (zh) * | 2019-03-18 | 2020-12-29 | 华中科技大学 | 一种基于均值化处理的红外图像非均匀性校正方法及*** |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5631466A (en) * | 1995-06-16 | 1997-05-20 | Hughes Electronics | Apparatus and methods of closed loop calibration of infrared focal plane arrays |
CN102538973A (zh) * | 2011-12-31 | 2012-07-04 | 南京理工大学 | 一种快速收敛的基于场景非均匀性校正方法 |
CN102589707A (zh) * | 2012-02-29 | 2012-07-18 | 华中光电技术研究所中国船舶重工集团公司第七一七研究所 | 红外焦平面阵列探测器非均匀性校正残差的实时补偿方法 |
CN103308178A (zh) * | 2013-06-04 | 2013-09-18 | 电子科技大学 | 一种非制冷红外焦平面阵列的非均匀性校正方法 |
CN103491318A (zh) * | 2013-09-25 | 2014-01-01 | 海视英科光电(苏州)有限公司 | 一种红外焦平面探测器图像校正方法及*** |
-
2014
- 2014-03-26 CN CN201410116715.3A patent/CN103868601B/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5631466A (en) * | 1995-06-16 | 1997-05-20 | Hughes Electronics | Apparatus and methods of closed loop calibration of infrared focal plane arrays |
CN102538973A (zh) * | 2011-12-31 | 2012-07-04 | 南京理工大学 | 一种快速收敛的基于场景非均匀性校正方法 |
CN102589707A (zh) * | 2012-02-29 | 2012-07-18 | 华中光电技术研究所中国船舶重工集团公司第七一七研究所 | 红外焦平面阵列探测器非均匀性校正残差的实时补偿方法 |
CN103308178A (zh) * | 2013-06-04 | 2013-09-18 | 电子科技大学 | 一种非制冷红外焦平面阵列的非均匀性校正方法 |
CN103491318A (zh) * | 2013-09-25 | 2014-01-01 | 海视英科光电(苏州)有限公司 | 一种红外焦平面探测器图像校正方法及*** |
Also Published As
Publication number | Publication date |
---|---|
CN103868601A (zh) | 2014-06-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106127688B (zh) | 一种超分辨率图像重建方法及其*** | |
CN103335724B (zh) | 基于定标的场景自适应irfpa非均匀性校正方法 | |
CN107255521B (zh) | 一种红外图像非均匀性校正方法及*** | |
CN109741267B (zh) | 基于三边滤波和神经网络的红外图像非均匀性校正方法 | |
CN104599254B (zh) | 基于组合模糊核结构先验的单透镜计算成像方法 | |
CN103868601B (zh) | Irfpa探测器非均匀响应的双边全变分正则化校正方法 | |
CN102778296B (zh) | 基于总变分的自适应红外焦平面非均匀性校正方法 | |
CN102968765B (zh) | 一种基于sigma滤波器的红外焦平面非均匀性校正方法 | |
CN102968776B (zh) | 一种线性与非线性滤波器相结合的非均匀性校正方法 | |
CN106017695B (zh) | 基于运动状态估计的自适应红外非均匀性校正方法 | |
CN101776486A (zh) | 一种基于红外焦平面非均匀性指纹模式的校正方法 | |
CN108230249B (zh) | 基于各向异性的l1范数全变分正则化非均匀性校正方法 | |
CN107680040A (zh) | 一种基于贝叶斯准则的多帧盲卷积超分辨重建方法及装置 | |
CN102063716A (zh) | 一种基于各向异性约束的多帧迭代盲解卷积图像复原方法 | |
CN101540043B (zh) | 单幅图像复原的解析迭代快速频谱外推方法 | |
CN109272520A (zh) | 一种联合运动指导与边缘检测的自适应红外焦平面非均匀校正方法 | |
CN102521797A (zh) | 扫描型红外成像***的场景非均匀校正方法 | |
CN103164846A (zh) | 一种基于转向核估计的红外条纹非均匀性校正方法 | |
CN108665425A (zh) | 基于帧间配准和自适应步长的红外图像非均匀性校正方法 | |
CN110782403A (zh) | 一种红外图像非均匀性校正方法 | |
CN103076096A (zh) | 基于中值直方图均衡的红外非均匀性校正算法 | |
Rong et al. | An improved non-uniformity correction algorithm and its hardware implementation on FPGA | |
CN104599261A (zh) | 短波红外焦平面的均值直方图神经网络非均匀性校正算法 | |
WO2019183843A1 (zh) | 基于帧间配准和自适应步长的红外图像非均匀性校正方法 | |
Zhao et al. | Fast iterative adaptive nonuniformity correction with gradient minimization for infrared focal plane arrays |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160817 |
|
CF01 | Termination of patent right due to non-payment of annual fee |