CN114543998A - Temperature field space-time distribution measuring device based on staring type snapshot type spectral imaging - Google Patents
Temperature field space-time distribution measuring device based on staring type snapshot type spectral imaging Download PDFInfo
- Publication number
- CN114543998A CN114543998A CN202210147293.0A CN202210147293A CN114543998A CN 114543998 A CN114543998 A CN 114543998A CN 202210147293 A CN202210147293 A CN 202210147293A CN 114543998 A CN114543998 A CN 114543998A
- Authority
- CN
- China
- Prior art keywords
- spectral
- measurement
- imaging
- temperature
- detector array
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000000701 chemical imaging Methods 0.000 title claims abstract description 60
- 230000003595 spectral effect Effects 0.000 claims abstract description 84
- 238000005259 measurement Methods 0.000 claims abstract description 61
- 238000003384 imaging method Methods 0.000 claims description 40
- 238000001228 spectrum Methods 0.000 claims description 30
- 230000005855 radiation Effects 0.000 claims description 25
- 239000011159 matrix material Substances 0.000 claims description 18
- 230000006835 compression Effects 0.000 claims description 14
- 238000007906 compression Methods 0.000 claims description 14
- 238000002834 transmittance Methods 0.000 claims description 14
- 238000005070 sampling Methods 0.000 claims description 13
- 238000000034 method Methods 0.000 claims description 11
- 239000004038 photonic crystal Substances 0.000 claims description 11
- XUIMIQQOPSSXEZ-UHFFFAOYSA-N Silicon Chemical compound [Si] XUIMIQQOPSSXEZ-UHFFFAOYSA-N 0.000 claims description 8
- 229910052710 silicon Inorganic materials 0.000 claims description 8
- 239000010703 silicon Substances 0.000 claims description 8
- 238000000691 measurement method Methods 0.000 claims description 6
- 230000003287 optical effect Effects 0.000 claims description 6
- 230000008859 change Effects 0.000 claims description 4
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 238000010276 construction Methods 0.000 claims description 3
- 238000001514 detection method Methods 0.000 claims description 3
- 230000002093 peripheral effect Effects 0.000 claims description 3
- 238000012545 processing Methods 0.000 claims description 3
- 230000004044 response Effects 0.000 claims description 3
- 230000035945 sensitivity Effects 0.000 claims description 3
- 239000011148 porous material Substances 0.000 claims description 2
- 230000002123 temporal effect Effects 0.000 claims 5
- 238000005516 engineering process Methods 0.000 abstract description 5
- 238000009529 body temperature measurement Methods 0.000 description 28
- 238000002485 combustion reaction Methods 0.000 description 2
- 238000005260 corrosion Methods 0.000 description 2
- 230000007797 corrosion Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 238000004880 explosion Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000012634 optical imaging Methods 0.000 description 2
- 230000002411 adverse Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007123 defense Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 239000000446 fuel Substances 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000004861 thermometry Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J5/00—Radiation pyrometry, e.g. infrared or optical thermometry
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J5/00—Radiation pyrometry, e.g. infrared or optical thermometry
- G01J2005/0077—Imaging
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Radiation Pyrometers (AREA)
Abstract
The invention belongs to the technical field of temperature field space-time distribution measuring devices, and particularly relates to a temperature field space-time distribution measuring device based on staring type snapshot spectral imaging. The invention is realized by adopting a spectral imaging technology, the spectral imaging can simultaneously acquire the image information and the spectral information of the target to be measured, the image information can realize the measurement of the spatial distribution of the temperature field, and the spectral information can realize the accurate measurement of the temperature of the target corresponding to each pixel.
Description
Technical Field
The invention belongs to the technical field of temperature field space-time distribution measuring devices, and particularly relates to a temperature field space-time distribution measuring device based on staring type snapshot spectral imaging.
Background
The temperature field measurement can obtain fuel combustion and explosion information and efficiency evaluation, can also obtain the health state and stealth performance evaluation of aviation and aerospace engines, and has very important effect in the fields of civil safety and national defense equipment. However, temperature measurement, in particular, time distribution and space distribution measurement of a temperature field, has not been a good technical problem to solve. Currently, the temperature measurement methods can be roughly divided into two categories: contact and contactless. The contact type temperature measurement utilizes good contact between sensors such as a thermocouple and a measured object, realizes temperature measurement by utilizing the property change of temperature-sensitive materials caused by temperature, hardly changes the temperature of the object when the sensors are added, but requires that the measured temperature does not exceed the upper limit temperature which can be borne by the sensors, the upper limit of the temperature measurement of the thermocouple is generally not more than 2000 ℃, and in addition, the contact temperature measurement generally limits the temperature measurement of severe environments such as corrosion, impact and the like. The non-contact type temperature measurement comprises infrared temperature measurement, colorimetric temperature measurement and multispectral radiation temperature measurement. The infrared temperature measurement adopts a thermal imager imaging mode to realize temperature measurement, can realize distribution measurement of a temperature field, takes products of FLIRE company in America as a representative, the upper limit of the temperature measurement is not more than 3500 ℃, the imaging time is in the magnitude of ms, the measurement application of the technology in an ultrahigh speed and ultrahigh temperature field is limited, the infrared temperature measurement requires to know the spectral emissivity of a measured object, and the spectral emissivity of the measured object is not only related to the material of the measured object, but also related to the temperature, so the infrared temperature measurement precision is limited in practical application. The colorimetric temperature measurement realizes temperature measurement by measuring spectral radiation of two wavelengths, eliminates the influence of emissivity through adjacent wavelength radiation energy ratio by utilizing the approximately equal relation of adjacent wavelength emissivity of a measured object, and realizes temperature measurement with higher precision in a certain temperature range, but the technology is not suitable for temperature measurement with a large dynamic range.
The multispectral temperature measurement realizes the measurement of a plurality of spectral radiations simultaneously, and the simultaneous measurement of the temperature and the spectral emissivity of a target to be measured is realized by combining the Planck radiation law according to the measured spectral information. However, the existing multispectral measurement technology mostly adopts a grating spectrometer and a Fourier transform spectrometer, so that the measurement speed is limited, and the imaging measurement cannot be realized. Therefore, multispectral thermometry has limited application in temperature field distribution measurements.
Disclosure of Invention
Aiming at the technical problems, the invention provides the temperature field space-time distribution measuring device based on the staring type snapshot spectral imaging, which has the advantages of accurate measurement, wide application range and high safety.
In order to solve the technical problems, the invention adopts the technical scheme that:
temperature field space-time distribution measuring device based on gazing type snapshot spectral imaging, including light imaging collection unit, spectrum imaging unit, signal sampling unit and computer PC end, be provided with the spectrum imaging unit on the light path direction of light imaging collection unit, spectrum imaging unit electric connection has signal sampling unit, signal sampling unit electric connection has computer PC end.
The spectral imaging unit comprises a spectral imaging coding module and a detector array, and the detector array is arranged in the light path direction of the spectral imaging coding module.
The light imaging collection unit adopts a reflection type telescope structure, a Cassegrain telescope or a Newton telescope.
The spectrum imaging coding module adopts a silicon-based flat photonic crystal, and the silicon-based flat photonic crystal performs target radiation spectrum pseudorandom transmittance coding modulation through cross holes, round holes or triangular holes of the flat photonic crystal air holes.
The detector array adopts a CMOS detector array, the wavelength range of the detector array is 400-900nm, the pixel size of the detector array is 5.6 mu m, the pixels of the detector array are not less than 2048 multiplied by 1024, and the frame frequency of the detector array is not less than 100 fps.
The spectral resolution of the spectral imaging unit is larger than 5nm, the number of spectral channels of the spectral imaging unit is N which is larger than or equal to 100, the actual compression measurement frequency M of the spectral imaging unit is smaller than or equal to 25, the compression ratio N of the spectral imaging unit is larger than or equal to 4:1, and each spectral measurement pixel of the spectral imaging unit is prepared by adopting a 5 x 5 silicon-based flat photonic crystal.
A measuring method of a temperature field space-time distribution measuring device based on staring type snapshot spectral imaging comprises the following steps:
s1, converging and imaging the radiation light of the target to be detected on the spectral imaging unit through the light imaging collection unit;
s2, the spectrum imaging unit converts the optical signal into an electric signal, and the detector array performs pseudo-random imaging detection on the imaging optical signal subjected to pseudo-random transmittance coding by the spectrum imaging coding module;
s3, the electric signal is collected and transmitted to a PC end of a computer through a signal sampling unit;
and S4, processing the spectral data and the image data in the PC terminal, and further performing inversion to obtain the spatial distribution and the time distribution information of the temperature field.
The compressed measurement signals obtained by the spectral imaging unit are restored at the PC end to obtain spectral signals of the target to be measured, at least 100 equations related to spectral emissivity and temperature are sequentially obtained according to the spectral signals obtained on each spectral measurement pixel by combining the Planck's radiation law, and linear fitting calculation is performed on the measurement data by using a linear least square method to obtain the spectral emissivity and the temperature of the target to be measured.
And the temperature obtained by each spectral measurement pixel is two-dimensionally arranged according to the position of the spectral measurement pixel to obtain the spatial distribution measurement of the temperature field, and the temperature field obtained at each moment is restored in the time dimension to master the temperature change rule of the target to be measured so as to obtain the time distribution measurement of the temperature field.
The method for sequentially obtaining at least 100 correlation equations with the spectral emissivity and the temperature comprises the following steps: light intensity V of target to be measured(λ)When the light is detected by the detector array, the transmittance regulation matrix Q of the spectral imaging coding modulei(λ)The regulation is carried out, the regulation spectrum signal of the ith time is detected by a detector array to form a compression measurement signal siIs composed of
si=∫V(λ)Qi(λ)η(λ)dλ i=1,2,…,M
Eta of(λ)Representing the spectral response conversion coefficient of the detector, eta(λ)The method is determined by the quantum efficiency of the detector array and the integral coefficient of a peripheral circuit of the detector array, and is obtained by performing wavelength radiation calibration on the detector array;
compressed measurement signal s realized by the spectral imaging unitiIs sampled by signalThe meta-sampling being a discrete digital quantity SiDispersing the light intensity signal of the target to be measured into N-dimensional vector representation, setting the number of corresponding spectral channels as N, andsatisfy the requirement ofMultiple compression of the measurement signal SiSatisfies Si∈RM×1SaidThe above-mentionedRepresenting a transmittance regulation compression measurement matrix;
the spectrum signal sparse representation is I ═ Ψ alpha, alpha is a K-sparse vector after the spectrum signal sparse representation, the K-sparse vector comprises K nonzero elements, and K is<<N; the Ψ is a sparse matrix; according to the theoretical framework of compressed sensing, willRewriting into S ═ QV ═ Q Ψ α ═ a α, where a ═ Q Ψ represents a sensing matrix, and is determined by both the compressed measurement matrix Q and the sparse matrix Ψ; then, the inverse problem of S ═ QV ═ Q Ψ α ═ a α is solved:
by finding optimal spectral sparse signalsThen the spectral signal to be measured is correctly recovered from the sparse signalTemperature measurement is carried out according to the spectral signals obtained on each spectral measurement pixel based on Planck's radiation lawInversion, according to planck's law of radiation, the formula for the spectral radiation temperature at absolute temperature T is:
l (lambda, T) is the radiance of the object, lambda is the wavelength, and T is the absolute temperature; the epsilon (lambda, T) is the spectral emissivity of the object; said C is1=3.7415×108W·μm4·m-2(ii) a Said C is2=1.43879×104μm4·K;
In the short wave band C2/λT>>1, now approximately replacing the planck formula with the wien formula:
the spectral radiant intensity output signal of the ith channel of the detector array is:
Vi=τ(λi)S(λi)L(λit), said τ (λ)i) Is spectral transmittance, S (λ)i) Is the detector sensitivity;
the output signal of the ith channel is then expressed as:
the output signal at the ith channel of the short band is represented as:
and the solution of the target true temperature requires the construction of the radiance temperature TiThe relationship with the target absolute temperature T is expressed by the formula:
establishing a 100-wavelength spectral emissivity model expressed as: epsilon (lambda, T) is exp (alpha)0+α1λ+α2λ2+…+αmλm) Said α is0,α1,α2,…,αmAnd if the measured data is constant, obtaining 100 equations in total, and performing linear fitting calculation on the measured data by using a linear least square method to obtain the spectral emissivity and the temperature of the target to be measured.
Compared with the prior art, the invention has the following beneficial effects:
1. the invention is realized by adopting a spectral imaging technology, the spectral imaging can simultaneously acquire the image information and the spectral information of the target to be measured, the image information can realize the measurement of the spatial distribution of the temperature field, and the spectral information can realize the accurate measurement of the temperature of the target corresponding to each pixel.
2. The invention adopts a multispectral radiation temperature measurement mode to realize temperature measurement, can simultaneously realize the spectral emissivity and the temperature measurement of a target to be measured, does not need to know the spectral emissivity, has accurate temperature measurement, and can realize temperature measurement in a larger dynamic range by acquiring spectral information in a wider range.
3. The invention adopts a staring type snapshot spectrum imaging mode, has high spectrum image information acquisition speed, can realize real-time continuous temperature measurement of the target to be measured, and obtains the time-varying rule of each temperature field on the basis of the spatial distribution measurement of the temperature fields.
4. The invention adopts staring type snapshot spectral imaging to realize temperature space-time distribution measurement, is a non-contact spectral imaging measurement method, adopts a telescopic imaging lens to image a temperature field on a spectral imaging acquisition module, is a non-contact remote measurement method, and can avoid the adverse environmental influences of explosion, combustion field impact, corrosion, vibration and the like.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described below. It should be apparent that the drawings in the following description are merely exemplary, and that other embodiments can be derived from the drawings provided by those of ordinary skill in the art without inventive effort.
The structures, ratios, sizes, and the like shown in the present specification are only used for matching with the contents disclosed in the specification, so that those skilled in the art will understand and read the present invention, and do not limit the conditions for implementing the present invention, so that the present invention has no technical essence, and any modifications of the structures, changes of the ratio relationships, or adjustments of the sizes, should still fall within the scope covered by the technical contents disclosed in the present invention without affecting the efficacy and the achievable purpose of the present invention.
FIG. 1 is a schematic structural view of the present invention;
FIG. 2 is a schematic diagram of a spectral imaging unit according to the present invention.
Wherein: the system comprises a light imaging collection unit 1, a spectrum imaging unit 2, a signal sampling unit 3, a computer PC (personal computer) terminal 4, a spectrum imaging coding module 2-1 and a detector array 2-2.
Detailed Description
To make the objects, technical solutions and advantages of the embodiments of the present invention clearer and more complete, the technical solutions in the embodiments of the present invention will be described below, obviously, the described embodiments are only a part of the embodiments of the present application, but not all embodiments, and the description is only for further explaining the features and advantages of the present invention, and not for limiting the claims of the present invention; all other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present application.
The following detailed description of embodiments of the present invention is provided in connection with the accompanying drawings and examples. The following examples are intended to illustrate the invention but are not intended to limit the scope of the invention.
Throughout the description of the present application, it is to be noted that, unless otherwise expressly specified or limited, the terms "mounted," "connected," and "connected" are to be construed broadly, e.g., as meaning either a fixed connection, a removable connection, or an integral connection; can be mechanically or electrically connected; they may be connected directly or indirectly through intervening media, or they may be interconnected between two elements. The specific meaning of the above terms in this application will be understood to be a specific case for those of ordinary skill in the art.
As shown in fig. 1, the radiation light of the target to be measured is converged and imaged on the spectral imaging unit 2 through the light imaging collecting unit 1, the spectral imaging unit 2 converts the light signal into an electrical signal, the electrical signal is collected and transmitted to the PC terminal 4 of the computer through the signal sampling unit 3, and the data processing is completed in the computer.
In this embodiment, the optical imaging collection unit 1 works in the visible/near infrared band, non-contact telemetry is realized by adopting a reflective telescope structure, and the optical imaging collection unit 1 is a cassegrain telescope or a newton telescope. The light radiated by the object to be measured is collected and focused to be imaged on the spectral imaging unit 2.
As shown in the schematic diagram of the spectral imaging unit in FIG. 2, the imaging unit includes a spectral imaging encoding module 2-1 and a detector array 2-2. The spectrum coding module 2-1 adopts a silicon-based flat photonic crystal design and adopts flat photonic crystal air pores such as cross holes, round holes, triangular holes and the like to perform target radiation spectrum pseudorandom transmittance coding modulation; the detector array 2-2 selects a CMOS detector array to carry out pseudo-random imaging detection on the imaging optical signal which is subjected to pseudo-random transmittance coding by the spectrum coding module 2-1.
The spectral imaging shown in fig. 2 is that the spectral imaging unit 2 has two-dimensional spatial information (x, y) and one-dimensional spectral information on the object to be measured. The incident light of each imaging unit of the target is incident to the detector array 2-2 through the spectral coding module 2-1 of the spectral imaging unit 2 in the whole spectral range. Light intensity V of target to be measured(λ)When detected by the detector, the transmissivity of the coding plate regulates and controls the matrix Qi(λ)The regulation is carried out, the regulation spectral signal of the ith time is detected by a detector to form a compression measurement signal siCan be described as
si=∫V(λ)Qi(λ)η(λ)dλ i=1,2,…,M
Wherein eta(λ)The spectral response conversion coefficient of the detector is mainly determined by the quantum efficiency of the detector and the integral coefficient of the peripheral circuit of the detector, and is generally obtained by carrying out wavelength radiometric calibration on the detector.
The compressed measuring signal realized by the spectral imaging unit 2 is sampled into discrete digital quantity S by the signal sampling unit 3i. Dispersing the light intensity signal of the target to be measured into N-dimensional vector representation, setting the number of corresponding spectral channels as N, andsatisfy the requirement ofMultiple compression measurement of spectral signal SiSatisfies Si∈RM×1. Using these matrices, the measurement process can be described as:
representing a transmittance adjustment compression measurement matrix. In the embodiment, the effective working wavelength of the detector array 2-2 is selected from the CMOS detector array and works at 400-900nm, the pixel size is 5.6 μm, the pixel is not less than 2048 × 1024, and the frame frequency is not less than 100 fps. The spectral resolution of the whole spectral imaging unit 2 is better than 5nm, the number N of spectral channels is set to be more than or equal to 100, the compression ratio N: M is more than or equal to 4:1, and each spectral measurement pixel of the spectral coding module 2-1 is prepared by adopting a 5 multiplied by 5 silicon-based flat plate photonic crystal. The spectral signal sparseness is represented as I ═ Ψ α, where α is a K-sparse vector (containing K non-zero elements, where K is the K-sparse vector after the spectral signal sparseness<<N), Ψ is a sparse matrix. According to the theoretical framework of compressed sensing, the above formula can be rewritten as
S=QV=QΨα=Aα
In the above formula, a ═ Q Ψ represents a sensing matrix, and is determined by both the compressed measurement matrix Q and the sparse matrix Ψ. The inverse problem of the above equation is then solved.
By finding optimal spectral sparse signalsThen the spectral signal to be measured is correctly recovered from the sparse signal
And carrying out temperature inversion based on the Planck's radiation law according to the spectrum signals obtained on each spectrum measurement pixel. Planck's law describes the distribution of the radiance of an object at different temperatures as a function of wavelength.
According to planck's radiation law, the spectral radiation temperature of the absolute temperature T can be represented by the formula:
wherein L (λ, T) is the radiance (W.m) of the object-2·μm-1sr-1) λ is the wavelength (μm), and T is the absolute temperature (K); ε (λ, T) is the spectral emissivity of the object; c1=3.7415×108W·μm4·m-2;C2=1.43879×104μm4·K。
In the short wave band C2/λT>>1, at this time, the Venn formula can be used to approximate the Planck formula
In the embodiment, the effective working wavelength of the CMOS detector array is selected to work at 400-900nm for the detector array 2-2, the spectral resolution of the whole spectral imaging unit 2 is better than 5nm, and the number N of spectral channels is set to be more than or equal to 100. Wherein the spectral radiant intensity output signal of the ith channel is:
Vi=τ(λi)S(λi)L(λi,T)
in the formula, tau (lambda)i) Is spectral transmittance, S (λ)i) For detector sensitivity
The output signal of the ith channel is then expressed as:
in the short band:
and the solving of the target true temperature requires the construction of the brightness temperature TiThe relationship with the target absolute temperature T,
expressed by the formula:
establishing a 100-wavelength spectral emissivity model expressed as: epsilon (lambda, T) is exp (alpha)0+α1λ+α2λ2+…+αmλm) In the formula, wherein alpha0,α1,α2,…,αmIs a constant
At the moment, 100 equations can be obtained in total, and linear fitting calculation is carried out on the measurement data by using a linear least square method, so that the spectral emissivity and the temperature of the target to be measured can be obtained. The temperature obtained by each spectrum pixel is arranged in two dimensions according to the position of the imaging pixel, so that the temperature field space distribution measurement can be completed, the temperature field obtained at each moment is restored in the time dimension, so that the temperature change rule of the target to be measured can be mastered, and the time distribution measurement of the temperature field can be obtained.
Although only the preferred embodiments of the present invention have been described in detail, the present invention is not limited to the above embodiments, and various changes can be made without departing from the spirit of the present invention within the knowledge of those skilled in the art, and all changes are encompassed in the scope of the present invention.
Claims (10)
1. Temperature field space-time distribution measuring device based on gazing type snapshot spectral imaging is characterized in that: including light imaging collection unit (1), spectrum imaging unit (2), signal sampling unit (3) and computer PC end (4), be provided with spectrum imaging unit (2) on the light path direction of light imaging collection unit (1), spectrum imaging unit (2) electric connection has signal sampling unit (3), signal sampling unit (3) electric connection has computer PC end (4).
2. The device for measuring the temporal and spatial distribution of the temperature field based on the staring type snapshot type spectral imaging as claimed in claim 1, wherein: the spectral imaging unit (2) comprises a spectral imaging coding module (2-1) and a detector array (2-2), and the detector array (2-2) is arranged in the optical path direction of the spectral imaging coding module (2-1).
3. The device for measuring the temporal and spatial distribution of the temperature field based on the staring type snapshot type spectral imaging as claimed in claim 1, wherein: the light imaging collection unit (1) adopts a reflection type telescope structure, a Cassegrain telescope or a Newton telescope.
4. The device for measuring the temporal and spatial distribution of the temperature field based on the staring type snapshot type spectral imaging as claimed in claim 2, wherein: the spectrum imaging coding module (2-1) adopts a silicon-based flat photonic crystal, and the silicon-based flat photonic crystal performs target radiation spectrum pseudorandom transmittance coding modulation through cross hole, round hole or triangular hole flat photonic crystal air pores.
5. The device for measuring the temporal and spatial distribution of the temperature field based on the staring type snapshot type spectral imaging as claimed in claim 2, wherein: the detector array (2-2) adopts a CMOS detector array, the wavelength range of the detector array (2-2) is 400-900nm, the pixel size of the detector array (2-2) is 5.6 mu m, the pixels of the detector array (2-2) are not less than 2048 multiplied by 1024, and the frame rate of the detector array (2-2) is not less than 100 fps.
6. The device for measuring the temporal and spatial distribution of the temperature field based on the staring type snapshot type spectral imaging as claimed in claim 1, wherein: the spectral resolution of the spectral imaging unit (2) is larger than 5nm, the number of spectral channels of the spectral imaging unit (2) is more than or equal to 100, the actual compression measurement frequency M of the spectral imaging unit (2) is less than or equal to 25, the compression ratio N: M of the spectral imaging unit (2) is more than or equal to 4:1, and each spectral measurement pixel of the spectral imaging unit (2) is prepared by adopting 5 x 5 silicon-based flat photonic crystals.
7. The measurement method of the temperature field space-time distribution measurement device based on the staring type snapshot type spectral imaging according to any one of claims 1 to 6, characterized in that: comprises the following steps:
s1, converging and imaging the radiation light of the target to be detected on the spectral imaging unit (2) through the light imaging collection unit (1);
s2, the spectrum imaging unit (2) converts the optical signals into electric signals, and the detector array (2-2) carries out pseudo-random imaging detection on the imaging optical signals subjected to pseudo-random transmittance coding by the spectrum imaging coding module (2-1);
s3, the electric signals are collected and transmitted to a PC (personal computer) end (4) through the signal sampling unit (3);
s4, processing the spectral data and the image data in the PC (4), and further performing inversion to obtain the spatial distribution and the time distribution information of the temperature field.
8. The measurement method of the temperature field space-time distribution measurement device based on the staring type snapshot type spectral imaging as claimed in claim 7, characterized in that: the compressed measurement signals obtained by the spectral imaging unit (2) are restored at a computer PC (personal computer) end (4) to obtain spectral signals of the target to be measured, at least 100 equations related to spectral emissivity and temperature are sequentially obtained according to the spectral signals obtained on each spectral measurement pixel by combining the Planck's radiation law, and linear fitting calculation is performed on the measurement data by using a linear least square method to obtain the spectral emissivity and the temperature of the target to be measured.
9. The measurement method of the temperature field space-time distribution measurement device based on the staring type snapshot type spectral imaging as claimed in claim 8, characterized in that: and the temperature obtained by each spectral measurement pixel is two-dimensionally arranged according to the position of the spectral measurement pixel to obtain the spatial distribution measurement of the temperature field, and the temperature field obtained at each moment is restored in the time dimension to master the temperature change rule of the target to be measured so as to obtain the time distribution measurement of the temperature field.
10. The measurement method of the temperature field space-time distribution measurement device based on the gaze-type snapshot spectral imaging as claimed in claim 9, wherein: the method for sequentially obtaining at least 100 correlation equations with the spectral emissivity and the temperature comprises the following steps: light intensity V of target to be measured(λ)When the light is detected by the detector array (2-2), the transmittance regulation matrix Q of the spectral imaging coding module (2-1)i(λ)The regulation is carried out, the regulation spectrum signal of the ith time is detected by a detector array (2-2) to form a compression measurement signal siIs composed of
si=∫V(λ)Qi(λ)η(λ)dλ i=1,2,…,M
Eta of(λ)Representing the spectral response conversion coefficient of the detector, eta(λ)The quantum efficiency of the detector array (2-2) and the integral coefficient of the peripheral circuit of the detector array (2-2) are used for determining, and the wavelength radiometric calibration is carried out on the detector array (2-2);
compressed measurement signal s realized by the spectral imaging unit (2)iIs sampled into discrete digital quantity S by a signal sampling unit (3)iDispersing the light intensity signal of the target to be measured into N-dimensional vector representation,the number of corresponding spectral channels is set to N, andsatisfy the requirement ofMultiple compression of the measurement signal SiSatisfies Si∈RM×1SaidThe above-mentionedRepresenting a transmittance regulation compression measurement matrix;
the spectrum signal sparse representation is I ═ Ψ alpha, alpha is a K-sparse vector after the spectrum signal sparse representation, the K-sparse vector comprises K nonzero elements, and K is<<N; the Ψ is a sparse matrix; according to the theoretical framework of compressed sensing, willRewriting into S ═ QV ═ Q Ψ α ═ a α, where a ═ Q Ψ represents a sensing matrix, and is determined by both the compressed measurement matrix Q and the sparse matrix Ψ; then, the inverse problem of S ═ QV ═ Q Ψ α ═ a α is solved:
by finding optimal spectral sparse signalsThen the spectral signal to be measured is correctly recovered from the sparse signalThe spectral signals obtained on each spectral measurement pixel are processed based on the Planck's radiation lawTemperature inversion, according to planck's radiation law, the formula of the spectral radiation temperature of absolute temperature T is:
l (lambda, T) is the radiance of the object, lambda is the wavelength, and T is the absolute temperature; the epsilon (lambda, T) is the spectral emissivity of the object; said C is1=3.7415×108W·μm4·m-2(ii) a Said C is2=1.43879×104μm4·K;
In the short wave band C2/λT>>1, now approximately replacing the planck formula with the wien formula:
the spectral radiation intensity output signal of the ith channel of the detector array (2-2) is as follows:
Vi=τ(λi)S(λi)L(λit), said τ (λ)i) Is spectral transmittance, S (λ)i) Is the detector sensitivity;
the output signal of the ith channel is then expressed as:
the output signal at the ith channel of the short band is represented as:
and the solution of the target true temperature requires the construction of the radiance temperature TiThe relationship with the target absolute temperature T is expressed by the formula:
establishing a 100-wavelength spectral emissivity model expressed as: epsilon (lambda, T) is exp (alpha)0+α1λ+α2λ2+…+αmλm) Said α is0,α1,α2,…,αmAnd if the data is constant, 100 equations are obtained in total, and linear fitting calculation is carried out on the measured data by using a linear least square method to obtain the spectral emissivity and the temperature of the target to be measured.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210147293.0A CN114543998A (en) | 2022-02-17 | 2022-02-17 | Temperature field space-time distribution measuring device based on staring type snapshot type spectral imaging |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210147293.0A CN114543998A (en) | 2022-02-17 | 2022-02-17 | Temperature field space-time distribution measuring device based on staring type snapshot type spectral imaging |
Publications (1)
Publication Number | Publication Date |
---|---|
CN114543998A true CN114543998A (en) | 2022-05-27 |
Family
ID=81675903
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210147293.0A Pending CN114543998A (en) | 2022-02-17 | 2022-02-17 | Temperature field space-time distribution measuring device based on staring type snapshot type spectral imaging |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114543998A (en) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110186566A (en) * | 2019-05-30 | 2019-08-30 | 上海交通大学 | Two-dimentional true temperature field imaging method and system based on the multispectral thermometric of light-field camera |
CN111458044A (en) * | 2020-05-28 | 2020-07-28 | 中国工程物理研究院流体物理研究所 | Transient temperature measuring device and method based on snapshot spectral imaging technology |
CN111664941A (en) * | 2020-07-08 | 2020-09-15 | 中北大学 | Compressed sensing spectrum measuring device and method for electro-optically regulating transmittance coding |
CN112345074A (en) * | 2020-11-11 | 2021-02-09 | 中北大学 | Chip-level satellite-borne hyperspectral imaging detector and spectral imaging method thereof |
-
2022
- 2022-02-17 CN CN202210147293.0A patent/CN114543998A/en active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110186566A (en) * | 2019-05-30 | 2019-08-30 | 上海交通大学 | Two-dimentional true temperature field imaging method and system based on the multispectral thermometric of light-field camera |
CN111458044A (en) * | 2020-05-28 | 2020-07-28 | 中国工程物理研究院流体物理研究所 | Transient temperature measuring device and method based on snapshot spectral imaging technology |
CN111664941A (en) * | 2020-07-08 | 2020-09-15 | 中北大学 | Compressed sensing spectrum measuring device and method for electro-optically regulating transmittance coding |
CN112345074A (en) * | 2020-11-11 | 2021-02-09 | 中北大学 | Chip-level satellite-borne hyperspectral imaging detector and spectral imaging method thereof |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102384841B (en) | Spectral responsivity test method for plane array detector | |
CN110487134B (en) | Device and method for measuring combustion speed and temperature of explosion flame | |
Chyzh et al. | Energy resolution of dual-channel opto-electronic surveillance system | |
CN110793632B (en) | High-speed high-precision spectrum video system and method for flame shooting | |
US4788428A (en) | Thermodynamics infrared imaging sensor | |
Runciman | Thermal imaging | |
Anderson et al. | The physical basis of current infrared remote-sensing techniques and the interpretation of data from aerial surveys | |
Bowers et al. | Unpolarized calibration and nonuniformity correction for long-wave infrared microgrid imaging polarimeters | |
CN114449079A (en) | High-temperature measuring device and method based on mobile phone camera | |
CN105806491A (en) | Three-wavelength two-dimensional temperature field measuring device and method | |
CN114543998A (en) | Temperature field space-time distribution measuring device based on staring type snapshot type spectral imaging | |
CN114894737B (en) | Spectral reflectivity reconstruction method based on infrared image | |
Méndez-Rial et al. | A high-speed MWIR uncooled multi-aperture snapshot spectral imager for IR surveillance and monitoring | |
RU2727349C1 (en) | Method of thermography of a remote object | |
Shaw et al. | Infrared cloud imager deployment at the north slope of Alaska during early 2002 | |
CN113418613A (en) | High-temperature transient measurement system and method based on multispectral colorimetry | |
CN201680912U (en) | High resolution thermal imaging system | |
Hinnrichs et al. | Comparison of QWIP to HgCdTe detectors for gas imaging | |
Soldani | Infrared signature: Theory and example of practical measurement methods | |
Goldberg et al. | Laboratory and field performance of megapixel QWIP focal plane arrays | |
Farley et al. | Radiometric calibration stability of the FIRST: a longwave infrared hyperspectral imaging sensor | |
RU2736094C1 (en) | True temperature pyrometer | |
Zhang et al. | Thermography of surface temperature of explosion fireball by visible and infrared band | |
US20220099491A1 (en) | Methods for spectral mapping | |
Yang et al. | Screening and testing methods of short-wave infrared area array detectors for synchronous monitoring atmospheric corrector |
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 |