Multispectral High Temperature Thermography

The paper considers the issues of creating high-temperature digital thermographs based on RGB photodetector arrays. It has been shown that increasing the reliability of temperature measurement of bodies with unknown spectral coefficient of thermal radiation can be ensured by optimal selection of the used spectral range and registration of the observed thermal radiation fields in three spectral ranges. The registration of thermal radiation in four or more spectral ranges was found to be inefficient due to the increasing error in temperature determination. This paper presents a method for forming three overlapping spectral regions in the NIR spectral range, which is based on the use of an external spectral filter and a combination of the spectral characteristics of an RGB photodetector array. It is shown that it is necessary to ensure the stability of the solution of the system of three nonlinear equations with respect to the influence of noise. For this purpose, the use of a priori information about the slope factor of the spectral dependence of the thermal radiation coefficient in the selected spectral range for the controlled bodies is proposed. The theoretical results are confirmed by examples of their application in a thermograph based on an array of CMOS RGB photodetectors.


Introduction
Today's product quality requirements impose a mandatory control of the conditions under which manufacturing processes are carried out. Process temperature control is one of the main requirements of modern standards. Temperature is a quantitative parameter of the internal energy of bodies, so recording and analysis of temperature fields gives information on the state of objects and the course of various processes [1,2]. The non-contact methods based on the registration of thermal radiation emitted by the surface of the heated product seem to be the most convenient ones for controlling the temperature of complex thermal processes. In many cases, it is effective to use an infrared imaging technique that enables the presentation of a partial image of the radiation temperature field T r (x,y) of the supervised body on the monitor screen. Extensive information on its capabilities and application can be found on the websites of leading manufacturers of thermal imaging and thermographic equipment. However, the issues concerning the physical basis of temperature measurements, the formation of differential radiative heat transfer fluxes and the corresponding digital signals, the methods of temperature calculation with the uncertainty of the thermal emission coefficient, as well as the sources of uncertainty in temperature measurements are considered official information that is not widely disseminated. For users of thermographic equipment, such absence causes difficulties in the interpretation of measurement results in complex cases [3], e.g., with non-uniform heating, dynamic change Only low-cost thermographic techniques with stable methods of digital processing of registered thermal fields can have practical application in the control of high-temperature thermal processes. It can have two designs: (1) monoblock, when the thermograph chamber together with processing and indication unit are assembled in one inseparable case, and (2) combined, consisting of two units: the thermograph chamber and the processing, indication, recording, and storage device [14]. Communication between them is performed through a wired (USB, GIG-E) or wireless (Wi-Fi) interface. A serious advantage of the mobile two-unit design is the possibility of using modern small-size computing equipment as a processing device-laptops and tablets, which are convenient for recording, storing, and subsequent analysis of recorded thermograms and video Thermal radiation of metals is generated by photons resulting from quantum processes in a thin surface layer of about 4 nm thick. Therefore, when the metal is heated to high temperatures, its surface roughness degree, its oxidation by atmospheric oxygen, presence of impurities, etc. can lead to a substantial change of spectral dependence of thermal radiation coefficient ε (λ) in comparison with the tabular values given in reference books. In addition, other factors related to the observation conditions of the heated surface and the parameters of the photodetector path of thermographic technology will also contribute to error in measuring temperature [8,9].
In large enterprises, when controlling processes with stably repeating parameters (for example, rolling mills), the deviation of the partial radiation temperature T r from the specified values is usually used, since the parameters of thermal processes are stable. In the case of new complex thermal processes, however, it is desirable to determine the real values of T. Therefore, when the ε eff , is uncertain, some authors suggest a more complex thermographic technique [10], using the measurement of thermal radiation in three [11], four [12], or even eight spectral bands [13].
Only low-cost thermographic techniques with stable methods of digital processing of registered thermal fields can have practical application in the control of high-temperature thermal processes. It can have two designs: (1) monoblock, when the thermograph chamber together with processing and indication unit are assembled in one inseparable case, and (2) combined, consisting of two units: the thermograph chamber and the processing, indication, recording, and storage device [14]. Communication between them is performed through a wired (USB, GIG-E) or wireless (Wi-Fi) interface. A serious advantage of the mobile two-unit design is the possibility of using modern small-size computing equipment as a processing device-laptops and tablets, which are convenient for recording, storing, and subsequent analysis of recorded thermograms and video recordings of complex dynamic processes. Thus, non-standard algorithms of processing the received thermal images optimized for a concrete task of temperature field control can be used.
In contrast to the thermographic technique of the LWIR range of the spectrum in thermography of the visible and NIR ranges, the spectral density of the brightness of the recorded thermal radiation L (x,y) is changed by several orders of magnitude [7][8][9][10][11]. Adjusting the illumination of the matrix of photodetectors by changing the f-number of the lens of the thermograph camera is unacceptable, since the solid angle of incidence of radiation on the light-sensitive elements changes, and there are deviations of the recorded signals from the calibration dependence. Therefore, it is necessary to introduce algorithms for the automatic selection of the exposure time of the frame τ ex , which can vary over a wide range, into the software of the photoreceiving path. At high temperatures and, accordingly, short exposure times, a noticeable effect of defects in the photodetector matrix used in the thermograph is manifested, which must be excluded. This requires the use of appropriate methods for carrying out calibration work, which makes it possible to take into account the peculiarities of the photodetector path and the noise properties of the used photodetector matrix.
When analyzing possible ways to solve the problems noted above, it is necessary to be guided by the basic principles of reducing the uncertainty of temperature measurement by thermographic technology in the high-temperature range. The following should be provided [15]: • the ability to measure the true temperature T and estimate the effective emissivity ε eff ; • minimization of the uncertainty of the temperature measurement results achieved through the optimal choice of spectral sections for registration of thermal radiation and the elimination of the influence of the deviation of the calibration curves from the calculated dependences; • determination of the maximum body temperature T max and its dependence on time, as well as the possibility of video recording of the temperature field and its subsequent frame-by-frame viewing, which is necessary when monitoring complex heat engineering processes; • invariance of the results of determining the maximum temperature T max to changes in the size of the image of the monitored bodies; • invariance of the measured values of T max to nonstationarity of the noise dispersion of the used photodetector matrix, i.e., the dependence of its noise on the value of the incident thermal radiation flux.
To introduce these principles into the design and algorithms of software for hightemperature thermographic technology, it is necessary to carefully analyze the mechanisms of the influence of several factors on the uncertainty of the measured values of T [16] and the existing capabilities of photodetector matrices, optoelectronic technologies and algorithms for digital processing of recorded heat fluxes.
Therefore, the authors have tried to cover the principles of multispectral high-temperature thermographs functioning, using all available (or existing) arsenal of modern knowledge in character analysis and modelling of optoelectronic systems and digital technology.

Design and Quality Indicators of Thermographic Equipment
An idea of the composition of modern thermographs based on uncooled photodetector matrices is given by the diagram shown in Figure 2. The absence of moving mechanical elements, miniaturization of electronic components up to the inclusion of an ADC in the single-chip matrix of photodetectors, and mass production of compact lenses with high resolution makes it possible to create small-sized digital thermographs with the ability to quickly record the received thermograms and their primary analysis. allows ensuring uniformity of the thermograph sensitivity over the field of view; a matrix of photodetectors with a unit for reading signals, amplifying them and converting them into digital form; a display device or monitor; a control unit for the operation of the thermograph with the necessary controls; a block for processing the received images and displaying thermograms on the monitor screen; an interface for communication with external devices and a storage device for received data, which are usually stored in the form of files of a certain structure. Since the assessment of the quality of thermographs based on the visually observed image is subjective, in the process of development of thermal imaging technology, a number of standard indicators (parameters and characteristics) have been developed, the measurement of which is carried out in the laboratory conditions [17][18][19]. They allow experts to predict how the camera will perform under real-world conditions. Standard indicators, a detailed description of which is given in [18,19], are usually divided into 8 groups. Of these, high-temperature thermography is usually used: • accuracy parameters the root-mean-square value of the error with which the temperature can be measured; • instantaneous field of view; • frame formation time or frequency of their repetition; • resolution.

Formation of a Digital Image of the Temperature Field
Thermographic equipment does not allow directly determining the distance to the object and therefore operates with signals that depend on two variables x, y, i.e., only with the dependence of the illumination E (x, y) of the object image and the background in the focal plane of the receiving lens. For incoherent radiation, the illumination E (x, y) in the image space is proportional to the brightness L (x1, y1) of conjugate points in the object space. When calculating the spectral illumination density E (λ) of the The currently available set of commercially available photodetector matrices of various spectral ranges makes it possible to create thermographic equipment with optimal parameters for measuring low, medium, and high temperatures. They usually include: a lens with a multilayer antireflection coating; an optical filter that forms the required spectral range and can be applied to the input window of the photodetector matrix; a built-in controllable shutter (in MWIR and LWIR ranges) to block the field of view in order to quickly correct the sensitivity of the thermograph elements, which allows ensuring uniformity of the thermograph sensitivity over the field of view; a matrix of photodetectors with a unit for reading signals, amplifying them and converting them into digital form; a display device or monitor; a control unit for the operation of the thermograph with the necessary controls; a block for processing the received images and displaying thermograms on the monitor screen; an interface for communication with external devices and a storage device for received data, which are usually stored in the form of files of a certain structure.
Since the assessment of the quality of thermographs based on the visually observed image is subjective, in the process of development of thermal imaging technology, a number of standard indicators (parameters and characteristics) have been developed, the measurement of which is carried out in the laboratory conditions [17][18][19]. They allow experts to predict how the camera will perform under real-world conditions. Standard indicators, a detailed description of which is given in [18,19], are usually divided into 8 groups. Of these, high-temperature thermography is usually used: • accuracy parameters the root-mean-square value of the error with which the temperature can be measured; • instantaneous field of view; • frame formation time or frequency of their repetition; • resolution.

Formation of a Digital Image of the Temperature Field
Thermographic equipment does not allow directly determining the distance to the object and therefore operates with signals that depend on two variables x, y, i.e., only with the dependence of the illumination E (x, y) of the object image and the background in the focal plane of the receiving lens. For incoherent radiation, the illumination E (x, y) in the image space is proportional to the brightness L (x 1 , y 1 ) of conjugate points in the object space. When calculating the spectral illumination density E (λ) of the photosensitive element of the matrix with the area s ae , the spectral radiation flux density at the entrance pupil of the objective is usually determined where s so ae is the area of the projection of this matrix element onto the surface of the observed object in the space of objects, Ω = πD 2 /(4z 2 ) is the solid angle in which the radiation of the body from this area falls on the lens, f and D are the focal length and diameter of the entrance pupil of the lens, accordingly, and z is the distance from the surface of the observed object to the main optical plane of the lens. For Lambert emitters, their brightness L is proportional to their surface luminosity M, i.e., they are related by the expression L = M/π, which leads to a simple expression for the spectral illumination density of the photodetector matrix where the K coefficient takes into account the aperture of the lens and the distance z to the object and is a dimensionless quantity, which can be called the coefficient of connection between the illumination of the thermograph matrix and the luminosity of the object at a unit transmittance of the lens τ l = 1. Then, the spectral flux entering the photosensitive element of the photodetector matrix through the lens is described by the expression Averaged spectral slope of current conversion of photosensitive elements of photonic detector arrays takes into account the external quantum yield η ex (λ) of the photosensitive element, the power of the difference flux of radiation heat transfer ∆Φ (λ) carried by photons, and physical constants: the magnitude of the electron charge e, Planck's constant h, and the speed of light in vacuum c. Multiplying S I (λ) by Φ ae (λ), we obtain the current I ph generated by photons in the photosensitive elements of the matrix, which forms a charge Q = I ph ·τ e . during the exposure of the frame τ e . When calculating the signal D ij at the output of the ADC of the thermograph, coming from the ij elements of the photonic detector matrix, it is necessary to integrate the fluxes of monochromatic radiation falling on the photosensitive elements of the matrix within the used spectral region λ sw -λ lw D ij (T, τ e ) = τ l Ks ae k ij gk QD   e hc where τ of (λ) is the transmission coefficient of the optical filter, which, in the form of multilayer coatings, can be applied to the optical surfaces of the lens or the input window of the photodetector array, k QD is the conversion factor of the charge entering the input of the matrix reading unit into a digital code, and g is the gain of the amplifier built in before the ADC. Modern arrays of photonic detectors have a sufficiently wide range of signal linearity and a system for subtracting the dark current; therefore, in further calculations, it is more convenient to use the rate of rise of digital counts obtained from (6) by dividing by the exposure time of the frame τ e , which is proportional to the magnitude of the photocurrents I ph generated by the photosensitive elements of the matrix, and, consequently, to the brightness of the observed objects. Since modern lenses have good spatial resolution, the effect of the spatial frequency transfer characteristic of the lens can be ignored. The program for processing the registered digital signals with the introduction of the effective value of the thermal radiation coefficient ε eff forms a palette image of the temperature field. The body temperature T is determined by comparing the rate of rise of digital signals V ji (7) divided by the effective thermal emissivity of the body ε eff with the calibration dependence F (T) (see Figure 3 [15]) obtained using a reference emitter-the black body model. When using the value ε eff = 1, the temperature of the partial radiation T r will be the result of measuring the temperature of real bodies. Modern arrays of photonic detectors have a sufficiently wide range of signal linearity and a system for subtracting the dark current; therefore, in further calculations, it is more convenient to use the rate of rise of digital counts obtained from (6) by dividing by the exposure time of the frame τe, which is proportional to the magnitude of the photocurrents Iph generated by the photosensitive elements of the matrix, and, consequently, to the brightness of the observed objects. Since modern lenses have good spatial resolution, the effect of the spatial frequency transfer characteristic of the lens can be ignored. The program for processing the registered digital signals with the introduction of the effective value of the thermal radiation coefficient εeff forms a palette image of the temperature field. The body temperature T is determined by comparing the rate of rise of digital signals Vji (7) divided by the effective thermal emissivity of the body εeff with the calibration dependence F (T) (see Figure 3 [15]) obtained using a reference emitter-the black body model. When using the value εeff = 1, the temperature of the partial radiation Tr will be the result of measuring the temperature of real bodies. To exclude the influence of the values of g and the distance z from the main plane of the lens of the thermograph to the body on the results of calculations of T, it is necessary to use the corrected measured values of the slew rates of digital signals where gс, Kс and gm, Km are the values of the corresponding parameters when calibrating the thermograph and during temperature measurements, respectively. To reduce the measurement error of T, it is necessary to use the optimal for a given range of measured temperatures Tmin-Tmax spectral region of the thermal radiation registration. Its position on the wavelength axis λ is selected from the condition of the need to provide a sufficient signal-to-noise ratio at the maximum possible slope of the V (T) dependence, which leads to the dependences shown in Figure 4 [11]. It follows from them that, if it is necessary to measure temperatures T > 800 °С, the spectral range from 0.7 to 0.8 μm should be used. To exclude the influence of the values of g and the distance z from the main plane of the lens of the thermograph to the body on the results of calculations of T, it is necessary to use the corrected measured values of the slew rates of digital signals where g c , K c and g m , K m are the values of the corresponding parameters when calibrating the thermograph and during temperature measurements, respectively. To reduce the measurement error of T, it is necessary to use the optimal for a given range of measured temperatures T min -T max spectral region of the thermal radiation registration. Its position on the wavelength axis λ is selected from the condition of the need to provide a sufficient signal-to-noise ratio at the maximum possible slope of the V (T) dependence, which leads to the dependences shown in Figure 4 [11]. It follows from them that, if it is necessary to measure temperatures T > 800 • C, the spectral range from 0.7 to 0.8 µm should be used.

Determining Temperature When Using Multiple Spectral Regions
When the thermograph uses one part of the spectrum, the automatic determination of the unknown effective coefficient of thermal radiation εeff is impossible, and to calculate the correct value of the temperature T, you need to know its value.
Registration of thermal radiation in two parts of the spectrum makes it possible to exclude the influence of only the emissivity of "gray" bodies on the results of determining the temperature, since for an unambiguous solution the system of two equations must contain only two unknown variables T and ε = const.
In the case of using three spectral regions, it becomes possible, when imposing an additional condition on the behavior of ε (λ), for example, a linear dependence ε (λ) = a + bλ with the obvious constraint 0 < ε(λ) < 1, to compose a system of three nonlinear equations and find T when solving it [20]. This approach is applicable in a narrow wavelength range, where the ε (λ) dependence is well described by a linear, exponential, hyperbolic, or other suitable dependence with two parameters.
A further increase in the number of spectral regions used requires registration of heat fluxes with high accuracy [21], which is achievable only in laboratory conditions. The capabilities of modern multispectral photodetector arrays for measuring T from thermal radiation in several spectral regions are limited due to a noticeable overlap of spectral characteristics. For example, the Imec corporation advertises hyperspectral photodetector arrays that register radiation in 16 narrow areas in the 460-630 nm range [22]. However, their application in multispectral thermography is complicated by the presence of residual radiation transmission of the used optical filters on the wings of their spectral characteristics. Practical application in monitoring the temperature of bodies with an unknown spectral coefficient of thermal radiation ε (λ) can have simpler thermographs that record thermal radiation in three spectral regions.

Determination of T Using Three Narrow Spectral Regions
When analyzing the capabilities of high-temperature thermography with registration of thermal radiation in three narrow spectral regions with the wavelengths of centers λ1, λ2 and λ3, the Wien approximation can be used to simplify the calculated expressions. When imposing restrictions on the behavior of ε (λ), the exponential dependence is used

Determining Temperature When Using Multiple Spectral Regions
When the thermograph uses one part of the spectrum, the automatic determination of the unknown effective coefficient of thermal radiation ε eff is impossible, and to calculate the correct value of the temperature T, you need to know its value.
Registration of thermal radiation in two parts of the spectrum makes it possible to exclude the influence of only the emissivity of "gray" bodies on the results of determining the temperature, since for an unambiguous solution the system of two equations must contain only two unknown variables T and ε = const.
In the case of using three spectral regions, it becomes possible, when imposing an additional condition on the behavior of ε (λ), for example, a linear dependence ε (λ) = a + bλ with the obvious constraint 0 < ε(λ) < 1, to compose a system of three nonlinear equations and find T when solving it [20]. This approach is applicable in a narrow wavelength range, where the ε (λ) dependence is well described by a linear, exponential, hyperbolic, or other suitable dependence with two parameters.
A further increase in the number of spectral regions used requires registration of heat fluxes with high accuracy [21], which is achievable only in laboratory conditions. The capabilities of modern multispectral photodetector arrays for measuring T from thermal radiation in several spectral regions are limited due to a noticeable overlap of spectral characteristics. For example, the Imec corporation advertises hyperspectral photodetector arrays that register radiation in 16 narrow areas in the 460-630 nm range [22]. However, their application in multispectral thermography is complicated by the presence of residual radiation transmission of the used optical filters on the wings of their spectral characteristics. Practical application in monitoring the temperature of bodies with an unknown spectral coefficient of thermal radiation ε (λ) can have simpler thermographs that record thermal radiation in three spectral regions.

Determination of T Using Three Narrow Spectral Regions
When analyzing the capabilities of high-temperature thermography with registration of thermal radiation in three narrow spectral regions with the wavelengths of centers λ 1 , λ 2 and λ 3 , the Wien approximation can be used to simplify the calculated expressions. When imposing restrictions on the behavior of ε (λ), the exponential dependence is used With a symmetric arrangement of λ 1 and λ 3 relative to λ 2 and the presence of calibration dependences F i (T), the system of corresponding equations can be represented in the form Multiplying the first equation by the third and excluding b, and then ε 2 , one can obtain a nonlinear equation for finding T: Since the product f 1 f 3 does not depend on b, the slope of the intersection of the abscissa of the curve ξ (T) for the exponential dependence ε (λ) depends only on the value of ε 2 , which is clearly seen from Figure 5 and confirms the analytical expression for the slope of the change in the function ξ (T) in point T = T real Thus, the solution (10) of nonlinear Equation (9) with a sufficient distance ΔΛ between the used narrow spectral regions of registration of thermal radiation makes it possible to provide an acceptable error in determining the temperature T. In practice, materials for the formation of narrow spectral characteristics of optical filters with high attenuation on the wings, which could be applied on the surface of the photodetector matrix, are absent. The use of RGB matrices when measuring T is ineffective since ε (λ) of many materials in the visible region has dependences far from exponential. Therefore, one should consider the possibility of forming three parts of the spectrum, shifted in the NIR spectrum range when using commercially available RGB photodetector matrices.

Formation of Three Spectrum Regions, Shifted in the NIR Range
Analysis of the spectral quantum efficiency ηex (λ) of digital single-chip color CMOS matrices shows the possibility of forming three spectrum regions shifted in the NIR range when using an external optical filter. For example, the MT9V034C12STC matrix with a built-in 10-bit ADC has a good ηex (λ) value in the NIR spectrum range (see Figure 6), Using the expression for estimation the error in indirect measurements, one can find the mean-root-square deviation of the quantity ξ from the relative errors in measuring the fluxes δf 1 , δf 2, δf 3 : transforming which, it is possible to obtain the dependence of the relative error in determining the temperature δT on δf : When using wavelengths λ 1 = 0.6, λ 2 = 0.75, and λ 3 = 0.9 µm, we obtain the transfer coefficient of the relative error [21] when measuring a temperature of 1000 • C, equal to 1.95 or about 2. Then, at δf 2 = 0.002, we obtain δT ≈ 0.004. When narrowing the used spectral region, i.e., at ∆Λ = 0.1 µm, we obtain the transfer coefficient of the random component of the error, equal to about 4.5, and δT = 0.009 or about one percent.
Thus, the solution (10) of nonlinear Equation (9) with a sufficient distance ∆Λ between the used narrow spectral regions of registration of thermal radiation makes it possible to provide an acceptable error in determining the temperature T. In practice, materials for the formation of narrow spectral characteristics of optical filters with high attenuation on the wings, which could be applied on the surface of the photodetector matrix, are absent. The use of RGB matrices when measuring T is ineffective since ε (λ) of many materials in the visible region has dependences far from exponential. Therefore, one should consider the possibility of forming three parts of the spectrum, shifted in the NIR spectrum range when using commercially available RGB photodetector matrices.

Formation of Three Spectrum Regions, Shifted in the NIR Range
Analysis of the spectral quantum efficiency η ex (λ) of digital single-chip color CMOS matrices shows the possibility of forming three spectrum regions shifted in the NIR range when using an external optical filter. For example, the MT9V034C12STC matrix with a built-in 10-bit ADC has a good η ex (λ) value in the NIR spectrum range (see Figure 6), since it is also designed to work in low light at night. Since the freedom to choose the spectral regions for registering thermal radiation for RGB matrices is limited, in addition to an external optical filter, a combination of the received RGB signals must also be used. It should be understood that due to the small separation of the effective wavelengths of these sections, high accuracy of measurements of heat fluxes and the imposition of additional restrictions on the behavior of ε (λ) will be required in order to form a stable method for determining the temperature T. To ensure the required accuracy, it is necessary to take into account possible deviations of the calibration dependences from theoretically calculated. It will also be necessary to take into account the dependence of the statistical parameters of the noise of the photodetector channel of the three-spectral thermograph on the value of the incident fluxes of thermal radiation. By using an external bandpass filter with the transmission of radiation in the NIR range, where an increase in the transmission of G and B filters of the Bayer mosaic is observed, and by calculating the difference between the signals formed by R and G by the light-sensitive elements of the matrix, it is possible to form the resulting quantum efficiencies η k (λ) in three spectral regions. For example, using a bandpass filter with a transmission in the 630-830 nm range and an MT9V034C12STC matrix, three dependences η (λ) can be formed, which are indicated in Figure 7 as 1st, 2nd, and 3rd sections. They completely overlap, but their effective wavelengths differ from each other, which makes it possible to use them to determine the temperature T, the effective value of the thermal emissivity ε eff2 , and the conditional temperatures of partial radiation T rk in the formed spectral regions. It should be understood that due to the small separation of the effective wavelengths of these sections, high accuracy of measurements of heat fluxes and the imposition of additional restrictions on the behavior of ε (λ) will be required in order to form a stable method for determining the temperature T. To ensure the required accuracy, it is necessary to take into account possible deviations of the calibration dependences from theoretically calculated. It will also be necessary to take into account the dependence of the statistical parameters of the noise of the photodetector channel of the three-spectral thermograph on the value of the incident fluxes of thermal radiation.

Temperature Determination Method
Proceeding similarly to the approach (9) already considered above, using narrow spectral regions, it is possible to compose a system of three nonlinear equations. The use of the exponential dependence ε (λ) = ε2•exp [b (λ − λ2)] allows one to correctly describe the behavior ε (λ) of many bodies in a limited spectral region. For small values of the product b(λ − λ2) < 0.05, which is the case for many structural materials in practice, the exponential expression can be simplified by expanding the exponent in a power series and discarding terms with powers greater than unity, which gives ε (λ) = ε2 + b (λ − λ2). Then, It should be understood that due to the small separation of the effective wavelengths of these sections, high accuracy of measurements of heat fluxes and the imposition of additional restrictions on the behavior of ε (λ) will be required in order to form a stable method for determining the temperature T. To ensure the required accuracy, it is necessary to take into account possible deviations of the calibration dependences from theoretically calculated. It will also be necessary to take into account the dependence of the statistical parameters of the noise of the photodetector channel of the three-spectral thermograph on the value of the incident fluxes of thermal radiation.

Temperature Determination Method
Proceeding similarly to the approach (9) already considered above, using narrow spectral regions, it is possible to compose a system of three nonlinear equations. The use of the exponential dependence ε (λ) = ε 2 ·exp [b (λ − λ 2 )] allows one to correctly describe the behavior ε (λ) of many bodies in a limited spectral region. For small values of the product b(λ − λ 2 ) < 0.05, which is the case for many structural materials in practice, the exponential expression can be simplified by expanding the exponent in a power series and discarding terms with powers greater than unity, which gives ε (λ) = ε 2 + b (λ − λ 2 ). Then, assuming that g = 1 and k ij ≈ 1 and taking into account (7), we can compose the following system of equations where M (λ, T) is the spectral density of the surface luminosity of an absolutely black body described by the corresponding Planck's formula.
After integration, (14) can be represented in a more convenient form for further calculations where F RF-GF (T), F RF (T), and F BF (T) correspond to the calibration dependences of the thermograph, which are determined experimentally using the black body model, q RF-GF (T), q RF (T), and q BF (T) are the integrals obtained by calculation and determining the dependences of these integrals on the temperature at a unit value of the coefficient b of the used dependence ε (λ) in the corresponding spectral region, and the dependences µ RF-GF (T), µ RF (T), and µ BF (T) are integrals q RF-GF (T), q RF (T), and q BF (T), normalized to the corresponding calibration dependences.
Due to the overlap of the used spectral regions and the inevitable spread of the recorded values of V RF -GF , V RF , and V BF due to the presence of noise, the numerical solution of system (15) gives an unacceptable spread of the obtained values of T. To obtain a more stable solution, an additional constraint on the behavior ε (λ) in the form of a functional connection b = γ·p (ε 2 ) between the slope coefficient b of the spectral dependence ε (λ) of the controlled surface, the value of the coefficient γ and the value of ε 2 in the form of a polynomial of the second degree with given values of the coefficients p 2 , p 1 , and p 0 . The feasibility of introducing this functional bond is due to an increase in the values of ε (λ) when the controlled surface is heated, which occurs due to its oxidation with atmospheric oxygen. As an example, Figure 8 shows the dependences ε (λ) of clean (black lines) and oxidized (gray lines) steel surface. Oxidation of the surface leads to a decrease in the slope coefficient b from the maximum possible value to close to zero. provides an adequate calculation of b when the value of ε2 changes during oxidation with atmospheric oxygen. For tungsten, the value of γ is γ = −1.07, and for ferrous metals, γ ≈ −0.5. Comparing Figures 8 and 9, it can be seen that the slope coefficient b during steel oxidation tends to zero as ε2 grows.  Several sets of coefficients p 2 , p 1 , and p 0 of the used polynomial can be proposed. For example, when describing the relationship between the values of b and ε 2 for bodies with possible ranges of variation of the thermal radiation coefficient 0.05 < ε 2 < 1, one can use a power polynomial in the form which describes a smooth decrease in the absolute value of |b| with an increase in ε 2 in the process of surface oxidation with atmospheric oxygen. The corresponding graphical dependencies are shown in Figure 9. The curves in the lower half-plane, where γ < 0, refer to metals, and in the upper, to dielectrics. Since the spectral dependences of the thermal radiation coefficients of most structural materials are known [23], the value of γ is calculated in advance from the value of b and ε 2 of pure unoxidized material, which provides an adequate calculation of b when the value of ε 2 changes during oxidation with atmospheric oxygen. For tungsten, the value of γ is γ = −1.07, and for ferrous metals, γ ≈ −0.5. Comparing Figures 8 and 9, it can be seen that the slope coefficient b during steel oxidation tends to zero as ε 2 grows.  The use of dependence (16), i.e., the involvement of a priori information on the possible behavior ε (λ) of the surface of the controlled body, allows one to significantly reduce the requirements for the accuracy of measuring heat fluxes and to use overlapping spectral regions in determining the temperature T. It should be noted that an increase in the distance between the used spectrum regions, on the one hand, reduces the requirements for the magnitude of measurement errors, but, on the other hand, due to the deviation of the behavior ε (λ) of bodies in the visible region from the linear law leads to violation of the condition ε (λ) = ε2 + b (λ − λ2), which, in turn, leads to an increase in the The use of dependence (16), i.e., the involvement of a priori information on the possible behavior ε (λ) of the surface of the controlled body, allows one to significantly reduce the requirements for the accuracy of measuring heat fluxes and to use overlapping spectral regions in determining the temperature T. It should be noted that an increase in the distance between the used spectrum regions, on the one hand, reduces the requirements for the magnitude of measurement errors, but, on the other hand, due to the deviation of the behavior ε (λ) of bodies in the visible region from the linear law leads to violation of the condition ε (λ) = ε 2 + b (λ − λ 2 ), which, in turn, leads to an increase in the error in determining T. This is the main reason that prevents the effective use of the R, G, and B regions of the visible spectrum in high-temperature thermography.
When using (16) and smoothing V RF-GF (T), V RF , and V BF using a median digital filter that suppresses impulse emissions of defective elements of the used photodetector matrix, the system of Equation (15) takes the form Subtracting the first equation from the third equation of system (16), we can obtain a quadratic equation for finding the dependence of ε 2 on the estimate of the temperature T *: Choosing a positive solution to this equation, and substituting it into the second equation of system (17), we obtain a nonlinear equation for finding the desired temperature value The disadvantage of the obtained solution is the uncertainty of the calculated values of ε 2 at values of γ close to zero. To eliminate it for small values of γ, one can use the smallness of the product bµ RF (T*) in comparison with the value of ε 2 (T*). This makes it possible to approximate the dependence ε 2 (T*) as a ratio V RF /F RF (T * ) or, more precisely, V RF /F RF (T * ) − bµ 2 (T * ), which makes it possible to obtain an approximate expression for Substituting (21) into the first and second equations of system (17) and forming their difference, we obtain the resulting equation for determining T To linearize Equation (15), we can multiply it by the ratio V RF /F RF (T * ), which greatly simplifies the search for a numerical solution T:  In the numerical solution of Equation (19) or (22), the range in which the solution is sought is set, i.e., the estimate of the temperature T*. As can be clearly seen from the illustration of the process of finding a solution shown in Figure 9, the low-temperature limit of this range is determined by the minimum of the three temperatures of partial radiation Tr min, recorded by the thermograph. The high-temperature limit for metals (γ < 0) is determined by the temperature of the spectral ratio for the first and third spectral illustration of the process of finding a solution shown in Figure 9, the low-temperature limit of this range is determined by the minimum of the three temperatures of partial radiation Tr min, recorded by the thermograph. The high-temperature limit for metals (γ < 0) is determined by the temperature of the spectral ratio for the first and third spectral regions Tsr13, which is calculated by the thermograph when solving the equations. For dielectrics (γ > 0), the Tsr13 value is usually less than the true surface temperature, so for them, the high-temperature boundary is chosen equal to 1.15•Tsr13. Figure 11. Illustration of the process of finding a solution to Equation (22) when determining the temperature of the region of maximum heating of the tungsten tape of the standard incandescent lamp SI10-300. In the numerical solution of Equation (19) or (22), the range in which the solution is sought is set, i.e., the estimate of the temperature T*. As can be clearly seen from the illustration of the process of finding a solution shown in Figure 9, the low-temperature limit of this range is determined by the minimum of the three temperatures of partial radiation T r min , recorded by the thermograph. The high-temperature limit for metals (γ < 0) is determined by the temperature of the spectral ratio for the first and third spectral regions T sr13 , which is calculated by the thermograph when solving the equations. For dielectrics (γ > 0), the T sr13 value is usually less than the true surface temperature, so for them, the high-temperature boundary is chosen equal to 1.15·T sr13 . Figure 11 shows an illustration of the process of finding a solution to Equation (22) using signals recorded by a three-spectral thermograph based on an MT9V034C12STC array of photodetectors. The spectral regions in which the thermograph records the thermal radiation correspond to those shown in Figure 7. The thermograph was calibrated according to the emission of the Mikron M390 black body model. The relative error in measuring T was 0.5% of the measured value T. This makes it possible to measure the temperature of bodies with a known value of ε eff with a relative error of about 0.5%. When measuring the temperature of the tungsten tape of the standard incandescent lamp SI10-300 in advance for the second spectral region used by the thermograph, the dependence of the effective thermal radiation coefficient of tungsten ε eff2 (T) was calculated taking into account the transmission of the quartz window of this lamp. By introducing the obtained value of ε eff2 (T m ) into the thermograph and adjusting the current through the tungsten tape of the lamp, we achieved that the temperature of the tungsten tape in the region of its maximum heating corresponded to T m . For the case shown in Figure 11, the T m value measured in the second spectral region was 2019 • C.
The measured signals V RF -GF , V RF , and V BF , due to the presence of noise and the deviation of the thermograph characteristics from the calculated values, have fluctuations, which lead to a scatter of the obtained values of the measured temperature T * m when solving Equation (22). For the example shown in Figure 11, the obtained estimate of T * m when using the value γ = −1.07 µm −1 turned out to be approximately 3 • C lower than T m . The root-mean-square deviation of the temperature estimation of the most heated region of the tungsten tape during the averaging was about 1.5% of the measured value T. The dashed line in Figure 11 shows how the value of ε 2 (T*) changes in the process of finding a solution.
The error in determining the temperature in the described way is also affected by the error in setting the value of the coefficient γ. The performed modeling shows that, for metals, the relative errors in determining the true temperature ∆T/T with deviations of the used value γ* from its true value γ by ±40% in practice, as shown in Figure 12, do not exceed 0.015. This is due to the weak influence of the slope coefficient b of the spectral dependence ε (λ) on the obtained values of ε eff . Such errors are permissible when solving the overwhelming majority of practical problems of contactless temperature measurement of various heat engineering processes. the tungsten tape in the region of its maximum heating corresponded to Tm. For the case shown in Figure 11, the Tm value measured in the second spectral region was 2019 °C.
The measured signals VRF-GF, VRF, and VBF, due to the presence of noise and the deviation of the thermograph characteristics from the calculated values, have fluctuations, which lead to a scatter of the obtained values of the measured temperature T * m when solving Equation (22). For the example shown in Figure 11, the obtained estimate of T * m when using the value γ = −1.07 μm −1 turned out to be approximately 3 °C lower than Tm. The root-mean-square deviation of the temperature estimation of the most heated region of the tungsten tape during the averaging was about 1.5% of the measured value T. The dashed line in Figure 11 shows how the value of ε2(T*) changes in the process of finding a solution.
The error in determining the temperature in the described way is also affected by the error in setting the value of the coefficient γ. The performed modeling shows that, for metals, the relative errors in determining the true temperature ΔT/T with deviations of the used value γ* from its true value γ by ±40% in practice, as shown in Figure 12, do not exceed 0.015. This is due to the weak influence of the slope coefficient b of the spectral dependence ε (λ) on the obtained values of εeff. Such errors are permissible when solving the overwhelming majority of practical problems of contactless temperature measurement of various heat engineering processes.  The described method for determining the temperature with an unknown coefficient of thermal emissivity makes it possible to determine ε eff2 = M(T r2 )/M(T * m ) from the obtained values of T r2 and T * m . Using the obtained value of ε eff2 , one can switch to the usual mode of temperature measurement taking into account ε eff2 , which reduces the variation in readings. As follows from the dependencies shown in Figure 4, the relative error in determining the temperature in the range from 0.7 to 0.8 µm will be an order of magnitude smaller than the relative error in determining ε eff2 . This indicates the possibility of using three-spectral thermographs when setting up and periodically monitoring complex high-temperature technological processes. However, in order to reduce the error in determining ε eff2 , it is necessary to use algorithms for averaging the recorded fields of thermal radiation, both in time and over the region of maximum heating.

Tehnique of Averaging Signals over the Area of Maximum Heating
An unpleasant feature of MT9V034C12STC CMOS sensors with a global exposure mode is a rapid increase in the root mean square deviation of the total noise at high thermal radiation fluxes and, accordingly, small frame exposure times τ e , which is illustrated in Figure 13. At high fluxes, dots begin to appear on the resulting image, the brightness of which is noticeably different from the brightness of neighboring pixels. Moreover, their number increases with decreasing τ e . Matrix elements corresponding to these points are usually called "faulty", since an increase in the currents generated by them is caused by the presence of defects in the structures of the insulator and the corresponding keys of the multiplexer. Therefore, manufacturers of high-temperature thermographs based on CMOS structures divide the range of measured temperatures into a number of subranges that differ in different values of the lens aperture D/f l . Partially fluctuating emissions of defective elements are excluded by median filtration. However, if it is necessary to control the heat engineering processes simultaneously in one wide temperature range, it is necessary to take into account the influence of the nonstationarity of the matrix noise and carefully take into account the existing deviations of the real characteristics of the thermograph from the calculated ones.
An unpleasant feature of MT9V034C12STC CMOS sensors with a global exposure mode is a rapid increase in the root mean square deviation of the total noise at high thermal radiation fluxes and, accordingly, small frame exposure times τe, which is illustrated in Figure 13. At high fluxes, dots begin to appear on the resulting image, the brightness of which is noticeably different from the brightness of neighboring pixels. Moreover, their number increases with decreasing τe. Matrix elements corresponding to these points are usually called "faulty", since an increase in the currents generated by them is caused by the presence of defects in the structures of the insulator and the corresponding keys of the multiplexer. Therefore, manufacturers of high-temperature thermographs based on CMOS structures divide the range of measured temperatures into a number of subranges that differ in different values of the lens aperture D/fl. Partially fluctuating emissions of defective elements are excluded by median filtration. However, if it is necessary to control the heat engineering processes simultaneously in one wide temperature range, it is necessary to take into account the influence of the nonstationarity of the matrix noise and carefully take into account the existing deviations of the real characteristics of the thermograph from the calculated ones.

Ensuring the Invariance of the Measured Values of the Maximum Temperature T max to the Image Size and the Nonstationarity of the Matrix Noise
In high-temperature thermography, the sizes of objects whose temperature must be controlled can be varied, and the temperature distribution over their surface is almost always inhomogeneous. In the process of heating objects, their temperature field can change quite quickly. Therefore, it is practically impossible to conduct a visual assessment of the temperature values from the image of the obtained temperature field. We have to implement the function of determining the maximum body temperature T max into the software of hightemperature thermographs. The algorithms used in it must ensure that the above principles are followed to increase the reliability of temperature measurements. In particular, they must ensure the invariance of the results of determining the maximum temperature T max to a change in the size of the image of the monitored bodies and the nonstationarity of the noise dispersion of the used photodetector matrix, i.e., the dependence of its noise on the value of the incident thermal radiation flux.
Neglect of these principles, when observing a sequence of frames, leads to a scatter of the obtained T max readings and an abrupt displacement of the point with the maximum temperature value over the zone of greatest heating of the generated temperature field. The reason for this phenomenon is due to the statistics of the number of registered photons, which obeys the Poisson distribution, the intrinsic noise of the thermograph chamber, and fluctuations in the temperature of the surface of the heated body. Fluctuations in surface temperature are caused by its interaction with colder convection air currents, which are formed due to the temperature difference between the heated body and the air temperature.
The presence of inhomogeneity of heating leads to distortion of the high-temperature mode of the histogram of the thermal image, which is illustrated in Figure 14. Even the radiating cavity of black body models can have a temperature gradient of several degrees, which distorts the corresponding histogram, as shown in Figure 15. The discreteness of the obtained digital images and the presence of fluctuations require averaging the obtained readings, which are used to calculate the maximum temperature. This raises the problem of optimal selection of image elements or pixels involved in averaging. temperature value over the zone of greatest heating of the generated temperature field. The reason for this phenomenon is due to the statistics of the number of registered photons, which obeys the Poisson distribution, the intrinsic noise of the thermograph chamber, and fluctuations in the temperature of the surface of the heated body. Fluctuations in surface temperature are caused by its interaction with colder convection air currents, which are formed due to the temperature difference between the heated body and the air temperature.
The presence of inhomogeneity of heating leads to distortion of the high-temperature mode of the histogram of the thermal image, which is illustrated in Figure 14. Even the radiating cavity of black body models can have a temperature gradient of several degrees, which distorts the corresponding histogram, as shown in Figure 15. The discreteness of the obtained digital images and the presence of fluctuations require averaging the obtained readings, which are used to calculate the maximum temperature. This raises the problem of optimal selection of image elements or pixels involved in averaging.
where the first term of the sum describes the variance of the radiation and shot noises, the second describes the geometric noises (caused by the inhomogeneity of the sensitivity of the matrix elements) with 2 n σ variance, and the third describes the noises of the readout unit of the matrix with variance 2 Dread σ . Expression (24) determines the algorithm for choosing the required averaging interval. Figure 15. The image of the output diaphragm of the M390 black body model recorded by the thermograph in the second spectral region at a temperature of its radiating core of 800 °C and an exposure time of 130 ms, as well as an image of its central zone with a range of D2 signals from 760 to 820 ADC counts and the corresponding histogram approximated using normal distribution.

Method for Obtaining Averaged Values k V * on the High-Temperature Slope of Histograms
Since the Poisson distribution used to describe radiation and shot noises is close to the normal distribution with a sufficiently large number of registered photons, the averaging interval can be chosen to be equal to ( ) Additionally, the algorithm should include a method for clearing the histograms Figure 15. The image of the output diaphragm of the M390 black body model recorded by the thermograph in the second spectral region at a temperature of its radiating core of 800 • C and an exposure time of 130 ms, as well as an image of its central zone with a range of D 2 signals from 760 to 820 ADC counts and the corresponding histogram approximated using normal distribution.
As seen in Figure 14d, when determining the maximum temperature T max , it is necessary to calculate the averaged values of digital readings D * k (T max ) on the high-temperature slopes of the right wing of the thermal image histograms in the used spectral regions. The determination technique of D * k (T max ) can be based on the existing relationship between the average values of the recorded signals and their variance where the first term of the sum describes the variance of the radiation and shot noises, the second describes the geometric noises (caused by the inhomogeneity of the sensitivity of the matrix elements) with σ 2 n variance, and the third describes the noises of the readout unit of the matrix with variance σ 2 Dread . Expression (24) determines the algorithm for choosing the required averaging interval.

Method for Obtaining Averaged Values V *
k on the High-Temperature Slope of Histograms Since the Poisson distribution used to describe radiation and shot noises is close to the normal distribution with a sufficiently large number of registered photons, the averaging interval can be chosen to be equal to 6σ ΣDk V * k in the determination of V * k . Additionally, the algorithm should include a method for clearing the histograms y hk (n) = H k [x hk (n)] from single outliers on its right wing, taking into account the dependence σ ΣDk V k .
Taking into account the described factors, the method of determination V * k consists of the following sequence of actions: (1) After receiving the frame, histograms are calculated for all three of its layers y hk (n) = H k [x hk (n)] with the width of pockets or bins equal to one or one level of the ADC conversion range.
(5) Normalize dependencies S xWk (n) by dividing them by S Fk (n), which corresponds As seen in Figure 16, the values of the elements of the series S xW (n)/S F (n) gradually decrease to the desired average value D when the point x moves from the maximum value D to the left, since the value of the integral is in the denominator of the ratio ∞ D−3σ ΣD W(x)dx ≈ 1. When moving to the left, the ratio S xW (n)/S F (n) has a small steepness, determined by the value of σ, which makes it possible to obtain a slightly biased estimate D with a small number of pixels involved in averaging. to the relationship As seen in Figure 16 Figure 17. The thermal image of the spiral of this lighting lamp in the second spectral region when viewed from the side is shown in the inset of Figure 17a. Uneven heating of the surface leads to stretching of the histograms towards lower signal values, which is illustrated by the dependencies in Figure 14d. In contrast to the dependencies shown in Figure 16, with uneven heating, the sum S xW (n) after crossing the dependence S xW (n)/S F (n) continues to increase, as shown in Figure 17b. However, the estimated value D * k at the point of intersection of the dependencies changes only slightly due to the low steepness of the ratio S xW (n)/S F (n). The described method for determining the average value of digital signals V * k on the right slope of the histogram was specially developed to obtain stable readings of the maximum temperature T max when inspecting bodies with uneven heating. Only 18 brightest pixels fall into the averaging interval in the second part of the spectrum. It is clearly seen that the shape of the high-temperature sections of the histograms has a random form with the presence of gaps between their bins or pockets.

Correction of the Deviation of the Real Characteristics of the Photodetector Array from the Calculated Ones
As shown above, algorithms for determining the temperature T in high-temperature thermography with an uncertainty in the value of the thermal emissivity must contain an option for calculating T and εeff2. It is based on solving a system of three nonlinear Equation (14) with the imposition of additional conditions on the behavior of the spectral The estimation of the current variance can be carried out by calculating the second moment of the distribution W (x), similarly to calculating the average value. Verification of the effectiveness of using this approach for evaluation σ ΣDk (V k ) showed that it leads to underestimated values in comparison with the dependences obtained in the course of calibration, shown in Figure 13. Therefore, the Gaussian approximations of the resulting histograms shown in Figure 17a, with lighter curves with dots, become narrower due to the limited sample, which increases the scatter of the determined values D * k . The use of dependencies in the calculation algorithm predetermined when calibrating the thermograph σ ΣDk (V k ) makes it possible to reduce the scatter of the estimate D * k , which is well illustrated in Figure 17c. Figure 17b illustrates the process of automatic search for an estimate D * k at the intersection of the dependencies S xWk (n)/S Fk (n) and S xWk (n). The root-mean-square deviation σ D k of the obtained estimates is small and ranges from two to four ADC readings, which ensures the stability of the calculated values of the maximum temperature of partial radiation T rmax even with small sizes of the thermal image of an object with uneven heating.

Correction of the Deviation of the Real Characteristics of the Photodetector Array from the Calculated Ones
As shown above, algorithms for determining the temperature T in high-temperature thermography with an uncertainty in the value of the thermal emissivity must contain an option for calculating T and ε eff2 . It is based on solving a system of three nonlinear Equation (14) with the imposition of additional conditions on the behavior of the spectral coefficient of thermal radiation ε (λ). In the final expressions for calculating T and ε eff2 , three estimates V * k of the recorded values are used, averaged over the brightest areas of the thermal image of the object. Usually, in these areas, the thermal conductivity is the best, that is, the surface is the least oxidized. Additionally, two more groups of previously obtained temperature dependences are used. The first of them is the calibration dependences F 1 (T) = F RF-GF (T), F 2 (T) = F RF (T), and F 3 (T) = F BF (T), which are determined during calibration according to the black body model. The second group is the calculated temperature dependences µ 2 (T) = µ RF (T) and µ 3 (T) − µ 1 (T) = µ BF (T) − µ RF-GF (T).
With knowledge of the spectral quantum efficiencies of the photosensitive elements of the matrix η k (ε) and taking into account the parameters of the calibration scheme, the calculated or hypothetical dependences F hk (T) and the calibration dependences F k (T) determined experimentally during calibration should coincide. However, the inevitable deviations η k (λ) of the elements of photodetectors matrices, lens parameters, and transmission coefficients τ of (λ) of optical filters during their manufacture from the averaged characteristics given in the documentation lead to differences in real dependences from those calculated. This is the basis for calibrating thermographs using standard emitters.
If it is necessary to measure T in one wide range of measured temperatures, it is necessary to take into account the deviations of the real characteristics of the photodetector array from the calculated ones. Comparison of the calibration dependences F k (T), obtained during the calibration according to the black body model M390, with the calculated ones F hk (T) shows that the values of F k (T) at high temperatures exceed the calculated values. The dependencies in Figure 18, obtained in the temperature range from 800 to 1700 • C with a discreteness of 100 • C, clearly demonstrate that, at high temperatures, i.e., at low frame exposure times τ e , a noticeable nonlinearity of the light characteristic of the MT9V034C12STC matrix in the global exposure mode appears. It is difficult to identify the reasons for this nonlinearity, since the technical documentation does not contain information about the used scheme for accumulating and reading signals generated by photodiodes. Only the dynamic range of illumination change is given, which is more than 55 dB in linear mode.
(Tm) are combined with the experimental ones Fk (Tm) at the calibration point T0, where the light characteristic of the matrix is almost linear. For example, for the MT9V034C12STC matrix, this point, as follows from Figure 18b, corresponds to the temperature of the radiating core of the black body model 1000 °С. Then, the initial values of the one-dimensional array of the discrete correcting function kFk0 (Tm) will be equal to the ratio Fk (T0)/Fhk (T0). This operation is allowable, since in nonlinear Equations (20) and (23), only the ratio s k V * (T*)/Fk (T*) are used. The method for calculating kF k(Tm) is as follows.  This behavior of the calibration curves does not cause problems in determining the temperatures of partial radiation T rk by comparing the measured values of the rise rates of digital readings V * k with the calibration curves (see Figure 3). However, when calculating the temperature of the spectral ratio T sr13 , the true temperature and the effective coefficient of thermal radiation ε eff2 , the nonlinearity of the light characteristic can lead to ambiguity of solutions. The foregoing is well illustrated by the temperature dependence of the ratio F 1 (T)/F 3 (T) for the available sample of the thermograph, shown in Figure 19. If the calculated ratios F h1 (T)/F h3 (T) are increased monotonically with increasing temperature, then a maximum is observed for the ratios of the obtained calibration curves F 1 (T)/F 3 (T), after which this ratio decreases. Obviously, starting from some values of T, it will not be possible to unambiguously determine the temperature of the spectral ratio T sr13 from the obtained relation. Similarly, corresponding distortions can arise when solving the system of nonlinear equations used to determine T. Changes in the values of the experimental calibration dependences Fk (Tm) are mainly caused by the inaccuracy of setting the lens f-number, which leads to a deviation of the actual D/fl ratio from the value used in theoretical calculations. These ratios, even for lenses of the same series, slightly differs from each other, which leads to a parallel shift of the experimental calibration dependences Fk (Tm) up or down due to small differences in the coefficient K (3). In some cases, small shifts of the calibration dependences Fk (Tm) relative to each other are observed, which is associated with insignificant displacements of the microlens arrays relative to the photosensitive elements of CMOS matrices at the production stage. Further, for the calibration points located to the right of T0, i.e., for Tm > T0, new values of the correcting function are calculated and for points with Tm < T0, the old values of kFk0 (Tm) remain. In this case, we obtain an array of values of the correcting function kFk(Tm), for which the first values at Tm ≤ T0 are equal to Fk (T0)/Fhk (T0), and the rest are additionally multiplied by the ratios Fhk (Tm)/Fk (Tm). The corrected calibration curves Fck (Tm) are close to the calculated ones Fhk (Tm). Correction of the obtained values k V * is somewhat more complicated, since it is necessary to take into account the weakening of heat fluxes due to the influence of the thermal radiation coefficient, which for real bodies is less than unity. When using the simplified method, it is advisable to pre-form the arrays of differences ΔFck (Tm) = Fk (Tm) − Fhk (Tm). Then, the correction of the recorded values of k V * Vпk can be carried out using the spline interpolation log2[ k V * (Trk)] between the calibration points log2[Fk(Tm)] with further subtraction of the value of the interpolated difference ΔFck (Trk) multiplied by the used effective thermal radiation coefficient Figure 19. The temperature dependence of the F 1 (T)/F 3 (T) ratio obtained by calibrating the available thermograph sample (solid line) and the calculated F h1 (T)/F h3 (T) (dashed line), as well as an illustration of the process of their use in determining the temperature of the spectral ratio T sr13 of a tungsten tape of a standard incandescent lamp SI10-300 at a temperature of partial radiation of its most heated section T r2 = 1206 • C.
In the case of nonlinearity of the light characteristic, it is necessary to use the correction of the registered calibration dependences F k (T) and measured values when activating the option for calculating ε eff2 and the true temperature T, as well as the temperature of the spectral ratio T sr13 . This simple correction consists of modifying the obtained discrete calibration dependences F k (T m ) by multiplying by discrete values of the correcting functions k Fk (T m ), the values of which are determined after calibrating the thermograph chamber according to the black body model. The method for calculating k Fk (T m ) is as follows. First, the calculated dependences F hk (T m ) are combined with the experimental ones F k (T m ) at the calibration point T 0 , where the light characteristic of the matrix is almost linear. For example, for the MT9V034C12STC matrix, this point, as follows from Figure 18b, corresponds to the temperature of the radiating core of the black body model 1000 • C. Then, the initial values of the one-dimensional array of the discrete correcting function k Fk0 (T m ) will be equal to the ratio F k (T 0 )/F hk (T 0 ). This operation is allowable, since in nonlinear Equations (20) and (23), only the ratio s V * k (T*)/F k (T*) are used. The method for calculating k F k (T m ) is as follows.
Changes in the values of the experimental calibration dependences F k (T m ) are mainly caused by the inaccuracy of setting the lens f-number, which leads to a deviation of the actual D/f l ratio from the value used in theoretical calculations. These ratios, even for lenses of the same series, slightly differs from each other, which leads to a parallel shift of the experimental calibration dependences F k (T m ) up or down due to small differences in the coefficient K (3). In some cases, small shifts of the calibration dependences F k (T m ) relative to each other are observed, which is associated with insignificant displacements of the microlens arrays relative to the photosensitive elements of CMOS matrices at the production stage.
Further, for the calibration points located to the right of T 0 , i.e., for T m > T 0 , new values of the correcting function are calculated and for points with T m < T 0 , the old values of k Fk0 (T m ) remain. In this case, we obtain an array of values of the correcting function k Fk (T m ), for which the first values at T m ≤ T 0 are equal to F k (T 0 )/F hk (T 0 ), and the rest are additionally multiplied by the ratios F hk (T m )/F k (T m ). The corrected calibration curves F ck (T m ) are close to the calculated ones F hk (T m ). Correction of the obtained values V * k is somewhat more complicated, since it is necessary to take into account the weakening of heat fluxes due to the influence of the thermal radiation coefficient, which for real bodies is less than unity.
When using the simplified method, it is advisable to pre-form the arrays of differences ∆F ck (T m ) = F k (T m ) − F hk (T m ). Then, the correction of the recorded values of V * k V IIk can be carried out using the spline interpolation log 2 [V * k (T rk )] between the calibration points log 2 [F k (T m )] with further subtraction of the value of the interpolated difference ∆F ck (T rk ) multiplied by the used effective thermal radiation coefficient where T rk -calculated from the registered values V * k of the temperature of partial radiation.

Discussion
The development of modern technologies for the production of complex products leads to the creation of new methods for their control. These methods and the corresponding equipment should provide the necessary reliability of control and ease of use. Non-contact measurement of high temperatures has features associated with a change in the coefficient of thermal radiation of the monitored objects during their heating and fluctuations of the recorded heat fluxes [8][9][10]24,25]. These fluctuations are caused by changes in the surface temperature of the controlled objects due to the effect of convection air flows, as well as the statistical nature of the emission of photons from the heated surface. Further development of high-temperature thermography requires the introduction of algorithms that allow one to take into account the influence of the unknown thermal radiation coefficient ε and the noted fluctuations of the recorded heat fluxes.
Currently used high-temperature thermographs use monochrome photodetector arrays and registration of thermal radiation in one part of the spectrum. They are used to control processes with consistently repeating parameters. This makes it possible, in the absence of information about the value of the thermal emissivity of the monitored objects, to use only the conditional temperature, i.e., the temperature of partial radiation T r . In this case, temperature deviations of technological processes are monitored by the deviation of T r from the specified limits of the range of permissible values.
If it is necessary to measure the temperature close to the true temperature during the adjustment of the thermal process, it is necessary to use methods and thermographic equipment that take into account the influence of the unknown coefficient of thermal radiation and the inevitable presence of fluctuations in the intensity of the recorded fluxes of thermal radiation. The methods used in LWIR thermal imaging to reduce the influence of ε and environmental radiation are not applicable in high temperature thermography. Therefore, researchers from different countries show interest in multispectral thermography and create appropriate laboratory samples [10][11][12][13]. Numerical modeling and practical testing of multispectral thermography methods convince that to solve practical problems of monitoring thermal processes, it is sufficient to use the registration of thermal in three parts of the spectrum. Although in difficult cases, for example, when monitoring laser hardening or welding processes, it is possible to use methods based on processing the spectrum of thermal radiation formed in the area of exposure to a laser beam [26].
In the most demanded temperature range from 800 to 1600 • C, typical for technological processes of heat treatment of ferrous metals, thermographs based on RGB photodetectors matrices can be used to register thermal radiation at the junction of the visible and NIR spectral ranges. The use of an external bandpass filter and a combination of signals generated by an RGB matrix makes it possible to register thermal radiation simultaneously in three overlapping spectral regions. To increase the stability of the solution of a system of three nonlinear equations, it is possible to use a priori information about the value of the slope b of the dependence ε (λ) of the clean unoxidized surface of the processed materials, which can be found in reference books. To take into account the oxidation of the surface by atmospheric oxygen when heated to high temperatures, one can use the relationship between the values of b and ε, since during oxidation, ε increases, and b falls almost to zero.
The existing inhomogeneity of the heating of bodies and the dynamic change in the image of the temperature field formed by the thermograph necessitates the introduction of the option for calculating the maximum temperature T max into the algorithms. Moreover, the algorithms for calculating T max should ensure the invariance of the determined values of the maximum temperature to the size of the image of the monitored bodies and take into account the variance of fluctuations of the recorded signals.
If it is necessary to form a wide range of measured temperatures, it is necessary to take into account the nonlinearity of the light characteristics of the photodetector matrix and the dependence of the dispersion of its noise on the intensity of the recorded heat fluxes, using the proposed methods for correcting the calibration dependences and recorded signals.
Compliance with the basic principles of reducing the uncertainty of temperature measurement and the introduction of the proposed solutions into the algorithms of operation allows creating, on the basis of modern RGB matrices of silicon photodetectors, rather cheap thermographic equipment for setting up and subsequent control of complex thermal processes.