CN108955936A - A method of non-uniform temperature distribution is carried out using multi-angle histogram information and is rebuild - Google Patents
A method of non-uniform temperature distribution is carried out using multi-angle histogram information and is rebuild Download PDFInfo
- Publication number
- CN108955936A CN108955936A CN201710354725.4A CN201710354725A CN108955936A CN 108955936 A CN108955936 A CN 108955936A CN 201710354725 A CN201710354725 A CN 201710354725A CN 108955936 A CN108955936 A CN 108955936A
- Authority
- CN
- China
- Prior art keywords
- matrix
- temperature
- histogram
- value
- rebuild
- 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.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01K—MEASURING TEMPERATURE; MEASURING QUANTITY OF HEAT; THERMALLY-SENSITIVE ELEMENTS NOT OTHERWISE PROVIDED FOR
- G01K11/00—Measuring temperature based upon physical or chemical changes not covered by groups G01K3/00, G01K5/00, G01K7/00 or G01K9/00
- G01K11/32—Measuring temperature based upon physical or chemical changes not covered by groups G01K3/00, G01K5/00, G01K7/00 or G01K9/00 using changes in transmittance, scattering or luminescence in optical fibres
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/001—Texturing; Colouring; Generation of texture or colour
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/20—Drawing from basic elements, e.g. lines or circles
- G06T11/206—Drawing of charts or graphs
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Radiation Pyrometers (AREA)
Abstract
The present invention provides a kind of method for carrying out non-uniform temperature distribution reconstruction using multi-angle histogram information.The following steps are included: arranging a plurality of measuring route on horizontal, vertical, left diagonal line and right diagonal line four direction around dut temperature region, and the histogram information in each direction each path is successively obtained, forms four histogram matrix；Four histogram matrix are combined using logic judgment and the algorithm of traversal, generate the probability matrix comprising multi-angle histogram all information；Temperature matrices are rebuild by analysis probability matrix, it is final to realize that two-dimension temperature distribution is rebuild.Effect of the invention is in the case where measurement angle is limited, realize that the distribution of non-uniform temperature field is rebuild using the measurement data of limited angle, the complexity of calculating can be effectively reduced while guaranteeing reconstruction quality in this method, it is easy to use, it is with a wide range of applications.
Description
Technical field
The present invention relates to a kind of non-uniform temperatures to be distributed method for reconstructing, and in particular to a kind of to use multi-angle histogram information
Carry out the method that non-uniform temperature distribution is rebuild.
Background technique
In many practical applications, due to effects such as heat transfer, cold boundary layer and gas diffusions, the distribution in temperature field is often
There are apparent inhomogeneities.For the measurement of non-uniform temperature field distribution, in addition to by light radiation thermometric and induced with laser
Fluorescence etc. directly obtains flame image by high speed camera and realizes that two-dimension temperature measurement is outer, can also be by the temperature field CT point
Cloth is rebuild, such as acoustics CT imaging technique and laser absorption spectrum CT imaging technique.
Laser absorption spectrum because of it without pretreatment, measurement fast response time, data it is accurate, can multi-parameter detect simultaneously etc.
Advantage is widely used in combustion process in the non-contact measurement of multi-parameter.Using TDLAS technology measurement Temperature Distribution without
Sample gas need to be acquired, laser and detector are directly installed on measure field, and measuring device itself and tested gas are completely isolated,
With fabulous environmental suitability.Generally, it can only obtain the average value of temperature in the measuring route using single-measurement angle
And the space distribution information of parameter can not be obtained；Few angle measurement data not can guarantee the precision of image reconstruction.In recent years, it utilizes
Multi-angle laser absorption light technology combination CT imaging technique realize non-uniform temperature field distribution measuring study achieve centainly into
Exhibition, such as Liu Chang in 2015 are published in the phase of volume 23 17 of " optics letter " (Optics Express), 22494-
Page 22511, entitled " gas temperature, concentration fast imaging sensor research based on fladellum TDLAS " (A Development
of a fan-beam TDLAS-based tomographic sensor for rapid imaging of temperature
And gas concentration) article；The above method needs to install multi-angle, multiple groups sensor around tested region,
Distribution reconstruction is carried out using the projection data of multi-angle.However, in practical applications, the factors such as site environment are limited to, it can not
Enough laser optical path and detector are arranged around tested region, meanwhile, a large amount of experimental data can also reduce temperature field
The operation burden and testing cost rebuild speed, increase system.Therefore, how temperature field is carried out using less measurement angle to divide
Quick, the Exact Reconstruction of cloth are the research emphasis of reconstruction of temperature field method.
In EO-1 hyperion laser absorption spectroscopy, it is disposed with a plurality of absorption line on single light path, can obtain more
Parameter information abundant, such as the histogram information of Temperature Distribution.For example, Trevor York in 2016 is published in " Applied Physics B "
The phase of volume 122 1 of (Applied Physics B), the 3rd (1-9) page, entitled " the methane air flame based on single channel TDLAS is non-
Uniform temperature field measurement " (Characterization of temperature non-uniformity over a
Premixed CH4-air flame based on line-of-sight TDLAS) article in, using EO-1 hyperion chromatography at
As technology, a plurality of absorption spectrum is arranged on single light path, obtains the temperature histogram information on the path.Single angle
For temperature histogram information is compared with temperature averages, richer data are contained, increase the information content of single angle.With
Toward the work for carrying out histogram reconstruction using EO-1 hyperion, it is intended to which inverting obtains the temperature histogram information on path, is not directed to
Two dimensional Distribution to temperature is rebuild.Present invention combination EO-1 hyperion laser absorption spectroscopy, using few angular histogram information into
Quick, the Exact Reconstruction of trip temperature field distribution.
Non-uniform temperature distribution, the water around dut temperature region are rebuild using the histogram information of four angles
A plurality of measuring route is arranged on flat, vertical, left diagonal line and right diagonal line four direction, and successively obtains every, each direction road
Histogram information on diameter forms four histogram matrix；Generate the probability square comprising four angular histogram all information
Battle array；Temperature matrices are rebuild by analysis probability matrix, it is final to realize that two-dimension temperature distribution is rebuild.
Summary of the invention
In order to realize point of non-uniform temperature field using the measurement data of limited angle in the case where measurement angle is limited
Cloth is rebuild, and the present invention provides a kind of method for carrying out non-uniform temperature distribution reconstruction using multi-angle histogram information, this method
The demand to measurement angle and data volume can be reduced, the complexity of calculating is effectively reduced while guaranteeing reconstruction quality.This
Technical solution used by inventing is as follows:
1. a kind of carry out the method that non-uniform temperature distribution is rebuild using multi-angle histogram information, this method includes following
Step:
Step 1 is arranged on horizontal, vertical, the left diagonal line and right diagonal line four direction around dut temperature region
A plurality of measuring route successively obtains the histogram information in each direction each path, forms four histogram matrix；Four sides
To temperature histogram matrix respectively by H, V, L and R indicate.
Step 2, the four direction histogram matrix H obtained by step 1, V, L, R generating probability matrix (PM,
Probability matrix), which shows that the element of temperature matrices corresponding position is the probability of a certain value.In PM
Element calculating process is as follows: pass through four row vectors of certain same position point in four direction histogram matrix H, V, L and R, into
Row Hadamard product, and regulation calculated value be greater than 1 when take 1, take zero when equal to zero.
Assuming that H ', V ', L ', R ' it is the fixed histogram matrix of temperature value, one temperature value of every determination, the histogram square
Battle array just will be updated once, and PM matrix updates therewith, and renewal process is as follows: four direction histogram matrix H-H ', V-V ', L-L '
And pass through four row vectors of certain same position point in R-R ', Hadamard product is carried out, calculated value takes 1 when being greater than 1, is equal to zero
When take zero.
Step 3, by H ', V ', L ', R ' is initialized as null matrix, generates PM matrix by formula (3).If in PM including monodrome
Row, then the element on temperature matrices corresponding position can determine, update H ', V ', L ', R '.Otherwise, Ergodic Matrices find row and most
Small row, and carried out branch operation and be converted to monodrome row.So circulation determines the temperature value of all positions undetermined, and will
Reconstructed value maps back original temperature range, constitutes final temperature matrices, completes non-homogeneous two-dimension temperature distribution and rebuilds.
Detailed description of the invention
Fig. 1 is the measuring route arrangement schematic diagram of horizontal, vertical, left diagonal line and right diagonal line four direction
Fig. 2 is the distribution algorithm for reconstructing flow chart based on histogram information
(a) original distribution and (b) reconstructed results of Fig. 3 temperature matrices 1
(a) original distribution and (b) reconstructed results of Fig. 4 temperature matrices 2
Specific embodiment
In the present embodiment, using the measurement of multiplex physicotherapy laser absorption spectroscopy techniques tested region temperature histogram, with quilt
It is divided into for the rectangular temperature region of 2 × 2 mesh point, the implementation process of this method is described in detail.Then with two components
It is not divided into for the Temperature Distribution of 14 × 14 and 16 × 16 mesh point, demonstrates the validity of algorithm.Utilize described one
Kind carries out the method that non-uniform temperature distribution is rebuild using multi-angle histogram information and rebuilds the Temperature Distribution in the region.With biography
The image reconstruction algorithm of system is compared, described that the method that non-uniform temperature distribution is rebuild is carried out using multi-angle histogram information, can
To reduce the demand to measurement angle and data volume, and reduce the complexity of calculating while guaranteeing reconstruction quality.
The present invention will be further described below with reference to the drawings:
Step 1 arranges the laser road of horizontal, vertical, left diagonal line and right diagonal line four direction around tested region
Diameter, the measurement of tested region temperature histogram is realized using multiplex physicotherapy laser absorption spectroscopy techniques, and measuring system includes 1 adjustable
The laser of humorous laser, laser output is divided into multiple laser by optical fiber splitter, includes a plurality of on single laser light beam
Absorption line, the collimated mirror of laser pass through tested region and are received by a detector after collimating.Fig. 1 is how to obtain histogram matrix
Signal.In figure, temperature region is divided into 2 × 2 mesh point, and temperature value is discretized as 3 equal portions, and is mapped to 0,1 He
2；Actual temperature matrix T is expressed as follows,
It is 0,1 and 2 institute that temperature value in four direction each path, which can be obtained, by multiplex physicotherapy laser absorption spectroscopy techniques
The normalization percentage accounted for, i.e. temperature histogram information.After obtaining histogram information, histogram matrix is obtained by following steps:
Arrow represents measuring route and direction in figure, and the number of arrow is equal to the line number i of histogram matrix, and the number of temperature is equal to straight
Square figure matrix column number j.Histogram matrix under four angles is indicated by H, V, L and R respectively.Element h in matrix HijValue etc.
In the number of temperature of horizontal histogram the i-th paths upper value equal to (j-1), similarly for V, L and R.The histogram of four direction
Figure matrix can be write:
Step 2, the four direction histogram matrix H obtained by step 1, V, L, R generating probability matrix (PM,
Probability matrix), which shows that the element of temperature matrices corresponding position is the probability of a certain value.
Element p in PM matrixijElement t in homography Tij, and indicate its possible value.If pijM is classified as 1,
Then indicate tijValue be equally likely to the corresponding temperature value of m column；If it is 0, tijValue can not for m column corresponding to
Temperature value.Vector pijBy passing point t in four histogram matrixijRow vector calculate Hadamard product obtain, such as formula
(10) shown in, and regulation calculated value be greater than 1 when take 1, take zero when equal to zero.
p11=h1.*v1.*l1.*r2=[h11v11l11r21 h12v12l12r22 h13v13l13r23] (7)
For example, if pijEqual to [0 0 1], then tijOnly possible value is 2.Only have a column to be equal to 1 row in matrix PM
Vector is referred to as monodrome row.
Assuming that H ', V ', L ', R ' it is the fixed histogram matrix of temperature value, one temperature value of every determination, the histogram square
Battle array just will be updated primary.For example, if having determined that t12Value be 2, t21Value be 0, then
Histogram matrix is updated to H-H', V-V', L-L', R-R', to wherein by four rows of certain same position point to
Amount, carry out Hadamard product, and regulation calculated value be greater than 1 when take 1, zero is taken when equal to zero, then PM matrix update formula is as follows,
Wherein, full zero row indicates that corresponding position element has determined that.
Step 3, by H ', V ', L ', R ' is initialized as null matrix, generates PM matrix by formula (3).If in PM including monodrome
Row, then the element on temperature matrices corresponding position can determine, update H ', V ', L ', R '.Otherwise, Ergodic Matrices, which are found, removes complete zero
Row layman and the smallest row, and carried out branch operation and be converted to monodrome row.For example, the value of certain row is [0 1 1], goes and be
2 and be the minimum value in all rows, then the row is chosen, [0 1 0] and [0 0 1] progress branch operation is set to.Such as
This circulation, determines the temperature value of all positions undetermined, and reconstructed value is mapped back original temperature range, constitutes final temperature square
Battle array is completed non-homogeneous two-dimension temperature distribution and is rebuild, and the flow chart of algorithm is as shown in Figure 2.
This example carrys out the validity of verification algorithm using two groups of different Temperature Distributions.Two groups of temperature datas are respectively 14 ×
14 and 16 × 16 temperature matrices, temperature range are mapped to 1-10 after being divided into 10 parts, and thus generate the straight of four direction
Square figure matrix.Fig. 3 is the original distribution (a) and (b) reconstructed results of temperature matrices 1；Fig. 4 is the original distribution (a) of temperature matrices 2
(b) reconstructed results.Reconstructed results to be assessed using PSNR index, calculation formula is as follows,
Wherein, N is the dimension of matrix, and R (x, y) and O (x, y) are the reconstructed value and original value of position (x, y), R respectivelymax
It is the maximum value of R, the PSNR that temperature matrices 1 and 2 rebuild front and back is respectively 51.5dB and 49.9dB.PSNR reflects quantization error
Size, value depend on temperature range and quantization number, distortion occur in departure process, and by histogram carry out
Reconstruction process is lossless.
In conclusion compared with traditional method for reconstructing, it is non-homogeneous using being carried out using multi-angle histogram information
The method that Temperature Distribution is rebuild can reduce the demand to measurement angle and data volume, effectively while guaranteeing reconstruction quality
Reduce the complexity calculated.Above to the description of the present invention and embodiments thereof, it is not limited to which this, is only this shown in attached drawing
One of embodiment of invention.Without departing from the spirit of the invention, it without creating designs and the technical side
Case similar structure or embodiment, the category scope of the present invention.
Claims (1)
1. a kind of carry out the method that non-uniform temperature distribution is rebuild using multi-angle histogram information, this method includes following step
It is rapid:
Step 1 is arranged on horizontal, vertical, the left diagonal line and right diagonal line four direction around dut temperature region a plurality of
Measuring route successively obtains the histogram information in each direction each path, forms four histogram matrix；Four direction
Temperature histogram matrix is indicated by H, V, L and R respectively.
Step 2, the four direction histogram matrix H obtained by step 1, V, L, R generating probability matrix (PM, probability
Matrix), which shows that the element of temperature matrices corresponding position is the probability of a certain value.Element calculating process in PM
It is as follows: to pass through four row vectors of certain same position point in four direction histogram matrix H, V, L and R, carry out Hadamard and multiply
Product, calculated value take 1 when being greater than 1, take zero when equal to zero.
Assuming that H ', V ', L ', R ' it is the fixed histogram matrix of temperature value, as soon as every determination temperature value, the histogram matrix
Will be updated once, PM matrix updates therewith, and renewal process is as follows: four direction histogram matrix H-H ', V-V ', L-L ' and R-R '
In by four row vectors of certain same position point, carry out Hadamard product, and when regulation calculated value is greater than 1 takes 1, is equal to zero
When take zero.
Step 3, by H ', V ', L ', R ' is initialized as null matrix, generates PM matrix by formula (3).If in PM including monodrome row,
Then the element on temperature matrices corresponding position can determine, update H ', V ', L ', R '.Otherwise, Ergodic Matrices are found capable and minimum
Row, and carried out branch operation and be converted to monodrome row.So circulation, determines the temperature value of all positions undetermined, and will weigh
Built-in value maps back original temperature range, constitutes final temperature matrices, completes non-homogeneous two-dimension temperature distribution and rebuilds.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710354725.4A CN108955936B (en) | 2017-05-19 | 2017-05-19 | Method for carrying out non-uniform temperature distribution reconstruction by adopting multi-angle histogram information |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710354725.4A CN108955936B (en) | 2017-05-19 | 2017-05-19 | Method for carrying out non-uniform temperature distribution reconstruction by adopting multi-angle histogram information |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108955936A true CN108955936A (en) | 2018-12-07 |
| CN108955936B CN108955936B (en) | 2019-12-24 |
Family
ID=64462032
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710354725.4A Active CN108955936B (en) | 2017-05-19 | 2017-05-19 | Method for carrying out non-uniform temperature distribution reconstruction by adopting multi-angle histogram information |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108955936B (en) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20080118783A1 (en) * | 2006-05-30 | 2008-05-22 | Cetegen Baki M | Fiber optic based in-situ diagnostics for PEM fuel cells |
| KR101159215B1 (en) * | 2011-12-22 | 2012-06-25 | 한국생산기술연구원 | Optics device for measuring gas temperature and density |
| CN203216634U (en) * | 2013-03-17 | 2013-09-25 | 浙江师范大学 | High temperature measurement and temperature field reconstruction device based on turnable diode laser absorption spectroscopy |
| CN103884663A (en) * | 2014-03-19 | 2014-06-25 | 中国人民解放军装备学院 | Two-dimensional reconstruction light distribution optimization method based on laser absorption spectrum technology |
-
2017
- 2017-05-19 CN CN201710354725.4A patent/CN108955936B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20080118783A1 (en) * | 2006-05-30 | 2008-05-22 | Cetegen Baki M | Fiber optic based in-situ diagnostics for PEM fuel cells |
| KR101159215B1 (en) * | 2011-12-22 | 2012-06-25 | 한국생산기술연구원 | Optics device for measuring gas temperature and density |
| CN203216634U (en) * | 2013-03-17 | 2013-09-25 | 浙江师范大学 | High temperature measurement and temperature field reconstruction device based on turnable diode laser absorption spectroscopy |
| CN103884663A (en) * | 2014-03-19 | 2014-06-25 | 中国人民解放军装备学院 | Two-dimensional reconstruction light distribution optimization method based on laser absorption spectrum technology |
Non-Patent Citations (2)
| Title |
|---|
| DOO-WON CHOI等: "Performance Improvements in Temperature Reconstructions of 2-D Tunable Diode Laser Absorption Spectroscopy (TDLAS)", 《JOURNAL OF THERMAL SCIENCE》 * |
| PEI-JIN LIU 等: "TDLAS for measurement of temperature in combustion environment", 《PROC. OF SPIE》 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108955936B (en) | 2019-12-24 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Nicolas et al. | A direct approach for instantaneous 3D density field reconstruction from background-oriented schlieren (BOS) measurements | |
| Venkatakrishnan et al. | Density measurements using the background oriented schlieren technique | |
| Roberti et al. | Microwave radiative transfer through horizontally inhomogeneous precipitating clouds | |
| Song et al. | Algebraic tomographic reconstruction of two-dimensional gas temperature based on tunable diode laser absorption spectroscopy | |
| CN108195419B (en) | Novel fluid multi-field observation system and method | |
| Ni et al. | Study on the detection of three-dimensional soot temperature and volume fraction fields of a laminar flame by multispectral imaging system | |
| CN106815878A (en) | A kind of image rebuilding method based on the multispectral tomography of various visual angles | |
| CN101680833A (en) | Binning and tomography for high spatial resolution temperature and species concentration measurements | |
| CN108627272A (en) | A kind of two-dimension temperature distribution method for reconstructing based on four angle laser absorption spectrums | |
| CN110514622B (en) | Two-dimensional temperature and concentration reconstruction system and method based on histogram information | |
| CN109061220B (en) | Airflow two-dimensional velocity distribution measuring method based on laser absorption spectrum tomography technology | |
| Altenhoff et al. | Soot aggregate sizing in an extended premixed flame by high-resolution two-dimensional multi-angle light scattering (2D-MALS) | |
| CN109100044A (en) | Method for reconstructing is fitted based on the multispectral gas temperature probability density distribution in monochromatic light road | |
| Hao et al. | Displacement field denoising for high-temperature digital image correlation using principal component analysis | |
| Ghandehari et al. | Surface temperatures in New York City: Geospatial data enables the accurate prediction of radiative heat transfer | |
| Démoulin et al. | The 3D geometry of active region upflows deduced from their limb-to-limb evolution | |
| KR101886540B1 (en) | Remote Sensed Image Simulation System and Method for feasibility test of observed image data | |
| Mariani et al. | A comparative study on the use of calibrated and rainbow schlieren techniques in axisymmetric supersonic jets | |
| Bichal | Development of 3D background oriented schlieren with a plenoptic camera | |
| CN107247035B (en) | A method of based on the height of observation rapid survey Liquid Diffusion Coefficient such as wick-containing cylindrical lens | |
| CN108955936A (en) | A method of non-uniform temperature distribution is carried out using multi-angle histogram information and is rebuild | |
| Li et al. | A real-time aircraft infrared imaging simulation platform | |
| CN107870159B (en) | Gas concentration two-dimensional reconstruction method for tunable semiconductor laser absorption spectrum | |
| CN109100045A (en) | Based on the multispectral gas temperature probability density distribution method for reconstructing in monochromatic light road | |
| Zou et al. | Developmental Trends in the Application and Measurement of the Bidirectional Reflection Distribution Function |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |