CN117232661B - Multichannel infrared radiation measurement system and multi-wavelength real-time temperature measurement method - Google Patents
Multichannel infrared radiation measurement system and multi-wavelength real-time temperature measurement method Download PDFInfo
- Publication number
- CN117232661B CN117232661B CN202311522856.0A CN202311522856A CN117232661B CN 117232661 B CN117232661 B CN 117232661B CN 202311522856 A CN202311522856 A CN 202311522856A CN 117232661 B CN117232661 B CN 117232661B
- Authority
- CN
- China
- Prior art keywords
- target
- measured
- temperature
- representing
- channel
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 230000005855 radiation Effects 0.000 title claims abstract description 134
- 238000005259 measurement Methods 0.000 title claims abstract description 101
- 238000000034 method Methods 0.000 title claims abstract description 44
- 238000009529 body temperature measurement Methods 0.000 title claims abstract description 20
- 238000005457 optimization Methods 0.000 claims abstract description 11
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 10
- 238000002834 transmittance Methods 0.000 claims description 9
- 238000012545 processing Methods 0.000 claims description 5
- 230000031700 light absorption Effects 0.000 claims description 4
- 230000006870 function Effects 0.000 description 28
- 230000003287 optical effect Effects 0.000 description 5
- 238000004590 computer program Methods 0.000 description 3
- 238000000691 measurement method Methods 0.000 description 3
- 230000008878 coupling Effects 0.000 description 2
- 238000010168 coupling process Methods 0.000 description 2
- 238000005859 coupling reaction Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000010354 integration Effects 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 238000001228 spectrum Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 238000004861 thermometry Methods 0.000 description 2
- 230000002238 attenuated effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000005316 response function Methods 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 239000013077 target material Substances 0.000 description 1
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E30/00—Energy generation of nuclear origin
- Y02E30/30—Nuclear fission reactors
Landscapes
- Radiation Pyrometers (AREA)
Abstract
A multichannel infrared radiation measurement system and a multi-wavelength real-time temperature measurement method relate to the field of aerospace. The method comprises the following steps: determining a first calibration coefficient and a second calibration coefficient according to system calibration of the multichannel infrared radiation measurement system; obtaining a target image of a target to be measured to determine an output image gray value of the target to be measured, and constructing a multichannel infrared radiation measurement formula of the target to be measured; using a lambertian formula to express the radiance of the target to be measured; acquiring a target image of the infrared standard star to determine an output image gray value of the infrared standard star; constructing a multichannel infrared radiation measurement formula of the infrared standard star by using the gray value of the output image of the infrared standard star; determining a temperature measurement equation according to a multi-channel infrared radiation measurement formula of the infrared standard star, a multi-channel infrared radiation measurement formula of the target to be measured and the radiance of the target to be measured expressed by a lambert formula; and solving a temperature measurement equation by utilizing a multi-target optimization algorithm to determine the temperature of the target to be measured.
Description
Technical Field
The invention relates to the field of aerospace, in particular to a multichannel infrared radiation measurement system and a multi-wavelength real-time temperature measurement method.
Background
The temperature is one of the most basic physical quantities of the space target, and in the field of aerospace, the temperature characteristics and the change rule thereof are important information for deducing the properties of the space target and the on-orbit working state thereof, so that an important basis can be provided for target identification.
In the measurement of the spatial target characteristics, a radiation thermometry is a common method for measuring the current temperature, and the spectrum parameters according to the measurement mainly comprise full radiation thermometry, dual-band temperature measurement and multispectral temperature measurement. The traditional total radiation temperature measurement method is to measure the total radiation energy of a measured target in a full wavelength range within a certain wave band, and is limited by the spatial resolution of an optical measurement system, so that a target signal and background noise cannot be separated fundamentally, and the physical characteristics of the target cannot be reflected accurately. The dual-band temperature measurement and multispectral temperature measurement method is to calculate the temperature of a target on the premise that the emissivity of the target is approximately equal at adjacent wavelengths or the emissivity model is known. In practical applications, on the one hand, the emissivity of the target material is often not a fixed value, and on the other hand, the emissivity of the target is not known a priori. At present, a universal emissivity model is not found to be suitable for all materials, so that the temperature obtained by inversion by the method is not the real temperature of the target, the temperature measurement accuracy is greatly influenced by the assumed emissivity model, and the reliable evaluation of the ground infrared radiation measurement equipment on the effect of identifying the space target and the like is seriously influenced.
Disclosure of Invention
The invention mainly solves the technical problems that: a method for processing general data independent of a spectral emissivity model is provided.
According to a first aspect, in one embodiment, a method for measuring a multi-wavelength real-time temperature is provided, including:
determining a first calibration coefficient and a second calibration coefficient according to system calibration of the multichannel infrared radiation measurement system;
obtaining a target image of a target to be detected so as to determine an output image gray value of the target to be detected;
constructing a multichannel infrared radiation measurement formula of the target to be measured by utilizing the gray value of the output image of the target to be measured, the first calibration coefficient, the second calibration coefficient and the radiance of the target to be measured; the radiance of the target to be measured is expressed by using a lambertian formula;
acquiring a target image of an infrared standard star so as to determine an output image gray value of the infrared standard star;
constructing a multichannel infrared radiation measurement formula of the infrared standard star by utilizing the gray value of the output image of the infrared standard star, the first calibration coefficient, the second calibration coefficient and the radiance of the infrared standard star;
determining a measurement equation of temperature according to a multi-channel infrared radiation measurement formula of the infrared standard star, a multi-channel infrared radiation measurement formula of the target to be measured and the radiance of the target to be measured expressed by a lambertian formula;
and solving the temperature measurement equation by utilizing a multi-target optimization algorithm to determine the temperature of the target to be measured.
In one embodiment, the determining the first calibration factor and the second calibration factor according to the system calibration of the multi-channel infrared radiation measurement system includes:
acquiring target images of a set number of known targets and corresponding radiance of the known targets, and determining an output image gray value of the known targets according to the target images of the known targets;
a linear fit is performed using a least squares method to determine a first calibration coefficient and a second calibration coefficient for each channel in the multi-channel infrared radiation measurement system in the following formula:
；
wherein,an output image gray value representing a known target; />Representing the irradiance of a known target; />Representing the wavelength of the ith channel; t represents temperature; />A first calibration factor representing an ith channel; />Representing a second calibration factor for the ith channel.
In one embodiment, the constructing a multi-channel infrared radiation measurement formula of the target to be measured by using the gray value of the output image of the target to be measured, the first calibration coefficient, the second calibration coefficient and the radiance of the target to be measured includes:
the multichannel infrared radiation measurement formula of the target to be measured is as follows:
；
rewriting a multichannel infrared radiation measurement formula of the target to be measured as follows:
；
wherein,representing the gray value of an output image of the object to be measured; />A first calibration factor representing an ith channel; />The atmospheric transmittance of the target to be measured of the ith channel and the multichannel infrared radiation measuring system is represented;representing the emissivity of the target to be measured of the ith channel; />Representing the wavelength of the ith channel; t represents temperature;the radiation degree of the target to be detected with the temperature T of the ith channel is represented; />Indicating that the temperature of the ith channel is +.>Is a background environment of the radiation; />Indicating that the i-th channel target to be measured and the temperature areAtmospheric irradiance of the multi-channel infrared radiation measurement system; />Representing a second calibration factor for the ith channel.
In one embodiment, the method for expressing the radiance of the target to be measured by using a lambertian formula includes:
；
wherein,the radiation degree of the target to be detected with the temperature T of the ith channel is represented; />Representing the wavelength of the ith channel; t represents temperature; k represents the light absorption coefficient; c (C) 1 Representing a first radiation constant; c (C) 2 Representing a second radiation constant.
In one embodiment, the constructing a multi-channel infrared radiation measurement formula of the infrared calibration star by using the gray value of the output image of the infrared calibration star, the first calibration coefficient, the second calibration coefficient and the radiance of the infrared calibration star includes:
the multichannel infrared radiation measurement formula of the infrared standard star is as follows:
；
the multichannel infrared radiation measurement formula of the infrared standard star is rewritten as follows:
；
wherein,an output image gray value representing an infrared standard star; />A second calibration factor representing an ith channel; />A first calibration factor representing an ith channel; />Indicating that the i-th channel target to be measured and the temperature are +.>Atmospheric irradiance of the multi-channel infrared radiation measurement system; />Representing the wavelength of the ith channel; />A first calibration factor representing an ith channel; />Representing the test of the ith channelAtmospheric transmittance of the target and multi-channel infrared radiation measurement system; />Indicating that the temperature of the ith channel is +.>Is provided.
In one embodiment, the measurement equation for determining the temperature according to the multi-channel infrared radiation measurement formula of the infrared standard star, the multi-channel infrared radiation measurement formula of the target to be measured and the radiance of the target to be measured expressed by the lambertian formula includes:
dividing the rewritten multichannel infrared radiation measurement formula of the target to be measured by the rewritten multichannel infrared radiation measurement formula of the infrared standard star:
；
and (3) making:
；
and determining a measurement equation of temperature by combining the radiance of the target to be measured expressed by the lambertian formula:
；
wherein T is i Representing the temperature of the object to be measured of the ith channel;representing the gray value of an output image of the object to be measured; />An output image gray value representing an infrared standard star; />First label representing ith channelDetermining coefficients;the atmospheric transmittance of the target to be measured of the ith channel and the multichannel infrared radiation measuring system is represented;representing the emissivity of the target to be measured of the ith channel; />Representing the wavelength of the ith channel; t represents temperature;the radiation degree of the target to be detected with the temperature T of the ith channel is represented; />Indicating that the temperature of the ith channel is +.>Is a background environment of the radiation; />Indicating that the i-th channel target to be measured and the temperature areAtmospheric irradiance of the multi-channel infrared radiation measurement system; />A second calibration factor representing an ith channel;indicating that the temperature of the ith channel is +.>Is provided.
In one embodiment, the method for solving the temperature measurement equation by using the multi-objective optimization algorithm to determine the temperature of the object to be measured includes:
constructing an objective function of the temperature of the object to be measured;
constructing a nonlinear equation constraint condition and a linear inequality constraint condition of the temperature of the target to be measured;
and determining the temperature of the target to be detected according to the objective function, the nonlinear equation constraint condition and the linear inequality constraint condition.
In one embodiment, the objective function of the temperature of the object to be measured includes a first objective function and a second objective function;
the first objective function is:
；
wherein,、/>…/>a function value representing a first objective function, T 1 、T 2 、T 3 、T 4 …T n-1 、T n Representing the temperature of the object to be measured of each channel; n represents the total number of channels;
the second objective function is:
；
wherein F is s A function value representing a second objective function; t (T) i Representing the temperature of the object to be measured of the ith channel;representing the average temperature of the targets to be measured of all channels; n represents the total number of channels.
In one embodiment, the nonlinear equation constraint of the temperature of the target to be measured is:
；
wherein Z is i Representing a temperature residual;representing the temperature of the object to be measured of the ith channel; />Representing the average temperature of the targets to be measured of all channels;
the linear inequality constraint condition of the temperature of the target to be measured is as follows:
；
wherein,representing the emissivity of the target to be measured of the ith channel; />Representing the wavelength of the ith channel; representing the emissivity of the target to be measured of the (i+1) th channel; />Representing the wavelength of the (i+1) th channel; t represents temperature.
According to a second aspect, there is provided in one embodiment a multi-channel infrared radiation measurement system comprising:
the observation device is used for acquiring a target image of a target to be detected and a target image of an infrared standard star;
a data processing system that performs the method of measuring a multi-wavelength real-time temperature described in any of the embodiments above.
According to the multi-channel infrared radiation measurement system and the multi-wavelength real-time temperature measurement method, after the system calibration of the multi-channel infrared radiation measurement system is completed, the infrared standard star is used as a reference target, target images of a target to be measured and the infrared standard star are respectively obtained, a multi-channel infrared radiation measurement formula of the target to be measured and a multi-channel infrared radiation measurement formula of the infrared standard star are constructed, then a measurement equation of the temperature of the target to be measured is determined according to the multi-channel infrared radiation measurement formula of the target to be measured and the multi-channel infrared radiation measurement formula of the infrared standard star, and a multi-target optimization algorithm is adopted for solving. According to the method and the device, an emissivity model of the target to be measured is not required to be established, temperature inversion of the target to be measured can be achieved without any treatment on brightness, temperature and wavelength, the principle of true temperature inversion is simplified, inversion speed is greatly improved compared with that of a traditional secondary measurement method, inversion accuracy is basically the same, and therefore complexity of a system is greatly reduced.
Drawings
FIG. 1 is a flow chart of a method for measuring multi-wavelength real-time temperature according to one embodiment;
FIG. 2 is a flowchart of step S200 in a method for measuring a multi-wavelength real-time temperature according to an embodiment;
FIG. 3 is a schematic illustration of radiation received by an i-th channel of a multi-channel infrared radiation measurement system in a multi-wavelength real-time temperature measurement method according to an embodiment;
FIG. 4 is a flowchart of step S300 in a method for measuring a multi-wavelength real-time temperature according to an embodiment;
FIG. 5 is a flowchart of step S500 in a method for measuring a multi-wavelength real-time temperature according to an embodiment;
FIG. 6 is a schematic diagram of a multi-channel infrared radiation measurement system according to another embodiment.
Detailed Description
The invention will be described in further detail below with reference to the drawings by means of specific embodiments. Wherein like elements in different embodiments are numbered alike in association. In the following embodiments, numerous specific details are set forth in order to provide a better understanding of the present application. However, one skilled in the art will readily recognize that some of the features may be omitted, or replaced by other elements, materials, or methods in different situations. In some instances, some operations associated with the present application have not been shown or described in the specification to avoid obscuring the core portions of the present application, and may not be necessary for a person skilled in the art to describe in detail the relevant operations based on the description herein and the general knowledge of one skilled in the art.
Furthermore, the described features, operations, or characteristics of the description may be combined in any suitable manner in various embodiments. Also, various steps or acts in the method descriptions may be interchanged or modified in a manner apparent to those of ordinary skill in the art. Thus, the various orders in the description and drawings are for clarity of description of only certain embodiments, and are not meant to be required orders unless otherwise indicated.
The numbering of the components itself, e.g. "first", "second", etc., is used herein merely to distinguish between the described objects and does not have any sequential or technical meaning. The terms "coupled" and "connected," as used herein, are intended to encompass both direct and indirect coupling (coupling), unless otherwise indicated.
Aiming at the problem of temperature inversion of a space target with unknown emissivity model, the application provides a universal data processing method independent of a spectrum emissivity model based on a space target multichannel infrared radiation measuring device with n narrow-band wavelength radiation measuring channels, a multi-target equation is established, and the temperature real-time inversion of the target with unknown emissivity model is realized by adopting a multi-target optimization algorithm. The specific conception of the application is as follows: the space target multichannel infrared radiation measurement telescope system is divided into n narrow-band wavelength radiation measurement channels through a light splitting structure, after the system calibration is completed for each measurement channel, the radiance information of a target to be measured under the n narrow-band wavelengths is respectively obtained, and n channel equations are established. And taking a certain determined infrared standard star as a reference target, acquiring the radiance information of the infrared standard star under n narrow-band wavelengths, and establishing n channel equations based on the infrared standard star. The multi-objective function and the constraint condition are preset by combining the two equations, and the extremum of the multi-objective function is solved by adopting a multi-objective optimization algorithm, so that the temperature of the inversion target to be measured is regarded as the temperature, and the specific description is given below.
Referring to fig. 1, in one embodiment, a method for measuring a multi-wavelength real-time temperature is provided, which includes the following steps:
step S100: and determining a first calibration coefficient and a second calibration coefficient according to the system calibration of the multichannel infrared radiation measurement system.
In one embodiment, the multi-channel infrared radiation measurement system is an apparatus for measuring radiation characteristics of an object in the infrared band. The system generally includes a multi-channel infrared radiation source, an optical assembly, a detector, a data acquisition system, and the like. According to the response wavelength range of the multichannel infrared radiation measurement system detector, the whole corresponding wave band is subdivided into n narrow wave bands. The system calibration test is carried out through the measured target with known radiation intensity and temperature, and specifically comprises the following steps:
and acquiring target images of the known targets of the set number group and the corresponding radiance of the known targets, and determining the gray value of the output image of the known targets according to the target images of the known targets.
Obtaining a relation model between output quantity and input quantity of different wavelength channel systems under different integration time through calibration:
；
wherein,an output image gray value representing a known target; />Representing the irradiance of a known target; />Representing the wavelength of the ith channel; t represents temperature; />A first calibration factor representing an ith channel; />Representing a second calibration factor for the ith channel. It should be noted that the first calibration coefficient and the second calibration coefficient are different under different integration time conditions of different optical systems, different attenuation sheets and devices.
And obtaining target images of the known targets and the corresponding radiance of the known targets through experiments of a set number of groups, and performing linear fitting on the data by using a least square method to obtain a radiance response function and a curve, so as to obtain a first calibration coefficient and a second calibration coefficient of each channel of the multichannel infrared radiation measurement system.
Step S200: and constructing a multichannel infrared radiation measurement formula of the target to be measured.
Referring to fig. 2, in executing step S200 to construct a multi-channel infrared radiation measurement formula of a target to be measured, an embodiment further includes the following steps:
step S210: and obtaining a target image of the target to be detected so as to determine the gray value of the output image of the target to be detected.
In one embodiment, the multi-channel infrared radiation measurement system can directly acquire the target image of the target to be measured, and the gray value of the output image of the target to be measured can be directly determined according to the target image of the target to be measured.
Step S220: and constructing a multichannel infrared radiation measurement formula of the target to be measured by using the gray value of the output image of the target to be measured, the first calibration coefficient, the second calibration coefficient and the radiance of the target to be measured.
In one embodiment, when radiation measurement is performed on a target to be measured, a radiation diagram received by an ith channel of the multi-channel infrared radiation measurement system is shown in fig. 3, and is mainly affected by three factors: (1) Atmospheric attenuated radiance of the target itselfThe method comprises the steps of carrying out a first treatment on the surface of the (2) The radiance of the object to be measured to the surrounding environment>The method comprises the steps of carrying out a first treatment on the surface of the (3) Atmospheric radiance of object to be measured and multichannel infrared radiation measurement system>. The atmosphere transmission parameters in each sub-band channel are obtained through a multichannel infrared radiation atmosphere correction system: atmospheric transmittance->And atmospheric emittance->The calculation of (2) and (3) can be performed directly or indirectly.
In one embodiment, for a multi-channel infrared radiation measurement system with n channels (n is generally an even number, and n channels are n narrow bands), the multi-channel infrared radiation measurement formula of the object to be measured of the ith channel is:
；
rewriting a multichannel infrared radiation measurement formula of a target to be measured as follows:
；
wherein,representing the gray value of an output image of the object to be measured; />A first calibration factor representing an ith channel;the atmospheric transmittance of the target to be measured of the ith channel and the multichannel infrared radiation measuring system is represented;represent the firstEmissivity of the target to be measured of i channels; />Representing the wavelength of the ith channel; t represents temperature;the radiation degree of the target to be detected with the temperature T of the ith channel is represented; />Indicating that the temperature of the ith channel is +.>The background environment radiance of (2) is generally obtained by integrating the environment radiation as a blackbody of the environment temperature in combination with a blackbody Planck formula in a response band of the detector; />Indicating that the i-th channel target to be measured and the temperature are +.>Atmospheric irradiance of the multi-channel infrared radiation measurement system; />Representing a second calibration factor for the ith channel.
Step S230: and (5) using a lambertian formula to express the radiance of the target to be measured.
In one embodiment, the radiance of the target to be measured expressed by using the lambertian body as a reference is shown as follows:
；
wherein,the radiation degree of the target to be detected with the temperature T of the ith channel is represented; />Representing the wavelength of the ith channel; t represents temperature; k represents the light absorption coefficient; c (C) 1 Representing a first radiation constant; c (C) 2 Representing a second radiation constant.
Step S300: and constructing a multi-channel infrared radiation measurement formula of the infrared standard star.
Referring to fig. 4, in executing step S300 to construct a multi-channel infrared radiation measurement formula of an infrared standard star, an embodiment further includes the following steps:
step S310: and acquiring a target image of the infrared standard star to determine an output image gray value of the infrared standard star.
In one embodiment, the multi-channel infrared radiation measurement system can directly acquire the target image of the infrared standard star, and the gray value of the output image of the target to be measured can be directly determined according to the target image of the infrared standard star.
Step S320: and constructing a multichannel infrared radiation measurement formula of the infrared standard star by using the gray value of the output image of the infrared standard star, the first calibration coefficient, the second calibration coefficient and the radiance of the infrared standard star.
In one embodiment, for a multi-channel infrared radiation measurement system with n channels (n is generally even, and n channels are n narrow bands), the multi-channel infrared radiation measurement formula of the infrared standard star of the ith channel is
；
The multichannel infrared radiation measurement formula of the infrared standard star is rewritten as follows:
；
wherein,an output image gray value representing an infrared standard star; />Represents the ith passA second calibration factor for the track; />A first calibration factor representing an ith channel; />Indicating that the i-th channel target to be measured and the temperature are +.>Atmospheric irradiance of the multi-channel infrared radiation measurement system; />Representing the wavelength of the ith channel; />A first calibration factor representing an ith channel; />The atmospheric transmittance of the target to be measured of the ith channel and the multichannel infrared radiation measuring system is represented; />Indicating that the temperature of the ith channel is +.>The radiance of the infrared standard star of the (c) can be known by inquiring the infrared standard star library.
Step S400: and determining a measurement equation of the temperature according to a multi-channel infrared radiation measurement formula of the infrared standard star, a multi-channel infrared radiation measurement formula of the target to be measured and the radiance of the target to be measured expressed by a lambert formula.
In one embodiment, the multi-channel infrared radiation measurement formula of the target to be measured after being rewritten in step S220 is divided by the multi-channel infrared radiation measurement formula of the infrared standard star after being rewritten in step S230, so as to obtain the following formula:
；
and (3) making:
；
and combining the radiance of the target to be measured expressed by using the lambertian formula in step S230, thereby obtaining a measurement equation of the radiance determination temperature of the target to be measured:
；
in this formula, n measurement equations containing n+1 unknowns, the n+1 unknowns including n, can be obtained for a n-channel multi-channel infrared radiation measurement systemAnd 1 temperature T.
Wherein T is i Representing the temperature of the object to be measured of the ith channel;representing the gray value of an output image of the object to be measured; />An output image gray value representing an infrared standard star; />A first calibration factor representing an ith channel;the atmospheric transmittance of the target to be measured of the ith channel and the multichannel infrared radiation measuring system is represented;representing the emissivity of the target to be measured of the ith channel; />Representing the wavelength of the ith channel; t represents temperature;the radiation degree of the target to be detected with the temperature T of the ith channel is represented; />Indicating that the temperature of the ith channel is +.>Is a background environment of the radiation; />Indicating that the i-th channel target to be measured and the temperature areAtmospheric irradiance of the multi-channel infrared radiation measurement system; />A second calibration factor representing an ith channel;indicating that the temperature of the ith channel is +.>Is provided.
Step S500: and solving a temperature measurement equation by utilizing a multi-target optimization algorithm to determine the temperature of the target to be measured.
Referring to fig. 5, in executing step S500, when solving the temperature measurement equation by using the multi-objective optimization algorithm to determine the temperature of the object to be measured, an embodiment further includes the following steps:
step S510: and constructing an objective function of the temperature of the object to be measured.
In one embodiment, a first objective function is constructed from the minimum temperature difference calculated for two adjacent channels:
；
wherein,、/>…/>a function value representing a first objective function; t (T) 1 、T 2 、T 3 、T 4 …T n-1 、T n Representing the temperature of the object to be measured of each channel; n represents the total number of channels.
In one embodiment, the standard deviation is expressed by using a Bessel formula with the goal of minimizing the standard deviation of the temperature obtained for the different channels, and the square of the standard deviation is minimized as a second objective function:
；
；
wherein F is s A function value representing a second objective function; t (T) i Representing the temperature of the object to be measured of the ith channel;representing the average temperature of the targets to be measured of all channels; n represents the total number of channels.
Step S520: and constructing a nonlinear equation constraint condition and a linear inequality constraint condition of the temperature of the target to be measured.
In one embodiment, the residual values of the individual channel measurements are obtained by subtracting the temperature average value from each channel equation to construct a nonlinear multivariate inequality constraint:
；
wherein Z is i Representing a temperature residual; t (T) i Represent the firsti temperatures of targets to be measured of the channels;representing the average temperature of the object to be measured for all channels. It should be noted that, ideally, the temperature residual of each channel is 0.
In one embodiment, a linear multivariate inequality constraint is constructed from the basic theory of radiation:
；
wherein,representing the emissivity of the target to be measured of the ith channel; />Representing the wavelength of the ith channel;representing the emissivity of the target to be measured of the (i+1) th channel; />Representing the wavelength of the (i+1) th channel; t represents temperature.
Step S530: and determining the temperature of the target to be measured according to the objective function, the nonlinear equation constraint condition and the linear inequality constraint condition.
In one embodiment, the computer system uses a multi-objective optimization method based on objective functions and constraints, under the constraint of nonlinear multivariable constraints and inequality constraints, when the temperature T is calculated i When the approximate solution or weak effective solution of the objective function is adopted, the temperature T is calculated i Consider the temperature of the object to be measured of the i-th channel.
According to the multi-wavelength real-time temperature measurement method, after system calibration is completed, the infrared standard star is used as a reference target, the target to be measured and the standard reference star are measured respectively, a multi-target extremum optimization algorithm is adopted, an emissivity model of the target to be measured is not required to be established, brightness, temperature and wavelength are not required to be processed, temperature inversion of the target to be measured can be achieved, the principle of true temperature inversion is simplified, inversion speed is greatly improved compared with that of a traditional secondary measurement method, inversion accuracy is basically the same as that of the traditional secondary measurement method, and complexity of the system is greatly reduced.
Referring to fig. 6, in another embodiment, the present application further provides a multi-channel infrared radiation measurement system 10, including an observation device 11 and a data processing system 12, where the observation device 11 is configured to obtain a target image of a target to be measured and a target image of an infrared standard star; the data processing system 12 performs the method for measuring the multi-wavelength real-time temperature described in any of the above embodiments, and since the method for measuring the multi-wavelength real-time temperature is described in the above embodiments, it is not described in detail here.
Those skilled in the art will appreciate that all or part of the functions of the various methods in the above embodiments may be implemented by hardware, or may be implemented by a computer program. When all or part of the functions in the above embodiments are implemented by means of a computer program, the program may be stored in a computer readable storage medium, and the storage medium may include: read-only memory, random access memory, magnetic disk, optical disk, hard disk, etc., and the program is executed by a computer to realize the above-mentioned functions. For example, the program is stored in the memory of the device, and when the program in the memory is executed by the processor, all or part of the functions described above can be realized. In addition, when all or part of the functions in the above embodiments are implemented by means of a computer program, the program may be stored in a storage medium such as a server, another computer, a magnetic disk, an optical disk, a flash disk, or a removable hard disk, and the program in the above embodiments may be implemented by downloading or copying the program into a memory of a local device or updating a version of a system of the local device, and when the program in the memory is executed by a processor.
The foregoing description of the invention has been presented for purposes of illustration and description, and is not intended to be limiting. Several simple deductions, modifications or substitutions may also be made by a person skilled in the art to which the invention pertains, based on the idea of the invention.
Claims (8)
1. A method for measuring temperature in real time at multiple wavelengths, comprising:
determining a first calibration coefficient and a second calibration coefficient according to system calibration of the multichannel infrared radiation measurement system;
obtaining a target image of a target to be detected so as to determine an output image gray value of the target to be detected;
constructing a multichannel infrared radiation measurement formula of the target to be measured by utilizing the gray value of the output image of the target to be measured, the first calibration coefficient, the second calibration coefficient and the radiance of the target to be measured; the radiance of the target to be measured is expressed by using a lambertian formula;
acquiring a target image of an infrared standard star so as to determine an output image gray value of the infrared standard star;
constructing a multichannel infrared radiation measurement formula of the infrared standard star by utilizing the gray value of the output image of the infrared standard star, the first calibration coefficient, the second calibration coefficient and the radiance of the infrared standard star;
dividing the rewritten multi-channel infrared radiation measurement formula of the target to be measured by the rewritten multi-channel infrared radiation measurement formula of the infrared standard star:
；
and (3) making:
；
and determining a measurement equation of temperature by combining the radiance of the target to be measured expressed by the lambertian formula:
；
wherein the method comprises the steps of，T i Representing the temperature of the object to be measured of the ith channel;representing the gray value of an output image of the object to be measured;an output image gray value representing an infrared standard star; />A first calibration factor representing an ith channel;the atmospheric transmittance of the target to be measured of the ith channel and the multichannel infrared radiation measuring system is represented;representing the emissivity of the target to be measured of the ith channel; />Representing the wavelength of the ith channel; t represents temperature;the radiation degree of the target to be detected with the temperature T of the ith channel is represented; />Indicating that the temperature of the ith channel is +.>Is a background environment of the radiation; />Indicating that the i-th channel target to be measured and the temperature areAtmosphere radiation of multi-channel infrared radiation measuring systemA emittance; />A second calibration factor representing an ith channel;indicating that the temperature of the ith channel is +.>The emittance of the infrared standard star of (2); k represents the light absorption coefficient; c (C) 1 Representing a first radiation constant; c (C) 2 Representing a second radiation constant;
solving the temperature measurement equation by using a multi-target optimization algorithm to determine the temperature of the target to be measured, wherein the method comprises the following steps:
constructing an objective function of the temperature of the object to be measured;
constructing a nonlinear equation constraint condition and a linear inequality constraint condition of the temperature of the target to be measured;
and determining the temperature of the target to be detected according to the objective function, the nonlinear equation constraint condition and the linear inequality constraint condition.
2. The method for measuring the temperature in real time at multiple wavelengths according to claim 1, wherein determining the first calibration factor and the second calibration factor according to the system calibration of the multi-channel infrared radiation measurement system comprises:
acquiring target images of a set number of known targets and corresponding radiance of the known targets, and determining an output image gray value of the known targets according to the target images of the known targets;
a linear fit is performed using a least squares method to determine a first calibration coefficient and a second calibration coefficient for each channel in the multi-channel infrared radiation measurement system in the following formula:
；
wherein,an output image gray value representing a known target; />Representing the irradiance of a known target;representing the wavelength of the ith channel; t represents temperature; />A first calibration factor representing an ith channel; />Representing a second calibration factor for the ith channel.
3. The method for measuring multi-wavelength real-time temperature according to claim 1, wherein constructing a multi-channel infrared radiation measurement formula of the target to be measured by using the gray value of the output image of the target to be measured, the first calibration coefficient, the second calibration coefficient and the radiance of the target to be measured, comprises:
the multichannel infrared radiation measurement formula of the target to be measured is as follows:
；
rewriting a multichannel infrared radiation measurement formula of the target to be measured as follows:
；
wherein,representing the gray value of an output image of the object to be measured; />A first calibration factor representing an ith channel;the atmospheric transmittance of the target to be measured of the ith channel and the multichannel infrared radiation measuring system is represented;representing the emissivity of the target to be measured of the ith channel; />Representing the wavelength of the ith channel; t represents temperature;the radiation degree of the target to be detected with the temperature T of the ith channel is represented; />Indicating that the temperature of the ith channel is +.>Is a background environment of the radiation; />Indicating that the i-th channel target to be measured and the temperature areAtmospheric irradiance of the multi-channel infrared radiation measurement system; />Representing a second calibration factor for the ith channel.
4. The method for measuring a multi-wavelength real-time temperature according to claim 3, wherein the expressing the radiance of the object to be measured using a lambertian formula comprises:
；
wherein,the radiation degree of the target to be detected with the temperature T of the ith channel is represented; />Representing the wavelength of the ith channel; t represents temperature; k represents the light absorption coefficient; c (C) 1 Representing a first radiation constant; c (C) 2 Representing a second radiation constant.
5. The method for measuring the temperature in real time with multiple wavelengths according to claim 4, wherein constructing a multi-channel infrared radiation measurement formula of the infrared calibration star by using the gray value of the output image of the infrared calibration star, the first calibration coefficient, the second calibration coefficient and the radiance of the infrared calibration star comprises:
the multichannel infrared radiation measurement formula of the infrared standard star is as follows:
；
the multichannel infrared radiation measurement formula of the infrared standard star is rewritten as follows:
；
wherein,an output image gray value representing an infrared standard star; />A second calibration factor representing an ith channel; />A first calibration factor representing an ith channel; />Indicating that the i-th channel target to be measured and the temperature are +.>Atmospheric irradiance of the multi-channel infrared radiation measurement system; />Representing the wavelength of the ith channel;the atmospheric transmittance of the target to be measured of the ith channel and the multichannel infrared radiation measuring system is represented;indicating that the temperature of the ith channel is +.>Is provided.
6. The method for measuring the temperature in real time at multiple wavelengths according to claim 1, wherein the objective function of the temperature of the object to be measured comprises a first objective function and a second objective function;
the first objective function is:
；
wherein,、/>…/>a function value representing a first objective function, T 1 、T 2 、T 3 、T 4 …T n-1 、T n Representing the temperature of the object to be measured of each channel; n represents the total number of channels;
the second objective function is:
；
wherein F is s A function value representing a second objective function; t (T) i Representing the temperature of the object to be measured of the ith channel;representing the average temperature of the targets to be measured of all channels; n represents the total number of channels.
7. The method for measuring the temperature in real time at multiple wavelengths according to claim 1, wherein the constraint condition of the nonlinear equation of the temperature of the object to be measured is:
；
wherein Z is i Representing a temperature residual; t (T) i Representing the temperature of the object to be measured of the ith channel;representing the average temperature of the targets to be measured of all channels;
the linear inequality constraint condition of the temperature of the target to be measured is as follows:
；
wherein,representing the emissivity of the target to be measured of the ith channel; />Representing the wavelength of the ith channel;representing the emissivity of the target to be measured of the (i+1) th channel; />Representing the wavelength of the (i+1) th channel; t represents temperature.
8. A multi-channel infrared radiation measurement system, comprising:
the observation device is used for acquiring a target image of a target to be detected and a target image of an infrared standard star;
a data processing system performing the method of measuring a multi-wavelength real-time temperature as claimed in any one of claims 1 to 7.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311522856.0A CN117232661B (en) | 2023-11-16 | 2023-11-16 | Multichannel infrared radiation measurement system and multi-wavelength real-time temperature measurement method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311522856.0A CN117232661B (en) | 2023-11-16 | 2023-11-16 | Multichannel infrared radiation measurement system and multi-wavelength real-time temperature measurement method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN117232661A CN117232661A (en) | 2023-12-15 |
| CN117232661B true CN117232661B (en) | 2024-02-23 |
Family
ID=89083001
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202311522856.0A Active CN117232661B (en) | 2023-11-16 | 2023-11-16 | Multichannel infrared radiation measurement system and multi-wavelength real-time temperature measurement method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN117232661B (en) |
Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4881823A (en) * | 1988-03-29 | 1989-11-21 | Purdue Research Foundation | Radiation thermometry |
| CN102080990A (en) * | 2010-12-01 | 2011-06-01 | 中国航空工业集团公司北京长城计量测试技术研究所 | Four-waveband high temperature measuring device and method |
| CN103674237A (en) * | 2012-09-25 | 2014-03-26 | 中国航天科工集团第二研究院二〇七所 | A method for calibrating cross radiance of an infrared fixed star and a sky background |
| CN107741276A (en) * | 2017-09-05 | 2018-02-27 | 中国航空工业集团公司洛阳电光设备研究所 | A kind of method for improving technics of temperature measurement precision with infrared thermal imager |
| CN111721423A (en) * | 2020-06-19 | 2020-09-29 | 中国人民解放军63660部队 | Three-band target surface temperature inversion method |
| CN112798110A (en) * | 2020-12-29 | 2021-05-14 | 杭州晨安科技股份有限公司 | Calibration fitting-based temperature detection method for infrared thermal imaging equipment |
| WO2022079680A1 (en) * | 2020-10-16 | 2022-04-21 | Arcelormittal | Method for estimating the temperature and the oxide thickness of a steel strip |
| CN114719996A (en) * | 2022-04-12 | 2022-07-08 | 中国科学院云南天文台 | High-precision spectral band radiance measuring system and method |
| CN115524012A (en) * | 2021-06-24 | 2022-12-27 | 哈尔滨工业大学 | Multispectral radiation temperature measuring device based on thermal infrared imager and measuring method thereof |
| CN116625519A (en) * | 2023-05-24 | 2023-08-22 | 河南师范大学 | Multispectral radiation temperature measurement inversion calculation method based on ISMA algorithm |
-
2023
- 2023-11-16 CN CN202311522856.0A patent/CN117232661B/en active Active
Patent Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4881823A (en) * | 1988-03-29 | 1989-11-21 | Purdue Research Foundation | Radiation thermometry |
| CN102080990A (en) * | 2010-12-01 | 2011-06-01 | 中国航空工业集团公司北京长城计量测试技术研究所 | Four-waveband high temperature measuring device and method |
| CN103674237A (en) * | 2012-09-25 | 2014-03-26 | 中国航天科工集团第二研究院二〇七所 | A method for calibrating cross radiance of an infrared fixed star and a sky background |
| CN107741276A (en) * | 2017-09-05 | 2018-02-27 | 中国航空工业集团公司洛阳电光设备研究所 | A kind of method for improving technics of temperature measurement precision with infrared thermal imager |
| CN111721423A (en) * | 2020-06-19 | 2020-09-29 | 中国人民解放军63660部队 | Three-band target surface temperature inversion method |
| WO2022079680A1 (en) * | 2020-10-16 | 2022-04-21 | Arcelormittal | Method for estimating the temperature and the oxide thickness of a steel strip |
| CN112798110A (en) * | 2020-12-29 | 2021-05-14 | 杭州晨安科技股份有限公司 | Calibration fitting-based temperature detection method for infrared thermal imaging equipment |
| CN115524012A (en) * | 2021-06-24 | 2022-12-27 | 哈尔滨工业大学 | Multispectral radiation temperature measuring device based on thermal infrared imager and measuring method thereof |
| CN114719996A (en) * | 2022-04-12 | 2022-07-08 | 中国科学院云南天文台 | High-precision spectral band radiance measuring system and method |
| CN116625519A (en) * | 2023-05-24 | 2023-08-22 | 河南师范大学 | Multispectral radiation temperature measurement inversion calculation method based on ISMA algorithm |
Also Published As
| Publication number | Publication date |
|---|---|
| CN117232661A (en) | 2023-12-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Marinetti et al. | Emissivity estimation for accurate quantitative thermography | |
| Kantzas et al. | Protocols for UV camera volcanic SO2 measurements | |
| US10337986B2 (en) | Method for optimisation of detection wavelengths for multi-gas detection | |
| Pereira et al. | A comparative study of calibration transfer methods for determination of gasoline quality parameters in three different near infrared spectrometers | |
| US8101905B2 (en) | Infrared sensor calibration system and method | |
| CN103913237B (en) | The three accurate temp measuring methods of wave band infra-red radiation | |
| Madura et al. | Automatic compensation of emissivity in three-wavelength pyrometers | |
| Zibordi et al. | In situ optical radiometry in the visible and near infrared | |
| Brouckaert et al. | Calibration transfer of a Raman spectroscopic quantification method for the assessment of liquid detergent compositions from at-line laboratory to in-line industrial scale | |
| CN108398191A (en) | High-precision colour temperature measurement method and system based on fiber spectrometer | |
| JP2021067611A5 (en) | ||
| Chrzanowski | Problem of determination of effective emissivity of some materials in MIR range | |
| CN117232661B (en) | Multichannel infrared radiation measurement system and multi-wavelength real-time temperature measurement method | |
| Huang et al. | A new wavelength selection criterion for two-color pyrometer interfered with participating media | |
| CN110132416B (en) | Broadband remote sensor on-orbit spectrum calibration method and device | |
| Barducci et al. | Maximum entropy temperature–emissivity separation in the TIR spectral range using the MaxEnTES algorithm | |
| Pinto | Uncertainty evaluation for in-flight radiometric calibration of Earth observation sensors | |
| Pohl et al. | Absolute calibration of the spectral responsivity of thermal detectors in the near-infrared (NIR) and mid-infrared (MIR) regions by using blackbody radiation | |
| Pereira et al. | Metrology of solar spectral irradiance at the top of the atmosphere in the near infrared measured at Mauna Loa Observatory: the PYR-ILIOS campaign | |
| Parr et al. | The use of filtered radiometers for radiance measurement | |
| Kuusk et al. | Estimation of variable atmospheric parameters for the atmospheric correction of satellite images | |
| Woolliams et al. | Intermediate uncertainty analysis for Earth Observation (Instrument Calibration Module): Training course textbook | |
| Bobrovnikov et al. | Comparison of signal processing methods in remote temperature measurements by pure rotational Raman spectra | |
| Krutikov et al. | The traceability of modern pyrometers to the primary standard of the unit of temperature and the classification of methods of pyrometry | |
| Wei et al. | Thermal infrared radiometer calibration and experimental measurements |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |