Movatterモバイル変換


[0]ホーム

URL:


CN118392324A - Infrared radiation measurement method for subtracting internal stray radiation based on background plate - Google Patents

Infrared radiation measurement method for subtracting internal stray radiation based on background plate
Download PDF

Info

Publication number
CN118392324A
CN118392324ACN202410865082.XACN202410865082ACN118392324ACN 118392324 ACN118392324 ACN 118392324ACN 202410865082 ACN202410865082 ACN 202410865082ACN 118392324 ACN118392324 ACN 118392324A
Authority
CN
China
Prior art keywords
background
calibration
background plate
infrared radiation
radiation
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.)
Granted
Application number
CN202410865082.XA
Other languages
Chinese (zh)
Other versions
CN118392324B (en
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.)
Institute of Optics and Electronics of CAS
Original Assignee
Institute of Optics and Electronics 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 Institute of Optics and Electronics of CASfiledCriticalInstitute of Optics and Electronics of CAS
Priority to CN202410865082.XApriorityCriticalpatent/CN118392324B/en
Publication of CN118392324ApublicationCriticalpatent/CN118392324A/en
Application grantedgrantedCritical
Publication of CN118392324BpublicationCriticalpatent/CN118392324B/en
Activelegal-statusCriticalCurrent
Anticipated expirationlegal-statusCritical

Links

Classifications

Landscapes

Abstract

The invention discloses an infrared radiation measurement method for subtracting internal stray radiation based on a background plate, which relates to the field of infrared system radiation measurement, and comprises the following steps: step 1, background plate and blackbody image acquisition strategies: after the temperature of the black body is stable, converting the attenuation sheet, sequentially converting the integral time gear, and collecting the black body image; cutting the black body image to a background plate after the black body image is acquired, and acquiring a background image; step 2, establishing an infrared radiation calibration model with background deduction; and 3, executing an infrared radiation measurement strategy of background subtraction according to the infrared radiation calibration model of background subtraction established in the step 2. According to the technical scheme, the influence of internal stray radiation on the target gray level is reduced by a background subtraction method, so that the measurement accuracy is higher.

Description

Translated fromChinese
一种基于背景板扣除内部杂散辐射的红外辐射测量方法A method for measuring infrared radiation based on background plate subtraction of internal stray radiation

技术领域Technical Field

本发明属于红外系统辐射测量领域,具体涉及一种基于背景板扣除内部杂散辐射的红外辐射测量方法。The invention belongs to the field of infrared system radiation measurement, and in particular relates to an infrared radiation measurement method based on background plate deduction of internal stray radiation.

背景技术Background technique

通过红外辐射特性测量,可获得目标的辐强度、辐照度、辐亮度等辐射量信息。利用这些信息,一方面可建立各类目标的辐射特性数据库,并分析它们的红外辐射特性规律,这对目标的侦察识别有极大帮助;另一方面可对目标进行定量的红外特征提取,特别是精确的红外辐射测量可获取目标准确的温度特征。因此目标的红外辐射测量已变得越来越重要。By measuring the infrared radiation characteristics, we can obtain the radiation information of the target, such as the radiation intensity, irradiance, and radiance. Using this information, on the one hand, we can establish a radiation characteristics database of various targets and analyze their infrared radiation characteristics, which is of great help to the reconnaissance and identification of the target; on the other hand, we can extract the infrared characteristics of the target quantitatively, especially the accurate infrared radiation measurement can obtain the accurate temperature characteristics of the target. Therefore, the infrared radiation measurement of the target has become more and more important.

红外辐射测量系统需要观测各种各样的目标,它们的辐射量变化范围十分大,若测量系统的动态范围不能覆盖目标辐射量变化范围,则会导致低辐射量目标和高辐射量目标测量失败。现有的红外辐射测量设备常采用加装衰减片的方法以实现宽动态范围的测量。红外辐射测量设备的自身热辐射很容易受环境温度的变化而变化,这很容易导致常规辐射定标方法精度降低甚至失效。目前针对该问题所采用的办法需建立红外辐射测量系统输出灰度和环境温度的经验关系,其可靠性和可行性较差。Infrared radiation measurement systems need to observe a variety of targets, and their radiation levels vary greatly. If the dynamic range of the measurement system cannot cover the target radiation level variation range, it will lead to the failure of measuring low-radiation targets and high-radiation targets. Existing infrared radiation measurement equipment often uses the method of adding attenuation plates to achieve wide dynamic range measurement. The thermal radiation of infrared radiation measurement equipment itself is easily affected by changes in ambient temperature, which can easily lead to reduced accuracy or even failure of conventional radiation calibration methods. The current method used to solve this problem requires the establishment of an empirical relationship between the output grayscale of the infrared radiation measurement system and the ambient temperature, which has poor reliability and feasibility.

发明内容Summary of the invention

本发明的目的在于提供一种基于背景板扣除内部杂散辐射的红外辐射测量方法,以力图解决或至少缓解上面提到的问题。The object of the present invention is to provide an infrared radiation measurement method based on background plate deduction of internal stray radiation, in an effort to solve or at least alleviate the above-mentioned problems.

本发明所采用的技术方案为:一种基于背景板扣除内部杂散辐射的红外辐射测量方法,该方法包括:The technical solution adopted by the present invention is: an infrared radiation measurement method based on background plate deduction of internal stray radiation, the method comprising:

步骤1,进行背景板设计;Step 1, design the background board;

步骤2,设计背景板和黑体图像采集策略:待黑体温度稳定后,变换衰减片,依次变换积分时间档位,采集黑体图像;各衰减片采集完后,切至背景板,采集背景板图像;Step 2, design the background plate and black body image acquisition strategy: after the black body temperature is stable, change the attenuation plate, change the integration time gear in sequence, and collect the black body image; after each attenuation plate is collected, switch to the background plate and collect the background plate image;

步骤3,建立背景扣除的红外辐射定标模型;Step 3, establishing a background-subtracted infrared radiation calibration model;

步骤4,设计背景扣除的红外辐射定标算法:用采集的黑体图像和背景板图像处理后的数据对步骤3的红外辐射定标模型进行最优化参数估计,以获得定标增益和定标偏置,定标增益是背景扣除的红外辐射定标模型的斜率,定标偏置是背景扣除的红外辐射定标模型的截距,计算定标误差并判断定标误差是否大于误差阈值,如果大于误差阈值,则红外辐射定标算法继续迭代,直到所有温度点的定标误差都小于误差阈值;Step 4, designing an infrared radiation calibration algorithm with background subtraction: using the collected black body image and the processed data of the background plate image to optimize the parameter estimation of the infrared radiation calibration model in step 3 to obtain the calibration gain and calibration bias, the calibration gain is the slope of the infrared radiation calibration model with background subtraction, the calibration bias is the intercept of the infrared radiation calibration model with background subtraction, calculating the calibration error and judging whether the calibration error is greater than the error threshold, if it is greater than the error threshold, the infrared radiation calibration algorithm continues to iterate until the calibration errors of all temperature points are less than the error threshold;

步骤5,执行背景扣除的红外辐射测量策略:使用步骤4获得的定标增益,获得目标表观辐亮度。Step 5, execute the infrared radiation measurement strategy with background subtraction: use the calibration gain obtained in step 4 to obtain the target apparent radiance.

本发明带来的有益效果:The beneficial effects brought by the present invention are:

(1)方法简单易行:现在的红外辐射测量系统大多都会装滤光轮和镜头盖,根据本发明的技术方案只需在滤光轮上多装一个背景板或对镜头盖做涂层处理即可,无需给系统增加额外的机械结构,其方法简单并且可行性高;(1) The method is simple and easy: Most of the current infrared radiation measurement systems are equipped with a filter wheel and a lens cover. According to the technical solution of the present invention, it is only necessary to install a background plate on the filter wheel or to coat the lens cover. There is no need to add additional mechanical structures to the system. The method is simple and has high feasibility.

(2)测量精度更高:通过扣除背景的方法,降低了内部杂散辐射对目标灰度的影响,使得测量精度更高;(2) Higher measurement accuracy: By subtracting the background, the influence of internal stray radiation on the target grayscale is reduced, making the measurement accuracy higher;

(3)适应环境能力更强:通过扣除背景,降低了环境温度对系统灰度值的影响,使其可以实现一次辐射定标,也可在更广泛的环境温度下使用,甚至只需在探测器本身发生老化后才需重新定标。(3) Stronger adaptability to the environment: By subtracting the background, the influence of ambient temperature on the grayscale value of the system is reduced, so that it can be calibrated once for radiation and can be used in a wider range of ambient temperatures. It even only needs to be recalibrated after the detector itself ages.

附图说明BRIEF DESCRIPTION OF THE DRAWINGS

图1示出根据本发明实施方式的基于背景板扣除内部杂散辐射的红外辐射测量方法的流程图;FIG1 shows a flow chart of an infrared radiation measurement method based on background plate subtraction of internal stray radiation according to an embodiment of the present invention;

图2(a)示出根据本发明实施方式的透10%衰减片、积分时间250us采集的140度黑体图像;FIG2 (a) shows a 140-degree black body image acquired through a 10% attenuation plate with an integration time of 250 us according to an embodiment of the present invention;

图2(b)示出根据本发明实施方式的透10%衰减片、积分时间250us采集的滤光轮背景图像;FIG2( b ) shows a filter wheel background image acquired through a 10% attenuation plate with an integration time of 250 μs according to an embodiment of the present invention;

图2(c)示出根据本发明实施方式的透10%衰减片、积分时间250us采集的镜头盖背景图像;FIG2( c ) shows a lens cover background image acquired through a 10% attenuation plate with an integration time of 250 μs according to an embodiment of the present invention;

图3(a)示出根据本发明实施方式的扣除了滤光轮背景的黑体图像;FIG3( a ) shows a black body image with the filter wheel background subtracted according to an embodiment of the present invention;

图3(b)示出扣除了镜头盖背景的黑体图像。FIG3( b ) shows a black body image with the lens cover background subtracted.

具体实施方式Detailed ways

下面将参照附图更详细地描述本公开的示例性实施方式。虽然附图中显示了本公开的示例性实施方式,然而应当理解,可以以各种形式实现本公开,且本公开不应被这里阐述的实施方式所限制。相反,提供这些实施方式是为了能够更透彻地理解本公开,并且能够将本公开的范围完整地传达给本领域的技术人员。The exemplary embodiments of the present disclosure will be described in more detail below with reference to the accompanying drawings. Although the exemplary embodiments of the present disclosure are shown in the accompanying drawings, it should be understood that the present disclosure can be implemented in various forms, and the present disclosure should not be limited by the embodiments described herein. On the contrary, these embodiments are provided to enable a more thorough understanding of the present disclosure and to fully convey the scope of the present disclosure to those skilled in the art.

本发明技术方案首先涉及到背景板的设计,背景板的设计有滤光轮背景板和镜头盖背景板两种设计方法;再者是背景板和黑体图像的采集策略:在每个黑体温度点,各衰减片和积分时间的图像采集完成后,立刻切至背景板,采集所需背景板图像。然后根据辐射测量理论建立扣除背景的红外辐射定标模型;采用辐射定标算法:根据扣除背景的红外辐射定标模型,由最小二乘法建立目标函数,由最优化算法求得定标增益k和定标偏置b;最后是背景扣除红外辐射测量策略,在测量目标前,将所需档位的背景板图像采集进行存储,在目标测量时调用,根据定标系数、目标测量系统之间的大气透过率和程辐射,反演出目标的表观辐亮度。The technical solution of the present invention firstly involves the design of the background plate. The background plate design includes two design methods: the filter wheel background plate and the lens cover background plate. The second is the acquisition strategy of the background plate and the black body image: after the image acquisition of each attenuation plate and the integration time is completed at each black body temperature point, the background plate is immediately switched to acquire the required background plate image. Then, an infrared radiation calibration model with background subtraction is established according to the radiation measurement theory; a radiation calibration algorithm is adopted: according to the infrared radiation calibration model with background subtraction, the objective function is established by the least square method, and the calibration gain k and the calibration bias b are obtained by the optimization algorithm; finally, the background subtraction infrared radiation measurement strategy is adopted. Before measuring the target, the background plate image of the required gear is collected and stored, and it is called when the target is measured. According to the calibration coefficient, the atmospheric transmittance and the path radiation between the target measurement system, the apparent radiance of the target is inverted.

下面参考图1描述根据本发明实施方式的基于背景板扣除内部杂散辐射的红外辐射测量方法的具体步骤:The specific steps of the infrared radiation measurement method based on background plate subtraction of internal stray radiation according to an embodiment of the present invention are described below with reference to FIG. 1 :

步骤1,进行背景板设计,该设计可以包括:在滤光轮上放置背景板,该背景板尺寸、材料与衰减片相同,表面镀有与衰减片相同的膜系,该背景板被称为滤光轮背景板;或者,在镜头盖上涂覆均匀的低发射率涂层,该涂层被称为镜头盖背景板。Step 1, designing a background plate, which may include: placing a background plate on the filter wheel, the background plate having the same size and material as the attenuation plate, and the surface being coated with the same film system as the attenuation plate, and the background plate is called the filter wheel background plate; or, applying a uniform low-emissivity coating on the lens cover, and the coating is called the lens cover background plate.

步骤2,设计背景板和黑体图像采集策略:待黑体温度稳定后,变换衰减片,依次变换积分时间档位,采集黑体图像;各衰减片采集完后,切至背景板,采集背景板图像,采集背景板图像包括:若是滤光轮背景板则采集不同积分时间的背景板图像,若是镜头盖背景板则采集不同积分时间和衰减片档位的背景板图像;直到设置的所有黑体温度点采集完成。Step 2, design the background plate and black body image acquisition strategy: after the black body temperature stabilizes, change the attenuation plate, change the integration time gear in turn, and collect the black body image; after each attenuation plate is collected, switch to the background plate and collect the background plate image. Collecting the background plate image includes: if it is a filter wheel background plate, collect background plate images with different integration times; if it is a lens cover background plate, collect background plate images with different integration times and attenuation plate gears; until all the set black body temperature points are collected.

判断黑体温度是否稳定的方法为根据黑体控制器反馈温度的稳定性,变换衰减片指的是变换不同透过率的光学衰减片,衰减片档位包括透100%、10%和5%三档衰减片,针对各衰减片有5个积分时间档位,设置的黑体温度点可以包括黑体温度范围的任意温度。The method for judging whether the blackbody temperature is stable is to judge whether the blackbody temperature is stable based on the stability of the feedback temperature of the blackbody controller. The conversion attenuation plate refers to the optical attenuation plate with different transmittances. The attenuation plate gears include 100%, 10% and 5% attenuation plates. There are 5 integration time gears for each attenuation plate. The set blackbody temperature point can include any temperature in the blackbody temperature range.

图2(a),图2(b),图2(c)示出根据本发明实施方式的一采集示例。图2所示的示例的采集背景为:在2022年9月1日,将大面源黑体放置在测量系统镜筒前,黑体温度设置为从50度开始间隔15度升至290度,采集各衰减片和积分时间的图像后,再采集相应背景板图像。图2(a)为透10%衰减片、积分时间250us采集的140度黑体图像,图2(b)为滤光轮背景板图像,图2(c)为镜头盖背景图像。Figure 2 (a), Figure 2 (b), and Figure 2 (c) show an acquisition example according to an embodiment of the present invention. The acquisition background of the example shown in Figure 2 is: on September 1, 2022, a large surface source black body was placed in front of the measurement system barrel, and the black body temperature was set to increase from 50 degrees to 290 degrees at intervals of 15 degrees. After collecting images of each attenuation plate and integration time, the corresponding background plate image was collected. Figure 2 (a) is a 140-degree black body image collected through a 10% attenuation plate and an integration time of 250us, Figure 2 (b) is a filter wheel background plate image, and Figure 2 (c) is a lens cover background image.

步骤3,建立背景扣除的红外辐射定标模型。Step 3: Establish a background-subtracted infrared radiation calibration model.

红外系统引起的辐射,分为衰减片之前的系统引起的杂散辐射和衰减片引起的背景辐射,根据辐射传输理论可得到采集黑体的数学模型和采集背景板的数学模型,从而得到背景扣除的红外辐射定标模型。The radiation caused by the infrared system is divided into stray radiation caused by the system before the attenuation plate and background radiation caused by the attenuation plate. According to the radiation transmission theory, the mathematical model of collecting black bodies and the mathematical model of collecting background plates can be obtained, thereby obtaining the infrared radiation calibration model with background subtraction.

采集黑体的数学表达式为The mathematical expression for collecting black bodies is:

(1) (1)

其中为采集黑体图像系统的灰度总量,为系统响应,为黑体辐射量,为衰减片之前的系统引起的杂散辐射引起的灰度,为背景板辐射引起的灰度,为探测器暗电流引起的灰度。in To collect the total grayscale of the blackbody image system, For system response, is the blackbody radiation, is the grayscale caused by stray radiation caused by the system before the attenuator, is the grayscale caused by the radiation of the background plate, is the grayscale caused by the dark current of the detector.

采集滤光轮处背景板的数学表达式为:The mathematical expression of the background plate at the acquisition filter wheel is:

(2) (2)

其中,为红外辐射测量设备的灰度总量,为衰减片及其之后系统辐射引起的灰度,为探测器暗电流引起的灰度。in, is the total amount of grayscale of the infrared radiation measurement device, is the grayscale caused by the attenuator and the subsequent system radiation, is the grayscale caused by the dark current of the detector.

采集镜头盖背景板的数学表达式为:The mathematical expression of the acquisition lens cover background plate is:

(3) (3)

其中,为红外辐射测量设备的灰度总量,为镜头盖辐射产生的灰度,为衰减片之前的系统引起的杂散辐射引起的灰度,为衰减片及其之后系统辐射引起的灰度,为探测器暗电流引起的灰度。in, is the total amount of grayscale of the infrared radiation measuring device, is the grayscale generated by the lens cover radiation, is the grayscale caused by stray radiation caused by the system before the attenuator, is the grayscale caused by the attenuation film and the subsequent system radiation, is the grayscale caused by the dark current of the detector.

减滤光轮背景板的辐射定标数学模型为:The mathematical model of radiation calibration of the filter wheel background plate is:

(4) (4)

减镜头盖背景板的辐射定标数学模型为:The mathematical model of radiation calibration of the background plate of the lens cover is:

(5) (5)

对式(4)进行简化可得减滤光轮背景板辐射定标数学模型:Simplifying equation (4) yields the mathematical model for radiation calibration of the filter wheel background plate:

(6) (6)

对式(5)进行简化可得减镜头盖背景板辐射定标数学模型:Simplifying equation (5) yields the mathematical model for radiation calibration with lens cover and background plate subtracted:

(7) (7)

步骤4,设计背景扣除的红外辐射定标算法:用采集的黑体图像和背景板图像处理后的数据对步骤3的定标模型进行最优化参数估计,以获得定标增益和定标偏置,定标增益是背景扣除的红外辐射定标模型的斜率,定标偏置是背景扣除的红外辐射定标模型的截距,计算定标误差并判断定标误差是否大于误差阈值,如果大于误差阈值,则红外辐射定标算法继续迭代,直到所有温度点的定标误差都小于误差阈值,该步骤包括首先从黑体图像减去背景图像,然后再去除不在线性区的图像,最后针对以上处理后的图像及其对应的黑体能量进行最优化参数估计,具体包括:Step 4, designing a background-subtracted infrared radiation calibration algorithm: using the collected black body image and the processed data of the background plate image to optimize the parameters of the calibration model in step 3 to obtain the calibration gain and calibration bias, the calibration gain is the slope of the background-subtracted infrared radiation calibration model, the calibration bias is the intercept of the background-subtracted infrared radiation calibration model, calculate the calibration error and determine whether the calibration error is greater than the error threshold, if it is greater than the error threshold, the infrared radiation calibration algorithm continues to iterate until the calibration errors of all temperature points are less than the error threshold, this step includes first subtracting the background image from the black body image, and then removing the image that is not in the linear region, and finally optimizing the parameters for the above processed image and its corresponding black body energy, specifically including:

步骤4.1,所有的黑体图像减相对应的背景图像,得到扣除了背景的图像;图3(a)示出根据本发明实施方式的扣除了滤光轮背景的黑体图像;图3(b)示出扣除了镜头盖背景的黑体图像;Step 4.1, all black body images are subtracted from the corresponding background images to obtain an image with the background subtracted; FIG3 (a) shows a black body image with the filter wheel background subtracted according to an embodiment of the present invention; FIG3 (b) shows a black body image with the lens cover background subtracted;

步骤4.2,根据探测器的线性区,去除不在线性区的黑体图像;其中,探测器是制冷型红外探测器,位于红外辐射测量设备内部;Step 4.2, according to the linear region of the detector, removing the black body image that is not in the linear region; wherein the detector is a cooling type infrared detector located inside the infrared radiation measuring device;

步骤4.3,根据式(8)计算以上(b)处理后的图像对应的不同温度黑体辐射到探测器像元的能量N,根据式(9)依据背景扣除的红外辐射定标模型由最小二乘法针对计算的能量和以上(b)处理后的图像的图像灰度建立目标函数,由最优化算法求得定标增益k和定标偏置b,定标增益k是式(6)、式(7)的斜率,定标偏置b是式(6)、式(7)的截距;最优化算法包括最速梯度法、拟牛顿法和高斯牛顿法等;Step 4.3, according to formula (8), calculate the energy N of the black body radiation of different temperatures corresponding to the image processed by (b) above to the detector pixel, according to formula (9), according to the infrared radiation calibration model after background subtraction, establish the objective function by the least square method for the calculated energy and the image grayscale of the image processed by (b) above, and obtain the calibration gain k and calibration bias b by the optimization algorithm, the calibration gain k is the slope of formula (6) and formula (7), and the calibration bias b is the intercept of formula (6) and formula (7); the optimization algorithm includes the steepest gradient method, the quasi-Newton method and the Gauss-Newton method, etc.;

(8) (8)

其中,M(T)为普朗克函数,F为测量系统的F数,为探测器像元面积,是波长,N为定标黑体的标准辐射量。在一个示例中,测量系统F=1.93、像元大小为15um,响应波长为7.7-9.3um,衰减片第一档为全透,衰减片第二档为透10%,衰减片第三档为透5%。Where M(T) is the Planck function, F is the F number of the measurement system, is the detector pixel area, is the wavelength, and N is the standard radiation of the calibration black body. In an example, the measurement system F=1.93, the pixel size is 15um, the response wavelength is 7.7-9.3um, the first level of the attenuation film is fully transparent, the second level of the attenuation film is 10% transparent, and the third level of the attenuation film is 5% transparent.

(9) (9)

其中,J是目标函数,DN为黑体图像灰度减背景板图像灰度,k为定标增益,b为定标偏置。Among them, J is the objective function, DN is the grayscale of the black body image minus the grayscale of the background plate image, k is the calibration gain, and b is the calibration bias.

步骤4.4,计算各温度点的定标误差:Step 4.4, calculate the calibration error at each temperature point:

(10) (10)

其中,δ为需计算的误差,DN为黑体图像灰度减背景板图像灰度,k为定标增益,b为定标偏置,N为定标黑体的标准辐射量。Among them, δ is the error to be calculated, DN is the grayscale of the blackbody image minus the grayscale of the background plate image, k is the calibration gain, b is the calibration bias, and N is the standard radiation of the calibration blackbody.

步骤4.5,判断定标误差是否大于误差阈值,如果大于误差阈值,则去掉该温度点图像,重新执行步骤4.3,直到所有温度点的定标误差都小于误差阈值,得到最终的定标增益k和定标偏置b;阈值H的设置可以为8%。Step 4.5, determine whether the calibration error is greater than the error threshold. If it is greater than the error threshold, remove the temperature point image and re-execute step 4.3 until the calibration errors of all temperature points are less than the error threshold, and obtain the final calibration gain k and calibration bias b; the threshold H can be set to 8%.

下面的表1至表3示出了250us各衰减片的定标系数。Tables 1 to 3 below show the calibration coefficients of each attenuator at 250us.

表1 250us各衰减片的常规方法定标系数Table 1 Conventional calibration coefficients of various attenuators at 250us

,

表2 减滤光轮背景板,250us各衰减片的定标系数Table 2: Calibration coefficients of filter wheel background plate, 250us attenuation plate

,

表3 减镜头盖背景板,250us各衰减片的定标系数Table 3. Lens cover background plate, calibration coefficients of each attenuation sheet at 250us

,

步骤5,执行背景扣除的红外辐射测量策略:使用步骤4获得的定标增益,获得目标表观辐亮度,该过程可以如下式(11)所示:Step 5, perform the infrared radiation measurement strategy of background subtraction: use the calibration gain obtained in step 4 to obtain the target apparent radiance. The process can be shown in the following equation (11):

(11) (11)

其中,为目标表观辐亮度,为目标灰度值,为对应背景板的灰度值,为目标和测量系统之前的大气程辐射,为目标和测量系统之前的大气透过率,F为测量系统的F数,为探测器像元面积。可以在测量目标前,存储所需档位的背景板图像采集,在目标测量时调用。in, is the target apparent radiance, is the target grayscale value, is the gray value of the corresponding background plate, is the atmospheric path radiation between the target and the measurement system, is the atmospheric transmittance between the target and the measurement system, F is the F number of the measurement system, is the pixel area of the detector. Before measuring the target, the background plate image acquisition of the required gear can be stored and called when measuring the target.

下面描述该步骤的一个示例。在2022年10月13日,将大面源黑体重新放置在测量系统镜筒前,设置几个黑体温度点,采集测量档位的黑体图像,并采集相应的背景板图像;此时大气透过率可认为1,大气程辐射为0,可计算出黑体140度、积分时间250us、各衰减片下的辐亮度,其与理论辐亮度和常规定标方法反演辐亮度的比较如表4所示。An example of this step is described below. On October 13, 2022, the large surface source blackbody was re-placed in front of the measurement system lens barrel, several blackbody temperature points were set, the blackbody image of the measurement gear was collected, and the corresponding background plate image was collected; at this time, the atmospheric transmittance can be considered to be 1, and the atmospheric path radiation is 0. The radiance of the blackbody at 140 degrees, integration time 250us, and each attenuation sheet can be calculated, and its comparison with the theoretical radiance and the radiance inverted by the conventional calibration method is shown in Table 4.

表4 250us本方法反演辐亮度与理论辐亮度和常规定标方法反演辐亮度的对比Table 4 Comparison of the radiance inverted by this method at 250us with the theoretical radiance and the radiance inverted by the conventional calibration method

,

在此处所提供的说明书中,说明了大量具体细节。然而,能够理解,本发明的实施方式可以在没有这些具体细节的情况下被实践。在一些实例中,并未详细示出公知的方法、结构和技术,以便不模糊对本说明书的理解。In the description provided herein, a large number of specific details are described. However, it is understood that embodiments of the present invention can be practiced without these specific details. In some instances, well-known methods, structures and techniques are not shown in detail so as not to obscure the understanding of this description.

尽管根据有限数量的实施方式描述了本发明,但是受益于上面的描述,本技术领域内的技术人员明白,在由此描述的本发明的范围内,可以设想其它实施方式。此外,应当注意,本说明书中使用的语言主要是为了可读性和教导的目的而选择的,而不是为了解释或者限定本发明的主题而选择的。Although the present invention has been described according to a limited number of embodiments, it will be apparent to those skilled in the art, with the benefit of the above description, that other embodiments may be envisioned within the scope of the invention thus described. In addition, it should be noted that the language used in this specification is primarily selected for readability and instructional purposes, rather than for explaining or limiting the subject matter of the present invention.

Claims (5)

Translated fromChinese
1.一种基于背景板扣除内部杂散辐射的红外辐射测量方法,其特征在于,该方法包括如下步骤:1. A method for measuring infrared radiation based on background plate subtraction of internal stray radiation, characterized in that the method comprises the following steps:步骤1,进行背景板设计;Step 1, design the background board;步骤2,设计背景板和黑体图像采集策略:待黑体温度稳定后,变换衰减片,依次变换积分时间档位,采集黑体图像;各衰减片采集完后,切至背景板,采集背景板图像;Step 2, design the background plate and black body image acquisition strategy: after the black body temperature is stable, change the attenuation plate, change the integration time gear in sequence, and collect the black body image; after each attenuation plate is collected, switch to the background plate and collect the background plate image;步骤3,建立背景扣除的红外辐射定标模型;Step 3, establishing a background-subtracted infrared radiation calibration model;步骤4,设计背景扣除的红外辐射定标算法:用采集的黑体图像和背景板图像处理后的数据对步骤3的红外辐射定标模型进行最优化参数估计,以获得定标增益和定标偏置,定标增益是背景扣除的红外辐射定标模型的斜率,定标偏置是背景扣除的红外辐射定标模型的截距,计算定标误差并判断定标误差是否大于误差阈值,如果大于误差阈值,则红外辐射定标算法继续迭代,直到所有温度点的定标误差都小于误差阈值;Step 4, designing an infrared radiation calibration algorithm with background subtraction: using the collected black body image and the processed data of the background plate image to optimize the parameter estimation of the infrared radiation calibration model in step 3 to obtain the calibration gain and calibration bias, the calibration gain is the slope of the infrared radiation calibration model with background subtraction, the calibration bias is the intercept of the infrared radiation calibration model with background subtraction, calculating the calibration error and judging whether the calibration error is greater than the error threshold, if it is greater than the error threshold, the infrared radiation calibration algorithm continues to iterate until the calibration errors of all temperature points are less than the error threshold;步骤5,执行背景扣除的红外辐射测量策略:使用步骤4获得的定标增益,获得目标表观辐亮度。Step 5, execute the infrared radiation measurement strategy with background subtraction: use the calibration gain obtained in step 4 to obtain the target apparent radiance.2.根据权利要求1所述的基于背景板扣除内部杂散辐射的红外辐射测量方法,其特征在于,步骤1包括:在滤光轮上放置背景板,该背景板的尺寸、材料与衰减片相同,表面镀有与衰减片相同的膜系,该背景板被称为滤光轮背景板;或者,在镜头盖上涂覆均匀的低发射率涂层,该涂层被称为镜头盖背景板;2. The infrared radiation measurement method based on background plate subtraction of internal stray radiation according to claim 1 is characterized in that step 1 comprises: placing a background plate on the filter wheel, the size and material of the background plate being the same as those of the attenuation plate, and the surface of the background plate being coated with the same film system as that of the attenuation plate, and the background plate is referred to as the filter wheel background plate; or, coating a uniform low-emissivity coating on the lens cover, and the coating is referred to as the lens cover background plate;步骤2包括:针对滤光轮背景板,采集不同积分时间的滤光轮背景板图像,针对镜头盖背景板,采集不同积分时间和衰减片档位的镜头盖背景板图像;直到设置的所有黑体温度点采集完成。Step 2 includes: for the filter wheel background plate, collecting filter wheel background plate images with different integration times, and for the lens cover background plate, collecting lens cover background plate images with different integration times and attenuation plate positions; until all the set black body temperature points are collected.3.根据权利要求2所述的基于背景板扣除内部杂散辐射的红外辐射测量方法,其特征在于,步骤3包括建立减滤光轮背景板的辐射定标数学模型和减镜头盖背景板的辐射定标数学模型:3. The infrared radiation measurement method based on background plate subtraction of internal stray radiation according to claim 2 is characterized in that step 3 comprises establishing a radiation calibration mathematical model of the filter wheel background plate and a radiation calibration mathematical model of the lens cover background plate:减滤光轮背景板的辐射定标数学模型为:The mathematical model of radiation calibration of the filter wheel background plate is:(6) (6)减镜头盖背景板的辐射定标数学模型为:The mathematical model of radiation calibration of the background plate of the lens cover is: (7) (7)其中,为采集黑体图像系统的灰度总量,为采集滤光轮背景板系统的灰度总量,为系统响应,为黑体辐射量,为衰减片之前的系统引起的杂散辐射引起的灰度,为镜头盖辐射产生的灰度,为采集镜头盖背景板系统的灰度总量。in, To collect the total grayscale of the blackbody image system, To collect the total grayscale of the filter wheel background system, For system response, is the blackbody radiation, is the grayscale caused by stray radiation caused by the system before the attenuator, is the grayscale generated by the lens cover radiation, To collect the total grayscale of the lens cover background system.4.根据权利要求3所述的基于背景板扣除内部杂散辐射的红外辐射测量方法,其特征在于,步骤4包括:4. The infrared radiation measurement method based on background plate subtraction of internal stray radiation according to claim 3, characterized in that step 4 comprises:步骤4.1,将所有的黑体图像减相对应的背景图像,得到扣除了背景的图像;Step 4.1, subtract the corresponding background image from all black body images to obtain an image with the background subtracted;步骤4.2,根据探测器的线性区,去除不在线性区的黑体图像;其中,探测器是制冷型红外探测器,位于红外辐射测量设备内部;Step 4.2, according to the linear region of the detector, removing the black body image that is not in the linear region; wherein the detector is a cooling type infrared detector located inside the infrared radiation measuring device;步骤4.3,根据式(8)计算不同温度黑体辐射到像元的能量N,根据式(9)依据背景扣除的红外辐射定标模型由最小二乘法针对计算的能量和测量灰度建立目标函数,由最优化算法求得定标增益k和定标偏置b,定标增益k是式(6)、式(7)的斜率,偏置b是式(6)、式(7)的截距;Step 4.3, calculate the energy N radiated to the pixel by the black body at different temperatures according to formula (8), establish the objective function for the calculated energy and the measured grayscale by the least square method according to the infrared radiation calibration model after background subtraction according to formula (9), and obtain the calibration gain k and calibration bias b by the optimization algorithm. The calibration gain k is the slope of formula (6) and formula (7), and the bias b is the intercept of formula (6) and formula (7);(8) (8)其中,M(T)为普朗克函数,F为测量系统的F数,为探测器像元面积,是波长,N为不同温度黑体辐射到像元的能量;Where M(T) is the Planck function, F is the F number of the measurement system, is the detector pixel area, is the wavelength, N is the energy radiated to the pixel by a black body of different temperatures;(9) (9)其中,J是目标函数,DN为黑体图像灰度减背景板图像灰度,k为定标增益,b为定标偏置;Where J is the objective function, DN is the grayscale of the black body image minus the grayscale of the background plate image, k is the calibration gain, and b is the calibration bias;步骤4.4,计算各温度点的定标误差:Step 4.4, calculate the calibration error at each temperature point:(10) (10)其中,δ为需计算的误差;Among them, δ is the error to be calculated;步骤4.5,判断定标误差是否大于误差阈值,如果大于误差阈值,则重新执行步骤4.3,直到所有温度点的定标误差都小于误差阈值。Step 4.5, determine whether the calibration error is greater than the error threshold. If it is greater than the error threshold, re-execute step 4.3 until the calibration errors of all temperature points are less than the error threshold.5.根据权利要求4所述的基于背景板扣除内部杂散辐射的红外辐射测量方法,其特征在于,步骤5通过下式执行背景扣除的红外辐射测量策略:5. The infrared radiation measurement method based on background plate subtraction of internal stray radiation according to claim 4 is characterized in that step 5 implements the infrared radiation measurement strategy of background subtraction by the following formula:(11) (11)其中,为目标表观辐亮度,为目标灰度值,为对应背景板的灰度值,为目标和测量系统之前的大气程辐射,为目标和测量系统之前的大气透过率,F为测量系统的F数,为探测器像元面积。in, is the target apparent radiance, is the target grayscale value, is the gray value of the corresponding background plate, is the atmospheric path radiation between the target and the measurement system, is the atmospheric transmittance between the target and the measurement system, F is the F number of the measurement system, is the detector pixel area.
CN202410865082.XA2024-07-012024-07-01Infrared radiation measurement method for subtracting internal stray radiation based on background plateActiveCN118392324B (en)

Priority Applications (1)

Application NumberPriority DateFiling DateTitle
CN202410865082.XACN118392324B (en)2024-07-012024-07-01Infrared radiation measurement method for subtracting internal stray radiation based on background plate

Applications Claiming Priority (1)

Application NumberPriority DateFiling DateTitle
CN202410865082.XACN118392324B (en)2024-07-012024-07-01Infrared radiation measurement method for subtracting internal stray radiation based on background plate

Publications (2)

Publication NumberPublication Date
CN118392324Atrue CN118392324A (en)2024-07-26
CN118392324B CN118392324B (en)2024-10-15

Family

ID=91999951

Family Applications (1)

Application NumberTitlePriority DateFiling Date
CN202410865082.XAActiveCN118392324B (en)2024-07-012024-07-01Infrared radiation measurement method for subtracting internal stray radiation based on background plate

Country Status (1)

CountryLink
CN (1)CN118392324B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication numberPriority datePublication dateAssigneeTitle
CN119374736A (en)*2024-12-262025-01-28中国科学院光电技术研究所 An efficient calibration method for wide dynamic range infrared radiation measurement system

Citations (13)

* Cited by examiner, † Cited by third party
Publication numberPriority datePublication dateAssigneeTitle
JPS577525A (en)*1980-06-161982-01-14Jeol LtdOptical system of infrared ray camera
JP2007251888A (en)*2006-03-202007-09-27Fujitsu Ltd Infrared imaging device
CN101470026A (en)*2007-12-242009-07-01南京理工大学Ununiformity emendation real-time calibration apparatus for staring type thermal imaging system
CN101776486A (en)*2009-12-312010-07-14华中科技大学Method for correcting non-uniformity fingerprint pattern on basis of infrared focal plane
US20130258313A1 (en)*2012-03-312013-10-03Daniel OrbandImage analysis system and methods for ir optics
US20150326799A1 (en)*2014-05-072015-11-12Microsoft CorporationReducing camera interference using image analysis
CN108230286A (en)*2016-12-162018-06-29北京机电工程研究所The influence rapid analysis method of environment and period factor to Infrared Targets detectivity
CN109282901A (en)*2018-11-162019-01-29北京遥感设备研究所 A split-block correction method for an uncooled infrared system
CN109737987A (en)*2018-12-292019-05-10中国科学院长春光学精密机械与物理研究所 A multi-photosynthetic large-aperture space camera infrared radiation calibration system on orbiting stars
US20200202569A1 (en)*2018-12-202020-06-25Flir Systems AbInfrared camera ambient temperature calibration systems and methods
CN115560853A (en)*2022-10-142023-01-03中国科学院长春光学精密机械与物理研究所Long-wave infrared spectrum band background radiation elimination method and system based on long-wave infrared detector
CN115615560A (en)*2022-10-262023-01-17中国科学院光电技术研究所 A Wide Dynamic Range High Precision Infrared Radiation Measurement Method
CN116989901A (en)*2023-07-242023-11-03中国人民解放军95859部队Air target infrared radiation brightness calculation method based on pixel-by-pixel calibration

Patent Citations (13)

* Cited by examiner, † Cited by third party
Publication numberPriority datePublication dateAssigneeTitle
JPS577525A (en)*1980-06-161982-01-14Jeol LtdOptical system of infrared ray camera
JP2007251888A (en)*2006-03-202007-09-27Fujitsu Ltd Infrared imaging device
CN101470026A (en)*2007-12-242009-07-01南京理工大学Ununiformity emendation real-time calibration apparatus for staring type thermal imaging system
CN101776486A (en)*2009-12-312010-07-14华中科技大学Method for correcting non-uniformity fingerprint pattern on basis of infrared focal plane
US20130258313A1 (en)*2012-03-312013-10-03Daniel OrbandImage analysis system and methods for ir optics
US20150326799A1 (en)*2014-05-072015-11-12Microsoft CorporationReducing camera interference using image analysis
CN108230286A (en)*2016-12-162018-06-29北京机电工程研究所The influence rapid analysis method of environment and period factor to Infrared Targets detectivity
CN109282901A (en)*2018-11-162019-01-29北京遥感设备研究所 A split-block correction method for an uncooled infrared system
US20200202569A1 (en)*2018-12-202020-06-25Flir Systems AbInfrared camera ambient temperature calibration systems and methods
CN109737987A (en)*2018-12-292019-05-10中国科学院长春光学精密机械与物理研究所 A multi-photosynthetic large-aperture space camera infrared radiation calibration system on orbiting stars
CN115560853A (en)*2022-10-142023-01-03中国科学院长春光学精密机械与物理研究所Long-wave infrared spectrum band background radiation elimination method and system based on long-wave infrared detector
CN115615560A (en)*2022-10-262023-01-17中国科学院光电技术研究所 A Wide Dynamic Range High Precision Infrared Radiation Measurement Method
CN116989901A (en)*2023-07-242023-11-03中国人民解放军95859部队Air target infrared radiation brightness calculation method based on pixel-by-pixel calibration

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
LI, T等: "Global optimization radiometric calibration method for an infrared system with a wide dynamic range", 《APPLIED OPTICS》, vol. 63, no. 2, 10 January 2024 (2024-01-10), pages 406 - 414*
MU GU等: "A new infrared radiation calibration method based on lens cover deduction of internal stray radiation", 《INFRARED PHYSICS AND TECHNOLOGY》, vol. 138, 24 March 2024 (2024-03-24), pages 1 - 5*
刘品伟等: "国产氧化钒红外探测器新型无挡板校正方法", 《光学与光电技术》, vol. 20, no. 01, 31 December 2022 (2022-12-31), pages 128 - 133*
张涛;: "内外定标结合的辐射特性测量方法", 智能计算机与应用, no. 06, 1 December 2014 (2014-12-01)*

Cited By (1)

* Cited by examiner, † Cited by third party
Publication numberPriority datePublication dateAssigneeTitle
CN119374736A (en)*2024-12-262025-01-28中国科学院光电技术研究所 An efficient calibration method for wide dynamic range infrared radiation measurement system

Also Published As

Publication numberPublication date
CN118392324B (en)2024-10-15

Similar Documents

PublicationPublication DateTitle
CN107843982B (en)Wave front-free detection self-adaptive optical system based on real-time phase difference technology
CN111751003B (en)Thermal imager temperature correction system and method and thermal imager
CN107677375A (en) Calibration device and calibration method for infrared radiation measurement system
CN118392324A (en)Infrared radiation measurement method for subtracting internal stray radiation based on background plate
CN105160631B (en)A kind of method for seeking radiant correction coefficient
CN106873152B (en) A high-speed aberration correction method based on machine learning
CN103267533B (en)A kind of practical high-spectrum remote sensing air automatic correcting method
CN107092752A (en)A kind of optical camera simulation imaging method and system based on ray tracing
CN105654430B (en)Contrast constrained pneumatic thermal radiation correction method
CN106679817A (en)Method for calibrating thermal infrared imager
CN110646099A (en)Method and device for inverting target infrared radiation image based on measured data
CN113865717A (en)Transient high-temperature colorimetric temperature measuring device based on high-speed camera
CN104296882A (en)Large-caliber and wide-dynamic-range infrared system radiometric calibration method
CN118624033A (en) A high-precision infrared radiation measurement method with wide dynamic range adapting to ambient temperature changes
CN108226059B (en)Satellite hyperspectral CO2On-orbit radiation calibration method for detector
CN117097997A (en)Noise image synthesis method for reverse image signal processing
CN114152352B (en)Method and system for experimental measurement of stray radiation of infrared optical system
CN115712202A (en)Zernike polynomial-based method for compensating image quality degradation caused by temperature
CN104360478A (en)Nested double-adaptive optical system
CN118822914A (en) An efficient non-uniformity correction method considering the influence of attenuation sheet and integration time
CN105866939A (en)Method for introducing wavefront distortion on basis of micro deformation mirror
CN113484255A (en)Modified visibility estimation method based on near-infrared double-spectrum real-time imaging and distance
CN113970374A (en) A calibration method for a defocusing plane polarization detection system
CN118794547A (en) A method for measuring target infrared radiation brightness considering ambient temperature correction
CN115931142A (en)Colorimetric temperature measuring black body calibration method

Legal Events

DateCodeTitleDescription
PB01Publication
PB01Publication
SE01Entry into force of request for substantive examination
SE01Entry into force of request for substantive examination
GR01Patent grant
GR01Patent grant

[8]ページ先頭

©2009-2025 Movatter.jp