CN106352987A - Calibration method, device and system for infrared thermal imager - Google Patents
Calibration method, device and system for infrared thermal imager Download PDFInfo
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 46
- 238000012887 quadratic function Methods 0.000 claims abstract description 282
- 230000003287 optical effect Effects 0.000 claims description 26
- 239000011159 matrix material Substances 0.000 claims description 15
- 230000005855 radiation Effects 0.000 claims description 14
- GNPVGFCGXDBREM-UHFFFAOYSA-N germanium atom Chemical compound [Ge] GNPVGFCGXDBREM-UHFFFAOYSA-N 0.000 claims description 6
- 241000872198 Serjania polyphylla Species 0.000 claims 1
- 238000005516 engineering process Methods 0.000 abstract description 4
- 238000001514 detection method Methods 0.000 abstract description 3
- 238000013213 extrapolation Methods 0.000 description 15
- 238000011088 calibration curve Methods 0.000 description 7
- 238000001931 thermography Methods 0.000 description 7
- 238000012360 testing method Methods 0.000 description 6
- 230000000694 effects Effects 0.000 description 5
- 238000002474 experimental method Methods 0.000 description 5
- 230000004313 glare Effects 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 230000005457 Black-body radiation Effects 0.000 description 3
- 238000012888 cubic function Methods 0.000 description 3
- 238000005259 measurement Methods 0.000 description 3
- 239000003595 mist Substances 0.000 description 3
- 230000009897 systematic effect Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 230000000007 visual effect Effects 0.000 description 2
- 241000208340 Araliaceae Species 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 230000036760 body temperature Effects 0.000 description 1
- 238000009529 body temperature measurement Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000002250 progressing effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000007789 sealing Methods 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 239000000779 smoke Substances 0.000 description 1
- 230000006641 stabilisation Effects 0.000 description 1
- 238000011105 stabilization Methods 0.000 description 1
- 230000002277 temperature effect Effects 0.000 description 1
- 238000013519 translation Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J5/00—Radiation pyrometry, e.g. infrared or optical thermometry
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J5/00—Radiation pyrometry, e.g. infrared or optical thermometry
- G01J5/80—Calibration
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
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,
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,
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,
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,
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:
The inverse function seeking d can get true temperature:
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,
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:
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,
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,
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.
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)
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)
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 |
-
2016
- 2016-08-26 CN CN201610742769.XA patent/CN106352987A/en active Pending
Patent Citations (6)
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)
Title |
---|
李云红 等: ""红外热像仪外场测温的大气透过率二次标定"", 《光学 精密工程》 * |
马丽娜 等: ""红外热像仪标定方法"", 《火力与指挥控制》 * |
Cited By (12)
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 |