CN110702231A - 处理图像的方法 - Google Patents

处理图像的方法 Download PDF

Info

Publication number
CN110702231A
CN110702231A CN201910618740.4A CN201910618740A CN110702231A CN 110702231 A CN110702231 A CN 110702231A CN 201910618740 A CN201910618740 A CN 201910618740A CN 110702231 A CN110702231 A CN 110702231A
Authority
CN
China
Prior art keywords
mask
bolometer
pix
masked
bolometers
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
Application number
CN201910618740.4A
Other languages
English (en)
Inventor
L.奇西
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.)
Schneider Electric SE
Schneider Electric Industries SAS
Original Assignee
Schneider Electric SE
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 Schneider Electric SE filed Critical Schneider Electric SE
Publication of CN110702231A publication Critical patent/CN110702231A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01JMEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J5/00Radiation pyrometry, e.g. infrared or optical thermometry
    • G01J5/10Radiation pyrometry, e.g. infrared or optical thermometry using electric radiation detectors
    • G01J5/20Radiation pyrometry, e.g. infrared or optical thermometry using electric radiation detectors using resistors, thermistors or semiconductors sensitive to radiation, e.g. photoconductive devices
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01JMEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J5/00Radiation pyrometry, e.g. infrared or optical thermometry
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01JMEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J5/00Radiation pyrometry, e.g. infrared or optical thermometry
    • G01J5/80Calibration
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/70Denoising; Smoothing
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01JMEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J5/00Radiation pyrometry, e.g. infrared or optical thermometry
    • G01J2005/0077Imaging
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01JMEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J5/00Radiation pyrometry, e.g. infrared or optical thermometry
    • G01J5/10Radiation pyrometry, e.g. infrared or optical thermometry using electric radiation detectors
    • G01J5/20Radiation pyrometry, e.g. infrared or optical thermometry using electric radiation detectors using resistors, thermistors or semiconductors sensitive to radiation, e.g. photoconductive devices
    • G01J2005/202Arrays
    • 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)
  • Spectroscopy & Molecular Physics (AREA)
  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Photometry And Measurement Of Optical Pulse Characteristics (AREA)
  • Radiation Pyrometers (AREA)
  • Transforming Light Signals Into Electric Signals (AREA)

Abstract

本发明涉及一种处理由检测器的测辐射热计收集的原始测量的方法,该方法包括:a)计算掩蔽项,该计算包括从给定列的每个掩蔽的测辐射热计的原始测量中减去所述列的掩蔽的测辐射热计的原始测量的平均值;b)使掩蔽项与掩蔽的测辐射热计的校准分量之间相关,根据以下关系确定:
Figure DDA0002124807580000011
所述项是在参考温度下并用掩模来掩蔽检测器获得的掩蔽的测辐射热计参考测量;c)校正原始图像的步骤,该步骤包括基于相关步骤b)的结果来计算每个测辐射热计的校正图像的校正测量。

Description

处理图像的方法
技术领域
本发明涉及一种处理由具有测辐射热计矩阵的检测器收集的原始图像 的方法。特别地,根据本发明的处理方法旨在校正由于检测器的测辐射热计 的特性的分散而导致的不均匀性。
本发明还涉及一种能够实现根据本发明的方法的所有步骤的计算机程 序。
本发明最后涉及一种包括测辐射热计和实现该计算机程序的计算单元 的检测器。
背景技术
通常,现有技术中已知的红外检测器包括根据n行和m列以矩阵方式布 置的测辐射热计。
当它们通过视域暴露于场景以获取图像时,这些对场景温度敏感的测辐 射热计会经历其电阻的变化。换句话说,在每个测辐射热计中流动的电流取 决于场景温度,并且还取决于环境温度。
特别地,测辐射热计矩阵的测辐射热计Pix_(i,j)的测量S_(i,j)根据以下 定律演变:
Figure BDA0002124807560000011
其中:
-Tamb是环境温度,更特别的是成像仪的温度;
-Tscene是测辐射热计观测到的场景的温度;
-Resp(Tamb)是测辐射热计的响应度,其取决于环境温度;
-S0,Tamb是对于场景温度等于环境温度的测辐射热计的输出值。
项Resp(Tamb)取决于用于测辐射热计的设计的材料及其架构。
场景温度的推导要求环境温度Tamb是已知的,因此检测器通常还设置有 温度传感器。
检测器还可以设置有额外的测辐射热计,即所谓的盲测辐射热计,它们 不暴露于场景。在其中流动的电流则仅取决于环境温度。
因此,根据这样的配置,暴露的测辐射热计的电阻变化的确定基于在所 述暴露的测辐射热计和盲测辐射热计中流动的电流之间的差分测量。
通常,测辐射热计矩阵的每列与盲测辐射热计相关,该盲测辐射热计在 差分测量期间针对所述列的每个测辐射热计实施。然而,可以设想其他配置, 尤其是用于若干列测辐射热计的单个盲测辐射热计的公共汇集。
能够用这种装置获得的场景的原始图像(图1)通常是不可利用的,并 且需要额外的处理。
特别地,图1中所示的图像揭示了检测器的测辐射热计的定位,更具体 地说是非均匀性(“像素化效应”)。这种效应起源于在各测辐射热计中的电 阻的显著分散。
该图像还呈现出列状外观,这是由于盲测辐射热计之间的电阻的分散。
为了减轻这些问题,可以设想各种解决方案。
特别地,可以提出在检测器上实施机械快门。特别地,机械快门被放置 在检测器的前面,以便收集与环境温度有关的参考图像,该参考图像随后从 场景的图像中减去。
然而,这种布置原则上相对简单,但并不令人满意。
实际上,快门的实施以及与之相关的机动化都会带来成本和体积问题。
此外,一旦环境温度变化,必须刷新参考图像。
可替代地,已经提出表征检测器的温度响应,特别是其每个测辐射热计 的温度响应。
该表征包括在各种温度下的参考测量,其中检测器的所有测辐射热计都 被快门遮住。
然后,参考测量使得可以确定每个测辐射热计的温度变化,从而构建保 持在检测器的存储空间中的校准表。
因此,当操作时,对于每个测辐射热计,检测器通过减去由基于校准表 的内插件获得的值来校正原始图像。
然而,使得可以减小各测辐射热计的图像的不均匀性效应的该解决方案 并不令人满意。
实际上,获取参考测量的过程很漫长,并且在检测器的制造中产生成本 超支。
此外,专用于保存校准表的存储空间(考虑到与其相关的成本)是不希 望的。
最后,在说明书末尾引用的文献[1]和[2]中提出了基于使得可以校正图 像的不均匀性的算法的第三种方法。
现有技术中已知的这些方法也不令人满意。
实际上,这些方法通常很难实施,并且它们的稳健性是有争议的。
而且,这些方法需要实施笨重的计算装置,这相应地以它们实施其中的 检测器的成本为代价。
本发明的一个目的是提出一种处理由检测器收集的图像的方法,该检测 器具有比现有技术中已知的技术更简单的测辐射热计,并且不需要实施机械 部件以及与之相关的机器自动化。
本发明的另一个目的是提出一种方法,使得可以校正在原始图像上观察 到的列状效果。
发明内容
本发明的目的至少部分地通过一种由计算单元执行的用于处理原始图 像的方法实现,其特征在于,原始测量S_(i,j)由检测器的并且根据n行和m 列以矩阵方式布置的测辐射热计Pix_(i,j)收集,所述测辐射热计Pix_(i,j)包括 暴露的测辐射热计和被遮住的掩蔽的测辐射热计,该方法包括:
a)计算与掩蔽的测辐射热计相关的掩蔽项S_mask_norm_(i,j),该计算 包括从给定列的每个掩蔽测辐射热计的原始测量S_(i,j)中减去所述列的掩蔽 的测辐射热计的原始测量S_(i,j)的平均值;
b)使掩蔽项S_mask_norm_(i,j)与掩蔽的测辐射热计的校准分量 S_mask_cal_(i,j)之间相关,根据以下关系确定:
Figure BDA0002124807560000031
所述项S_mask_ref_(i,j)是在参考温度下并且通过用也保持在参考温度 的掩蔽件来掩蔽检测器获得的掩蔽的测辐射热计参考测量;
c)校正原始图像的步骤,该步骤包括基于相关步骤b)的结果来计算每 个测辐射热计的校正图像的校正测量S_Cor_(i,j)。
根据一实施方式,检测器包括专用于保存校准分量S_mask_cal_(i,j)的存 储空间。
根据一实施方式,相关步骤b)包括线性回归,使得每个掩蔽项 S_mask_norm_(i,j)和校准分量S_mask_cal_(i,j)满足以下关系:
S_mask_norm_(i,j)=β.S_mask_cal_(i,j)+Doffset
其中β和Doffset是线性回归期间确定的项。
根据一实施方式,校正步骤c)包括根据以下关系计算每个测辐射热计 的校正测量S_Cor_(i,j):
S_Cor_(i,j)=S_(i,j)-β.S_cal_(i,j)-Doffset
其中,对于给定的辐射热计,S_cal_(i,j)是根据以下关系确定的校准分量:
Figure BDA0002124807560000041
项S_ref_(i,j)是在参考温度下并且通过用也保持在参考温度的掩蔽件来 掩蔽检测器的测辐射热计的参考测量。
根据一实施方式,校准分量S_cal_(i,j)保存在检测器的存储空间中。
根据一实施方式,该方法还包括以下步骤:
d)对于列中的每一列,计算对应于所考虑的列的掩蔽的测辐射热计的 测量S_(i,j)的平均值的列项Col_mask_(j);
e)使列项Col_mask_(j)和掩蔽的测辐射热计的校准坐标C_mask_cal_(j) 之间相关,根据以下关系确定:
f)将校正图像进行附加校正以形成最终图像,附加校正包括:
基于校正测量S_Cor_(i,j)并且基于相关步骤e)的结果来计算最终图像 的最终测量S_fin_(i,j)。
根据一实施方式,相关步骤e)包括线性回归,使得列项Col_mask_(j) 和校准坐标C_mask_cal_(j)满足以下关系:
Col_mask_(j)=α.C_mask_cal_(j)+Coloffset
其中α和Coloffset是线性回归期间确定的项。
根据一实施方式,附加校正步骤f)包括根据以下关系计算每个测辐射 热计的最终测量S_fin_(i,j):
S_fin_(i,j)=S_Cor_(i,j)-α.Col_cal_(j)-Coloffset
其中,对于给定列的测辐射热计(Pix_(i,j)),Col_cal_(i,j)是根据以下关系 确定的校准坐标:
Figure BDA0002124807560000051
根据一实施方式,检测器设置有安装在光圈上的透镜,所述光圈在检测 器的角部的水平处遮住掩蔽的测辐射热计。
根据一实施方式,检测器还包括盲测辐射热计,每个盲测辐射热计被实 施用于特定于此的至少一列测辐射热计中的测辐射热计的差分测量,有利 地,每个盲测辐射热计与单列测辐射热计相关。
本发明还涉及一种计算机程序,当其由计算机执行时,其实现本发明的 方法的实施。
本发明还涉及一种装置,包括:
-检测器,其具有根据n行和m列以矩阵方式布置的多个测辐射热计, 所述多个测辐射热计包括暴露于场景的暴露的测辐射热计和被遮住的掩蔽 的测辐射热计,
-计算单元,其配备有根据本发明的计算机程序。
本发明还涉及根据本发明的装置的实施,用于房间中的检测,尤其是人 检测。
附图说明
参考附图,其他特征和优点在通过非限制性示例给出的用于处理图像的 方法的以下描述中将变得显而易见,其中:
-图1是由检测器获得的场景的原始图像,该检测器具有根据80行×80 列的矩阵布置的测辐射热计;
-图2是检测器的示意图,该检测器具有安装在能够根据本发明实施的 光圈上的透镜;
-图3是由与检测器保持在参考温度的快门遮住的检测器收集的参考图 像;
-图4是基于图3的参考图像计算的矢量(所谓的列参考矢量)的分量 的图示;
-图5是基于图3的参考图像和图4的矢量确定的矩阵(所谓的列参考 矩阵)的分量的图示;
-图6是基于图3的参考图像计算的矢量(所谓的掩蔽的参考矢量)的 分量的图示;
-图7是基于图3的参考图像和图6的矢量确定的矩阵(所谓的掩蔽的 参考矩阵)的分量的图示;
-图8是根据掩蔽的校准分量S_mask_cal_(i,j)(水平轴)的掩蔽项 S_mask_norm_(i,j)(在垂直轴上)的值的图表图示;
-图9a和9b分别是根据本发明的原始图像和校正的均匀图像;
-图10是根据掩蔽的校准分量C_mask_cal_(i,j)的掩蔽项 C_mask_norm_(i,j)(在垂直轴上)的值的图表图示;
-图11表示最终图像;
-图12仅示出了图9a的原始图像的列状外观的校正。
具体实施方式
现在将结合图1至11描述本发明。
图2图示了具有根据n行(表示为“Li”)和m列(表示为“Cj”)以矩阵 方式布置的多个测辐射热计(表示为Pix_(i,j))的检测器。
编号i、j的测辐射热计对应于设置在行i与列j的交叉点处的测辐射热 计。
在测辐射热计Pix_(i,j)中,可以区分暴露于要成像的场景的所谓的暴露 的测辐射热计Pix_Exp_(i,j)和测辐射热计被遮住的所谓的掩蔽的测辐射热计 Pix_mask_(i,j)。
掩蔽的测辐射热计Pix_mask_(i,j)例如被光圈遮住,在光圈上安装有设置 在场景和检测器1之间的透镜2。根据这种配置,掩蔽的测辐射热计 Pix_mask_(i,j)设置在检测器的角部。然而,本发明不限于这种布置,并且掩 蔽的测辐射热计Pix_mask_(i,j)可以例如形成完整的测辐射热计列,特别是沿 测辐射热计矩阵的边的列。
图2的检测器还可以包括多个盲测辐射热计Bol_Blind_(j)。
检测器1还包括计算单元4,计算单元4配备有用于执行根据本发明的 方法的各个步骤的计算处理器。计算单元还可以包括用于保存原始、参考或 中间测量的存储空间。
最后,检测器1可包括用于评估所述检测器所处环境的温度的温度探针5。温度探针可以例如包括PN结。
根据本发明的方法提出校正由检测器1收集的原始图像的缺陷。
在这方面,“现有技术”部分中讨论的图1展示出均匀性和列状外观的缺 陷。
不均匀性主要是由于形成检测器的测辐射热计的电阻的分散。
也可以在图像角部的水平处观察到较暗的区域。这些暗区是由于掩蔽的 测辐射热计Pix_mask_(i,j)被其上安装透镜的光圈遮住,并将用于校正在原始 图像上观察到的不均匀性。掩蔽的测辐射热计Pix_mask_(i,j)的原始测量 S_(i,j)在后面提交文件中被表示为掩蔽原始测量S_M_(i,j)。
可以观察到对应于原始测量S_(i,j)(表示暴露的测辐射热计Pix_Exp_(i,j) 的暴露原始测量S_E_(i,j))的更清晰的中心区域。
最后,图1还揭示了列状外观,这是由于盲测辐射热计Bol_Blind_(j)的 电阻的分散。在这方面,每个盲测辐射热计Bol_Blind_(j)被实施用于特定于 此的至少一列测辐射热计中的测辐射热计的差分测量,有利地,每个盲测辐 射热计Bol_Blind_(j)与测辐射热计Pix_(i,j)的单列Cj相关。
要注意的是,如果所有列都与相同的盲测辐射热计相关,则不会观察到 列状效果。然而,这种配置是不希望的,因为该单盲测辐射热计的失效会使 检测器完全不可操作。
根据本发明的方法可以包括当测辐射热计Pix_(i,j)由与检测器保持在温 度(所谓的参考温度Tr,例如在20℃)的快门遮住时确定测辐射热计Pix_(i,j) 的响应S_(i,j)(表示为S_ref_(i,j))。
响应S_Ref_(i,j)可以通过检测器的简单测量获得,其中快门位于测辐射 热计Pix_(i,j)的前面。在这方面,图3图示了由检测器收集的参考图像。
基于这些参考测量确定矢量,即具有分量C_ref_col(j)(对于j范围从1 到m)的所谓的列参考矢量V_ref_col(在图4中图示)。
分量C_ref_col(j)由以下关系确定:
Figure BDA0002124807560000071
换句话说,每个分量C_ref_col(j)是列j的参考测量S_Ref_(i,j)的平均值。
还可以确定每个与测辐射热计Pix_(i,j)相关的校准分量S_col_cal_(i,j)的 矩阵,即所谓的列参考矩阵M_Pix_col_cal(在图5中图示)。
然后,以下列方式计算列j的每个测辐射热计Pix_(i,j)的分量 S_col_cal_(i,j):
S_col_cal_(i,j)=S_ref_(i,j)-C_ref_col_(j)
换句话说,检测器的列j的测辐射热计Pix_(i,j)的校准分量S_col_cal_(i,j) 对应于所述测辐射热计的参考测量,从该参考测量中减去列j的所有测辐射 热计的参考测量S_ref_(i,j)的平均值。
列校准矩阵M_Pix_col_cal,比如随后在提交文件中描述的,有利地实 施为校正由检测器收集的原始图像的列状外观。
补充地或可替代地,仍然基于参考测量S_ref_(i,j),可以确定具有分量 C_mask_cal_(j)的矢量,即所谓的掩蔽的参考矢量V_ref_pix(在图6中图示)。
分量C_mask_cal_(j)由以下关系确定:
Figure BDA0002124807560000081
其中,数据S_mask_ref_(i,j)仅对应于掩蔽的测辐射热计Pix_mask_(i,j) 的响应S_Ref_(i,j)。
换句话说,每个分量C_mask_cal_(j)是与列j的掩蔽的测辐射热计相关 的参考测量S_Ref_(i,j)的平均值。
还可以确定具有每个与掩蔽的测辐射热计Pix_mask_(i,j)相关的分量 S_mask_cal_(i,j)的矩阵,即所谓的掩蔽参考矩阵M_Pix_mask_cal(在图7中 表示)。
针对列j的每个掩蔽的测辐射热计Pix_mask_(i,j)的校准分量 S_mask_cal_(i,j)根据以下公式确定:
S_mask_cal_(i,j)=S_mask-ref_(i,j)-C_mask_cal_(j)
换句话说,检测器的列j的测辐射热计Pix_mask_(i,j)的校准分量 S_mask_cal_(i,j)对应于所述测辐射热计的参考测量,从该参考测量中减去列 j的所有掩蔽的测辐射热计的参考测量的平均值。
掩蔽校准矩阵M_Pix_mask_cal和列状校准矩阵M_Pix_col_cal的建立在 检测器的制造期间实施,并且它们的分量有利地保存在检测器的专用存储空 间中,例如保存在计算单元4的存储空间中。
根据本发明的方法包括步骤a)计算矩阵(所谓的掩蔽矩阵)的项(所 谓的掩蔽项S_mask_norm_(i,j))。
特别地,每个掩蔽项S_mask_norm_(i,j)与掩蔽的测辐射热计 Pix_mask_(i,j)相关。
确定与列Cj的掩蔽的测辐射热计Pix_mask_(i,j)有关的掩蔽项 S_mask_norm_(i,j)包括从掩蔽的原始测量S_M_(i,j)中减去列Cj的掩蔽的原 始测量(表示为S_M_(j))的平均值。
特别地,列j的每个掩蔽项S_mask_norm_(i,j)根据以下关系确定:
Figure BDA0002124807560000091
S_mask_norm_(i,j)=S_M_(i,j)-S_M_(j)
确定掩蔽项S_mask_norm_(i,j)的步骤之后是步骤b),即使所述掩蔽项 S_mask_norm_(i,j)与掩蔽的测辐射热计Pix_mask_(i,j)的掩蔽校准分量 S_mask_cal_(i,j)之间相关。
图8特别针对每个掩蔽的测辐射热计Pix_mask_(i,j)图示了根据掩蔽校准 分量S_mask_cal_(i,j)的掩蔽项S_mask_norm_(i,j)(在垂直轴上)的值。在该 图中,掩蔽项S_mask_norm_(i,j)根据掩蔽的校准分量S_mask_cal_(i,j)线性变 化。
然后,这两个项之间的相关性可以包括根据以下关系的线性回归:
S_mask_norm_(i,j)=β.S_mask_cal_(i,j)+Doffset
其中β和Doffset是线性回归期间确定的项。
然后,这些可以在校正步骤c)期间根据以下关系校正每个测辐射热计 的原始测量S_(i,j):
S_Cor_(i,j)=S_(i,j)-β.S_col_cal_(i,j)-Doffset
图9a和9b图示了上文所述的校正的效果。特别地,图9a是原始图像, 而图9b示出了图像的不均匀性校正的效果。尽管列状外观在图9b中仍是明 显的,但其提供了迄今为止在原始图像上无法检测到的人的闪现。
根据本发明的方法还可以实施列状外观的校正。列状外观的校正可以独 立于上文给出的均匀性校正执行。
此外,根据本发明的方法可以包括在步骤d)期间,对于列Cj中的每一 列,计算对应于所考虑的列的掩蔽的测辐射热计Pix_mask_(i,j)的测量S_(i,j) 的平均值的列项Col_mask_(j)。
然后,在步骤e)期间,可以将这些列项Col_mask_(j)与掩蔽的测辐射 热计Pix_mask_(i,j)的校准坐标C_mask_cal_(j)相关。
图10特别针对每个掩蔽的测辐射热计Pix_mask_(i,j)图示了根据列j的 掩蔽测辐射热计Pix_mask_(i,j)的校准坐标C_mask_cal_(j)的列项 Col_mask_(j)(在垂直轴上)的值。在该图中,项Col_mask_(j)根据校准的坐 标C_mask_cal_(j)线性变化。
相关性可以包括线性回归,使得列项Col_mask_(j)和校准的坐标 C_mask_cal_(j)满足以下关系:
Col_mask_(j)=α.C_mask_cal_(j)+Coloffset
其中α和Coloffset是线性回归期间确定的项。
然后,根据本发明的方法可以包括将校正图像进行附加校正(步骤f) 以形成最终图像。
附加校正尤其包括基于校正测量S_Cor_(i,j)并且基于相关步骤的结果来 计算最终图像的最终测量S_fin_(i,j)。
例如,附加校正包括根据以下关系计算每个测辐射热计Pix_(i,j)的最终 测量S_fin_(i,j):
S_fin_(i,j)=S_Cor_(i,j)-α.Col_cal_(j)-Coloffset
在图11中图示的如此校正的图像不再呈现列状外观。
应当理解,原始图像的列状和非均匀性外观具有不同的起源,因此可以 独立地校正。
特别地,可以反转这些校正的实施顺序。
在这方面,图12仅示出了图9a的原始图像的列状外观的校正。
本发明还涉及一种计算机程序,当该计算机程序由计算机或计算单元执 行时,实现根据本发明的方法的实施。
本发明使得可以以差异化的方式校正具有测辐射热计矩阵的检测器的 原始图像的列状和非均匀性外观。
特别地,本发明仅需要在参考温度下的参考测量以校正原始图像。
本发明可以有利地用于房间中的检测,尤其是人的检测。
参考文献
[1]EP2940991B1
[2]US2010237245A1

Claims (13)

1.一种由计算单元执行的用于处理原始图像的方法,其特征在于,原始测量S_(i,j)由检测器(1)的根据n行(Li)和m列(Cj)以矩阵方式布置的测辐射热计Pix_(i,j)收集,并且,所述测辐射热计Pix_(i,j)包括暴露的测辐射热计(Pix_exp_(i,j))和被遮住的掩蔽的测辐射热计(Pix_mask_(i,j)),该方法包括:
a)计算与掩蔽的测辐射热计(Pix_mask_(i,j))相关的掩蔽项S_mask_norm_(i,j),该计算包括从给定列(Cj)的每个掩蔽的测辐射热计(Pix_mask_(i,j))的原始测量S_(i,j)中减去所述列的掩蔽的测辐射热计的原始测量S_(i,j)的平均值;
b)使掩蔽项S_mask_norm_(i,j)与掩蔽的测辐射热计(Pix_mask_(i,j))的校准分量S_mask_cal_(i,j)之间相关,根据以下关系确定:
Figure FDA0002124807550000011
所述项S_mask_ref_(i,j)是在参考温度(Tr)下并且通过用也保持在参考温度(Tr)的掩蔽件来掩蔽检测器获得的掩蔽测辐射热计参考测量(Pix_mask_(i,j));
c)校正原始图像的步骤,该步骤包括基于相关步骤b)的结果来计算每个测辐射热计(Pix_(i,j)的校正图像的校正测量S_Cor_(i,j)。
2.根据权利要求1所述的方法,其中,所述检测器包括专用于保存校准分量S_mask_cal_(i,j)的存储空间。
3.根据权利要求1或2所述的方法,其中,所述相关步骤b)包括线性回归,使得每个掩蔽项S_mask_norm_(i,j)和校准分量S_mask_cal_(i,j)满足以下关系:
S_mask_norm_(i,j)=β.S_mask_cal_(i,j)+Doffset
其中β和Doffset是线性回归期间确定的项。
4.根据权利要求3所述的方法,其中,所述校正步骤c)包括根据以下关系计算每个测辐射热计(Pix_(i,j))的校正测量S_Cor_(i,j):
S_Cor_(i,j)=S_(i,j)-β.S_cal_(i,j)-Doffset
其中,对于给定的辐射热计(Pix_(i,j)),S_cal_(i,j)是根据以下关系确定的校准分量:
项S_ref_(i,j)是在参考温度(Tr)下并且通过用也保持在参考温度(Tr)的掩蔽件来掩蔽检测器的测辐射热计(Pix_(i,j))的参考测量。
5.根据权利要求4所述的方法,其中,所述校准分量S_cal_(i,j)保存在检测器的存储空间中。
6.根据权利要求4或5所述的方法,其中,所述方法还包括以下步骤:
d)对于列(Cj)中的每一列,计算对应于所考虑的列的掩蔽的测辐射热计Pix_mask_(i,j)的测量S_(i,j)的平均值的列项Col_mask_(j);
e)使列项Col_mask_(j)和掩蔽的测辐射热计(Pix_mask_(i,j))的校准坐标C_mask_cal_(j)之间相关,根据以下关系确定:
Figure FDA0002124807550000022
f)将校正图像附加校正以形成最终图像,附加校正包括:
基于校正测量S_Cor_(i,j)并且基于相关步骤e)的结果来计算最终图像的最终测量S_fin_(i,j)。
7.根据权利要求6所述的方法,其中,所述相关步骤e)包括线性回归,使得列项Col_mask_(j)和校准坐标C_mask_cal_(j)满足以下关系:
Col_mask_(j)=α.C_mask_cal_(j)+Coloffset
其中α和Coloffset是线性回归期间确定的项。
8.根据权利要求7所述的方法,其中,所述附加校正步骤f)包括根据以下关系计算每个测辐射热计(Pix_(i,j))的最终测量S_fin_(i,j):
S_fin_(i,j)=S_Cor_(i,j)-α.Col_cal_(j)-Coloffset
其中,对于给定列(Cj)的测辐射热计(Pix_(i,j)),Col_cal_(i,j)是根据以下关系确定的校准坐标:
9.根据权利要求1至8中任一项所述的方法,其中,所述检测器设置有安装在光圈上的透镜,所述光圈在检测器(1)的角部的水平处遮住掩蔽的测辐射热计(Pix_mask_(i,j))。
10.根据权利要求1至9中任一项所述的方法,其中,所述检测器还包括盲测辐射热计(Bol_Blind_(j)),每个盲测辐射热计(Bol_Blind_(j))被实施用于特定于此的至少一列测辐射热计中的测辐射热计的差分测量,有利地,每个盲测辐射热计(Bol_Blind_(j))与单列(Cj)测辐射热计(Pix_(i,j))相关。
11.一种计算机程序,当其由计算机执行时,实现根据权利要求1至10中任一项所述的方法的实施。
12.一种装置,包括:
-检测器,其具有根据n行(Li)和m列(Cj)以矩阵方式布置的多个测辐射热计,所述多个测辐射热计包括暴露于场景的暴露的测辐射热计(Pix_exp_(i,j))和被遮住的掩蔽的测辐射热计(Pix_mask_(i,j)),
-计算单元,其配备有根据权利要求11所述的计算机程序。
13.一种根据权利要求12所述的装置的实施,用于房间中的检测,尤其是人检测。
CN201910618740.4A 2018-07-10 2019-07-10 处理图像的方法 Pending CN110702231A (zh)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
FR1856355A FR3083901B1 (fr) 2018-07-10 2018-07-10 Methode de traitement d'une image
FR1856355 2018-07-10

Publications (1)

Publication Number Publication Date
CN110702231A true CN110702231A (zh) 2020-01-17

Family

ID=65031405

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910618740.4A Pending CN110702231A (zh) 2018-07-10 2019-07-10 处理图像的方法

Country Status (4)

Country Link
US (1) US20200018652A1 (zh)
EP (1) EP3594643A1 (zh)
CN (1) CN110702231A (zh)
FR (1) FR3083901B1 (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR3111699A1 (fr) 2020-06-23 2021-12-24 Schneider Electric Industries Sas Procede de traitement d’une image brute collectee par un detecteur a bolometres et dispositif associe
FR3134212A1 (fr) 2022-03-30 2023-10-06 Schneider Electric Industries Sas Procede de comptage de personnes et de determination de leur position dans un espace

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070120058A1 (en) * 1995-10-24 2007-05-31 Bae, Inc. Uncooled infrared sensor
CN101458125A (zh) * 2007-12-12 2009-06-17 Ulis股份公司 检测电磁辐射的器件和系统、读取成像测辐射热计的方法
CN102192791A (zh) * 2010-01-27 2011-09-21 索尼公司 图像拾取装置
US20140232875A1 (en) * 2011-06-10 2014-08-21 Flir Systems, Inc. Determination of an absolute radiometric value using blocked infrared sensors

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4795610B2 (ja) * 2000-05-01 2011-10-19 ビーエイイー・システムズ・インフォメーション・アンド・エレクトロニック・システムズ・インテグレイション・インコーポレーテッド 放射センサの温度変動を補償するための方法および装置
US7683321B1 (en) 2008-12-19 2010-03-23 Fluke Corporation Shutterless infrared imager algorithm with drift correction
FR3020735B1 (fr) 2014-04-30 2017-09-15 Ulis Procede de traitement d'une image infrarouge pour une correction des non uniformites

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070120058A1 (en) * 1995-10-24 2007-05-31 Bae, Inc. Uncooled infrared sensor
CN101458125A (zh) * 2007-12-12 2009-06-17 Ulis股份公司 检测电磁辐射的器件和系统、读取成像测辐射热计的方法
CN102192791A (zh) * 2010-01-27 2011-09-21 索尼公司 图像拾取装置
US20140232875A1 (en) * 2011-06-10 2014-08-21 Flir Systems, Inc. Determination of an absolute radiometric value using blocked infrared sensors

Also Published As

Publication number Publication date
FR3083901B1 (fr) 2021-10-08
US20200018652A1 (en) 2020-01-16
FR3083901A1 (fr) 2020-01-17
EP3594643A1 (fr) 2020-01-15

Similar Documents

Publication Publication Date Title
US9883178B2 (en) Method for measuring performance parameters and detecting bad pixels of an infrared focal plane array module
CN108871588B (zh) 一种红外成像系统多维度联合非均匀性校正方法
EP2693740B1 (en) Radiation imaging apparatus
KR20150125609A (ko) 비-균일성 정정을 위한 적외선 이미지 처리 방법
EP1678485B1 (en) Method and ir-camera for determining the risk of condensation
US9696210B2 (en) Extended temperature range mapping process of a furnace enclosure using various device settings
CN109525840B (zh) 一种成像芯片上微弱缺陷的检测方法
KR101158259B1 (ko) 적외선 센서의 신호 검출 회로 및 그 보정방법
CN110631706B (zh) 红外图像的校正方法、装置及存储介质
CN110702231A (zh) 处理图像的方法
CN101889869B (zh) 摄像设备及其控制方法
US8481918B2 (en) System and method for improving the quality of thermal images
CN105869129B (zh) 针对非均匀校正后的热红外图像剩余非均匀噪声去除方法
CN105890873A (zh) 一种热红外高光谱成像仪盲元检测的装置和方法
US9581500B2 (en) Diagnosis of the defective state of a bolometric detection array
CN113138026A (zh) 一种基于测温系统的测温方法、装置、计算机设备
CN110542482B (zh) 一种盲元检测方法、装置及电子设备
CN113096041A (zh) 图像校正方法及装置
CN113837954A (zh) 处理由辐射热计检测器和相关设备收集的原始图像的方法
US11262245B2 (en) Method for processing an image
CN113049118B (zh) 一种红外焦平面探测器闪元测试装置及测试方法
CN114518175B (zh) 一种红外热成像图像的温度矫正方法以及相关装置
Hakim et al. Analysis of radiometric calibration on matrix imager of LAPAN-A3 satellite payload
US20050194539A1 (en) Method for improving measurement accuracy of infrared imaging radiometers
EP4068759A1 (en) Methods for dark signal non-uniformity 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