CN106352987A - Calibration method, device and system for infrared thermal imager - Google Patents

Calibration method, device and system for infrared thermal imager Download PDF

Info

Publication number
CN106352987A
CN106352987A CN201610742769.XA CN201610742769A CN106352987A CN 106352987 A CN106352987 A CN 106352987A CN 201610742769 A CN201610742769 A CN 201610742769A CN 106352987 A CN106352987 A CN 106352987A
Authority
CN
China
Prior art keywords
equation
temperature
temperature drift
coefficient
quadratic function
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
CN201610742769.XA
Other languages
Chinese (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.)
Hunan University of Arts and Science
Original Assignee
Hunan University of Arts and Science
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 Hunan University of Arts and Science filed Critical Hunan University of Arts and Science
Priority to CN201610742769.XA priority Critical patent/CN106352987A/en
Publication of CN106352987A publication Critical patent/CN106352987A/en
Pending legal-status Critical Current

Links

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
    • 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

Abstract

The invention belongs to the field of infrared detection technology, and provides a calibration method, device and system for infrared thermal imager. The method comprises the following steps: acquiring the temperature of the standard blackbody, the grayscale of the image corresponding to the standard blackbody, and the detector temperature; working out coefficients of the quadratic function calibration equation according to the quadratic function, the temperature of the standard blackbody and the grayscale of the image corresponding to the standard blackbody; working out nine parameters of the quadratic function temperature excursion equation according to the quadratic function temperature excursion equation, coefficients of the quadratic function calibration equation, and the detector temperature; determining the corresponding temperature of the target object according to the nine parameters of the quadratic function temperature excursion equation, the quadratic function temperature excursion equation, the quadratic function calibration equation, and the pre-acquired image grayscale of the target object. The calibration method, device and system for infrared thermal imager can effectively calibrate the infrared thermal imager, with high accuracy and small error, improving the performance of the device.

Description

Thermal infrared imager scaling method, apparatus and system
Technical field
The present invention relates to infrared detection technology field is and in particular to a kind of thermal infrared imager scaling method, device and be System.
Background technology
At present, thermal infrared imager all must be demarcated before dispatching from the factory, also referred to as indoor standardization or Laboratory Calibration, dispatches from the factory Also must periodically be demarcated afterwards, with the performance of checking system.If site environment differs larger with demarcating environment, there is a need to Field calibration is carried out to equipment, to ensure the temperature measurement accuracy of tested target.So-called demarcate it is simply that with standard black body radiation source pair The work that infrared heat image instrument measuring result is corrected obtains the Digital output of thermal infrared imager and thermal infrared imager connects Relation between the spoke brightness receiving.
Easily worked by ambient temperature, distance, movement temperature or detector especially for uncooled IRFPA type Infrared Detectorss Temperature effects, camera lens impact, and existing scaling method complex disposal process, accuracy is poor, has a strong impact on thermal infrared imager Accuracy in temperature sensing and reliability.
How effectively thermal infrared imager to be demarcated, improves equipment precision and performance, be those skilled in the art urgently The problem that need to solve.
Content of the invention
For defect of the prior art, the present invention provides a kind of thermal infrared imager scaling method, apparatus and system, energy Enough accurately and efficiently thermal infrared imager is demarcated, and simplify processing procedure, precision high level error is little, improve the performance of equipment.
In a first aspect, the present invention provides a kind of thermal infrared imager scaling method, the method specifically comprises the following steps that
Obtain the temperature of standard blackbody and gradation of image corresponding with standard blackbody, and detector temperature, standard blackbody Number be multiple;
Temperature according to quadratic function calibration equation and standard blackbody and gradation of image corresponding with this standard blackbody, obtain Obtain each term coefficient of quadratic function calibration equation;
According to each term coefficient of quadratic function temperature drift equation and quadratic function calibration equation, and detector temperature, obtain Nine parameters of quadratic function temperature drift equation;
According to nine parameters of quadratic function temperature drift equation, quadratic function temperature drift equation and quadratic function calibration equation, with And the gradation of image of pre-acquiring target object, determine the corresponding temperature of target object.
Further, the temperature according to quadratic function calibration equation and standard blackbody and figure corresponding with this standard blackbody As gray scale, obtain each term coefficient of quadratic function calibration equation, specifically include:
According to following quadratic function calibration equation, and the temperature of standard blackbody and this standard blackbody corresponding image ash Degree, determines quadratic coefficients, the coefficient of first order of quadratic function calibration equation and the quadratic function demarcation side of quadratic function calibration equation The constant term of journey,
D=a*t2+b*t+c
Wherein, d is the corresponding gradation of image of standard blackbody, and t is the temperature of standard blackbody, and a is quadratic function calibration equation Quadratic coefficients, b be quadratic function calibration equation coefficient of first order, c be quadratic function calibration equation constant term, quadratic function Each term coefficient of calibration equation include the quadratic coefficients of quadratic function calibration equation, the coefficient of first order of quadratic function calibration equation, The constant term of quadratic function calibration equation.
Further, each term coefficient according to quadratic function temperature drift equation and quadratic function calibration equation, and detector Temperature, obtains nine parameters of quadratic function temperature drift equation, specifically includes:
Quadratic coefficients according to following quadratic coefficients temperature drift equation and quadratic function calibration equation and detector temperature, Obtain quadratic coefficients, the coefficient of first order of quadratic coefficients temperature drift equation and the quadratic coefficients temperature drift equation of quadratic coefficients temperature drift equation Constant term,
A=a2*tx 2+a1*tx+a0
Wherein, a is the quadratic coefficients of quadratic function calibration equation, txFor detector temperature, a2For quadratic coefficients temperature drift equation Quadratic coefficients, a1For the coefficient of first order of quadratic coefficients temperature drift equation, a0For the constant term of quadratic coefficients temperature drift equation,
Coefficient of first order according to following coefficient of first order temperature drift equation and quadratic function calibration equation and detector temperature, Obtain quadratic coefficients, the coefficient of first order of coefficient of first order temperature drift equation and the coefficient of first order temperature drift equation of coefficient of first order temperature drift equation Constant term,
B=b2*tx 2+b1*tx+b0
Wherein, b is the coefficient of first order of quadratic function calibration equation, b2For the quadratic coefficients of coefficient of first order temperature drift equation, b1For The coefficient of first order of coefficient of first order temperature drift equation, b0For the constant term of coefficient of first order temperature drift equation,
Constant term according to following constant term temperature drift equation and quadratic function calibration equation and detector temperature, obtain The constant term of the quadratic coefficients of constant term temperature drift equation, the coefficient of first order of constant term temperature drift equation and constant term temperature drift equation,
C=c2*tx 2+c1*tx+c0
Wherein, c is the constant term of quadratic function calibration equation, c2For the quadratic coefficients of constant term temperature drift equation, c1For constant The coefficient of first order of item temperature drift equation, c0For the constant term of constant term temperature drift equation,
Quadratic function temperature drift equation includes quadratic coefficients temperature drift equation, coefficient of first order temperature drift equation and constant term temperature drift side Journey, nine parameters of quadratic function temperature drift equation include the quadratic coefficients of quadratic coefficients temperature drift equation, quadratic coefficients temperature drift equation Coefficient of first order, the constant term of quadratic coefficients temperature drift equation, the quadratic coefficients of coefficient of first order temperature drift equation, coefficient of first order temperature drift side The coefficient of first order of journey, the constant term of coefficient of first order temperature drift equation, the quadratic coefficients of constant term temperature drift equation, constant term temperature drift equation Coefficient of first order and constant term temperature drift equation constant term.
Further, according to nine parameters of quadratic function temperature drift equation, quadratic function temperature drift equation and quadratic function mark Determine equation, and the gradation of image of pre-acquiring target object, determine the corresponding temperature of target object, specifically include:
According to nine parameters of quadratic function temperature drift equation, quadratic function temperature drift equation and quadratic function calibration equation, with And the gradation of image of pre-acquiring target object, by equation below, determine the corresponding temperature of target object,
t ( d 0 ) = - f b ( t x ) + f b ( t x ) 2 - 4 * f a ( t x ) * ( f c ( t x ) - d 0 ) 2 * f a ( t x ) f a ( t x ) = a 2 * t x 2 + a 1 * t x + a 0 f b ( t x ) = b 2 * t x 2 + b 1 * t x + b 0 f c ( t x ) = c 2 * t x 2 + c 1 * t x + c 0
Wherein, t (d0) it is the corresponding temperature of target object, d0Gradation of image for pre-acquiring target object.
Second aspect, the present invention provides a kind of thermal infrared imager caliberating device, and this device includes data acquisition module, mark Determine coefficient determination module, temperature drift coefficient determination module and target temperature acquisition module.Data acquisition module is used for obtaining standard black The temperature of body and gradation of image corresponding with standard blackbody, and detector temperature, the number of standard blackbody is multiple.Demarcate system Number determining module is used for the temperature according to quadratic function calibration equation and standard blackbody and image corresponding with this standard blackbody Gray scale, obtains each term coefficient of quadratic function calibration equation.Temperature drift coefficient determination module is used for according to quadratic function temperature drift equation With each term coefficient of quadratic function calibration equation, and detector temperature, obtain nine parameters of quadratic function temperature drift equation.Mesh Mark temperature acquisition module is used for according to nine parameters of quadratic function temperature drift equation, quadratic function temperature drift equation and quadratic function mark Determine equation, and the gradation of image of pre-acquiring target object, determine the corresponding temperature of target object.
Further, calibration coefficient determining module the temperature according to quadratic function calibration equation and standard blackbody and with The corresponding gradation of image of this standard blackbody, during each term coefficient of acquisition quadratic function calibration equation, specifically for: according to following two Secondary function calibration equation, and the temperature of standard blackbody and the corresponding gradation of image of this standard blackbody, determine that quadratic function is demarcated The constant term of the quadratic coefficients of equation, the coefficient of first order of quadratic function calibration equation and quadratic function calibration equation,
D=a*t2+b*t+c
Wherein, d is the corresponding gradation of image of standard blackbody, and t is the temperature of standard blackbody, and a is quadratic function calibration equation Quadratic coefficients, b be quadratic function calibration equation coefficient of first order, c be quadratic function calibration equation constant term, quadratic function Each term coefficient of calibration equation include the quadratic coefficients of quadratic function calibration equation, the coefficient of first order of quadratic function calibration equation, The constant term of quadratic function calibration equation.
Further, temperature drift coefficient determination module is each according to quadratic function temperature drift equation and quadratic function calibration equation Term coefficient, and detector temperature, during nine parameters of acquisition quadratic function temperature drift equation, specifically for: according to secondary as follows The quadratic coefficients of coefficient temperature drift equation and quadratic function calibration equation and detector temperature, obtain quadratic coefficients temperature drift equation Quadratic coefficients, the coefficient of first order of quadratic coefficients temperature drift equation and quadratic coefficients temperature drift equation constant term,
A=a2*tx 2+a1*tx+a0
Wherein, a is the quadratic coefficients of quadratic function calibration equation, txFor detector temperature, a2For quadratic coefficients temperature drift equation Quadratic coefficients, a1For the coefficient of first order of quadratic coefficients temperature drift equation, a0For the constant term of quadratic coefficients temperature drift equation,
Coefficient of first order according to following coefficient of first order temperature drift equation and quadratic function calibration equation and detector temperature, Obtain quadratic coefficients, the coefficient of first order of coefficient of first order temperature drift equation and the coefficient of first order temperature drift equation of coefficient of first order temperature drift equation Constant term,
B=b2*tx 2+b1*tx+b0
Wherein, b is the coefficient of first order of quadratic function calibration equation, txFor detector temperature, b2For coefficient of first order temperature drift equation Quadratic coefficients, b1For the coefficient of first order of coefficient of first order temperature drift equation, b0For the constant term of coefficient of first order temperature drift equation,
Constant term according to following constant term temperature drift equation and quadratic function calibration equation and detector temperature, obtain The constant term of the quadratic coefficients of constant term temperature drift equation, the coefficient of first order of constant term temperature drift equation and constant term temperature drift equation,
C=c2*tx 2+c1*tx+c0
Wherein, c is the constant term of quadratic function calibration equation, txFor detector temperature, c2For constant term temperature drift equation two Ordered coefficients, c1For the coefficient of first order of constant term temperature drift equation, c0For the constant term of constant term temperature drift equation,
Quadratic function temperature drift equation includes quadratic coefficients temperature drift equation, coefficient of first order temperature drift equation and constant term temperature drift side Journey, nine parameters of quadratic function temperature drift equation include the quadratic coefficients of quadratic coefficients temperature drift equation, quadratic coefficients temperature drift equation Coefficient of first order, the constant term of quadratic coefficients temperature drift equation, the quadratic coefficients of coefficient of first order temperature drift equation, coefficient of first order temperature drift side The coefficient of first order of journey, the constant term of coefficient of first order temperature drift equation, the quadratic coefficients of constant term temperature drift equation, constant term temperature drift equation Coefficient of first order and constant term temperature drift equation constant term.
Further, target temperature acquisition module is in nine parameters according to quadratic function temperature drift equation, quadratic function temperature Drift equation and quadratic function calibration equation, and the gradation of image of pre-acquiring target object, determine the corresponding temperature of target object When, specifically for: according to nine parameters of quadratic function temperature drift equation, quadratic function temperature drift equation and quadratic function demarcation side Journey, and the gradation of image of pre-acquiring target object, by equation below, determine the corresponding temperature of target object,
t ( d 0 ) = - f b ( t x ) + f b ( t x ) 2 - 4 * f a ( t x ) * ( f c ( t x ) - d 0 ) 2 * f a ( t x ) f a ( t x ) = a 2 * t x 2 + a 1 * t x + a 0 f b ( t x ) = b 2 * t x 2 + b 1 * t x + b 0 f c ( t x ) = c 2 * t x 2 + c 1 * t x + c 0
Wherein, t (d0) it is the corresponding temperature of target object, d0Gradation of image for pre-acquiring target object.
The third aspect, the present invention provides a kind of thermal infrared imager calibration system, and this system includes thermal infrared imager, multiple Standard blackbody and seal box, thermal infrared imager is located at seal box, and thermal infrared imager includes flir camera lens, the detector being sequentially connected And processor, filr camera lens is for by the infrared radiation collection of standard blackbody or target object to detector, detector is used for adopting The infra-red radiation of collection standard blackbody, and be converted to the infrared hybrid optical system of standard blackbody, or the infra-red radiation of collection target object, And be converted to the infrared hybrid optical system of target object, and the infrared hybrid optical system of standard blackbody and/or target object is red Outer gray level image sends to processor.Processor is used for obtaining the temperature of standard blackbody and infrared gray-scale maps corresponding with standard blackbody The gradation of image of picture, and detector temperature, the temperature according to quadratic function calibration equation and standard blackbody and with this standard The corresponding gradation of image of black matrix, obtains each term coefficient of secondary calibration equation, always according to quadratic function temperature drift equation and two deutero-albumoses Determine each term coefficient of equation, and detector temperature, obtain nine parameters of quadratic function temperature drift equation, and according to secondary letter Nine parameters of number temperature drift equation, quadratic function temperature drift equation and quadratic function calibration equation, and pre-acquiring target object is red The gradation of image of outer gray level image, determines the corresponding temperature of target object.The side wall of seal box is provided with opening, opening and flir mirror Head coupling, and it is provided with germanium wafer.
As shown from the above technical solution, thermal infrared imager scaling method of the present invention, apparatus and system, is intended using biquadratic Close nine parameter calibration methods, scaling method accuracy is high, error is little.Calibration coefficient determining module is according to the temperature of multiple standard blackbody The corresponding gradation of image with standard blackbody, using quadratic function calibration equation, is fitted, and in extrapolation measure, accuracy is high, Error in data is little.
Meanwhile, each term coefficient of quadratic function calibration equation can fully combine detector temperature, and temperature drift coefficient determines mould Block is fitted using quadratic function temperature drift equation, and determines nine parameters of temperature drift function, and mean error is little, and degree of accuracy is high. It is fitted using biquadratic function, and determines nine parameters of temperature drift equation, target temperature acquisition module can be effectively to red Outer thermal imaging system is demarcated, and improves equipment dependability, in conjunction with the gradation of image of target object, can effectively determine object The temperature of body.
Additionally, being corrected to the measurement result of thermal infrared imager using multiple standard black body radiation sources, accuracy is high.Close Joint sealing can be placed in refrigerator or high-temperature cabinet, conveniently provides a metastable experimental situation, germanium wafer for thermal infrared imager Can effectively prevent flir camera lens from producing solidifying mist.
Therefore, thermal infrared imager scaling method of the present invention, apparatus and system, can be carried out to thermal infrared imager effectively Demarcate, and simplify processing procedure, precision high level error is little, improves equipment performance.
Brief description
In order to be illustrated more clearly that the specific embodiment of the invention or technical scheme of the prior art, below will be to concrete In embodiment or description of the prior art, the accompanying drawing of required use is briefly described.In all of the figs, similar element Or partly typically identified by similar reference.In accompanying drawing, each element or part might not be drawn according to actual ratio.
Fig. 1 shows a kind of thermal infrared imager scaling method flow chart provided by the present invention;
Fig. 2 shows a kind of multiple matched curve figure with standard blackbody temperature as independent variable provided by the present invention;
Fig. 3 shows a kind of multiple matched curve figure with gradation of image as independent variable provided by the present invention;
Fig. 4 shows the relation between a kind of infrared hybrid optical system gradation of image provided by the present invention and target temperature Figure;
Fig. 5 shows that a kind of temperature spot of demarcating provided by the present invention is chosen and image gray levels choice schematic diagram;
Fig. 6 shows a kind of structural representation of thermal infrared imager calibration system provided by the present invention;
Fig. 7 shows the veiling glare effectiveness comparison figure of three kinds of camera lenses provided by the present invention;
Fig. 8 shows the thermometric matched curve figure under varying environment provided by the present invention;
Fig. 9 shows a kind of quadratic coefficients temperature drift equation curve figure provided by the present invention;
Figure 10 shows a kind of coefficient of first order temperature drift equation curve figure provided by the present invention;
Figure 11 shows a kind of constant term temperature drift equation curve figure provided by the present invention.
Specific embodiment
Below in conjunction with accompanying drawing, the embodiment of technical solution of the present invention is described in detail.Following examples are only used for Clearly technical scheme is described, is therefore intended only as example, and the protection of the present invention can not be limited with this Scope.
It should be noted that unless otherwise stated, technical term used in this application or scientific terminology should be this The ordinary meaning that bright one of ordinary skill in the art are understood.
In a first aspect, a kind of thermal infrared imager scaling method that first embodiment of the invention is provided, in conjunction with Fig. 1, should Method specifically comprises the following steps that
Step s1, obtains the temperature of standard blackbody and gradation of image corresponding with standard blackbody, and detector temperature, mark The number of quasi- black matrix is multiple;
Step s2, the temperature according to quadratic function calibration equation and standard blackbody and image corresponding with this standard blackbody Gray scale, obtains each term coefficient of quadratic function calibration equation;
Step s3, according to each term coefficient of quadratic function temperature drift equation and quadratic function calibration equation, and detector temperature Degree, obtains nine parameters of quadratic function temperature drift equation;
Step s4, demarcates according to nine parameters of quadratic function temperature drift equation, quadratic function temperature drift equation and quadratic function Equation, and the gradation of image of pre-acquiring target object, determine the corresponding temperature of target object.
In step s2, when temperature and gradation of image are fitted, number is gathered using Fu Jirui movement and flir camera lens According to, concrete gradation of image and blackbody temperature data, as shown in table 1:
Table 1 gradation of image blackbody temperature data
Wherein, overstriking font component is that the partial dot chosen does extrapolation measure, respectively with standard blackbody temperature and image ash Spend for independent variable, carry out 2 times, 3 times, 4 matchings, the relatively degree of fitting of each curve.In conjunction with Fig. 2, Fig. 2 shows one kind with standard Blackbody temperature is the multiple matched curve figure of independent variable, and transverse axis is standard blackbody temperature, and the longitudinal axis is gradation of image, wherein, solid line It is the multiple matched curve with standard blackbody temperature as independent variable, dotted line is the matching that selected point does extrapolation measure, curve l211 It is the quadratic function matched curve with standard blackbody temperature as independent variable, curve l212 is that the quadratic function matching of extrapolation measure is bent Line, curve l221 is the cubic function matched curve with standard blackbody temperature as independent variable, and curve l222 is the three of extrapolation measure Secondary Function Fitting curve, curve l231 is the biquadratic function matched curve with standard blackbody temperature as independent variable, and curve l232 is The biquadratic function matched curve of extrapolation measure is it is known that do the effect of quadratic function matching during extrapolation measure preferably, such as in figure is red Curve.In conjunction with Fig. 3, Fig. 3 shows a kind of multiple matched curve figure with gradation of image as independent variable, and the longitudinal axis is standard blackbody temperature Degree, transverse axis is gradation of image, and wherein, solid line is the multiple matched curve with gradation of image as independent variable, and dotted line is that selected point is done The matching of extrapolation measure, curve l311 is the quadratic function matched curve with gradation of image as independent variable, and curve l312 is extrapolation The quadratic function matched curve of test, curve l321 is the cubic function matched curve with gradation of image as independent variable, curve L322 is the cubic function matched curve of extrapolation measure, and curve l331 is that the biquadratic function matching with gradation of image as independent variable is bent Line, curve l332 is the biquadratic function matched curve of extrapolation measure it is known that doing the effect of quadratic function matching during extrapolation measure Good, but there is larger difference in high-temperature extrapolation.Therefore, the method is fitted from quadratic function, forms quadratic function mark Determine equation, error in data is little.Quadratic function calibration equation, specifically with target temperature as independent variable, simulates quadratic polynomial every Coefficient.
In step s3, detector temperature change does not change the type function of calibration curve, calibration curve only do stretching and Translation, impact be calibration curve function each term coefficient, each term coefficient of quadratic function calibration equation can be detector The function of temperature, and each term coefficient of quadratic function calibration equation is in monotonic relationshi with the relation of detector temperature, below by Relatively each term coefficient temperature drift equation once fitting and the difference of quadratic fit, to judge quality.According to 3 DEG C, 13.3 DEG C, 24.8 DEG C, 28.6 DEG C, 36.7 DEG C, the data fitting calibrating equation under 6 environment such as 49.5 DEG C, calculate 13 DEG C, the black matrix at 29.6 DEG C Temperature, and compare calculation error with actual temperature, it is shown in Table 2.
Table 2 different detector temperature difference blackbody temperature next, quadratic fit error
From table 2, poorer than the effect of quadratic linear matching using once linear matching, the method determines using secondary letter Number matching each coefficient temperature drift equation.
As shown from the above technical solution, the present embodiment thermal infrared imager scaling method, using biquadratic matching nine parameter Standardizition, scaling method accuracy is high, and error is little.The method is according to the temperature of multiple standard blackbody and the corresponding figure of standard blackbody As gray scale, using quadratic function calibration equation, it is fitted, in extrapolation measure, accuracy is high, and error in data is little.Meanwhile, two Each term coefficient of secondary function calibration equation can fully combine detector temperature, is fitted using quadratic function temperature drift equation, And determining nine parameters of temperature drift function, mean error is little, and degree of accuracy is high.It is fitted using biquadratic function, and determine temperature Nine parameters of drift equation, you can effectively thermal infrared imager is demarcated, improve equipment dependability and systematic function, in conjunction with The gradation of image of target object, can determine the temperature of target object effectively.
Therefore, the present embodiment thermal infrared imager scaling method, can demarcate to thermal infrared imager effectively, and letter Change processing procedure, improve equipment performance.
Specifically, the present embodiment thermal infrared imager scaling method is obtaining each term coefficient of quadratic function calibration equation When, according to following quadratic function calibration equation, and the temperature of standard blackbody and the corresponding gradation of image of this standard blackbody, determine The constant of the quadratic coefficients of quadratic function calibration equation, the coefficient of first order of quadratic function calibration equation and quadratic function calibration equation ,
D=a*t2+b*t+c
Wherein, d is the corresponding gradation of image of standard blackbody, and t is the temperature of standard blackbody, and a is quadratic function calibration equation Quadratic coefficients, b be quadratic function calibration equation coefficient of first order, c be quadratic function calibration equation constant term, quadratic function Each term coefficient of calibration equation include the quadratic coefficients of quadratic function calibration equation, the coefficient of first order of quadratic function calibration equation, The constant term of quadratic function calibration equation.Using above-mentioned quadratic function calibration equation, and the temperature of combined standard black matrix and standard The corresponding gradation of image of black matrix, you can obtain each term coefficient of quadratic function calibration equation, meet thermal infrared imager radiation calibration Ultimate principle.In actual testing temperature, obtain the gradation of image d of target object0, now, the equation that the method is formed is: d0 =a*t2+ b*t+c, then pass through equation below, you can obtain the temperature of target object,
t ( d 0 ) = - b + b 2 - 4 a ( c - d 0 ) 2 a
And, in conjunction with Fig. 4, Fig. 4 shows the relation between the gradation of image d of infrared hybrid optical system and target temperature t, two Person is in s type curve, is respectively present a saturation point in high/low temperature, i.e. high temperature saturation and low temperature saturation, because thermal imaging system receives Heat radiation is independent variable, and gradation of image is dependent variable, so, when quadratic function calibration equation is fitted, preferably adopt target Temperature is as independent variable, and carries out segmentation choice according to the saturation trend of gray value to relation curve, concrete as shown in figure 5, Fig. 5 shows that demarcation temperature spot is chosen and image gray levels accept or reject schematic diagram, when carrying out calibration experiment, can be in conjunction with Fig. 5 Gray value interval and blackbody temperature point, the temperature of setting standard blackbody.
Specifically, the present embodiment thermal infrared imager scaling method is determining nine parameters of quadratic function temperature drift equation When, the quadratic coefficients according to following quadratic coefficients temperature drift equation and quadratic function calibration equation and detector temperature, obtain two The constant of the quadratic coefficients of ordered coefficients temperature drift equation, the coefficient of first order of quadratic coefficients temperature drift equation and quadratic coefficients temperature drift equation ,
A=a2*tx 2+a1*tx+a0
Wherein, a is the quadratic coefficients of quadratic function calibration equation, txFor detector temperature, a2For quadratic coefficients temperature drift equation Quadratic coefficients, a1For the coefficient of first order of quadratic coefficients temperature drift equation, a0For the constant term of quadratic coefficients temperature drift equation,
Coefficient of first order according to following coefficient of first order temperature drift equation and quadratic function calibration equation and detector temperature, Obtain quadratic coefficients, the coefficient of first order of coefficient of first order temperature drift equation and the coefficient of first order temperature drift equation of coefficient of first order temperature drift equation Constant term,
B=b2*tx 2+b1*tx+b0
Wherein, b is the coefficient of first order of quadratic function calibration equation, b2For the quadratic coefficients of coefficient of first order temperature drift equation, b1For The coefficient of first order of coefficient of first order temperature drift equation, b0For the constant term of coefficient of first order temperature drift equation,
Constant term according to following constant term temperature drift equation and quadratic function calibration equation and detector temperature, obtain The constant term of the quadratic coefficients of constant term temperature drift equation, the coefficient of first order of constant term temperature drift equation and constant term temperature drift equation,
C=c2*tx 2+c1*tx+c0
Wherein, c is the constant term of quadratic function calibration equation, c2For the quadratic coefficients of constant term temperature drift equation, c1For constant The coefficient of first order of item temperature drift equation, c0For the constant term of constant term temperature drift equation, above-mentioned quadratic function temperature drift equation includes secondary Coefficient temperature drift equation, coefficient of first order temperature drift equation and constant term temperature drift equation, nine parameters of above-mentioned quadratic function temperature drift equation Including the quadratic coefficients of quadratic coefficients temperature drift equation, the coefficient of first order of quadratic coefficients temperature drift equation, quadratic coefficients temperature drift equations Constant term, the quadratic coefficients of coefficient of first order temperature drift equation, the coefficient of first order of coefficient of first order temperature drift equation, coefficient of first order temperature drift equation Constant term, the quadratic coefficients of constant term temperature drift equation, the coefficient of first order of constant term temperature drift equation and constant term temperature drift equation Constant term.
It is fitted between detector temperature and each term coefficient of quadratic function calibration equation using quadratic function temperature drift equation Relation, take into full account the impact that detector temperature is brought, improve accuracy and the reliability of the method, and error is little.
Specifically, the present embodiment thermal infrared imager scaling method, in the gradation of image of combining target object, determines target During the corresponding temperature of object, according to nine parameters of quadratic function temperature drift equation, quadratic function temperature drift equation and quadratic function mark Determine equation, and the gradation of image of pre-acquiring target object, by equation below, determine the corresponding temperature of target object,
t ( d 0 ) = - f b ( t x ) + f b ( t x ) 2 - 4 * f a ( t x ) * ( f c ( t x ) - d 0 ) 2 * f a ( t x ) f a ( t x ) = a 2 * t x 2 + a 1 * t x + a 0 f b ( t x ) = b 2 * t x 2 + b 1 * t x + b 0 f c ( t x ) = c 2 * t x 2 + c 1 * t x + c 0
Wherein, t (d0) it is the corresponding temperature of target object, d0Gradation of image for pre-acquiring target object.Thermal infrared imager Thermometric is carried out using above-mentioned equation, the accuracy of equipment is high, highly reliable.Thermal infrared imager systematic function after corrected Accurately, temperature test is accurately and reliably.t(d0) it is alternatively the corresponding temperature of standard blackbody, d0It is alternatively the image ash of standard blackbody Degree.Thermal infrared imager adopts said method, the temperature of the temperature in conjunction with multiple standard blackbody and calculating acquisition, checks red further The systematic function of outer thermal imaging system.
Here describes in detail demarcates flow process: set and step into room temperature to1, ambient temperature and detector temperature stable after, lead to Cross the infrared hybrid optical system that each thermal infrared imager of computer controls gathers series of standards black matrix, record detector temperature simultaneously t1.The number of thermal infrared imager is less than or equal to the number of black matrix.After the completion of collection, set temperature to2 stepping into room, and series The temperature of standard blackbody, same method gathers the infrared hybrid optical system of series of standards black matrix.
Table 3 gathers the temperature data of standard blackbody
The data such as table 4 of single thermal infrared imager collection, according to the standard blackbody temperature under each environment and detector temperature Degree matching, obtains quadratic function calibration equation.I.e. when detector temperature is for t1, the quadratic function calibration equation that matching obtains For: d=a*t2+ b*t+c, determines that each term coefficient is a (1), b (1), c (1);When detector temperature is for t2, matching obtain two Secondary function calibration equation is: d=a*t2+ b*t+c, determines that each term coefficient is a (2), b (2), c (2);It is t3 in detector temperature When, the quadratic function calibration equation that matching obtains is: d=a*t2+ b*t+c, determines that each term coefficient is a (3), b (3), c (3).
To be listed in Table 5 in each parameter set, be detector temperature t further according to achievement in research a, b before, cxLetter Number, can draw specific function expression in table 5 last column.
The single thermal infrared imager calibration curve matching of table 4
Table 5 calibration curve parameter fitting
Final detector temperature tx, blackbody temperature t, the relation of gradation of image d can represent with lower peg model:
d = f a ( t x ) * t 2 + f b ( t x ) * t + f c ( t x ) f a ( t x ) = a 2 * t x 2 + a 1 * t x + a 0 f b ( t x ) = b 2 * t x 2 + b 1 * t x + b 0 f c ( t x ) = c 2 * t x 2 + c 1 * t x + c 0
The inverse function seeking d can get true temperature:
t ( d 0 ) = - f b ( t x ) + f b ( t x ) 2 - 4 * f a ( t x ) * ( f c ( t x ) - d 0 ) 2 * f a ( t x )
Second aspect, a kind of thermal infrared imager caliberating device that second embodiment of the invention is provided, in conjunction with Fig. 6, should Device includes data acquisition module 1, calibration coefficient determining module 2, temperature drift coefficient determination module 3 and target temperature acquisition module 4. Data acquisition module 1 is used for obtaining the temperature of standard blackbody and gradation of image corresponding with standard blackbody, and detector temperature, The number of standard blackbody is multiple.Calibration coefficient determining module 2 is used for the temperature according to quadratic function calibration equation and standard blackbody Degree and gradation of image corresponding with this standard blackbody, obtain each term coefficient of quadratic function calibration equation.Temperature drift coefficient determines Module 3 is used for each term coefficient according to quadratic function temperature drift equation and quadratic function calibration equation, and detector temperature, obtains Nine parameters of quadratic function temperature drift equation.Target temperature acquisition module 4 is used for nine ginsengs according to quadratic function temperature drift equation Number, quadratic function temperature drift equation and quadratic function calibration equation, and the gradation of image of pre-acquiring target object, determine object The corresponding temperature of body.
As shown from the above technical solution, the present embodiment thermal infrared imager caliberating device, using biquadratic matching nine parameter Standardizition, scaling method accuracy is high, and error is little.Data acquisition module 1 obtains required data when demarcating, and calibration coefficient determines Module 2, according to the temperature of multiple standard blackbody and the corresponding gradation of image of standard blackbody, using quadratic function calibration equation, is carried out Matching, in extrapolation measure, accuracy is high, and error in data is little.
Meanwhile, each term coefficient of quadratic function calibration equation can fully combine detector temperature, and temperature drift coefficient determines mould Block 3 is fitted using quadratic function temperature drift equation, and determines nine parameters of temperature drift function, and mean error is little, and degree of accuracy is high. It is fitted using biquadratic function, and determines nine parameters of temperature drift equation, target temperature acquisition module 4 can be effectively to red Outer thermal imaging system is demarcated, and improves equipment dependability, in conjunction with the gradation of image of target object, can effectively determine object The temperature of body.
Therefore, thermal infrared imager caliberating device of the present invention, can accurately and efficiently demarcate to thermal infrared imager, And precision high level error is little, improve the performance of equipment.
In order to improve the accuracy of the present embodiment thermal infrared imager caliberating device, calibration coefficient determining module 2 further In the temperature according to quadratic function calibration equation and standard blackbody and gradation of image corresponding with this standard blackbody, obtain secondary During each term coefficient of function calibration equation, specifically for: according to following quadratic function calibration equation, and the temperature of standard blackbody The corresponding gradation of image with this standard blackbody, determines the quadratic coefficients of quadratic function calibration equation, quadratic function calibration equation Coefficient of first order and the constant term of quadratic function calibration equation,
D=a*t2+b*t+c
Wherein, d is the corresponding gradation of image of standard blackbody, and t is the temperature of standard blackbody, and a is quadratic function calibration equation Quadratic coefficients, b be quadratic function calibration equation coefficient of first order, c be quadratic function calibration equation constant term, quadratic function Each term coefficient of calibration equation include the quadratic coefficients of quadratic function calibration equation, the coefficient of first order of quadratic function calibration equation, The constant term of quadratic function calibration equation.Calibration coefficient determining module 2 adopts above-mentioned quadratic function calibration equation, and combined standard The temperature of black matrix and the corresponding gradation of image of standard blackbody, you can obtain each term coefficient of quadratic function calibration equation, meet red The ultimate principle of outer thermal imaging system radiation calibration.In actual testing temperature, obtain the gradation of image d of target object0, now obtain Equation is: d0=a*t2+ b*t+c, you can obtain the temperature of target object.
Temperature drift coefficient determination module 3 in each term coefficient according to quadratic function temperature drift equation and quadratic function calibration equation, And detector temperature, during nine parameters of acquisition quadratic function temperature drift equation, specifically for: according to following quadratic coefficients temperature drift The quadratic coefficients of equation and quadratic function calibration equation and detector temperature, obtain the secondary system of quadratic coefficients temperature drift equation The constant term of number, the coefficient of first order of quadratic coefficients temperature drift equation and quadratic coefficients temperature drift equation,
A=a2*tx 2+a1*tx+a0
Wherein, a is the quadratic coefficients of quadratic function calibration equation, txFor detector temperature, a2For quadratic coefficients temperature drift equation Quadratic coefficients, a1For the coefficient of first order of quadratic coefficients temperature drift equation, a0For the constant term of quadratic coefficients temperature drift equation,
Coefficient of first order according to following coefficient of first order temperature drift equation and quadratic function calibration equation and detector temperature, Obtain quadratic coefficients, the coefficient of first order of coefficient of first order temperature drift equation and the coefficient of first order temperature drift equation of coefficient of first order temperature drift equation Constant term,
B=b2*tx 2+b1*tx+b0
Wherein, b is the coefficient of first order of quadratic function calibration equation, txFor detector temperature, b2For coefficient of first order temperature drift equation Quadratic coefficients, b1For the coefficient of first order of coefficient of first order temperature drift equation, b0For the constant term of coefficient of first order temperature drift equation,
Constant term according to following constant term temperature drift equation and quadratic function calibration equation and detector temperature, obtain The constant term of the quadratic coefficients of constant term temperature drift equation, the coefficient of first order of constant term temperature drift equation and constant term temperature drift equation,
C=c2*tx 2+c1*tx+c0
Wherein, c is the constant term of quadratic function calibration equation, txFor detector temperature, c2For constant term temperature drift equation two Ordered coefficients, c1For the coefficient of first order of constant term temperature drift equation, c0For the constant term of constant term temperature drift equation,
Quadratic function temperature drift equation includes quadratic coefficients temperature drift equation, coefficient of first order temperature drift equation and constant term temperature drift side Journey, nine parameters of quadratic function temperature drift equation include the quadratic coefficients of quadratic coefficients temperature drift equation, quadratic coefficients temperature drift equation Coefficient of first order, the constant term of quadratic coefficients temperature drift equation, the quadratic coefficients of coefficient of first order temperature drift equation, coefficient of first order temperature drift side The coefficient of first order of journey, the constant term of coefficient of first order temperature drift equation, the quadratic coefficients of constant term temperature drift equation, constant term temperature drift equation Coefficient of first order and constant term temperature drift equation constant term.
Temperature drift coefficient determination module 3 is fitted detector temperature using quadratic function temperature drift equation and quadratic function is demarcated Relation between each term coefficient of equation, takes into full account the impact that detector temperature is brought, and improves accuracy and the reliability of the method Property, and error is little.
Target temperature acquisition module 4 nine parameters according to quadratic function temperature drift equation, quadratic function temperature drift equation and Quadratic function calibration equation, and the gradation of image of pre-acquiring target object, when determining the corresponding temperature of target object, concrete use In: according to nine parameters of quadratic function temperature drift equation, quadratic function temperature drift equation and quadratic function calibration equation, and obtain in advance Take the gradation of image of target object, by equation below, determine the corresponding temperature of target object,
t ( d 0 ) = - f b ( t x ) + f b ( t x ) 2 - 4 * f a ( t x ) * ( f c ( t x ) - d 0 ) 2 * f a ( t x ) f a ( t x ) = a 2 * t x 2 + a 1 * t x + a 0 f b ( t x ) = b 2 * t x 2 + b 1 * t x + b 0 f c ( t x ) = c 2 * t x 2 + c 1 * t x + c 0
Wherein, t (d0) it is the corresponding temperature of target object, d0Gradation of image for pre-acquiring target object.
Target temperature acquisition module 4 carries out thermometric using above-mentioned equation, and parameters are calibration coefficient determining module 2 He The data that temperature drift coefficient determination module 3 process refers to, when carrying out temperature test, accuracy is lifted, and reliability is high.
The third aspect, a kind of thermal infrared imager calibration system that third embodiment of the invention is provided, this system includes Thermal infrared imager, multiple standard blackbody and seal box.Thermal infrared imager is located at seal box, and thermal infrared imager includes being sequentially connected Flir camera lens, detector and processor, filr camera lens is used for the infrared radiation collection of standard blackbody or target object to detection Device, detector is used for gathering the infra-red radiation of standard blackbody, and is converted to the infrared hybrid optical system of standard blackbody, or collection target The infra-red radiation of object, and be converted to the infrared hybrid optical system of target object, and by the infrared hybrid optical system of standard blackbody and/ Or the infrared hybrid optical system of target object sends to processor.Processor is used for obtaining temperature and and the standard blackbody of standard blackbody The gradation of image of corresponding infrared hybrid optical system, and detector temperature, according to the temperature of quadratic function calibration equation and standard blackbody Degree and gradation of image corresponding with this standard blackbody, obtain each term coefficient of secondary calibration equation, always according to quadratic function temperature Drift equation and each term coefficient of secondary calibration equation, and detector temperature, obtain nine parameters of quadratic function temperature drift equation, And nine parameters according to quadratic function temperature drift equation, quadratic function temperature drift equation and quadratic function calibration equation, and in advance Obtain the gradation of image of target object infrared hybrid optical system, determine the corresponding temperature of target object.The side wall of seal box is provided with out Mouth, opening and flir shots match, and it is provided with germanium wafer.
As shown from the above technical solution, the present embodiment thermal infrared imager calibration system, using multiple standard black body radiation Source is corrected to the measurement result of thermal infrared imager, and accuracy is high.Seal box can be placed in refrigerator or high-temperature cabinet, is conveniently Thermal infrared imager provides a metastable experimental situation, and germanium wafer can effectively prevent flir camera lens from producing solidifying mist.Infrared Thermal imaging system accurately and effectively gathers the infrared hybrid optical system of standard blackbody or target object, determine quadratic function calibration equation, two Secondary function temperature drift equation and nine parameters of quadratic function temperature drift equation, can correct the systematicness of thermal infrared imager exactly Can, be conducive to improving the reliability of infrared thermal imagery instrument system.Meanwhile, the thermal infrared imager after corrected, combining target thing The gradation of image of the infrared hybrid optical system of body, can effectively obtain the temperature of target object.
Therefore, the present embodiment thermal infrared imager calibration system, accurately and efficiently can enter rower to thermal infrared imager Fixed, and simplify processing procedure, precision high level error is little, improves the performance of equipment.
From the angle avoiding and evading veiling glare, by comparing the degree of susceptibility to veiling glare for the camera lens, to select Camera lens.Here selects flir camera lens, Beijing Hong Yuan camera lens and Nanjing wavelength photocell lenses, the electric furnace that will be energized as substantial radiation source, It is placed on turntable, by the visual field area progressing into thermal infrared imager outside the area of visual field, observe and record the change of infrared hybrid optical system Change process.Infrared hybrid optical system according to three kinds of camera lens collections and the gray value of each infrared hybrid optical system, as shown in fig. 7, transverse axis For infrared hybrid optical system sequence number, the longitudinal axis is gray value, and as shown in Figure 7, Nanjing wavelength photocell lenses are better than Beijing Hong Yuan mirror Head, but also have different with flir camera lens, flir camera lens roguing light-dispersing effect is best, therefore, the present embodiment thermal infrared imager Adopt flir camera lens with calibration system, at utmost to reduce the impact to picture quality for the veiling glare.
When being demarcated using this system, standard blackbody temperature setting is: 8 DEG C, 30 DEG C, 50 DEG C, 90 DEG C, 130 DEG C, Its temperature of infrared heat image instrument measuring is respectively 8.5 DEG C, 31 DEG C, 49.2 DEG C, 92.6 DEG C, 116 DEG C.Thermal infrared imager and black matrix distance 3m, and keep optical center just right, enables thermal infrared imager obtain standard blackbody clearly as and this is apart from energy after focusing Guarantee to ignore the air radiation impact between thermal infrared imager and standard blackbody.
When obtaining low temperature environment, thermal infrared imager is energized, refrigerator integrally put into by seal box, takes after its temperature stabilization Go out.Due to have seal box isolate, in the short time after taking-up, thermal infrared imager environment will not suddenly change, have time enough Collection standard blackbody image.Slow with laboratory heat exchanging process, in seal box, environment is relatively stable, same in its temperature-rise period Sample can gather black matrix image.Solidifying mist is had on germanium wafer during low temperature, wipes collection image after water smoke with lens paper.
When obtaining hot environment, thermal infrared imager is energized, as in high-temperature cabinet, after stablizing, collection standard is black for seal box Body image.
When obtaining normal temperature environment, room air-conditioner temperature, under different temperatures environment, gathered data such as table 6 institute are demarcated in setting Show,
Corresponding detector temperature at a temperature of table 6 varying environment
It is assumed that detector temperature parameter and ambient temperature are linear after detector energising is stable, in follow-up matching Measure the temperature drift that detector is caused by ambient temperature with during measurement with detector temperature parameter.
According to previous experiments conclusion, in matching gray scale temperature curve, with temperature as independent variable, using quadratic function Matching.Each detector temperature txUnder simulate a quadratic function calibration equation:
D=fa (tx)*t2+fb(tx)*t+fc(tx)
Here, d is the gradation of image of infrared hybrid optical system, t is the temperature of standard blackbody, fa (tx), fb (tx), fc (tx) It is each term coefficient of quadratic function calibration equation, each term coefficient is detector temperature t againxFunction.By a series of detector Temperature can get following coefficient sequence:
fa(t1) fa(t2) fa(t3)
fb(t1) fb(t2) fb(t3)
fc(t1) fc(t2) fc(t3)
Temperature relation matching further according to each coefficient sequence and detector draws the function of each coefficient:
fa(tx)=a2*tx 2+a1*tx+a0Quadratic coefficients temperature drift equation
fb(tx)=b2*tx 2+b1*tx+b0Coefficient of first order temperature drift equation
fc(tx)=c2*tx 2+c1*tx+c0Constant term temperature drift equation
Wherein, txIt is detector temperature, can use three values t1, t2, t3, fa also has three values, it is possible to simulating a1, a2, a0, and then draw every fitting coefficient a2…c0:
a2a1a0
b2b1b0
c2c1c0
Simulating a2…c0Afterwards, you can realize removing temperature drift thermometric, obtain temperature computation after detector temperature and gradation of image Process is as follows:
According to detector temperature tx, calculate demarcate fit equation each term coefficient:
fa(tx), fb (tx), fc (tx)
According to gradation of image d0With quadratic function calibration equation, it is calculated target temperature:
t ( d 0 ) = - f b ( t x ) + f b ( t x ) 2 - 4 * f a ( t x ) * ( f c ( t x ) - d 0 ) 2 * f a ( t x )
Choose wherein 3 DEG C, 13.3 DEG C, 24.8 DEG C, 28.6 DEG C, 36.7 DEG C, in 49.5 DEG C of environment, its gradation of image is black As shown in figure 8, transverse axis is standard blackbody temperature, the longitudinal axis is gradation of image, wherein, l1 to body temperature degrees of data fitting calibrating equation curve For 3 DEG C of matched curve of environment, l2 is the matched curve of 13.3 DEG C of environment, and l3 is the matched curve of 24.8 DEG C of environment, and l4 is ring The matched curve that 28.6 DEG C of border, l5 is the matched curve of 36.7 DEG C of environment, and l6 is the matched curve of 49.5 DEG C of environment.Can by Fig. 8 Know, the regularity that quadratic function calibration equation changes with detector temperature is stronger, draw the demarcation matching under each detector temperature Equation coefficient, then matching each term coefficient temperature drift equation, specifically combine Fig. 9~Figure 11, from Fig. 9~Figure 11, fit equation is each During term coefficient, curve is close linear, and combines table 2, this system poorer than the effect of quadratic linear matching using once linear matching Determine and adopt quadratic function matching each coefficient temperature drift equation.
Here takes different ambient temperature points to carry out fitted calibration curve, and does thermometric test.
Experiment one: fixed point or ambient temperature: demarcate in this condition for 3 DEG C, 24.8 DEG C, 36.7 DEG C, 49.3 DEG C, obtain Calibration curve;
Point for measuring temperature or movement temperature: 13 DEG C, 29.6 DEG C;
Temperature-measuring results: as shown in table 7
Temperature-measuring results tested by table 7
Experiment two: fixed point or ambient temperature: 3 DEG C, 24.8 DEG C, 36.7 DEG C;
Point for measuring temperature or movement temperature: 13 DEG C, 29.6 DEG C, 49.5 DEG C
Temperature-measuring results: as shown in table 8
Two temperature-measuring results tested by table 8
From above-mentioned experimental result: using this calibration system, within 1 DEG C, this system disclosure satisfy that survey to mean error Temperature is required with thermal imaging, and the less precision of detector temperature is higher, can work below room temperature.
Should be respectively with standard blackbody temperature and detector temperature as independent variable during calibration experiment, calibration equation and each term coefficient Temperature drift equation all should take quadratic function matching.Black matrix quantity preferably selects 5,6, at least needs four ambient temperature points during demarcation, temperature Degree scope cover will width, example: less than 5 DEG C, 20 DEG C, 35 DEG C, be higher than 50 DEG C.
Finally it is noted that various embodiments above, only in order to technical scheme to be described, is not intended to limit;To the greatest extent Pipe has been described in detail to the present invention with reference to foregoing embodiments, it will be understood by those within the art that: its according to So the technical scheme described in foregoing embodiments can be modified, or wherein some or all of technical characteristic is entered Row equivalent;And these modifications or replacement, do not make the essence of appropriate technical solution depart from various embodiments of the present invention technology The scope of scheme, it all should be covered in the middle of the claim of the present invention and the scope of description.

Claims (9)

1. a kind of thermal infrared imager scaling method is it is characterised in that include:
Obtain the temperature of standard blackbody and gradation of image corresponding with described standard blackbody, and detector temperature, described standard The number of black matrix is multiple;
Temperature according to quadratic function calibration equation and described standard blackbody and gradation of image corresponding with this standard blackbody, obtain Obtain each term coefficient of described quadratic function calibration equation;
According to each term coefficient of quadratic function temperature drift equation and described quadratic function calibration equation, and described detector temperature, Obtain nine parameters of described quadratic function temperature drift equation;
Demarcated according to nine parameters of described quadratic function temperature drift equation, described quadratic function temperature drift equation and described quadratic function Equation, and the gradation of image of pre-acquiring target object, determine the corresponding temperature of described target object.
2. according to claim 1 thermal infrared imager scaling method it is characterised in that
The described temperature according to quadratic function calibration equation and standard blackbody and gradation of image corresponding with this standard blackbody, obtain Obtain each term coefficient of quadratic function calibration equation, specifically include:
According to following quadratic function calibration equation, and the temperature of described standard blackbody and this standard blackbody corresponding image ash Degree, determines quadratic coefficients, the coefficient of first order of described quadratic function calibration equation and described two of described quadratic function calibration equation The constant term of secondary function calibration equation,
D=a*t2+b*t+c
Wherein, d is the corresponding gradation of image of described standard blackbody, and t is the temperature of described standard blackbody, and a is described quadratic function The quadratic coefficients of calibration equation, b is the coefficient of first order of described quadratic function calibration equation, and c is described quadratic function calibration equation Constant term, each term coefficient of described quadratic function calibration equation includes the quadratic coefficients of described quadratic function calibration equation, institute State coefficient of first order, the constant term of described quadratic function calibration equation of quadratic function calibration equation.
3. according to claim 2 thermal infrared imager scaling method it is characterised in that
Described each term coefficient according to quadratic function temperature drift equation and quadratic function calibration equation, and detector temperature, obtain Nine parameters of quadratic function temperature drift equation, specifically include:
Quadratic coefficients according to following quadratic coefficients temperature drift equation and described quadratic function calibration equation and described detector temperature Degree, obtains quadratic coefficients, the coefficient of first order and described two of described quadratic coefficients temperature drift equation of described quadratic coefficients temperature drift equation The constant term of ordered coefficients temperature drift equation,
A=a2*tx 2+a1*tx+a0
Wherein, a is the quadratic coefficients of described quadratic function calibration equation, txFor described detector temperature, a2For described quadratic coefficients The quadratic coefficients of temperature drift equation, a1For the coefficient of first order of described quadratic coefficients temperature drift equation, a0For described quadratic coefficients temperature drift equation Constant term,
Coefficient of first order according to following coefficient of first order temperature drift equation and described quadratic function calibration equation and described detector temperature Degree, obtains quadratic coefficients, the coefficient of first order and described of described coefficient of first order temperature drift equation of described coefficient of first order temperature drift equation The constant term of ordered coefficients temperature drift equation,
B=b2*tx 2+b1*tx+b0
Wherein, b is the coefficient of first order of described quadratic function calibration equation, b2For the quadratic coefficients of described coefficient of first order temperature drift equation, b1For the coefficient of first order of described coefficient of first order temperature drift equation, b0For the constant term of described coefficient of first order temperature drift equation,
Constant term according to following constant term temperature drift equation and described quadratic function calibration equation and described detector temperature, Obtain quadratic coefficients, the coefficient of first order of described constant term temperature drift equation and the described constant term temperature drift of described constant term temperature drift equation The constant term of equation,
C=c2*tx 2+c1*tx+c0
Wherein, c is the constant term of described quadratic function calibration equation, c2For the quadratic coefficients of described constant term temperature drift equation, c1For The coefficient of first order of described constant term temperature drift equation, c0For the constant term of described constant term temperature drift equation,
Described quadratic function temperature drift equation include described quadratic coefficients temperature drift equation, described coefficient of first order temperature drift equation and described often Several temperature drift equations, nine parameters of described quadratic function temperature drift equation include the secondary system of described quadratic coefficients temperature drift equation The coefficient of first order of several, described quadratic coefficients temperature drift equation, the constant term of described quadratic coefficients temperature drift equation, described coefficient of first order temperature The drift quadratic coefficients of equation, the coefficient of first order of described coefficient of first order temperature drift equation, the constant term of described coefficient of first order temperature drift equation, The quadratic coefficients of described constant term temperature drift equation, the coefficient of first order of described constant term temperature drift equation and described constant term temperature drift equation Constant term.
4. according to claim 3 thermal infrared imager scaling method it is characterised in that
Described nine parameters according to quadratic function temperature drift equation, quadratic function temperature drift equation and quadratic function calibration equation, with And the gradation of image of pre-acquiring target object, determine the corresponding temperature of target object, specifically include:
Demarcated according to nine parameters of described quadratic function temperature drift equation, described quadratic function temperature drift equation and described quadratic function Equation, and the gradation of image of described pre-acquiring target object, by equation below, determine the corresponding temperature of described target object Degree,
t ( d 0 ) = - f b ( t x ) + f b ( t x ) 2 - 4 * f a ( t x ) * ( f c ( t x ) - d 0 ) 2 * f a ( t x ) f a ( t x ) = a 2 * t x 2 + a 1 * t x + a 0 f b ( t x ) = b 2 * t x 2 + b 1 * t x + b 0 f c ( t x ) = c 2 * t x 2 + c 1 * t x + c 0
Wherein, t (d0) it is the corresponding temperature of described target object, d0Gradation of image for described pre-acquiring target object.
5. a kind of thermal infrared imager caliberating device is it is characterised in that include:
Data acquisition module, for obtaining the temperature of standard blackbody and gradation of image corresponding with described standard blackbody, and visits Survey device temperature, the number of described standard blackbody is multiple;
Calibration coefficient determining module, for the temperature according to quadratic function calibration equation and described standard blackbody and with this standard The corresponding gradation of image of black matrix, obtains each term coefficient of described quadratic function calibration equation;
Temperature drift coefficient determination module, for each term system according to quadratic function temperature drift equation and described quadratic function calibration equation Number, and described detector temperature, obtain nine parameters of described quadratic function temperature drift equation;
Target temperature acquisition module, for according to nine parameters of described quadratic function temperature drift equation, described quadratic function temperature drift Equation and described quadratic function calibration equation, and the gradation of image of pre-acquiring target object, determine that described target object corresponds to Temperature.
6. according to claim 5 thermal infrared imager caliberating device it is characterised in that
Described calibration coefficient determining module is in the temperature according to quadratic function calibration equation and standard blackbody and black with this standard The corresponding gradation of image of body, during each term coefficient of acquisition quadratic function calibration equation, specifically for: according to following quadratic function mark Determine equation, and the temperature of described standard blackbody and the corresponding gradation of image of this standard blackbody, determine that described quadratic function is demarcated The constant term of the quadratic coefficients of equation, the coefficient of first order of described quadratic function calibration equation and described quadratic function calibration equation,
D=a*t2+b*t+c
Wherein, d is the corresponding gradation of image of described standard blackbody, and t is the temperature of described standard blackbody, and a is described quadratic function The quadratic coefficients of calibration equation, b is the coefficient of first order of described quadratic function calibration equation, and c is described quadratic function calibration equation Constant term, each term coefficient of described quadratic function calibration equation includes the quadratic coefficients of described quadratic function calibration equation, institute State coefficient of first order, the constant term of described quadratic function calibration equation of quadratic function calibration equation.
7. according to claim 6 thermal infrared imager caliberating device it is characterised in that
Described temperature drift coefficient determination module in each term coefficient according to quadratic function temperature drift equation and quadratic function calibration equation, with And detector temperature, during nine parameters of acquisition quadratic function temperature drift equation, specifically for: according to following quadratic coefficients temperature drift side The quadratic coefficients of journey and described quadratic function calibration equation and described detector temperature, obtain described quadratic coefficients temperature drift side The constant term of the quadratic coefficients of journey, the coefficient of first order of described quadratic coefficients temperature drift equation and described quadratic coefficients temperature drift equation,
A=a2*tx 2+a1*tx+a0
Wherein, a is the quadratic coefficients of described quadratic function calibration equation, txFor described detector temperature, a2For described quadratic coefficients The quadratic coefficients of temperature drift equation, a1For the coefficient of first order of described quadratic coefficients temperature drift equation, a0For described quadratic coefficients temperature drift equation Constant term,
Coefficient of first order according to following coefficient of first order temperature drift equation and described quadratic function calibration equation and described detector temperature Degree, obtains quadratic coefficients, the coefficient of first order and described of described coefficient of first order temperature drift equation of described coefficient of first order temperature drift equation The constant term of ordered coefficients temperature drift equation,
B=b2*tx 2+b1*tx+b0
Wherein, b is the coefficient of first order of described quadratic function calibration equation, txFor described detector temperature, b2For described coefficient of first order The quadratic coefficients of temperature drift equation, b1For the coefficient of first order of described coefficient of first order temperature drift equation, b0For described coefficient of first order temperature drift equation Constant term,
Constant term according to following constant term temperature drift equation and described quadratic function calibration equation and described detector temperature, Obtain quadratic coefficients, the coefficient of first order of described constant term temperature drift equation and the described constant term temperature drift of described constant term temperature drift equation The constant term of equation,
C=c2*tx 2+c1*tx+c0
Wherein, c is the constant term of described quadratic function calibration equation, txFor described detector temperature, c2For described constant term temperature drift The quadratic coefficients of equation, c1For the coefficient of first order of described constant term temperature drift equation, c0Constant for described constant term temperature drift equation ,
Described quadratic function temperature drift equation include described quadratic coefficients temperature drift equation, described coefficient of first order temperature drift equation and described often Several temperature drift equations, nine parameters of described quadratic function temperature drift equation include the secondary system of described quadratic coefficients temperature drift equation The coefficient of first order of several, described quadratic coefficients temperature drift equation, the constant term of described quadratic coefficients temperature drift equation, described coefficient of first order temperature The drift quadratic coefficients of equation, the coefficient of first order of described coefficient of first order temperature drift equation, the constant term of described coefficient of first order temperature drift equation, The quadratic coefficients of described constant term temperature drift equation, the coefficient of first order of described constant term temperature drift equation and described constant term temperature drift equation Constant term.
8. according to claim 7 thermal infrared imager caliberating device it is characterised in that
Described target temperature acquisition module is in nine parameters according to quadratic function temperature drift equation, quadratic function temperature drift equation and two Secondary function calibration equation, and the gradation of image of pre-acquiring target object, when determining the corresponding temperature of target object, concrete use In: demarcated according to nine parameters of described quadratic function temperature drift equation, described quadratic function temperature drift equation and described quadratic function Equation, and the gradation of image of described pre-acquiring target object, by equation below, determine the corresponding temperature of described target object Degree,
t ( d 0 ) = - f b ( t x ) + f b ( t x ) 2 - 4 * f a ( t x ) * ( f c ( t x ) - d 0 ) 2 * f a ( t x ) f a ( t x ) = a 2 * t x 2 + a 1 * t x + a 0 f b ( t x ) = b 2 * t x 2 + b 1 * t x + b 0 f c ( t x ) = c 2 * t x 2 + c 1 * t x + c 0
Wherein, t (d0) it is the corresponding temperature of described target object, d0Gradation of image for described pre-acquiring target object.
9. a kind of thermal infrared imager calibration system is it is characterised in that include:
Thermal infrared imager, multiple standard blackbody and seal box,
Described thermal infrared imager is located at described seal box, and described thermal infrared imager includes flir camera lens, the detector being sequentially connected And processor,
Described filr camera lens is used for the infrared radiation collection of described standard blackbody or target object to described detector,
Described detector is used for gathering the infra-red radiation of described standard blackbody, and is converted to the infrared gray-scale maps of described standard blackbody Picture, or gather the infra-red radiation of described target object, and be converted to the infrared hybrid optical system of described target object, and will be described The infrared hybrid optical system of the infrared hybrid optical system of standard blackbody and/or described target object sends to described processor,
Described processor is used for obtaining the temperature of described standard blackbody and the figure of infrared hybrid optical system corresponding with described standard blackbody As gray scale, and described detector temperature, the temperature according to quadratic function calibration equation and described standard blackbody and with this mark The corresponding gradation of image of quasi- black matrix, obtains each term coefficient of described secondary calibration equation, always according to quadratic function temperature drift equation and Each term coefficient of described secondary calibration equation, and described detector temperature, obtain nine of described quadratic function temperature drift equation Parameter, and nine parameters according to described quadratic function temperature drift equation, described quadratic function temperature drift equation and described secondary letter Number calibration equation, and the gradation of image of target object infrared hybrid optical system described in pre-acquiring, determine that described target object corresponds to Temperature,
The side wall of described seal box is provided with opening, described opening and described flir shots match, and is provided with germanium wafer.
CN201610742769.XA 2016-08-26 2016-08-26 Calibration method, device and system for infrared thermal imager Pending CN106352987A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610742769.XA CN106352987A (en) 2016-08-26 2016-08-26 Calibration method, device and system for infrared thermal imager

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610742769.XA CN106352987A (en) 2016-08-26 2016-08-26 Calibration method, device and system for infrared thermal imager

Publications (1)

Publication Number Publication Date
CN106352987A true CN106352987A (en) 2017-01-25

Family

ID=57855320

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610742769.XA Pending CN106352987A (en) 2016-08-26 2016-08-26 Calibration method, device and system for infrared thermal imager

Country Status (1)

Country Link
CN (1) CN106352987A (en)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107124600A (en) * 2017-04-26 2017-09-01 北京环境特性研究所 One kind is based on DMD infrared projections equipment and its quantitative output intent
CN107941348A (en) * 2017-11-24 2018-04-20 广东电网有限责任公司电力科学研究院 A kind of infrared measurement of temperature Calibration Method
CN110006535A (en) * 2019-04-15 2019-07-12 中科和光(天津)应用激光技术研究所有限公司 A kind of temp measuring method based on infrared thermal imaging technique
CN110793635A (en) * 2019-10-10 2020-02-14 深圳市景阳信息技术有限公司 Thermal imaging temperature measurement method and device and terminal equipment
CN112050948A (en) * 2020-08-05 2020-12-08 沈阳上博智像科技有限公司 Non-barrier infrared temperature measurement method based on detector temperature drift model
CN113465752A (en) * 2021-09-02 2021-10-01 北方夜视科技(南京)研究院有限公司 Multi-point calibration correction method and system for infrared self-adaptive selection of calibration point
CN114441043A (en) * 2022-01-21 2022-05-06 武汉新朗光电科技有限公司 Temperature correction method of uncooled infrared thermometer and thermometer
JP7318049B1 (en) 2022-03-31 2023-07-31 國家中山科學研究院 How to correct the temperature measurement of a black body under temperature fluctuations of the thermosensitive element

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101183026A (en) * 2007-11-16 2008-05-21 中华人民共和国中山出入境检验检疫局 Infrared temperature measurement online self-calibrating equipment and method thereof
CN102564598A (en) * 2012-01-06 2012-07-11 电子科技大学 Method for calibrating and correcting temperature measurement of infrared probe as well as corresponding temperature measuring method
CN102798475A (en) * 2012-08-02 2012-11-28 电子科技大学 Method for determining temperature of object according to standard infrared video image
CN102818636A (en) * 2012-09-03 2012-12-12 北京理工大学 Radiometric calibration method for low-temperature measurement of thermal infrared imager
US20130193326A1 (en) * 2012-01-27 2013-08-01 Bae Systems Information And Electronic Systems Integration Inc. Non-linear calibration of night vision microbolometer
CN104614076A (en) * 2015-01-27 2015-05-13 天津大学 Precision calibrating method of infrared thermometer with wide temperature range

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101183026A (en) * 2007-11-16 2008-05-21 中华人民共和国中山出入境检验检疫局 Infrared temperature measurement online self-calibrating equipment and method thereof
CN102564598A (en) * 2012-01-06 2012-07-11 电子科技大学 Method for calibrating and correcting temperature measurement of infrared probe as well as corresponding temperature measuring method
US20130193326A1 (en) * 2012-01-27 2013-08-01 Bae Systems Information And Electronic Systems Integration Inc. Non-linear calibration of night vision microbolometer
CN102798475A (en) * 2012-08-02 2012-11-28 电子科技大学 Method for determining temperature of object according to standard infrared video image
CN102818636A (en) * 2012-09-03 2012-12-12 北京理工大学 Radiometric calibration method for low-temperature measurement of thermal infrared imager
CN104614076A (en) * 2015-01-27 2015-05-13 天津大学 Precision calibrating method of infrared thermometer with wide temperature range

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
李云红 等: ""红外热像仪外场测温的大气透过率二次标定"", 《光学 精密工程》 *
马丽娜 等: ""红外热像仪标定方法"", 《火力与指挥控制》 *

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107124600A (en) * 2017-04-26 2017-09-01 北京环境特性研究所 One kind is based on DMD infrared projections equipment and its quantitative output intent
CN107941348A (en) * 2017-11-24 2018-04-20 广东电网有限责任公司电力科学研究院 A kind of infrared measurement of temperature Calibration Method
CN110006535A (en) * 2019-04-15 2019-07-12 中科和光(天津)应用激光技术研究所有限公司 A kind of temp measuring method based on infrared thermal imaging technique
CN110793635A (en) * 2019-10-10 2020-02-14 深圳市景阳信息技术有限公司 Thermal imaging temperature measurement method and device and terminal equipment
CN110793635B (en) * 2019-10-10 2020-10-23 深圳市景阳信息技术有限公司 Thermal imaging temperature measurement method and device and terminal equipment
CN112050948A (en) * 2020-08-05 2020-12-08 沈阳上博智像科技有限公司 Non-barrier infrared temperature measurement method based on detector temperature drift model
CN112050948B (en) * 2020-08-05 2021-12-03 沈阳上博智像科技有限公司 Non-barrier infrared temperature measurement method based on detector temperature drift model
CN113465752A (en) * 2021-09-02 2021-10-01 北方夜视科技(南京)研究院有限公司 Multi-point calibration correction method and system for infrared self-adaptive selection of calibration point
CN113465752B (en) * 2021-09-02 2021-11-05 北方夜视科技(南京)研究院有限公司 Multi-point calibration correction method and system for infrared self-adaptive selection of calibration point
CN114441043A (en) * 2022-01-21 2022-05-06 武汉新朗光电科技有限公司 Temperature correction method of uncooled infrared thermometer and thermometer
CN114441043B (en) * 2022-01-21 2022-10-04 武汉新朗光电科技有限公司 Temperature correction method of uncooled infrared thermometer and thermometer
JP7318049B1 (en) 2022-03-31 2023-07-31 國家中山科學研究院 How to correct the temperature measurement of a black body under temperature fluctuations of the thermosensitive element

Similar Documents

Publication Publication Date Title
CN106352987A (en) Calibration method, device and system for infrared thermal imager
US10416076B2 (en) Quantifying gas in passive optical gas imaging
Güttinger et al. Corrections to scaling in the susceptibility of xenon
CN106404352A (en) Method for measuring large field-of-view telescope optical system distortion and field curvature
CN112284687B (en) Imaging simulation system and method suitable for deep space exploration imaging spectrometer
CN110514392B (en) Temperature deformation measurement system and method based on wide-angle lens and image distortion correction
CN206146624U (en) Blind first detection device of thermal infrared hyperspectral imager appearance
US20160165152A1 (en) Drift Correction Method for Infrared Imaging Device
US10070076B2 (en) Drift correction method for infrared imaging device
US20210343005A1 (en) Thermal camera, and method thereof for early diagnosis of infectious diseases
CN111854964A (en) Method and device for measuring body temperature, temperature measuring tool and storage medium
CN107340064A (en) Thermal imaging system heterogeneity evaluation means for correcting based on scan rectangle black matrix
CN208366471U (en) It is a kind of for automatically determining the device of the minimum discernable temperature difference of infrared thermoviewer
Miklavec et al. A large aperture blackbody bath for calibration of thermal imagers
CN110455417B (en) Quantitative measurement error correction method for stray radiation of infrared optical system
CN106855433A (en) A kind of method and device of imaging spectrometer auto-scaling
Hoover Jr et al. What is MRT and how do I get one?
Guesalaga et al. An on-line turbulence profiler for the AOF: on-sky results
US20190154510A1 (en) Method for Determining a Temperature without Contact and Infrared Measuring System
KR102492065B1 (en) Apparatus and method for optical measuring
CN106597375B (en) A kind of imaging spectrometer calculates the method and device of target position
Mayen et al. Depletion curves in cluster lenses: simulations and application to the cluster MS1008-1224
CN103411757B (en) Measuring system based on optical transfer function of endoscope with special spatial frequency and calibrating method thereof
Burks et al. Signal intensity transfer function determination on thermal systems with stray light or scattering present
CN106404177A (en) Infrared scanning measured temperature correcting method

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
AD01 Patent right deemed abandoned

Effective date of abandoning: 20191108

AD01 Patent right deemed abandoned