CN109084899B - 一种红外焦平面探测器非均匀性本体输出校正方法及装置 - Google Patents

一种红外焦平面探测器非均匀性本体输出校正方法及装置 Download PDF

Info

Publication number
CN109084899B
CN109084899B CN201811102187.0A CN201811102187A CN109084899B CN 109084899 B CN109084899 B CN 109084899B CN 201811102187 A CN201811102187 A CN 201811102187A CN 109084899 B CN109084899 B CN 109084899B
Authority
CN
China
Prior art keywords
temperature
value
slope
working temperature
body output
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
CN201811102187.0A
Other languages
English (en)
Other versions
CN109084899A (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.)
Iray Technology Co Ltd
Original Assignee
Iray Technology Co Ltd
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 Iray Technology Co Ltd filed Critical Iray Technology Co Ltd
Priority to CN201811102187.0A priority Critical patent/CN109084899B/zh
Publication of CN109084899A publication Critical patent/CN109084899A/zh
Application granted granted Critical
Publication of CN109084899B publication Critical patent/CN109084899B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J5/00—Radiation pyrometry, e.g. infrared or optical thermometry
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J5/00—Radiation pyrometry, e.g. infrared or optical thermometry
    • G01J5/80—Calibration

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Radiation Pyrometers (AREA)
  • Photometry And Measurement Of Optical Pulse Characteristics (AREA)

Abstract

本发明实施例公开了一种非制冷红外焦平面探测器非均匀性本体输出校正方法、装置、设备及计算机可读存储介质。其中，方法包括首先获取非制冷红外焦平面探测器快门在三个连续相同间隔温度点触发更新时的三对温度‑本体输出值，根据这6个值分别计算第二温度区间和第三温度区间的直线斜率，并根据计算得到的两个直线斜率预测第一温度区间和第四温度区间的直线斜率，最后根据探测器当前工作环境温度所处的工作温度区间的直线斜率计算相对应的本体输出值。本申请提供的技术方案不依赖场景的运动，不需要进行标定，不限制非制冷红外焦平面探测器的工作环境温度，且有利于提高非制冷红外焦平面探测器的量产效率，普适性强。

Description

一种红外焦平面探测器非均匀性本体输出校正方法及装置
技术领域
本发明实施例涉及非制冷红外焦平面探测器图像处理技术领域，特别是涉及一种非制冷红外焦平面探测器非均匀性本体输出校正方法、装置、设备及计算机可读存储介质。
背景技术
随着红外技术的发展，红外焦平面阵列技术应运而生。该技术制备的探测器(包括制冷型红外焦平面探测器和非制冷红外焦平面探测器)为当今技术性能最先进的红外探测器之一，在军事和民用领域都有着广泛的应用。
制冷型红外焦平面探测器的优势在于灵敏度高，能够分辨更细微的温度差别，探测距离远，但是其结构复杂且成本高昂，主要应用于高端军事装备。非制冷红外焦平面探测器无需制冷装置，能够在室温状态下工作，具有启动快、功耗低、体积小、重量轻、寿命长、成本低等诸多优点。尽管非制冷红外焦平面探测器在灵敏度上与制冷器件有一定差距，但经过近十余年的发展，其在性价比上已经明显优于制冷型探测器，具有更加广阔的应用前景。
但是，由于受到制作器件的半导体材料和加工工艺等条件的限制，非制冷红外焦平面探测器的输出响应并不相同，导致了红外焦平面阵列响应的非均匀性。而响应的非均匀性会直接影响探测器最终输出图像的清晰度，成为进一步提高图像质量的瓶颈，并在一定程度上限制了红外成像系统的应用。
从本质上讲，要彻底解决非制冷红外焦平面阵列响应的非均匀性问题，必须从提高非制冷红外焦平面阵列的加工工艺水平上着手。但是，就目前各学科的发展现状很难保证非制冷红外焦平面阵列每个探测单元响应输出的均匀性。而通过非均匀性校正技术有效地减小或去除非均匀性，成为提高非制冷红外焦平面阵列成像质量的关键所在。
现有的非均匀性校正的方法包括基于标定技术的算法和基于场景技术的算法。
基于标定技术的算法指在实验室内利用均匀的高温和低温黑体对红外焦平面进行标定，从而计算出增益和偏移系数的方法，例如两点法、扩展两点法，多点标定查表法、多点标定拟合法。
两点法和扩展两点法通过标定不同输入能量的两点输出电压，根据标定数据计算出每个探测器单元一次函数的校正参数。使用时，每个探测器单元经过一次函数校正后输出。此类算法虽然结构简单，易于硬件实现，但是没有补偿各探测元差别性漂移的能力，在非制冷焦平面探测器上没有补偿探测器温度变化产生漂移的补偿能力，会导致探测器工作的温度范围较窄，为了降低算法的误差，一般会限制探测器工作环境的温度。
多点标定查表法为将非制冷红外探测器放入高低温控温箱中，每隔一段温度(例如温度小于5度)进行一次标定，记录探测器温度值及标定值。使用时，根据探测器温度值，查到标定值进行补偿计算。多点标定拟合法为将非制冷红外探测器放入高低温控温箱中，每隔一段温度(例如温度小于5度)进行一次标定，记录探测器温度值及标定值，然后对这些标定数据进行拟合，得到拟合曲线。使用时，根据探测器温度值和拟合曲线，算出当前温度探测器的补偿值，再进行补偿计算。由于每颗探测器都需要提前标定大量数据，不仅占用大量的外部存储器Flash和内部存储器DDR，而且标定测试也需要很长的时间及多台高低温控温箱及黑体设备。大容量的外部存储器Flash和内部存储器DDR会提高整个系统成本，而大量数据的标定使得探测器量产效率低，如果增加多套设备，不仅提高生产成本，而且效率也不会有质的提升。
基于场景技术的算法不需要黑体标定，而是利用场景的统计特性，获得每个像素的校正系数，例如时域高通滤波法，恒定统计法以及人工神经网络法等。由于校正系数的获取和成像过程是同时进行的，故基于场景技术的校正算法均有一定的漂移跟踪和补偿能力。但是，该类校正算法的计算量大，收敛速度慢，依赖场景的运动，难于硬件实现。
发明内容
本发明实施例的目的是提供一种非制冷红外焦平面探测器非均匀性本体输出校正方法、装置、设备及计算机可读存储介质，不依赖场景的运动，不需要进行标定，不限制非制冷红外焦平面探测器的工作环境温度，且有利于提高非制冷红外焦平面探测器的量产效率，普适性强。
为解决上述技术问题，本发明实施例提供以下技术方案：
本发明实施例一方面提供了一种非制冷红外焦平面探测器非均匀性本体输出校正方法，包括：
获取非制冷红外焦平面探测器的快门在三个连续相同间隔温度点触发更新时的第一本体输出值和相对应的第一温度值、第二本体输出值和对应的第二温度值、第三本体输出值和对应的第三温度值；
判断所述非制冷红外焦平面探测器的当前工作温度所处的工作温度区间，并根据该工作温度区间对应的斜率计算当前工作温度的本体输出值；
其中，第一工作温度区间小于所述第一温度值，第二工作温度区间大于所述第一温度值且小于所述第二温度值，第三工作温度区间大于所述第二温度值且小于所述第三温度值，第四工作温度区间大于所述第三温度值；
根据所述第一本体输出值、所述第一温度值、所述第二本体输出值、所述第二温度值计算所述第二温度区间的温度-本体输出值直线的第二斜率，根据所述第二本体输出值、所述第二温度值、所述第三本体输出值、所述第三温度值计算所述第三温度区间的第三斜率；根据所述第二斜率、所述第三斜率和预设斜率调整值预测所述第一工作温度区间和所述第四工作温度区间的斜率。
可选的，所述根据所述第二斜率、所述第三斜率和预设斜率调整值预测所述第一工作温度区间和所述第四工作温度区间的斜率包括：
根据下述公式计算第一预测斜率，以作为所述第一工作温度区间的斜率：
Bk(n-1)′＝Bk(n)-(Bk(n+1)-Bk(n))*α；
根据下述公式计算第四预测斜率，以作为所述第四工作温度区间的斜率：
Bk(n+2)′＝Bk(n+1)+(Bk(n+1)-Bk(n))*α；
其中，Bk(n-1)'为所述第一预测斜率，Bk(n)为所述第二斜率，Bk(n+1)为所述第三斜率，α为预设斜率调整值、Bk(n+2)'为所述第四预测斜率。
可选的，所述根据该工作温度区间对应的斜率计算当前工作温度的本体输出值包括：
若当前工作温度位于所述第一工作温度区间，根据下述公式计算当前工作温度对应的本体输出值：
Bc(n-1)＝Bs(n)-Bk(n-1)*(Tc-Tn)；
若当前工作温度位于所述第二工作温度区间，根据下述任意一个公式计算当前工作温度对应的本体输出值：
Tc＜Tn+1-Td，Bc(n)＝Bs(n+1)+Bk(n)*(Tc-Tn+1)；或
Tc＞Tn+Td，Bc(n)＝Bs(n)+Bk(n)*(Tc-Tn)；
若当前工作温度位于所述第三工作温度区间，根据下述任意一个公式计算当前工作温度对应的本体输出值：
Tc＜Tn+2-Td，Bc(n+1)＝Bs(n+2)+Bk(n+1)*(Tc-Tn+2)；或
Tc＞Tn+1+Td，Bc(n+1)＝Bs(n+1)+Bk(n+1)*(Tc-Tn+1)；
若当前工作温度位于所述第四工作温度区间，根据下述公式计算当前工作温度对应的本体输出值：
Bc(n+2)＝Bs(n+2)+Bk(n+2)*(Tc-Tn+2)；
其中，Bs(n)为所述第一本体输出值，Tn为所述第一温度值，Bs(n+1)为所述第二本体输出值，Tn+1为所述第二温度值，Bs(n+2)为所述第三本体输出值，Tn+2为所述第三温度值，Tc为当前工作温度，Bk(n-1)为所述第一工作温度区间的斜率，Bk(n)为所述第二斜率，Bk(n+1)为所述第三斜率，Bk(n+2)为所述第四工作温度区间的斜率，Bc(n-1)为当前工作温度位于所述第一工作温度区间的本体输出值，Bc(n)为当前工作温度位于所述第二工作温度区间的本体输出值，Bc(n+1)为当前工作温度位于所述第三工作温度区间的本体输出值，Bc(n+2)为当前工作温度位于所述第四工作温度区间的本体输出值；Td为触发温度阈值。
可选的，所述获取非制冷红外焦平面探测器的快门在连续三个时间段触发更新时的第一本体输出值和相对应的第一温度值、第二本体输出值和对应的第二温度值、第三本体输出值和对应的第三温度值包括：
所述非制冷红外焦平面探测器处于探测器启动阶段，获取其在快门第一次触发更新时的第一本体输出值和相对应的第一温度值；
所述当前工作温度处于第二工温度作区间，相应的，所述根据该工作温度区间对应的斜率计算当前工作温度的本体输出值包括：
根据下述公式计算所述非制冷红外焦平面探测器的当前工作温度处于所述第二工作温度区间的本体输出值：
B＝Bs(0)+Bk*(Tc-T0)；
式中，Bc为Tc对应的本体输出值，Bs(0)为所述第一本体输出值，T0为所述第一温度值，Bk为预设经验调整值，Tc为当前工作温度。
可选的，还包括：
获取所述非制冷红外焦平面探测器的快门在当前工作温度为第四温度值时，触发更新的第四本体输出值；所述第四温度值为所述第三温度值和触发温度阈值之和；
根据所述第三本体输出值、所述第三温度值、所述第四本体输出值及所述第四温度值计算第四斜率；
利用所述第四斜率更新利用所述第二斜率、所述第三斜率和预设斜率调整值预测得到的斜率，以作为所述第四工作温度区间的斜率；
其中，所述第四工作温度区间大于所述第三温度值且小于所述第四温度值。
可选的，在所述获取所述非制冷红外焦平面探测器的快门在当前工作温度为第四温度值时，触发更新的第四本体输出值之后，还包括：
根据所述第三斜率、所述第四斜率和所述预设斜率调整值预测第五工作温度区间的斜率；所述第五工作温度区间大于所述第四温度值；
相应的，根据下述公式计算第五预测斜率，以作为所述第五工作温度区间的斜率：
Bk(n+3)′＝Bk(n+2)+(Bk(n+2)-Bk(n+1))*α；
其中，Bk(n+3)'为所述第五预测斜率，Bk(n+1)为所述第三斜率，α为预设斜率调整值、Bk(n+2)为所述第四斜率。
本发明实施例另一方面提供了一种非制冷红外焦平面探测器非均匀性本体输出校正装置，包括：
快门更新信息获取模块，用于获取非制冷红外焦平面探测器的快门在三个连续相同间隔温度点触发更新时的第一本体输出值和相对应的第一温度值、第二本体输出值和对应的第二温度值、第三本体输出值和对应的第三温度值；
工作区间判断模块，用于判断所述非制冷红外焦平面探测器的当前工作温度所处的工作温度区间；第一工作温度区间小于所述第一温度值，第二工作温度区间大于所述第一温度值且小于所述第二温度值，第三工作温度区间大于所述第二温度值且小于所述第三温度值，第四工作温度区间大于所述第三温度值；
本体输出值计算模块，用于根据该工作温度区间对应的斜率计算当前工作温度的本体输出值；其中，根据所述第一本体输出值、所述第一温度值、所述第二本体输出值、所述第二温度值计算所述第二温度区间的第二斜率，根据所述第二本体输出值、所述第二温度值、所述第三本体输出值、所述第三温度值计算所述第三温度区间的第三斜率；根据所述第二斜率、所述第三斜率和预设斜率调整值预测所述第一工作温度区间和所述第四工作温度区间的斜率。
可选的，还包括斜率更新模块，所述斜率更新模块包括：
获取子模块，用于获取所述非制冷红外焦平面探测器的快门在当前工作温度为第四温度值时，触发更新的第四本体输出值；所述第四温度值为所述第三温度值和触发温度阈值之和；
计算子模块，用于根据所述第三本体输出值、所述第三温度值、所述第四本体输出值及所述第四温度值计算第四斜率；
更新子模块，用于利用所述第四斜率更新利用所述第二斜率、所述第三斜率和预设斜率调整值预测得到的斜率，以作为所述第四工作温度区间的斜率；其中，所述第四工作温度区间大于所述第三温度值且小于所述第四温度值。
本发明实施例还提供了一种非制冷红外焦平面探测器非均匀性本体输出校正设备，包括处理器，所述处理器用于执行存储器中存储的计算机程序时实现如前任一项所述非制冷红外焦平面探测器非均匀性本体输出校正方法的步骤。
本发明实施例最后还提供了一种计算机可读存储介质，所述计算机可读存储介质上存储有非制冷红外焦平面探测器非均匀性本体输出校正程序，所述非制冷红外焦平面探测器非均匀性本体输出校正程序被处理器执行时实现如前任一项所述非制冷红外焦平面探测器非均匀性本体输出校正方法的步骤。
本发明实施例提供了一种非制冷红外焦平面探测器非均匀性本体输出校正方法，首先获取非制冷红外焦平面探测器快门在三个连续相同间隔温度点触发更新时的三对温度-本体输出值，根据这6个值分别计算第二温度区间和第三温度区间的直线斜率，并根据计算得到的两个直线斜率预测第一温度区间和第四温度区间的直线斜率，最后根据探测器当前工作环境温度所处的工作温度区间的直线斜率计算相对应的本体输出值。
本申请提供的技术方案的优点在于，由于非制冷红外焦平面探测器的本体输出值和探测器工作温度在一个温度区域内近似为单调线性关系，故利用每个工作温度区间的温度-本体输出值直线的斜率计算处于该工作温度区间中任何一个温度点对应的本体输出值，从而实现对探测器像元本体非均匀性的精确校正。该方法在非控温条件下且不需要进行标定，可以节省大量时间和人力，在探测器正常工作的范围内不限制其工作环境温度，解决了传统标定校正方法温度受限的缺陷，有利于提升探测器的量产效率，节约生成成本，适用环境广，增强普适性；此外，不依赖场景运动，简单易于实现，解决了传统场景校正方法的计算量大、收敛速度慢的缺陷。
此外，本发明实施例还针对非制冷红外焦平面探测器非均匀性本体输出校正方法提供了相应的实现装置、设备及计算机可读存储介质，进一步使得所述方法更具有实用性，所述装置、设备及计算机可读存储介质具有相应的优点。
附图说明
为了更清楚的说明本发明实施例或现有技术的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单的介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的非制冷红外焦平面探测器工作温度-像元本体输出值之间的线性关系示意图；
图2为本发明实施例提供的非制冷红外焦平面成像系统的框架示意图；
图3为本发明实施例提供的一种非制冷红外焦平面探测器非均匀性本体输出校正方法的流程示意图；
图4为本发明实施例提供的一种探测器像元本体输出值随温度变化的示意图；
图5为本发明实施例提供的另一种探测器像元本体输出值随温度变化的示意图；
图6为本发明实施例提供的另一种非制冷红外焦平面探测器非均匀性本体输出校正方法的流程示意图；
图7为本发明实施例提供的一种非制冷红外焦平面探测器各工作温度区间的斜率计算过程的流程示意图；
图8为本发明实施例提供的另一种非制冷红外焦平面探测器各工作温度区间的本体输出值计算过程的流程示意图；
图9为本发明实施例提供的非制冷红外焦平面探测器非均匀性本体输出校正装置的一种具体实施方式结构图；
图10为本发明实施例提供的非制冷红外焦平面探测器非均匀性本体输出校正装置的另一种具体实施方式结构图。
具体实施方式
为了使本技术领域的人员更好地理解本发明方案，下面结合附图和具体实施方式对本发明作进一步的详细说明。显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
本申请的说明书和权利要求书及上述附图中的术语“第一”、“第二”、“第三”“第四”等是用于区别不同的对象，而不是用于描述特定的顺序。此外术语“包括”和“具有”以及他们任何变形，意图在于覆盖不排他的包含。例如包含了一系列步骤或单元的过程、方法、系统、产品或设备没有限定于已列出的步骤或单元，而是可包括没有列出的步骤或单元。
本申请的发明人经过研究发现，非制冷红外焦平面探测器的像元本体输出的非均匀性受温度影响较大，为了研究非制冷红外焦平面探测器输出与探测器工作温度之间的关系，将非制冷红外焦平面探测器放到高低温控温箱内，对着固定温度的黑体，调节控温箱的温度，测试得到图1的数据，探测器工作环境的温度和本体输出值都是从红外探测器读出的。根据图1中的数据可分析可以得到，非制冷红外焦平面探测器随着自身温度变化是单调变化的，且自身工作温度变化比较小时，可以近似为一条直线。
由此可知，可以通过已知的直线计算出小范围内的探测器输出值，而且由于其单调性还可以预测出部分直线。
鉴于此，本申请通过利用每个工作温度区间的温度-本体输出值直线的斜率计算处于该工作温度区间中任何一个温度点对应的本体输出值，从而实现对探测器像元本体非均匀性的精确校正。
基于上述本发明实施例的技术方案，下面首先结合图2对本发明实施例的技术方案涉及的一些可能的应用场景进行举例介绍，图2为本发明实施例提供的非制冷红外焦平面成像系统的框架示意图。
如图2所示，非制冷红外焦平面成像系统可包含光学镜头，非制冷红外焦平面探测器，图像处理芯片，存储器，快门及其电机和电机驱动。
其中，光学镜头用于收集目标的光信号并传输给非制冷红外焦平面探测器；快门关闭时用于给非制冷红外焦平面探测器均匀的光信号，不受外部光信号干扰；非制冷焦平面探测器用于将接收到的光信号转换为模拟电压信号，再通过数模转换将其转换为数字信号输出；图像处理芯片用于将接收到的图像数字信号，做图像效果算法处理(即本申请提供非制冷红外焦平面探测器非均匀性本体输出校正方法)及图像数据格式转换处理，得到想要的图像效果和图像数据格式并传输出去；非易失型存储器Flash用来存放响应率SiTF标定数据，易失型存储器DDR用于图像处理芯片算法的数据存放。
图像处理芯片首先获取非制冷红外焦平面探测器快门在三个连续相同间隔温度点触发更新时的三对温度-本体输出值，根据这6个值分别计算第二温度区间和第三温度区间的直线斜率，并根据计算得到的两个直线斜率预测第一温度区间和第四温度区间的直线斜率，最后根据探测器当前工作环境温度所处的工作温度区间的直线斜率计算相对应的本体输出值。
将得到的探测器像元的本体输出B[ij](i、j为像元的坐标)，再根据非制冷红外焦平面探测器非均匀性校正方法的成像公式就可以计算出经过非均匀性校正过后的图像。非制冷红外焦平面探测器非均匀性校正方法的成像公式：
Y[ij]＝K[ij](X[ij]-B[ij])+Bavr；
式中，X[ij]为非制冷红外焦平面探测器原始数据，B[ij]为非制冷红外焦平面探测器像元本体输出值，Bavr为非制冷红外焦平面探测器像元本体输出值的均值，K[ij]为非均匀性响应率校正系数、i、j为像素的坐标值。
需要注意的是，上述应用场景仅是为了便于理解本申请的思想和原理而示出，本申请的实施方式在此方面不受任何限制。相反，本申请的实施方式可以应用于适用的任何场景。
在介绍了本发明实施例的技术方案后，下面详细的说明本申请的各种非限制性实施方式。
首先参见图3，图3为本发明实施例提供的一种非制冷红外焦平面探测器非均匀性本体输出校正方法的流程示意图，本发明实施例可包括以下内容：
S301：获取非制冷红外焦平面探测器的快门在三个连续相同间隔温度点触发更新时的第一本体输出值和相对应的第一温度值、第二本体输出值和对应的第二温度值、第三本体输出值和对应的第三温度值。
S302：根据非制冷红外焦平面探测器的当前工作温度确定其所处的工作温度区间。
S303：根据该工作温度区间对应的斜率计算当前工作温度的本体输出值。
快门在当前工作环境温度达到预设温度时，便会自动触发进行像元本体输出值的更新，也即通过打快门的方式，可以得到当前温度下每个探测器单元的本体输出值，当温度超过快门阈值后，会再次触发打快门得到当前温度下每个探测器单元的本体输出值，三个连续相同间隔温度点即为快门连续三次更新时对应的温度。举例来说，快门第一次本体输出值更新的温度为T1，触发温度阈值为Tth，那么在温度为T1+Tth时，便会触发第二次更新，在T1+2Tth时，便会触发第三次更新，那么，三个连续相同间隔温度点可为T1、T1+Tth、T1+2Tth。
至于如何触发快门得到当前温度下每个探测器单元的本体输出值，可参见现有技术，本申请对此不做任何限定。
工作温度区间可包括四个，即第一工作温度区间、第二工作温度区间、第三工作温度区间及第四工作温度区间。第一工作温度区间小于第一温度值，第二工作温度区间大于第一温度值且小于第二温度值，第三工作温度区间大于第二温度值且小于第三温度值，第四工作温度区间大于第三温度值。第二温度值为第一温度值和快门的触发温度阈值(Tth)之和、第三温度值为第二温度值和快门的触发温度阈值之和、第四温度值为第三温度值和快门的触发温度阈值之和，请参阅图4所示。
举例来说以图4为例，n＝1，三个连续相同间隔温度点为T1、T2＝T1+Tth、T3＝T1+2Tth，那么工作温度区间为：
第一工作温度区间T＜T1，第二工作温度区间为T1＜T＜T2，第三工作温度区间为T2＜T＜T3，第四工作温度区间为T＞T3。
对应每段工作温度区间，均对应描述该段工作温度区间中温度与本体输出值的线性关系，也即可以利用一次函数(即温度-本体输出值直线)进行表示其关系。根据第一本体输出值、第一温度值、第二本体输出值、第二温度值计算第二温度区间的温度-本体输出值直线的第二斜率，根据第二本体输出值、第二温度值、第三本体输出值、第三温度值计算第三温度区间的第三斜率；根据第二斜率、第三斜率和预设斜率调整值预测第一工作温度区间和第四温度区间的斜率。具体的实现过程可如下所述：
在第二工作温度区间，该区间端点值为(T1，Bs(1))、(T2，Bs(2))，根据两点确定一条直线的数学基本知识，可知第二工作温度区间的直线表达式可为(横坐标为温度值、纵坐标为本体输出值)，相应的，该段工作温度区间的斜率也可通过基本数学知识换算得到，根据斜率和任何一个该工作温度区间的温度值，便可得到该温度值对应的本体输出值。
在第三工作温度区间，该区间端点值为(T2，Bs(2))、(T3，Bs(3))，根据两点确定一条直线的数学基本知识，可知第三工作温度区间的直线表达式可为
Figure BDA0001806954460000131
(横坐标为温度值、纵坐标为本体输出值)，相应的，该段工作温度区间的斜率也可通过基本数学知识换算得到，根据斜率和任何一个该工作温度区间的温度值，便可得到该温度值对应的本体输出值。
在第一工作温度区间，该区间的一个端点值为(T1，Bs(1))，由于相邻工作温度区间温度和本体输出值的变化关系相差不是太大，可预先设置斜率调整值(例如可为1)，该斜率调整值为经验值，可衡量相邻两个工作温度区间的斜率变化的预测趋势，故可根据相邻的工作温度区间(第二工作温度区间)的斜率(第二斜率)和斜率调整值可预测得到该工作温度区间的斜率，本领域技术人员可根据当前非制冷红外焦平面探测器的器件自身的结构参数和工作环境进行确定斜率调整值，本申请对此不做任何限定。利用已知点和预测的斜率可确定这个工作温度区间的温度-本体输出值直线，利用该直线可预测得到这个温度区间中任何一个温度值对应的本体输出值。
在第四工作温度区间，该区间的一个端点值为(T3，Bs(3))，由于相邻工作温度区间温度和本体输出值的变化关系相差不是太大，可预先设置斜率调整值(例如可为1)，该斜率调整值为经验值，可衡量相邻两个工作温度区间的斜率变化的预测趋势，故可根据相邻的工作温度区间(第三工作温度区间)的斜率(第三斜率)和斜率调整值可预测得到该工作温度区间的斜率，本领域技术人员可根据当前非制冷红外焦平面探测器的器件自身的结构参数和工作环境进行确定斜率调整值，本申请对此不做任何限定。利用已知点和预测的斜率可确定这个工作温度区间的温度-本体输出值直线，利用该直线可预测得到这个温度区间中任何一个温度值对应的本体输出值。
根据非制冷红外焦平面探测器当前的工作环境温度，确定其所属的工作温度区间，利用相应工作温度区间的斜率计算其相对应的像元本体输出值。
若当前工作温度不在各个工作温度区间之内，但为各工作温度区间的端点值，则不需要计算，直接获取相应的快门更新的本体输出值即可，具体可包括：
若当前工作温度为第一温度值，则其对应的本体输出值即为第一本体输出值；
若当前工作温度为第二温度值，则其对应的本体输出值即为第二本体输出值；
若当前工作温度为第三温度值，则其对应的本体输出值即为第三本体输出值。
需要说明的是，本发明实施例只是针对非制冷红外焦平面探测器的任意一个像元本体而言，非制冷红外焦平面探测器中每个像元本体的输出值校正均可采用本发明实施例提供的技术方案。且该发明实施例只是提供了在已知三对温度-快门更新本体输出值(也即确定两个已知工作环境温度区间)的情况下，对这两个工作温度区间及两个相邻的工作环境温度区间中任何一个温度下对应的本体输出值进行计算的方法，非制冷红外焦平面探测器在工作过程中，其工作环境温度会在正常工作温度范围内进行不断变化，而在当前温度与快门在上个触发更新温度点之差的绝对值大于预设的触发温度阈值时，快门便会自发进行更新本体输出值，这样随着时间的推移，便会得到其他温度点的本体输出值，根据上述实施例提供的技术方案，可实现根据连续相同间隔温度点计算整个正常工作温度范围内任何一个温度点对应的本体输出值。
在本发明实施例提供的技术方案中，由于非制冷红外焦平面探测器的本体输出值和探测器工作温度在一个温度区域内近似为单调线性关系，故利用每个工作温度区间的温度-本体输出值直线的斜率计算处于该工作温度区间中任何一个温度点对应的本体输出值，从而实现对探测器像元本体非均匀性的精确校正。该方法在非控温条件下且不需要进行标定，可以节省大量时间和人力，在探测器正常工作的范围内不限制其工作环境温度，解决了传统标定校正方法温度受限的缺陷，有利于提升探测器的量产效率，节约生成成本，适用环境广，增强普适性；此外，不依赖场景运动，简单易于实现，解决了传统场景校正方法的计算量大、收敛速度慢的缺陷。
可选的，根据第二斜率、第三斜率和预设斜率调整值预测第一工作温度区间和第四温度区间的斜率具体的实施方式可为：
根据下述公式计算第一预测斜率，以作为第一工作温度区间的斜率：
Bk(n-1)′＝Bk(n)-(Bk(n+1)-Bk(n))*α；
根据下述公式计算第四预测斜率，以作为第四工作温度区间的斜率：
Bk(n+2)′＝Bk(n+1)+(Bk(n+1)-Bk(n))*α；
其中，Bk(n-1)'为第一预测斜率，Bk(n)为第二斜率，Bk(n+1)为第三斜率，α为预设斜率调整值、Bk(n+2)'为第四预测斜率。
当然，也可通过其他方式计算预测的斜率，这均不影响本申请的实现。
根据各工作温度区间对应的斜率计算当前工作温度的本体输出值具体可包括：
若当前工作温度位于第一工作温度区间，可根据下述公式计算当前工作温度对应的本体输出值：
Bc(n-1)＝Bs(n)-Bk(n-1)*(Tc-Tn)；
若当前工作温度位于第二工作温度区间，可根据下述任意一个公式计算当前工作温度对应的本体输出值(也即可利用下边两个公式中的任何一个公式进行计算)：
Tc＜Tn+1-Td，Bc(n)＝Bs(n+1)+Bk(n)*(Tc-Tn+1)；或
Tc＞Tn+Td，Bc(n)＝Bs(n)+Bk(n)*(Tc-Tn)；
若当前工作温度位于第三工作温度区间，根据下述任意一个公式计算当前工作温度对应的本体输出值(也即可利用下边两个公式中的任何一个公式进行计算)：
Tc＜Tn+2-Td，Bc(n+1)＝Bs(n+2)+Bk(n+1)*(Tc-Tn+2)；或
Tc＞Tn+1+Td，Bc(n+1)＝Bs(n+1)+Bk(n+1)*(Tc-Tn+1)；
若当前工作温度位于第四工作温度区间，根据下述公式计算当前工作温度对应的本体输出值：
Bc(n+2)＝Bs(n+2)+Bk(n+2)*(Tc-Tn+2)；
其中，Bs(n)为第一本体输出值，Tn为第一温度值，Bs(n+1)为第二本体输出值，Tn+1为第二温度值，Bs(n+2)为第三本体输出值，Tn+2为第三温度值，Tc为当前工作温度，Bk(n-1)为第一工作温度区间的斜率，Bk(n)为第二斜率，Bk(n+1)为第三斜率，Bk(n+2)为第四工作温度区间的斜率，Bc(n-1)为当前工作温度位于第一工作温度区间的本体输出值，Bc(n)为当前工作温度位于第二工作温度区间的本体输出值，Bc(n+1)为当前工作温度位于第三工作温度区间的本体输出值，Bc(n+2)为当前工作温度位于第四工作温度区间的本体输出值；Td为触发温度阈值。
在计算各工作温度区间内工作温度的本体输出值时，如果该段工作温度区间的斜率为已知端点计算得到的，便可用该实际斜率进行计算；如若当前工作温度区间的斜率为根据已知斜率预测得到的斜率，则利用预测的斜率进行计算，后期可在快门得到更新的本体输出值后，对该预测斜率进行更新，利用更新后的斜率重新计算，从而可对先前得到的本体输出值进行更新。当前，为了不增加整个图像处理过程中的计算量，也可不进行更新本体输出值。
在一种具体的实施方式中，当还没有完全得到三个连续相同间隔温度点触发更新时的本体输出值时，也即在非制冷红外焦平面探测器启动阶段时，计算启动时间的本体输出值，可根据下述方法进行计算：
非制冷红外焦平面探测器处于探测器启动阶段，当快门在第一次触发更新时，得到第一本体输出值和相对应的第一温度值；也即当前工作温度处于第二工作温度区间，相应的，根据该工作温度区间对应的斜率计算当前工作温度的本体输出值包括：
根据下述公式计算非制冷红外焦平面探测器的当前工作温度处于第二工作温度区间(T0-T1)的本体输出值：
B＝Bs(0)+Bk*(Tc-T0)；
式中，Bc为Tc对应的本体输出值，Bs(0)为第一本体输出值，T0为第一温度值，Bk为预设经验调整值，Tc为当前工作温度。
本领域技术人员可根据当前非制冷红外焦平面探测器的器件自身的结构参数和工作环境确定Bk(例如可为一个量化后的值，-12872)，本申请对此不做任何限定。
开机后温度慢慢上升，请参阅图5所示，温度变化值从(Tc-T0)大于触发快门更新的触发温度阈值Tth，触发快门B值更新得到第二组(Bs(1)，T1)，可以算出第二工作温度区间的第二斜率Bk2＝(Bs(1)-Bs(0))/(T1–T0)。
当前工作温度在工作环境温度区间T0～T2时，当前本体输出值计算如下：
B＝Bs(1)+Bk2(Tc–T1)。
当前工作温度在工作环境温度区间T-1～T0时，当前本体输出值计算如下：
B＝Bs(0)+Bk2(Tc–T0)。
通过对非制冷红外焦平面探测器在启动阶段的本体输出值的校正，进一步的提升了整个非制冷红外焦平面探测器成像的准确性。
基于上述实施例，本申请还提供了另外一个实施例，请参阅图6，具体可包括：
S304：获取非制冷红外焦平面探测器的快门在当前工作温度为第四温度值时，触发更新的第四本体输出值。
第四温度值为第三温度值和触发温度阈值之和，也即在第三温度值和第四温度值的差值大于触发温度阈值时，快门便会自动触发进行像元本体输出值的更新，得到与S301的第三温度值相邻的下一个温度点的本体输出值。
S305：根据第三本体输出值、第三温度值、第四本体输出值及第四温度值计算第四斜率。
对第四工作温度区间而言，此时，便得到两个端点(T3，Bs(3))、(T4，Bs(4))，根据两点确定一条直线的数学基本知识，可知第四工作温度区间的直线表达式可为
Figure BDA0001806954460000181
(横坐标为温度值、纵坐标为本体输出值)，相应的，该段工作温度区间的斜率也可通过基本数学知识换算得到，根据斜率和任何一个该工作温度区间的温度值，便可得到该温度值对应的本体输出值。
S306：利用第四斜率更新利用第二斜率、第三斜率和预设斜率调整值预测得到的斜率，以作为第四工作温度区间的斜率。
第四工作温度区间大于第三温度值且小于第四温度值。在得到第四工作温度区间的斜率后，便可利用该值更新预测得到的斜率值，并可相应的更新下个工作温度区间的斜率和本工作温度区间的本体输出值。
S307：根据第三斜率、第四斜率和预设斜率调整值预测第五工作温度区间的斜率。
第五工作温度区间大于第四温度值；
相应的，可根据下述公式计算第五预测斜率，以作为第五工作温度区间的斜率：
Bk(n+3)′＝Bk(n+2)+(Bk(n+2)-Bk(n+1))*α；
其中，Bk(n+3)'为第五预测斜率，Bk(n+1)为第三斜率，α为预设斜率调整值、Bk(n+2)为第四斜率。
由上可知，可根据S301-S306对整个非制冷红外焦平面探测器在其正常工作温度范围区间(分为多个工作温度区间)内进行斜率预测-斜率更新，然后根据斜率计算本体输出值，从而实现对工作温度区间中的任意一个温度的本体输出值的校正，从而提升整个非制冷红外焦平面探测器成像质量。
为了本领域技术人员更加清楚明白本申请技术方案的原理和实现过程，本申请还提供了具体的实例，请参阅图7及图8所示，具体可包括：
首先，图像处理芯片出图开始时，算法进入启动阶段，开启出图时，就会做一次快门Bs值更新，得到第一组(Bs(0)，T0)，当前温度Tc对应的本体输出值计算如下：
B＝Bs(0)+Bk*(Tc-T0)。
开机后温度上升，温度变化值(Tc-T0)大于触发快门更新的触发温度阈值Tth，触发快门B值更新得到第二组(Bs(1)，T1)，可以算出该工作温度区间的斜率为：
Bk2＝(Bs(1)-Bs(0))/(T1–T0)
当前工作温度在工作环境温度区间T0～T2时，当前本体输出值计算如下：
B＝Bs(1)+Bk2(Tc–T1)。
当前工作温度在工作环境温度区间T-1～T0时，当前本体输出值计算如下：
B＝Bs(0)+Bk2(Tc–T0)。
当温度继续变化，经过多次快门更新后，得到三个连续相同间隔温度下生成的Bs[ij](快门时，得到的原始本体输出值)后，进入常规的算法计算阶段，在这个阶段包括斜率Bk值更新和本体输出值Bc值的计算两部分，可参阅图4。
1、Bk值更新
当温度变化ΔT(Tc-T[n])触发Tth(触发温度阈值)时，快门本体输出值Bs值更新。Tth初始值一般设定比较小，刚开始的时候，比较快的进行Bk更新，逐步加大到目标值。
本体输出值Bc值通过快门本体输出值Bs值更新，温度T通过读取探测器内部温度数据更新。
需要计算更新的是Bk值。
凡是触发Bs更新，其前后(相邻的两个工作温度区间)的Bk都需要更新，举例来说，触发Bs[n+2](为Tn+2时刻的本体输出值)更新，
先更新前一个Bk[n+1](Bk[n+1]为Tn+1-Tn+2工作温度区间的斜率)：
Bk[n+1]＝(Bs[n+2]–Bs[n+1])/(T[n+2]–T[n+1])；
然后更新Bk[n+2]`(Bk[n+2]`为Tn+2-Tn+3工作温度区间的斜率)：
Bk[n+2]`＝Bk[n+1]+α*(Bk[n+1]-Bk[n])(α为预设斜率调整值)。
2、非制冷红外焦平面探测器本体输出值输出计算：
探测器温度在T[n+2]～T[n+3]或者T[n]～T[n-1]：
用预测的斜率值计算本体输出值：
B＝Bs[n+2]+Bk[n+2]`*(Tc–T[n+2])
B＝Bs[n]+BK[n-1]`*(Tc–T[n])。
探测器温度在T[n+2]～T[n]：
用已生成的Bk和Bs进行计算
B＝Bs[n]+BK[n](Tc–T[n])。
由上可知，本发明实施例不依赖场景的运动，不需要进行标定，不限制非制冷红外焦平面探测器的工作环境温度，且有利于提高非制冷红外焦平面探测器的量产效率，普适性强。
本发明实施例还针对非制冷红外焦平面探测器非均匀性本体输出校正方法提供了相应的实现装置，进一步使得所述方法更具有实用性。下面对本发明实施例提供的非制冷红外焦平面探测器非均匀性本体输出校正装置进行介绍，下文描述的非制冷红外焦平面探测器非均匀性本体输出校正装置与上文描述的非制冷红外焦平面探测器非均匀性本体输出校正方法可相互对应参照。
参见图9，图9为本发明实施例提供的非制冷红外焦平面探测器非均匀性本体输出校正装置在一种具体实施方式下的结构图，该装置可包括：
快门更新信息获取模块901，用于获取非制冷红外焦平面探测器的快门在三个连续相同间隔温度点触发更新时的第一本体输出值和相对应的第一温度值、第二本体输出值和对应的第二温度值、第三本体输出值和对应的第三温度值；
工作区间判断模块902，用于判断非制冷红外焦平面探测器的当前工作温度所处的工作温度区间；第一工作温度区间小于第一温度值，第二工作温度区间大于第一温度值且小于第二温度值，第三工作温度区间大于第二温度值且小于第三温度值，第四工作温度区间大于第三温度值；
本体输出值计算模块903，用于根据该工作温度区间对应的斜率计算当前工作温度的本体输出值；其中，根据第一本体输出值、第一温度值、第二本体输出值、第二温度值计算第二温度区间的第二斜率，根据第二本体输出值、第二温度值、第三本体输出值、第三温度值计算第三温度区间的第三斜率；根据第二斜率、第三斜率和预设斜率调整值预测第一工作温度区间和第四工作温度区间的斜率。
可选的，在本实施例的一些实施方式中，请参阅图10，所述装置例如还可以包括斜率更新模块904，所述斜率更新模块904包括：
获取子模块，用于获取非制冷红外焦平面探测器的快门在当前工作温度为第四温度值时，触发更新的第四本体输出值；第四温度值为第三温度值和触发温度阈值之和；
计算子模块，用于根据第三本体输出值、第三温度值、第四本体输出值及第四温度值计算第四斜率；
更新子模块，用于利用第四斜率更新利用第二斜率、第三斜率和预设斜率调整值预测得到的斜率，以作为第四工作温度区间的斜率；其中，第四工作温度区间大于第三温度值且小于第四温度值。
具体的，还可包括斜率预测更新模块905，所述斜率预测更新模块905可包括：
计算子模块，用于根据第三斜率、第四斜率和预设斜率调整值预测第五工作温度区间的斜率；第五工作温度区间大于第四温度值；
相应的，根据下述公式计算第五预测斜率，以作为第五工作温度区间的斜率：
Bk(n+3)′＝Bk(n+2)+(Bk(n+2)-Bk(n+1))*α；
其中，Bk(n+3)'为第五预测斜率，Bk(n+1)为第三斜率，α为预设斜率调整值、Bk(n+2)为第四斜率。
可选的，在本实施例的另一些实施方式中，所述快门更新信息获取模块901例如还可以包括启动阶段计算子模块，启动阶段计算子模块具体可包括：
获取单元，非制冷红外焦平面探测器处于探测器启动阶段，获取其在快门第一次触发更新时的第一本体输出值和相对应的第一温度值；
计算单元，当前工作温度处于第二工温度作区间，相应的，根据该工作温度区间对应的斜率计算当前工作温度的本体输出值包括：
根据下述公式计算非制冷红外焦平面探测器的当前工作温度处于第二工作温度区间的本体输出值：
B＝Bs(0)+Bk*(Tc-T0)；
式中，Bc为Tc对应的本体输出值，Bs(0)为第一本体输出值，T0为第一温度值，Bk为预设经验调整值，Tc为当前工作温度。
在一些具体的实施方式中，所述本体输出值计算模块903具体可包括：
第一预测斜率计算子模块，用于根据下述公式计算第一预测斜率，以作为第一工作温度区间的斜率：
Bk(n-1)′＝Bk(n)-(Bk(n+1)-Bk(n))*α；
第四预测斜率计算子模块，用于根据下述公式计算第四预测斜率，以作为第四工作温度区间的斜率：
Bk(n+2)′＝Bk(n+1)+(Bk(n+1)-Bk(n))*α；
其中，Bk(n-1)'为第一预测斜率，Bk(n)为第二斜率，Bk(n+1)为第三斜率，α为预设斜率调整值、Bk(n+2)'为第四预测斜率。
第一工作温度区间本体输出值计算子模块，用于若当前工作温度位于第一工作温度区间，根据下述公式计算当前工作温度对应的本体输出值：
Bc(n-1)＝Bs(n)-Bk(n-1)*(Tc-Tn)；
第二工作温度区间本体输出值计算子模块，用于若当前工作温度位于第二工作温度区间，根据下述任意一个公式计算当前工作温度对应的本体输出值：
Tc＜Tn+1-Td，Bc(n)＝Bs(n+1)+Bk(n)*(Tc-Tn+1)；或
Tc＞Tn+Td，Bc(n)＝Bs(n)+Bk(n)*(Tc-Tn)；
第三工作温度区间本体输出值计算子模块，用于若当前工作温度位于第三工作温度区间，根据下述任意一个公式计算当前工作温度对应的本体输出值：
Tc＜Tn+2-Td，Bc(n+1)＝Bs(n+2)+Bk(n+1)*(Tc-Tn+2)；或
Tc＞Tn+1+Td，Bc(n+1)＝Bs(n+1)+Bk(n+1)*(Tc-Tn+1)；
第四工作温度区间本体输出值计算子模块，用于若当前工作温度位于第四工作温度区间，根据下述公式计算当前工作温度对应的本体输出值：
Bc(n+2)＝Bs(n+2)+Bk(n+2)*(Tc-Tn+2)；
其中，Bs(n)为第一本体输出值，Tn为第一温度值，Bs(n+1)为第二本体输出值，Tn+1为第二温度值，Bs(n+2)为第三本体输出值，Tn+2为第三温度值，Tc为当前工作温度，Bk(n-1)为第一工作温度区间的斜率，Bk(n)为第二斜率，Bk(n+1)为第三斜率，Bk(n+2)为第四工作温度区间的斜率，Bc(n-1)为当前工作温度位于第一工作温度区间的本体输出值，Bc(n)为当前工作温度位于第二工作温度区间的本体输出值，Bc(n+1)为当前工作温度位于第三工作温度区间的本体输出值，Bc(n+2)为当前工作温度位于第四工作温度区间的本体输出值；Td为触发温度阈值。
本发明实施例所述非制冷红外焦平面探测器非均匀性本体输出校正装置的各功能模块的功能可根据上述方法实施例中的方法具体实现，其具体实现过程可以参照上述方法实施例的相关描述，此处不再赘述。
由上可知，本发明实施例不依赖场景的运动，不需要进行标定，不限制非制冷红外焦平面探测器的工作环境温度，且有利于提高非制冷红外焦平面探测器的量产效率，普适性强。
本发明实施例还提供了一种非制冷红外焦平面探测器非均匀性本体输出校正设备，具体可包括：
存储器，用于存储计算机程序；
处理器，用于执行计算机程序以实现如上任意一实施例所述非制冷红外焦平面探测器非均匀性本体输出校正方法的步骤。
本发明实施例所述非制冷红外焦平面探测器非均匀性本体输出校正设备的各功能模块的功能可根据上述方法实施例中的方法具体实现，其具体实现过程可以参照上述方法实施例的相关描述，此处不再赘述。
由上可知，本发明实施例不依赖场景的运动，不需要进行标定，不限制非制冷红外焦平面探测器的工作环境温度，且有利于提高非制冷红外焦平面探测器的量产效率，普适性强。
本发明实施例还提供了一种计算机可读存储介质，存储有非制冷红外焦平面探测器非均匀性本体输出校正程序，所述非制冷红外焦平面探测器非均匀性本体输出校正程序被处理器执行时如上任意一实施例所述非制冷红外焦平面探测器非均匀性本体输出校正方法的步骤。
本发明实施例所述计算机可读存储介质的各功能模块的功能可根据上述方法实施例中的方法具体实现，其具体实现过程可以参照上述方法实施例的相关描述，此处不再赘述。
由上可知，本发明实施例不依赖场景的运动，不需要进行标定，不限制非制冷红外焦平面探测器的工作环境温度，且有利于提高非制冷红外焦平面探测器的量产效率，普适性强。
本说明书中各个实施例采用递进的方式描述，每个实施例重点说明的都是与其它实施例的不同之处，各个实施例之间相同或相似部分互相参见即可。对于实施例公开的装置而言，由于其与实施例公开的方法相对应，所以描述的比较简单，相关之处参见方法部分说明即可。
专业人员还可以进一步意识到，结合本文中所公开的实施例描述的各示例的单元及算法步骤，能够以电子硬件、计算机软件或者二者的结合来实现，为了清楚地说明硬件和软件的可互换性，在上述说明中已经按照功能一般性地描述了各示例的组成及步骤。这些功能究竟以硬件还是软件方式来执行，取决于技术方案的特定应用和设计约束条件。专业技术人员可以对每个特定的应用来使用不同方法来实现所描述的功能，但是这种实现不应认为超出本发明的范围。
结合本文中所公开的实施例描述的方法或算法的步骤可以直接用硬件、处理器执行的软件模块，或者二者的结合来实施。软件模块可以置于随机存储器(RAM)、内存、只读存储器(ROM)、电可编程ROM、电可擦除可编程ROM、寄存器、硬盘、可移动磁盘、CD-ROM、或技术领域内所公知的任意其它形式的存储介质中。
以上对本发明所提供的一种非制冷红外焦平面探测器非均匀性本体输出校正方法、装置、设备及计算机可读存储介质进行了详细介绍。本文中应用了具体个例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想。应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以对本发明进行若干改进和修饰，这些改进和修饰也落入本发明权利要求的保护范围内。

Claims (10)

1.一种非制冷红外焦平面探测器非均匀性本体输出校正方法，其特征在于，包括：
获取非制冷红外焦平面探测器的快门在三个连续相同间隔温度点触发更新时的第一本体输出值和相对应的第一温度值、第二本体输出值和对应的第二温度值、第三本体输出值和对应的第三温度值；
判断所述非制冷红外焦平面探测器的当前工作温度所处的工作温度区间，并根据该工作温度区间对应的斜率计算当前工作温度的本体输出值；
其中，第一工作温度区间小于所述第一温度值，第二工作温度区间大于所述第一温度值且小于所述第二温度值，第三工作温度区间大于所述第二温度值且小于所述第三温度值，第四工作温度区间大于所述第三温度值；
根据所述第一本体输出值、所述第一温度值、所述第二本体输出值、所述第二温度值计算所述第二温度区间的温度-本体输出值直线的第二斜率，根据所述第二本体输出值、所述第二温度值、所述第三本体输出值、所述第三温度值计算所述第三温度区间的第三斜率；根据所述第二斜率、所述第三斜率和预设斜率调整值预测所述第一工作温度区间和所述第四工作温度区间的斜率。
2.根据权利要求1所述的非制冷红外焦平面探测器非均匀性本体输出校正方法，其特征在于，所述根据所述第二斜率、所述第三斜率和预设斜率调整值预测所述第一工作温度区间和所述第四工作温度区间的斜率包括：
根据下述公式计算第一预测斜率，以作为所述第一工作温度区间的斜率：
Bk(n-1)'＝Bk(n)-(Bk(n+1)-Bk(n))*α；
根据下述公式计算第四预测斜率，以作为所述第四工作温度区间的斜率：
Bk(n+2)'＝Bk(n+1)+(Bk(n+1)-Bk(n))*α；
其中，Bk(n-1)'为所述第一预测斜率，Bk(n)为所述第二斜率，Bk(n+1)为所述第三斜率，α为预设斜率调整值，Bk(n+2)'为所述第四预测斜率。
3.根据权利要求1所述的非制冷红外焦平面探测器非均匀性本体输出校正方法，其特征在于，所述根据该工作温度区间对应的斜率计算当前工作温度的本体输出值包括：
若当前工作温度位于所述第一工作温度区间，根据下述公式计算当前工作温度对应的本体输出值：
Bc(n-1)＝Bs(n)-Bk(n-1)*(Tc-Tn)；
若当前工作温度位于所述第二工作温度区间，根据下述任意一个公式计算当前工作温度对应的本体输出值：
Tc＜Tn+1-Td，Bc(n)＝Bs(n+1)+Bk(n)*(Tc-Tn+1)；或
Tc＞Tn+Td，Bc(n)＝Bs(n)+Bk(n)*(Tc-Tn)；
若当前工作温度位于所述第三工作温度区间，根据下述任意一个公式计算当前工作温度对应的本体输出值：
Tc＜Tn+2-Td，Bc(n+1)＝Bs(n+2)+Bk(n+1)*(Tc-Tn+2)；或
Tc＞Tn+1+Td，Bc(n+1)＝Bs(n+1)+Bk(n+1)*(Tc-Tn+1)；
若当前工作温度位于所述第四工作温度区间，根据下述公式计算当前工作温度对应的本体输出值：
Bc(n+2)＝Bs(n+2)+Bk(n+2)*(Tc-Tn+2)；
其中，Bs(n)为所述第一本体输出值，Tn为所述第一温度值，Bs(n+1)为所述第二本体输出值，Tn+1为所述第二温度值，Bs(n+2)为所述第三本体输出值，Tn+2为所述第三温度值，Tc为当前工作温度，Bk(n-1)为所述第一工作温度区间的斜率，Bk(n)为所述第二斜率，Bk(n+1)为所述第三斜率，Bk(n+2)为所述第四工作温度区间的斜率，Bc(n-1)为当前工作温度位于所述第一工作温度区间的本体输出值，Bc(n)为当前工作温度位于所述第二工作温度区间的本体输出值，Bc(n+1)为当前工作温度位于所述第三工作温度区间的本体输出值，Bc(n+2)为当前工作温度位于所述第四工作温度区间的本体输出值，Td为触发温度阈值。
4.根据权利要求1所述的非制冷红外焦平面探测器非均匀性本体输出校正方法，其特征在于，所述获取非制冷红外焦平面探测器的快门在连续三个时间段触发更新时的第一本体输出值和相对应的第一温度值、第二本体输出值和对应的第二温度值、第三本体输出值和对应的第三温度值包括：
所述非制冷红外焦平面探测器处于探测器启动阶段，获取其在快门第一次触发更新时的第一本体输出值和相对应的第一温度值；
所述当前工作温度处于第二工温度作区间，相应的，所述根据该工作温度区间对应的斜率计算当前工作温度的本体输出值包括：
根据下述公式计算所述非制冷红外焦平面探测器的当前工作温度处于所述第二工作温度区间的本体输出值：
B＝Bs(0)+Bk*(Tc-T0)；
式中，B为Tc对应的本体输出值，Bs(0)为所述第一本体输出值，T0为所述第一温度值，Bk为预设经验调整值，Tc为当前工作温度。
5.根据权利要求1至4任意一项所述的非制冷红外焦平面探测器非均匀性本体输出校正方法，其特征在于，还包括：
获取所述非制冷红外焦平面探测器的快门在当前工作温度为第四温度值时，触发更新的第四本体输出值；所述第四温度值为所述第三温度值和触发温度阈值之和；
根据所述第三本体输出值、所述第三温度值、所述第四本体输出值及所述第四温度值计算第四斜率；
利用所述第四斜率更新利用所述第二斜率、所述第三斜率和预设斜率调整值预测得到的斜率，以作为所述第四工作温度区间的斜率；
其中，所述第四工作温度区间大于所述第三温度值且小于所述第四温度值。
6.根据权利要求5所述的非制冷红外焦平面探测器非均匀性本体输出校正方法，其特征在于，在所述获取所述非制冷红外焦平面探测器的快门在当前工作温度为第四温度值时，触发更新的第四本体输出值之后，还包括：
根据所述第三斜率、所述第四斜率和所述预设斜率调整值预测第五工作温度区间的斜率；所述第五工作温度区间大于所述第四温度值；
相应的，根据下述公式计算第五预测斜率，以作为所述第五工作温度区间的斜率：
Bk(n+3)'＝Bk(n+2)+(Bk(n+2)-Bk(n+1))*α；
其中，Bk(n+3)'为所述第五预测斜率，Bk(n+1)为所述第三斜率，α为预设斜率调整值，Bk(n+2)为所述第四斜率。
7.一种非制冷红外焦平面探测器非均匀性本体输出校正装置，其特征在于，包括：
快门更新信息获取模块，用于获取非制冷红外焦平面探测器的快门在三个连续相同间隔温度点触发更新时的第一本体输出值和相对应的第一温度值、第二本体输出值和对应的第二温度值、第三本体输出值和对应的第三温度值；
工作区间判断模块，用于判断所述非制冷红外焦平面探测器的当前工作温度所处的工作温度区间；第一工作温度区间小于所述第一温度值，第二工作温度区间大于所述第一温度值且小于所述第二温度值，第三工作温度区间大于所述第二温度值且小于所述第三温度值，第四工作温度区间大于所述第三温度值；
本体输出值计算模块，用于根据该工作温度区间对应的斜率计算当前工作温度的本体输出值；其中，根据所述第一本体输出值、所述第一温度值、所述第二本体输出值、所述第二温度值计算所述第二温度区间的第二斜率，根据所述第二本体输出值、所述第二温度值、所述第三本体输出值、所述第三温度值计算所述第三温度区间的第三斜率；根据所述第二斜率、所述第三斜率和预设斜率调整值预测所述第一工作温度区间和所述第四工作温度区间的斜率。
8.根据权利要求7所述的非制冷红外焦平面探测器非均匀性本体输出校正装置，其特征在于，还包括斜率更新模块，所述斜率更新模块包括：
获取子模块，用于获取所述非制冷红外焦平面探测器的快门在当前工作温度为第四温度值时，触发更新的第四本体输出值；所述第四温度值为所述第三温度值和触发温度阈值之和；
计算子模块，用于根据所述第三本体输出值、所述第三温度值、所述第四本体输出值及所述第四温度值计算第四斜率；
更新子模块，用于利用所述第四斜率更新利用所述第二斜率、所述第三斜率和预设斜率调整值预测得到的斜率，以作为所述第四工作温度区间的斜率；其中，所述第四工作温度区间大于所述第三温度值且小于所述第四温度值。
9.一种非制冷红外焦平面探测器非均匀性本体输出校正设备，其特征在于，包括处理器，所述处理器用于执行存储器中存储的计算机程序时实现如权利要求1至6任一项所述非制冷红外焦平面探测器非均匀性本体输出校正方法的步骤。
10.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质上存储有非制冷红外焦平面探测器非均匀性本体输出校正程序，所述非制冷红外焦平面探测器非均匀性本体输出校正程序被处理器执行时实现如权利要求1至6任一项所述非制冷红外焦平面探测器非均匀性本体输出校正方法的步骤。
CN201811102187.0A 2018-09-20 2018-09-20 一种红外焦平面探测器非均匀性本体输出校正方法及装置 Active CN109084899B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811102187.0A CN109084899B (zh) 2018-09-20 2018-09-20 一种红外焦平面探测器非均匀性本体输出校正方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811102187.0A CN109084899B (zh) 2018-09-20 2018-09-20 一种红外焦平面探测器非均匀性本体输出校正方法及装置

Publications (2)

Publication Number Publication Date
CN109084899A CN109084899A (zh) 2018-12-25
CN109084899B true CN109084899B (zh) 2020-01-10

Family

ID=64842052

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811102187.0A Active CN109084899B (zh) 2018-09-20 2018-09-20 一种红外焦平面探测器非均匀性本体输出校正方法及装置

Country Status (1)

Country Link
CN (1) CN109084899B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109839196A (zh) * 2019-02-22 2019-06-04 烟台艾睿光电科技有限公司 一种红外焦平面探测器的测试板及测试系统
CN111556252B (zh) * 2020-05-22 2021-12-21 烟台艾睿光电科技有限公司 一种快门校正的方法、装置及电子设备

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4209285B2 (ja) * 2003-08-07 2009-01-14 日本電信電話株式会社 温度検出方法及び温度検出器
CN103604503B (zh) * 2013-06-17 2016-07-20 中国航天科工集团第三研究院第八三五八研究所 一种积分时间动态调整的高精度非均匀性校正方法
CN104580894B (zh) * 2014-12-25 2017-10-31 烟台艾睿光电科技有限公司 一种红外焦平面的多点校正方法及系统

Also Published As

Publication number Publication date
CN109084899A (zh) 2018-12-25

Similar Documents

Publication Publication Date Title
US9930324B2 (en) Time based offset correction for imaging systems
Budzier et al. Calibration of uncooled thermal infrared cameras
CN104964746B (zh) 基于无快门的非制冷红外成像系统和非均匀性校正方法
US20160065848A1 (en) Thermography for a thermal imaging camera
US20160198102A1 (en) Time based offset correction for imaging systems and adaptive calibration control
US20170179912A1 (en) Gain calibration for an imaging system
CN103308178A (zh) 一种非制冷红外焦平面阵列的非均匀性校正方法
CN105841821B (zh) 基于定标的无挡片的非均匀性校正装置及其方法
CN108871590B (zh) 非制冷红外焦平面探测器非均匀性响应率校正方法及装置
CN109084899B (zh) 一种红外焦平面探测器非均匀性本体输出校正方法及装置
CN109540297B (zh) 基于fpa温度的红外热像仪标定方法
US8823800B2 (en) Method and apparatus for the correction of a non-uniform sensitivity of detector elements in thermographic cameras
CN102687496A (zh) 无需温度调节用于对检测器输出的图像进行校正的方法以及实现该方法的检测器
US20180276803A1 (en) Image adjustment based on locally flat scenes
CN114018416B (zh) 一种热红外相机辐射定标方法及系统
US8624178B2 (en) Method of correcting sensitivity and matrix image sensor for implementing this method
CN116244662B (zh) 一种多源高程数据融合方法、装置、计算机设备及介质
CN110536085B (zh) 一种读出电路及图像校正方法
Huawei et al. An adaptive two-point non-uniformity correction algorithm based on shutter and its implementation
KR101854355B1 (ko) 다중센서를 선택적으로 이용하는 영상보정 장치
CN114518175B (zh) 一种红外热成像图像的温度矫正方法以及相关装置
CN107782447B (zh) 成像光谱仪光谱定标中空间维自动识别方法及系统
CN115734089A (zh) 半导体光电传感器暗电流的校正方法和系统
CN114441043A (zh) 非制冷红外测温仪的温度修正方法、测温仪
Zhou et al. Realized multipoint method for real-time correction of nonuniformity of uncooled IRFPA

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
PE01 Entry into force of the registration of the contract for pledge of patent right

Denomination of invention: A non-uniformity body output correction method and device for infrared focal plane detector

Effective date of registration: 20211230

Granted publication date: 20200110

Pledgee: Yantai Branch of China Merchants Bank Co.,Ltd.

Pledgor: Yantai Airui Photo-Electric Technology Co.,Ltd.

Registration number: Y2021980017076

PE01 Entry into force of the registration of the contract for pledge of patent right
PC01 Cancellation of the registration of the contract for pledge of patent right

Date of cancellation: 20230106

Granted publication date: 20200110

Pledgee: Yantai Branch of China Merchants Bank Co.,Ltd.

Pledgor: Yantai Airui Photo-Electric Technology Co.,Ltd.

Registration number: Y2021980017076

PC01 Cancellation of the registration of the contract for pledge of patent right
PE01 Entry into force of the registration of the contract for pledge of patent right

Denomination of invention: A method and device for correcting nonuniformity bulk output of infrared focal plane detector

Effective date of registration: 20230113

Granted publication date: 20200110

Pledgee: Yantai Branch of China Merchants Bank Co.,Ltd.

Pledgor: Yantai Airui Photo-Electric Technology Co.,Ltd.

Registration number: Y2023980031040

PE01 Entry into force of the registration of the contract for pledge of patent right
PC01 Cancellation of the registration of the contract for pledge of patent right

Granted publication date: 20200110

Pledgee: Yantai Branch of China Merchants Bank Co.,Ltd.

Pledgor: Yantai Airui Photo-Electric Technology Co.,Ltd.

Registration number: Y2023980031040

PC01 Cancellation of the registration of the contract for pledge of patent right
EE01 Entry into force of recordation of patent licensing contract

Application publication date: 20181225

Assignee: INFIRAY TECHNOLOGIES CO.,LTD.

Assignor: Yantai Airui Photo-Electric Technology Co.,Ltd.

Contract record no.: X2024980006380

Denomination of invention: A non-uniform body output correction method and device for infrared focal plane detectors

Granted publication date: 20200110

License type: Common License

Record date: 20240530

EE01 Entry into force of recordation of patent licensing contract