CN109870424A - Hartmann's ray tracing method based on colored three steps transposition technology - Google Patents
Hartmann's ray tracing method based on colored three steps transposition technology Download PDFInfo
- Publication number
- CN109870424A CN109870424A CN201910163686.9A CN201910163686A CN109870424A CN 109870424 A CN109870424 A CN 109870424A CN 201910163686 A CN201910163686 A CN 201910163686A CN 109870424 A CN109870424 A CN 109870424A
- Authority
- CN
- China
- Prior art keywords
- formula
- light
- phase
- measurement
- colored
- 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
- 238000000034 method Methods 0.000 title claims abstract description 67
- 238000005516 engineering process Methods 0.000 title claims abstract description 14
- 230000017105 transposition Effects 0.000 title claims abstract description 12
- 238000005259 measurement Methods 0.000 claims abstract description 114
- 230000003287 optical effect Effects 0.000 claims abstract description 51
- 230000001052 transient effect Effects 0.000 claims abstract description 7
- 239000011159 matrix material Substances 0.000 claims description 22
- 238000006073 displacement reaction Methods 0.000 claims description 13
- 238000003384 imaging method Methods 0.000 claims description 10
- 238000006243 chemical reaction Methods 0.000 claims description 5
- 230000007704 transition Effects 0.000 claims description 3
- 238000013461 design Methods 0.000 abstract description 7
- 239000007788 liquid Substances 0.000 description 27
- 239000007789 gas Substances 0.000 description 13
- 238000009826 distribution Methods 0.000 description 8
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 8
- 238000002156 mixing Methods 0.000 description 8
- 239000003795 chemical substances by application Substances 0.000 description 7
- 230000008859 change Effects 0.000 description 5
- 238000002485 combustion reaction Methods 0.000 description 5
- 230000000694 effects Effects 0.000 description 5
- 238000005305 interferometry Methods 0.000 description 5
- 238000004364 calculation method Methods 0.000 description 4
- 238000010438 heat treatment Methods 0.000 description 4
- 239000003345 natural gas Substances 0.000 description 4
- 235000019476 oil-water mixture Nutrition 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 238000000691 measurement method Methods 0.000 description 3
- 235000019198 oils Nutrition 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 239000000126 substance Substances 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 2
- 230000001427 coherent effect Effects 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- 125000001475 halogen functional group Chemical group 0.000 description 2
- 238000005210 holographic interferometry Methods 0.000 description 2
- 230000033001 locomotion Effects 0.000 description 2
- 239000003921 oil Substances 0.000 description 2
- 239000003208 petroleum Substances 0.000 description 2
- 230000009290 primary effect Effects 0.000 description 2
- 230000001105 regulatory effect Effects 0.000 description 2
- 238000010008 shearing Methods 0.000 description 2
- 238000004611 spectroscopical analysis Methods 0.000 description 2
- 238000001228 spectrum Methods 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 101100117236 Drosophila melanogaster speck gene Proteins 0.000 description 1
- 238000000342 Monte Carlo simulation Methods 0.000 description 1
- 238000005452 bending Methods 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 239000000567 combustion gas Substances 0.000 description 1
- 230000021615 conjugation Effects 0.000 description 1
- 238000011217 control strategy Methods 0.000 description 1
- 230000001276 controlling effect Effects 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 239000013078 crystal Substances 0.000 description 1
- 238000005520 cutting process Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 210000003746 feather Anatomy 0.000 description 1
- 239000000446 fuel Substances 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000003754 machining Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000005096 rolling process Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000009738 saturating Methods 0.000 description 1
- 230000035939 shock Effects 0.000 description 1
- 238000005507 spraying Methods 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 238000012800 visualization Methods 0.000 description 1
Landscapes
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
The invention discloses a kind of Hartmann's ray tracing methods based on colored three steps transposition technology, including are calculated image using colored three step phase-shifting methods and obtained phase information；The present invention is using colored three step phase-shifting methods, in conjunction with the anti-integral curve three-dimensional reconstruction algorithm of three-dimensional non-uniform dielectric field, pass through the combination settings of projection screen and telecentric optical system, it realizes the accurate trace of the measurement light to non-uniform dielectric field and the transient state refracting characteristic of three-dimensional space refractive index is measured, substantially increase measuring accuracy and efficiency；And the measuring system whole design of the non-uniform dielectric field based on Hartmann's ray tracing is accurate, and measurement accuracy is high, and cost is relatively low, has a wide range of application, and has important theory significance and engineering application value, is suitble to promote and apply.
Description
Technical field
The present invention relates to a kind of measurement methods more particularly to a kind of Hartmann's light based on colored three steps transposition technology to chase after
Mark method, non-uniform dielectric field refer to since the spatial distribution of the factors such as pressure, speed, temperature, density, substance is unevenly formed
Space physics characteristic spatial field heterogeneous, including due to temperature, pressure, flowing formed non-uniform dielectric field, also include
The non-uniform dielectric field that liquid liquid blending agent, gas-liquid mixed media, gas medium are formed.
Background technique
Non-uniform dielectric field refers to due to the uneven shape of the spatial distribution of the factors such as pressure, speed, temperature, density, substance
At space physics characteristic spatial field heterogeneous, the objects such as temperature, pressure, air-flow are able to achieve by the measurement of non-uniform dielectric field
The reproduction of reason amount.Therefore, to the three-dimensional measurement of non-uniform dielectric field, the scientific researches such as hydrodynamics, thermal conduction study and the energy,
There is highly important scientific meaning and broad application prospect in many industries such as environment, power and chemical industry.
Chemical reaction flow, laminar flow, turbulent flow, thin stream, low temperature flow, high-temperature aerodynamics, supersonic and hypersonic
Flowing and wind tunnel experiment, water hole flow field are all typical non-uniform dielectric fields, are accurately realized the inspection of non-uniform dielectric field
It surveys, effective parameter can be provided for chemical reaction processes, the design of aerodynamics equipment and the controls such as burning etc., realized optimal
Control strategy, in Industrial Boiler, power plant is coal-fired, combustion gas turbine and aviation, automobile engine burning, marine propeller,
The fields such as spraying operation are all widely used.For example, the ignition temperature field for being applied to boiler monitors, hearth combustion operating condition
It is unreasonable be the principal element for leading to boiler operatiopn accident, combustion state it is reasonable with it is unreasonable, directly influence boiler hot
The height of efficiency and the ratio of power supply and fuel consumption.Therefore, it is the uniform and stable operation of combustion field for keeping burner hearth, improves burning
Efficiency realizes on-line monitoring to the Temperature Distribution of furnace chamber inner combustion field, is the requisite measure for regulating and controlling combustion state；It also can be air
The design of dynamics equipment and the measurement data that offer basis is provided, as rocket engine ejects under high, vacuum state
Shape like feather jet flow plume measurement, it is unstable that plume flow pattern defect will lead to the swing of aircraft, thrust decay, flight
Determine, cannot fly by planned orbit, final aircraft is likely to that scheduled height and region is not achieved in flight course, to plumage
The real-time and accurate measurement for flowing flow field is the basis for adjusting rocket blast type；It can also be provided in real time for petroleum, gas extraction
The data of mixing flow field, such as in petroleum, natural gas extraction Oil, Water, Gas ratio, this be related to oil and natural gas exploitation plan
It slightly, is the basis of oil and natural gas exploitation.
Currently, mainly passing through the change of the refractive index of measurement non-uniform dielectric field to the universal measurement method of non-uniform dielectric field
Change to realize.When light penetrates the refractive index field of inhomogeneous medium, it will two kinds of variations occur, first is that light deviates original side
To second is that the opposite variation phase that light generates position phase is moved.Index method is exactly to determine refractive index according to the two characteristics.
Index method has the advantages that other methods are irreplaceable, they are noiseless to former field, belong to non-intrusive measurement, temporal resolution
Height can be used to study very fast transient process, or even can track the process of spatial non-uniform medium field variation；Spatial discrimination
Rate is high, can be used to study temperature gradient minimax spatial non-uniform medium field, such as shock wave, flammule；Sky can be provided simultaneously
Between non-uniform dielectric field measurement, rather than put measurement, contain much information.Other main methods include schlieren method, shadowing method, light
Learn interferometry, holographic interferometry, More's deflection, laser speckle photography, laser speckle shear interference and platform uncle's interference etc..Shade
Method, schlieren method are measured according to the Diffraction dispersion of light；Through Optical Interference Spectra and holographic interference rule directly measure phase change, and
Assuming that light is without deviation；More's deflection, laser speckle photography, laser speckle shear interference and platform uncle's interference etc. directly measure light
Deflection angle.
Schlieren method can generate this phenomenon of deviation in the different Propagation of internal density using light, by using knife
Mouth device blocks the light of deviation, and then the throughput of light generates variation, on the receiving planes such as optical screen or camera
Light and shade variation can be observed and realize refractometry.In addition to traditional schlieren method, currently used schlieren method also focus schlieren method,
Background schlieren method, rainbow schlieren method etc..Shadowing method principle is as schlieren method, but it is the deviation bring by recording light
The measurement of refractive index is realized in the displacement put on projection screen.The experimental principle and system of schlieren method and shadowing method constitute all simpler
It is single, but the precision of shadowing method is relatively low, and traditional schlieren method is influenced vulnerable to background luminance, and background schlieren method cannot be determined accurately often
The displacement of one point limits the precision of measurement.
Common Through Optical Interference Spectra is usually using Mach-Zehnder interferometer, bimirror interferometer, Mai Keerxun interferometer, platform
Graceful interferometer and schlieren interferometer etc..The light that same light source issues is divided into two beams by they, a branch of after measured zone and another
A branch of reference light is relevant, by analyzing the phase change of the light through measured zone, realizes that the measurement of spatial refractive index, reaction are empty
Between non-uniform dielectric field information.Mach-Zehnder interferometer, bimirror interferometer, Mai Keerxun interferometer, the graceful interferometer of platform etc.
Amplitude optical spectroscopy is all used, and schlieren interferometer then uses polarization-nulling method.Holographic interferometry and ordinary optical interferometry principle
It is identical, but same light beam is recorded on same holographic dry plate in different times, is then made using time optical spectroscopy by it
These wavefront interfere again, by analyzing the phase change of the light through measured zone, realize the measurement of spatial refractive index, instead
Answer the non-uniform dielectric field information in space.Interferometry is widely applied in heat transfer, burning and aerodynamic studies at present, and
A kind of most basic and most classic method in Measurement in Thermophysics.
More's deflection method realizes the measurement of refractive index in non-uniform dielectric field by moire phenomenon, what light source issued
Light places grating by two front and backs, forms the Moire fringe of two grating fringes superposition, under normal circumstances, Moire fringe is flat
The orderly distribution of capable or by setting rule, when measuring space is non-uniform dielectric field, light can be because of the variation of refractive index
Deviation occurs, Moire fringe can be deformed or be displaced, and the displacement for analyzing the striped can be obtained the deflection angle of light, to count
Calculate space non-uniform dielectric field distribution.More's deflection method reduces interferometry to machinery, optical component and environmental Kuznets Curves
It is required that application range is wider, but the diffraction phenomena in More's skewing method reduces the contrast and spatial resolution of striped.
Platform uncle's interference is measured using the primary effect of platform.When the directional light vertical irradiation period is the grating of d, at it
A certain plane will appear the identical picture of the grating afterwards, and referred to as from being imaged, i.e. the primary effect of platform, imaging surface are known as the primary face of platform for this.
Another identical grating is placed in the primary face of the platform of previous grating, latter grating and previous parallel gratings and coaxial, and
Two gratings are identical in vertical axis positioning, and zero light field is observed after latter grating at this time.When being inserted into a phase between two gratings
The imaging certainly of object non-uniform dielectric field, light generation disturbance, first grating can generate distortion, this distortion is imaged together certainly
Second grating interaction, interference fringe can be observed on the film viewing screen after second grating, investigates the interference fringe energy
The information at light deflection angle is obtained, uncle interferes here it is platform.Platform uncle's interference is chiefly used in the measurement of optical instrument and parameter, has now also sent out
Exhibition is used for the measurement of non-uniform dielectric field.But the requirement on machining accuracy of above-mentioned these interferometry measuring devices manufacture is very high, if
Standby adjustment is also very time-consuming laborious, so that equipment prevailing price is higher.
Laser speckle photography generates countless small points in the diffusing reflection of body surface using the good laser beam of coherence
Light source, the coherent light that they are issued are concerned in spatial and each other, and different divide is generated on the projection device of different location
The speck and blackening of cloth, when it passes through spatial non-uniform medium field, the pattern generated on projection device can be according to certain
Rule movement, tracks its motion information, obtains the displacement of each speckle point, so that it may which Inversion Calculation goes out spatial non-uniform medium field folding
Penetrate the variation of rate.The technology is insensitive to vibrating, and does not need the gradient direction that non-uniform dielectric field is known in advance.But speckle
The point-to-point analysis of photograph is time-consuming and troublesome, and needs to correct the width of fringe generated due to Halo effect.
Speckle is imaged on the CCD of camera by laser speckle shear interference by diplopore, allow CCD crystal oscillator face slightly offset from
The picture plane of conjugation, so that this is sheared, they are relevant by the speckle of diplopore at the picture of a slightly defocus on CCD
Fringe pattern is formed, here it is speckle-shearing interference.It is the deflection angle of light, shape measured by laser speckle shearing interference method
At coherent fringe be light equal displacement gradient line, represent the distribution situation of spatial heterogeneity, therefore laser speckle is cut
Online visualization measurement can be carried out to spatial non-uniform medium field by cutting interference.Laser speckle shear interference is commonly used in strain, inclines
Tiltedly, bending and the measurement of other mechanical parameters, the measurement for non-uniform dielectric field are also developing.The method reduce laser to dissipate
The calculation amount of point-to-point analysis in spot photograph, but need the correction to Halo effect width of fringe.
The above-mentioned spatial measuring method to non-uniform dielectric field is made a general survey of, respectively there are advantage and disadvantage, interference class method processes equipment
Installation accuracy requires height, so that its higher cost, and schlieren method, the measurement accuracy of shadowing method are lower.Therefore, the invention patent
It is proposed that a kind of inexpensive, high accuracy three-dimensional refractive index measurement method, this method are based on Hartmann's ray tracing principle, pass through throwing
The combination settings of shadow screen and telecentric optical system realize the accurate trace of measurement light, to realize the survey of three-dimensional space refractive index
Amount.
Summary of the invention
Present invention solves the technical problem that being to overcome the deficiencies of existing technologies, provide a kind of by projection screen and telecentricity light
The combination settings of system realize the accurate trace of measurement light and realize the measurement of three-dimensional space refractive index based on colour three
Walk Hartmann's ray tracing method of transposition technology.
In order to solve the above-mentioned technical problem, the present invention is addressed by following technical proposals:
Based on Hartmann's ray tracing method of colored three steps transposition technology, including following particular content:
Using each luminous point on projection screen be point light source to measurement space project spherical surface incident ray, wherein have one with
Optical axis forms the non-uniform dielectric field refraction that deflection angle penetrates measurement space interior for the spherical surface incident ray of θ and is formed in parallel with light
The deviation light of axis, the deviation light directive lens for being parallel to optical axis projected from measurement space, lens are to being parallel to optical axis
Deviation light is reflected again, and the deviation light after reflecting again is screened and converged in by object space telecentric stop
In CCD camera, CCD camera shoots the deviation light of convergence and obtains image；Incident light is derived by by formula (A)
Angle corresponding relationship between line and refracted light is derived by deflection angle theta and shadow as shown in formula (B) according to formula (A)
The relationship of image position deviation, to obtain the folding of the non-uniform dielectric field as shown in formula (C) by formula (A) and formula (B)
Penetrate the relationship between the integral of rate gradient and blurring amount；Image is calculated using three step phase-shifting methods of colour and obtains phase
Position information determines the light of the image generated after the refraction that the luminous point on projection screen passes through non-uniform dielectric field by phase information
Refracting characteristic and corresponding relationship are learned, to realize ray tracing；Wherein, formula (A), formula (B), the following institute of formula (C) difference
Show:
Wherein, z is optical axis direction,Variable quantity respectively in X, Y, Z axis, n are spatial refractive index,For the variation of spatial refractive index, F is the focal length of lens, and δ is displacement of the pixel of image on the imaging screen of CCD camera
Amount, S1For the distance between the rear focus of object space telecentric stop and imaging screen of CCD camera, L is that projection screen and measurement are empty
Between the distance between ejecting end.
Preferably, the non-uniform dielectric of measurement space interior flowing generates temperature field and/or pressure field, realize to by
Measurement in the non-uniform dielectric field that temperature, pressure, flowing are formed.Wherein, non-uniform dielectric field includes by measurement space interior
The temperature field and/or pressure field that the liquid liquid blending agent or gas-liquid mixed media or air dielectric of flowing generate, measurement space are set
There is the temperature-adjusting device for adjusting temperature field existing for measurement space interior, realizes to due to temperature, pressure, flowing shape
At non-uniform dielectric field measurement, also achieve and formed not to by liquid liquid blending agent, gas-liquid mixed media, gas medium
The measurement of same non-uniform dielectric field.Wherein, liquid liquid blending agent is oil-water mixture, and gas-liquid mixed media is natural gas and oil
Mixed liquor, gas medium be air draught or other mixed gas air-flows, pressure field is by liquid liquid mixed flow or gas-liquid mixed
Flowing or gas flowing generate；Temperature-adjusting device includes heating temperature-control system, thermostatic control system, and stationary temperature field is by temperature
Regulating device is spent to provide.
Preferably, colored three step phase-shifting methods are based on three step phase-shifting methods, colored three step phase-shifting methods include following particular content:
By modulation intelligence I1、I2、I3、I4、I5、I6It is modulated to tri- differences of R, G, B of colored RGB respectively by brightness subregion
In channel, make tri- different channels R, G, B of the colored RGB on projection screen and CCD camera respectively at horizontal and vertical two
Three phase images are respectively exported and acquired on direction and obtain three phase informations, and CCD camera transmits three phase informations
To projection screen；By formula (D) calculate CCD camera acquisition RGB triple channel information M=[RGB] ', RGB triple channel information M with
Actual RGB triple channel information N=[rgb] ' by transition matrix A and background matrix M0=[R0G0B0] ' realized by formula (E)
Conversion；It is final to realize the measurement of transient state refracting characteristic according to the actual RGB triple channel information N being calculated；Wherein, projection screen
Show the cosine wave that three outs of phase are moved, the cosine wave that respective three outs of phase are moved in horizontal and vertical both direction
Modulation intelligence I1、I2、I3、I4、I5、I6It is calculated respectively by formula (F), formula (G)；Pass through formula (H), formula respectively
(I) phase information of each pixel is calculated and is obtained to three phase images of CCD camera acquisition, it is each to what is obtained
The cosine wave striped that three outs of phase of the phase information of a pixel and projection screen display are moved compares, so that it is determined that throwing
The light refractive properties and corresponding relationship of the image generated after the refraction that luminous point on shadow screen passes through inhomogeneous medium field, thus
Realize ray tracing；Formula (D), (E), (F), (G), (H), (I) are as follows respectively:
M=AN+M0(E),
Wherein, I'(x, y) it is background gray levels, I " (x, y) is modulation gray value；
Wherein, φ1(x,y)、φ2(x, y) is phase information.
In the present invention, non-homogeneous Jie existing for measurement space interior is calculated using anti-integral curve three-dimensional reconstruction algorithm
The delta data of the refractive index of matter field；Anti- integral curve three-dimensional reconstruction algorithm includes following particular content: being established with non-homogeneous Jie
The index distribution of matter field is the ill-condition matrix of variable, is then solved to ill-condition matrix, so that Problems of Reconstruction is converted to
The Solve problems of ill-condition matrix；Wherein, the disease as shown in formula (J) is obtained later with discrete by carrying out deformation to formula (C)
State matrix:
T=WDS2(J),
Wherein T is the displacement information matrix of image picture point, and measurement light of the characterization compared with uniform dielectric field passes through non-equal
The displacement that even medium field generates directly is obtained by the corresponding relationship of image；W is model coefficient matrix, passes through direct volume reconstruction
Model coefficient calculation method under strategy obtains, and proposed adoption monte carlo method regards it as and asks any in volume elements space
Any falls in the probability in light beam, then a large amount of stochastical sampling can be directly carried out in volume elements space, according to sampled point and light
The spatial relation of beam, and then calculate sampled point in sample space and fall in the probability in light beam, when sample size is sufficiently large,
Think that this probability value is model coefficient；D is variations in refractive index vector matrix to be solved；S2It is projection matrix, by three-dimensional track
It projects in the offset of distance；The calculation method of formula (J) can specifically be rebuild using addition algebraic reconstruction, multiplication algebra or filter
Wave backprojection reconstruction algorithm is solved, and realizes the reconstruction of three-dimensional non-uniform dielectric field.
Hartmann's ray tracing method based on colored three steps transposition technology of the invention is by based on Hartmann's light
The measuring system of the non-uniform dielectric field of trace realizes, the measuring system of the non-uniform dielectric field based on Hartmann's ray tracing
Space, the projection screen for projecting spherical surface incident ray to measurement space, use are measured there are non-uniform dielectric field including internal
In on projection screen luminous point project spherical surface incident ray penetrate measurement space interior non-uniform dielectric field reflect to be formed it is flat
Row is screened and is pooled the object space telecentric optical system of image in the deviation light of optical axis, measurement be spatially located at projection screen with
Between object space telecentric optical system；Phase information is calculated image and obtained using three step phase-shifting methods of colour, passes through phase
Information determines the light refractive properties and corresponding relationship of image；It is empty that measurement is calculated using anti-integral curve three-dimensional reconstruction algorithm
Between non-uniform dielectric field refractive index existing for inside delta data；Non-homogeneous Jie based on Hartmann's ray tracing of the invention
The fundamental measurement principle of the measuring system of matter field is as follows: projecting using each luminous point on projection screen as point light source to measurement space
Spherical surface incident ray, wherein having one to form the spherical surface incident ray that deflection angle is θ with optical axis penetrates the non-of measurement space interior
Uniform dielectric field reflects the deviation light for being formed in parallel with optical axis, and the deviation light for being parallel to optical axis projected from measurement space is penetrated
To object space telecentric optical system, object space telecentric optical system screens deviation light and pools image；Using colour three
Step phase-shifting method is calculated image and is obtained phase information, and the light refractive properties and correspondence of image are determined by phase information
Relationship, to realize ray tracing；Measurement space interior is calculated using anti-integral curve three-dimensional reconstruction algorithm and there is flowing
Non-uniform dielectric field refractive index delta data, thus realize to measurement space interior existing for non-uniform dielectric field three-dimensional
The real-time measurement of spatial refractive index；For the present invention using colored three step phase-shifting methods, the anti-integral in conjunction with three-dimensional non-uniform dielectric field is bent
Line three-dimensional reconstruction algorithm realizes the measurement to non-uniform dielectric field by the combination settings of projection screen and telecentric optical system
The accurate trace of light and the transient state refracting characteristic of three-dimensional space refractive index is measured, substantially increases measuring accuracy and effect
Rate；And the overall system design is accurate, and measurement accuracy is high, and cost is relatively low, has a wide range of application, with important theory significance and
Engineering application value is suitble to promote and apply.
In the present invention, object space telecentric optical system includes for being reflected to the deviation light projected from measurement space
Lens, for the deviation light of convergence is shot and is obtained the CCD camera of image, by lens project deviation light converge
Object space telecentric stop in CCD camera, lens are located between measurement space and object space telecentric stop, and object space telecentric stop is located at
Between lens and CCD camera, projection screen, measurement space, lens, object space telecentric stop, CCD camera are using optical axis as center axis
And successively face is arranged.In the present invention, spherical surface is projected to measurement space using each luminous point on projection screen as point light source
Incident ray, wherein having one to form the spherical surface incident ray that deflection angle is θ with optical axis penetrates the non-homogeneous of measurement space interior
Medium field reflects the deviation light for being formed in parallel with optical axis, and the deviation light directive for being parallel to optical axis projected from measurement space is saturating
Mirror, lens reflect the deviation light for being parallel to optical axis again, and the deviation light after reflecting again passes through object space telecentricity
Diaphragm is screened and is converged in CCD camera, and CCD camera shoots the deviation light of convergence and obtains image, it may be assumed that logical
The incident ray that the luminous point crossed on projection screen projects is successively by the refraction twice in measurement space, lens by object space telecentric stop pair
Deviation light is screened and is converged in CCD camera, so that CCD camera shoots the deviation light of convergence and obtains shadow
Picture, compared with traditional Hartmann's method, the present invention using projection screen, lens, object space telecentric stop, CCD camera combination light
Design is learned instead of Hartmann's screen, the track of imaging ray is limited, realizes the quick trace of light.
In the present invention, the measuring system of the non-uniform dielectric field based on Hartmann's ray tracing further includes measurement annulus, is thrown
Shadow screen, lens are each provided on the inner wall of measurement annulus, and measurement space is the border circular areas being arranged concentrically with measurement annulus, circle
The radius in domain is less than the radius of measurement annulus, and optical axis passes through the center of circle of measurement annulus and border circular areas, and projection screen and lens are with circle
It is symmetrical set centered on the heart.The measurement annulus and border circular areas being arranged concentrically further improve the accuracy and effect of measurement
Rate.
In the present invention, projection screen, lens, object space telecentric stop, CCD camera are made of 4-8 group, 4-8 group projection screen, thoroughly
Mirror, object space telecentric stop, CCD camera are uniformly distributed along 360 degree of circumferential directions of measurement annulus, form multiple groups back projection type Hart
Graceful system and its cross section for measuring light section, every group of back projection type Hartmann system obtain light and exist according to ray tracing, measurement
The linear displacement generated in the measurement direction, and according to multiple measurement direction glazed thread traces as a result, in conjunction with corresponding anti-
Integral curve three-dimensional reconstruction algorithm calculates the delta data for obtaining the refractive index of non-uniform dielectric field of tested three-dimensional space, into
One step improves the accuracy and efficiency of measurement.
In the present invention, the angle corresponding relationship between incident ray and refracted light, foundation are derived by by formula (A)
Formula (A) is derived by the relationship of deflection angle theta and image position deviation as shown in formula (B), thus by formula (A) and public
Formula (B) obtains the relationship between the integral and blurring amount of the gradient of non-uniform dielectric field refractive index as shown in formula (C).
In the present invention, CCD camera is connect with projection screen, and CCD camera is by the image transmission for shooting and obtaining to projection screen, shadow
As being shown on projection screen.
In the present invention, the non-uniform dielectric field of measurement space interior flowing generates temperature field and/or pressure field, measures space
Equipped with the temperature-adjusting device for adjusting temperature field existing for measurement space interior, realize to due to temperature, pressure, flowing
The measurement of the non-uniform dielectric field of formation.Wherein, non-uniform dielectric field is oil-water mixture, and pressure field is by non-uniform dielectric field flow
Movable property is raw；Temperature-adjusting device includes heating temperature-control system, thermostatic control system, and stationary temperature field is mentioned by temperature-adjusting device
For.
The present invention has significant technical effect: using colored three step phase-shifting methods, knot due to using above technical scheme
The anti-integral curve three-dimensional reconstruction algorithm for closing three-dimensional non-uniform dielectric field, is set by the combination of projection screen and telecentric optical system
It sets, realize the accurate trace of the measurement light to non-uniform dielectric field and the transient state refracting characteristic of three-dimensional space refractive index is surveyed
Amount, substantially increases measuring accuracy and efficiency；And the overall system design is accurate, and measurement accuracy is high, and cost is relatively low, application
Range is wide, has important theory significance and engineering application value, is suitble to promote and apply.
Detailed description of the invention
Fig. 1 is the structural schematic diagram using the measuring system of traditional Hartmann's method.
Fig. 2 is that the measuring principle of the measuring system of the non-uniform dielectric field of the invention based on Hartmann's ray tracing is implemented
The structural schematic diagram of example.
Fig. 3 is the back projection type Hartmann of the invention being made of 4 groups of projection screens, lens, object space telecentric stop, CCD camera
The structural schematic diagram of measuring system embodiment.
Fig. 4 is the phase images of -2 π/3 that projection screen of the invention is shown in a lateral direction.
Fig. 5 is 0 phase images that projection screen of the invention is shown in a lateral direction.
Fig. 6 is the phase images of 2 π/3 that projection screen of the invention is shown in a lateral direction.
Fig. 7 is the phase images of -2 π/3 that projection screen of the invention is shown in a longitudinal direction.
Fig. 8 is 0 phase images that projection screen of the invention is shown in a longitudinal direction.
Fig. 9 is the phase images of 2 π/3 that projection screen of the invention is shown in a longitudinal direction.
Figure 10 is the RGB color phase images synthesized by Fig. 4,5,6,7,8,9 that projection screen of the invention is shown.
Specific embodiment
Present invention is further described in detail with embodiment with reference to the accompanying drawing.
Based on Hartmann's ray tracing method of colored three steps transposition technology, as shown in Figure 2,3, including it is following specific interior
Hold:
Spherical surface incident ray is projected to measurement space 3 using each luminous point on projection screen 2 as point light source, wherein there is one
With optical axis 1 formed deflection angle be θ spherical surface incident ray penetrate measurement space 3 inside non-uniform dielectric field reflect to be formed it is parallel
In the deviation light of optical axis 1, the deviation light directive lens 4 for being parallel to optical axis 1 projected from measurement space 3, lens 4 are to parallel
It is reflected again in the deviation light of optical axis 1, the deviation light after reflecting again is screened by object space telecentric stop 5
And converge in CCD camera 6, the deviation light of 6 pairs of CCD camera convergences is shot and obtains image；It is derived by formula (A)
The angle corresponding relationship between incident ray and refracted light is obtained, is derived by as shown in formula (B) partially according to formula (A)
The relationship of rotational angle theta and image position deviation, to obtain non-homogeneous Jie as shown in formula (C) by formula (A) and formula (B)
Relationship between the integral and blurring amount of the refractive index gradient of matter field；Image is calculated using three step phase-shifting methods of colour
And phase information is obtained, it is generated after determining the refraction that the luminous point on projection screen 2 passes through non-uniform dielectric field by phase information
Image light refractive properties and corresponding relationship, to realize ray tracing；Wherein, formula (A), formula (B), formula (C)
It is as follows respectively:
Wherein, z is 1 direction of optical axis,Variable quantity respectively in X, Y, Z axis, n are space refraction
Rate,For the variation of spatial refractive index, F is 4 focal length of lens, and δ is the pixel of image on the imaging screen of CCD camera 6
Displacement, S1For object space telecentric stop 5 rear focus and CCD camera 6 the distance between imaging screen, L be projection screen 2 with
Measure the distance between the ejecting end in space 3.
In the present embodiment, the non-uniform dielectric of measurement 3 internal flow of space generates temperature field and/or pressure field, realizes
Measurement to the non-uniform dielectric field formed due to temperature, pressure, flowing.Wherein, non-uniform dielectric field includes by measuring space 3
The temperature field and/or pressure field that the liquid liquid blending agent or gas-liquid mixed media or air dielectric of internal flow generate, measurement are empty
Between 3 be equipped with for adjusts measure space interior existing for temperature field temperature-adjusting device, realize to due to temperature, pressure,
The measurement for flowing the non-uniform dielectric field formed, also achieves to by liquid liquid blending agent, gas-liquid mixed media, gas medium shape
At different non-uniform dielectric fields measurement.Wherein, liquid liquid blending agent is oil-water mixture, and gas-liquid mixed media is natural
The mixed liquor of gas and oil, gas medium are air draught or other mixed gas air-flows, and pressure field is by liquid liquid mixed flow or gas
Liquid mixed flow or gas flowing generate；Temperature-adjusting device includes heating temperature-control system, thermostatic control system, stationary temperature
Field is provided by temperature-adjusting device
In the present embodiment, colored three step phase-shifting methods include following particular content:
By modulation intelligence I1、I2、I3、I4、I5、I6It is modulated to tri- differences of R, G, B of colored RGB respectively by brightness subregion
In channel, make tri- different channels R, G, B of the colored RGB on projection screen 2 and CCD camera 6 respectively horizontal and vertical two
Three phase images are respectively exported and acquired on a direction and obtain three phase informations, and CCD camera 6 passes three phase informations
Transport to projection screen 2；The RGB triple channel information M=[RGB] ' that CCD camera 6 acquires, RGB triple channel information are calculated by formula (D)
M and actual RGB triple channel information N=[rgb] ' by transition matrix A and background matrix M0=[R0G0B0] ' real by formula (E)
Now convert；It is final to realize the measurement of transient state refracting characteristic according to the actual RGB triple channel information N being calculated；Wherein, it projects
The cosine wave that 2 three outs of phase of display of screen are moved, in horizontal and vertical both direction more than respective three outs of phase shifting
String wave modulation information I1、I2、I3、I4、I5、I6It is calculated respectively by formula (F), formula (G), at horizontal and vertical two
The cosine wave that respective three outs of phase are moved on direction as shown in Fig. 4,5,6,7,8,9, is synthesized respectively by Fig. 4,5,6,7,8,9
RGB color phase images be the spacing phase images of 2 π/3 as shown in Figure 10；Respectively by formula (H), public formula (I) to CCD
Three phase images that camera 6 acquires are calculated and are obtained the phase information of each pixel, to each obtained pixel
Phase information and projection screen 2 on the cosine wave striped that moves of three outs of phase showing compare, so that it is determined that projection screen 2
On luminous point by the light refractive properties and corresponding relationship of the image generated after the refraction of inhomogeneous medium field, to realize
Ray tracing；Formula (D), (E), (F), (G), (H), (I) are as follows respectively:
M=AN+M0(E),
Wherein, I'(x, y) it is background gray levels, I " (x, y) is modulation gray value；
Wherein, φ1(x,y)、φ2(x, y) is phase information.
In the present embodiment, it is calculated using anti-integral curve three-dimensional reconstruction algorithm non-equal existing for 3 inside of measurement space
The delta data of the refractive index of even medium field；Anti- integral curve three-dimensional reconstruction algorithm includes following particular content: being established with non-equal
The index distribution of even medium field is the ill-condition matrix of variable, is then solved to ill-condition matrix, so that Problems of Reconstruction be turned
It is changed to the Solve problems of ill-condition matrix；Wherein, it is obtained as shown in formula (J) later by carrying out deformation to formula (C) with discrete
Ill-condition matrix:
T=WDS2(J),
Wherein T be image picture point displacement information matrix, W be model coefficient matrix, D be variations in refractive index to be solved to
Moment matrix, S2It is projection matrix.
Hartmann's ray tracing method based on colored three steps transposition technology of the present embodiment is by based on hartmann diaphragm
The measuring system of the non-uniform dielectric field of line trace realizes, the measurement system of the non-uniform dielectric field based on Hartmann's ray tracing
System enters including internal there are the measurement space 3 of non-uniform dielectric field, for projecting spherical surface to measurement space 3 as shown in Figure 2,3
Penetrate light projection screen 2, for on projection screen 2 luminous point project spherical surface incident ray penetrate measurement space 3 inside it is non-
The deviation light that the refraction of uniform dielectric field is formed in parallel with optical axis 1 is screened and pools the object space telecentric optics system of image
System measures space 3 between projection screen 2 and object space telecentric optical system；Image is counted using three step phase-shifting methods of colour
Phase information is calculated and obtained, the light refractive properties and corresponding relationship of image are determined by phase information；Using anti-integral curve
The delta data of non-uniform dielectric field refractive index existing for 3 inside of measurement space is calculated in three-dimensional reconstruction algorithm.
In the present embodiment, object space telecentric optical system includes for rolling over to the deviation light projected from measurement space 3
The lens 4 penetrated, the CCD camera 6 for image to be shot and obtained to the deviation light of convergence, for what is projected to lens 4
Deviation light is screened and is converged in the object space telecentric stop 5 in CCD camera 6, and it is remote with object space that lens 4 are located at measurement space 3
Between heart diaphragm 5, object space telecentric stop 5 is between lens 4 and CCD camera 6, projection screen 2, measurement space 3, lens 4, object space
Telecentric stop 5, CCD camera 6 are center axis with optical axis 1 and successively face is arranged.In the present invention, on projection screen
Each luminous point is point light source to measurement space injection spherical surface incident ray, wherein having one to form deflection angle with optical axis is θ's
The non-uniform dielectric field refraction that spherical surface incident ray penetrates measurement space interior is formed in parallel with the deviation light of optical axis, from measurement
The deviation light directive lens for being parallel to optical axis that space is projected, lens roll over the deviation light for being parallel to optical axis again
It penetrates, the deviation light after reflecting again is screened by object space telecentric stop and converged in CCD camera, CCD camera pair
The deviation light of convergence is shot and obtains image, it may be assumed that the spherical surface incident ray projected by luminous point on projection screen is successively
Refraction twice by measuring space, lens is screened deviation light by object space telecentric stop and is converged in CCD camera,
So that CCD camera is shot and is obtained image to the deviation light of convergence, with traditional Hartmann's method phase shown in FIG. 1
Than, the present embodiment using projection screen, lens, object space telecentric stop, CCD camera combination optical design instead of Hartmann's screen 8,
The track of imaging ray is limited, the quick trace of light is realized.
In the present embodiment, the measuring system of the non-uniform dielectric field based on Hartmann's ray tracing further includes measurement annulus 7,
Projection screen 2, lens 4 are each provided on the inner wall of measurement annulus 7, and measurement space 3 is the circle being arranged concentrically with measurement annulus 7
Domain, the radius of border circular areas are less than the radius of measurement annulus 7, and optical axis 1 passes through the center of circle of measurement annulus 7 and border circular areas, projection
Screen 2 is symmetrical set centered on the center of circle with lens 4.
In the present embodiment, projection screen 2, lens 4, object space telecentric stop 5, CCD camera 6 are constituted by 4 groups, 4 groups of projection screens 2,
Lens 4, object space telecentric stop 5, CCD camera 6 are uniformly distributed along 360 degree of circumferential directions of measurement annulus 7.
In the present embodiment, the angle corresponding relationship between incident ray and refracted light is derived by by formula (A), according to
The relationship of deflection angle theta and image position deviation as shown in formula (B) is derived by according to formula (A), thus by formula (A) and
Formula (B) obtains the relationship between the integral and blurring amount of the non-uniform dielectric field refractive index gradient as shown in formula (C).
In the present embodiment, CCD camera 6 is connect with projection screen 2, and CCD camera 6 extremely projects the image transmission for shooting and obtaining
Screen 2, image is shown on projection screen 2.
In the present embodiment, the non-uniform dielectric field of measurement space interior flowing generates temperature field and/or pressure field, and measurement is empty
Between be equipped with for adjusts measure space interior existing for temperature field temperature-adjusting device, realize to due to temperature, pressure, stream
The measurement of the dynamic non-uniform dielectric field formed.Wherein, non-uniform dielectric field is oil-water mixture, and pressure field is by non-uniform dielectric field
Flowing generates；Temperature-adjusting device includes heating temperature-control system, thermostatic control system, and stationary temperature field is by temperature-adjusting device
It provides.
In short, the foregoing is merely presently preferred embodiments of the present invention, it is all according to equalization made by scope of the present invention patent
Variation and modification, shall all be covered by the patent of the invention.
Claims (3)
1. Hartmann's ray tracing method based on colored three steps transposition technology, it is characterised in that: including following particular content:
Spherical surface incident ray is projected to measurement space (3) using each luminous point on projection screen (2) as point light source, wherein there is one
Deflection angle is formed with optical axis (1) to reflect to be formed for non-uniform dielectric field that the spherical surface incident ray of θ penetrates measurement space (3) inside
It is parallel to the deviation light of optical axis (1), the deviation light directive lens for being parallel to optical axis (1) projected from measurement space (3)
(4), lens (4) reflect the deviation light for being parallel to optical axis (1) again, and the deviation light after reflecting again passes through object
Square telecentric stop (5) is screened and is converged on CCD camera (6), and CCD camera (6) shoots the deviation light of convergence
And obtain image；It is derived by the angle corresponding relationship between incident ray and refracted light by formula (A), according to formula
(A) it is derived by the relationship of deflection angle theta and image position deviation as shown in formula (B), to pass through formula (A) and formula
(B) relationship between the integral and blurring amount of the refractive index gradient of the inhomogeneous medium field as shown in formula (C) is obtained；It adopts
Phase information is calculated image and obtained with three step phase-shifting methods of colour, the light on projection screen (2) is determined by phase information
The light refractive properties and corresponding relationship of the image generated after the refraction that point passes through inhomogeneous medium field, to realize that light chases after
Mark；Wherein, formula (A), formula (B), formula (C) are as follows respectively:
Wherein, z is optical axis (1) direction,Variable quantity respectively in X, Y, Z axis, n are spatial refractive index,For the variation of spatial refractive index, F is lens (4) focal length, and δ is the pixel of image on the imaging screen of CCD camera (6)
Displacement, S1For the rear focus of object space telecentric stop (5) and the distance between the imaging screen of CCD camera (6), L is to throw
The distance between the ejecting end of shadow screen (2) and measurement space (3).
2. Hartmann's ray tracing method according to claim 1 based on colored three steps transposition technology, it is characterised in that:
The non-uniform dielectric for measuring space (3) internal flow generates temperature field and/or pressure field.
3. Hartmann's ray tracing method according to claim 1 or 2 based on colored three steps transposition technology, feature exist
In: colored three step phase-shifting methods include following particular content:
By modulation intelligence I1、I2、I3、I4、I5、I6It is modulated to R, G, B tri- different channels of colored RGB respectively by brightness subregion
In, the channel for keeping projection screen (2) different with R, G, B tri- of colored RGB on CCD camera (6) is respectively horizontal and vertical two
Three phase images are respectively exported and acquired on a direction and obtain three phase informations, and CCD camera (6) is by three phase informations
It is transmitted to projection screen (2)；The RGB triple channel information M=[RGB] ' of CCD camera (6) acquisition, RGB threeway are calculated by formula (D)
Road information M and actual RGB triple channel information N=[rgb] ' by transition matrix A and background matrix M0=[R0G0B0] ' pass through public affairs
Formula (E) realizes conversion；It is final to realize the measurement of transient state refracting characteristic according to the actual RGB triple channel information N being calculated；Its
In, projection screen (2) shows the cosine wave that three outs of phase are moved, the respective three different phases in horizontal and vertical both direction
The modulation intelligence I of the cosine wave of displacement1、I2、I3、I4、I5、I6It is calculated respectively by formula (F), formula (G)；Pass through respectively
Formula (H), public formula (I) calculate three phase images that CCD camera (6) acquires and obtain the phase letter of each pixel
It ceases, the cosine wave striped that three outs of phase shown in the phase information and projection screen (2) to each obtained pixel are moved
It compares, so that it is determined that the optics of the image generated after the refraction that the luminous point on projection screen (2) passes through inhomogeneous medium field
Refracting characteristic and corresponding relationship, to realize ray tracing；Formula (D), (E), (F), (G), (H), (I) are as follows respectively:
M=AN+M0(E),
Wherein, I'(x, y) it is background gray levels, I " (x, y) is modulation gray value；
Wherein, φ1(x,y)、φ2(x, y) is phase information.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910163686.9A CN109870424B (en) | 2019-03-05 | 2019-03-05 | Hartmann ray tracing method based on color three-step phase shifting technology |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910163686.9A CN109870424B (en) | 2019-03-05 | 2019-03-05 | Hartmann ray tracing method based on color three-step phase shifting technology |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109870424A true CN109870424A (en) | 2019-06-11 |
| CN109870424B CN109870424B (en) | 2021-05-11 |
Family
ID=66919786
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910163686.9A Expired - Fee Related CN109870424B (en) | 2019-03-05 | 2019-03-05 | Hartmann ray tracing method based on color three-step phase shifting technology |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109870424B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116580076A (en) * | 2023-07-10 | 2023-08-11 | 中国空气动力研究与发展中心计算空气动力研究所 | Device, method, apparatus and medium for acquiring particle graph containing model boundary |
Citations (26)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP0402691A1 (en) * | 1989-05-29 | 1990-12-19 | Rainer Dipl.-Ing. Thiessen | Laser-sensor with external resonator |
| EP0871164A2 (en) * | 1997-04-10 | 1998-10-14 | TDK Corporation | Optical recording medium |
| EP1284418A2 (en) * | 2001-08-17 | 2003-02-19 | Fraunhofer-Gesellschaft Zur Förderung Der Angewandten Forschung E.V. | Method for optical analysis and optical detector device |
| JP2004028843A (en) * | 2002-06-27 | 2004-01-29 | Nikon Corp | Ray tracing method |
| CN1181334C (en) * | 1999-10-11 | 2004-12-22 | 因诺瓦蒂斯股份公司 | Method and device for characterizing culture liquid |
| CN101082577A (en) * | 2007-07-06 | 2007-12-05 | 云南大学 | Method for accurate measuring trace quantity liquid refractivity |
| US20080019252A1 (en) * | 2006-07-21 | 2008-01-24 | Fujifilm Corporation | Optical recording method, optical recording apparatus, optical recording medium, optical reproducing method, and optical reproducing apparatus |
| US20080043248A1 (en) * | 2006-08-18 | 2008-02-21 | Meric Ozcan | Photonic crystal sensors using band edge and/or defect mode modulation |
| CN101285855A (en) * | 2008-01-31 | 2008-10-15 | 西北工业大学 | Self focusing -type optical current mutual inductor |
| CN101769821A (en) * | 2010-02-04 | 2010-07-07 | 北京理工大学 | Lens refractive index and thickness measuring method and device based on differential confocal technology |
| JP2010186060A (en) * | 2009-02-12 | 2010-08-26 | Shin-Etsu Chemical Co Ltd | Patterning process |
| CN101819317A (en) * | 2009-02-27 | 2010-09-01 | 索尼公司 | Zoomar and imaging device |
| US20110134438A1 (en) * | 2009-12-07 | 2011-06-09 | Canon Kabushiki Kaisha | Refractive index distribution measuring method and refractive index distribution measuring apparatus |
| CN102183307A (en) * | 2011-03-17 | 2011-09-14 | 北京航空航天大学 | Method for acquiring thermal radiation noise of gradient temperature gas medium in infrared imaging guidance |
| CN102621097A (en) * | 2012-03-30 | 2012-08-01 | 常熟微纳激光光子技术有限公司 | Method for high-sensitivity measurement of surface optical nonlinearity of material |
| CN102975383A (en) * | 2012-10-18 | 2013-03-20 | 西安交通大学 | Preparation method for microwave device with gradient refractive indexes based on mixed liquid medium |
| US20150204723A1 (en) * | 2011-05-20 | 2015-07-23 | Beihang University | Optical phase device, method and system |
| CN105138759A (en) * | 2015-08-17 | 2015-12-09 | 哈尔滨工业大学 | Method for transmitting light in variable-refractive index flat side window |
| CN106501214A (en) * | 2016-10-21 | 2017-03-15 | 厦门大学嘉庚学院 | Complex refractivity index measuring method based on the absorbing medium of real refraction horn cupping |
| CN107401977A (en) * | 2017-08-15 | 2017-11-28 | 合肥工业大学 | The imaging compensating method of refraction deviation is considered in the measurement of high temperature binocular stereo vision |
| CN107832532A (en) * | 2017-11-14 | 2018-03-23 | 哈尔滨工业大学 | A kind of Pneumatic light transmission effects are to high-speed aircraft image quality numerical computation method |
| CN108225737A (en) * | 2017-12-06 | 2018-06-29 | 天津大学 | A kind of novel optical lens measurement method |
| CN108444877A (en) * | 2018-06-11 | 2018-08-24 | 浙江大学 | A kind of phase particle interference imaging method and device for measuring drop |
| CN108507495A (en) * | 2018-03-19 | 2018-09-07 | 中国计量大学 | A kind of free form surface detection method based on reverse Hartmann detection |
| CN109029932A (en) * | 2018-08-03 | 2018-12-18 | 长春理工大学 | A kind of unified full filed ray tracing method of refraction-reflection type telescopic system |
| CN109307480A (en) * | 2018-09-30 | 2019-02-05 | 中国计量大学 | A kind of transmissive element multi-surface surface testing method |
-
2019
- 2019-03-05 CN CN201910163686.9A patent/CN109870424B/en not_active Expired - Fee Related
Patent Citations (26)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP0402691A1 (en) * | 1989-05-29 | 1990-12-19 | Rainer Dipl.-Ing. Thiessen | Laser-sensor with external resonator |
| EP0871164A2 (en) * | 1997-04-10 | 1998-10-14 | TDK Corporation | Optical recording medium |
| CN1181334C (en) * | 1999-10-11 | 2004-12-22 | 因诺瓦蒂斯股份公司 | Method and device for characterizing culture liquid |
| EP1284418A2 (en) * | 2001-08-17 | 2003-02-19 | Fraunhofer-Gesellschaft Zur Förderung Der Angewandten Forschung E.V. | Method for optical analysis and optical detector device |
| JP2004028843A (en) * | 2002-06-27 | 2004-01-29 | Nikon Corp | Ray tracing method |
| US20080019252A1 (en) * | 2006-07-21 | 2008-01-24 | Fujifilm Corporation | Optical recording method, optical recording apparatus, optical recording medium, optical reproducing method, and optical reproducing apparatus |
| US20080043248A1 (en) * | 2006-08-18 | 2008-02-21 | Meric Ozcan | Photonic crystal sensors using band edge and/or defect mode modulation |
| CN101082577A (en) * | 2007-07-06 | 2007-12-05 | 云南大学 | Method for accurate measuring trace quantity liquid refractivity |
| CN101285855A (en) * | 2008-01-31 | 2008-10-15 | 西北工业大学 | Self focusing -type optical current mutual inductor |
| JP2010186060A (en) * | 2009-02-12 | 2010-08-26 | Shin-Etsu Chemical Co Ltd | Patterning process |
| CN101819317A (en) * | 2009-02-27 | 2010-09-01 | 索尼公司 | Zoomar and imaging device |
| US20110134438A1 (en) * | 2009-12-07 | 2011-06-09 | Canon Kabushiki Kaisha | Refractive index distribution measuring method and refractive index distribution measuring apparatus |
| CN101769821A (en) * | 2010-02-04 | 2010-07-07 | 北京理工大学 | Lens refractive index and thickness measuring method and device based on differential confocal technology |
| CN102183307A (en) * | 2011-03-17 | 2011-09-14 | 北京航空航天大学 | Method for acquiring thermal radiation noise of gradient temperature gas medium in infrared imaging guidance |
| US20150204723A1 (en) * | 2011-05-20 | 2015-07-23 | Beihang University | Optical phase device, method and system |
| CN102621097A (en) * | 2012-03-30 | 2012-08-01 | 常熟微纳激光光子技术有限公司 | Method for high-sensitivity measurement of surface optical nonlinearity of material |
| CN102975383A (en) * | 2012-10-18 | 2013-03-20 | 西安交通大学 | Preparation method for microwave device with gradient refractive indexes based on mixed liquid medium |
| CN105138759A (en) * | 2015-08-17 | 2015-12-09 | 哈尔滨工业大学 | Method for transmitting light in variable-refractive index flat side window |
| CN106501214A (en) * | 2016-10-21 | 2017-03-15 | 厦门大学嘉庚学院 | Complex refractivity index measuring method based on the absorbing medium of real refraction horn cupping |
| CN107401977A (en) * | 2017-08-15 | 2017-11-28 | 合肥工业大学 | The imaging compensating method of refraction deviation is considered in the measurement of high temperature binocular stereo vision |
| CN107832532A (en) * | 2017-11-14 | 2018-03-23 | 哈尔滨工业大学 | A kind of Pneumatic light transmission effects are to high-speed aircraft image quality numerical computation method |
| CN108225737A (en) * | 2017-12-06 | 2018-06-29 | 天津大学 | A kind of novel optical lens measurement method |
| CN108507495A (en) * | 2018-03-19 | 2018-09-07 | 中国计量大学 | A kind of free form surface detection method based on reverse Hartmann detection |
| CN108444877A (en) * | 2018-06-11 | 2018-08-24 | 浙江大学 | A kind of phase particle interference imaging method and device for measuring drop |
| CN109029932A (en) * | 2018-08-03 | 2018-12-18 | 长春理工大学 | A kind of unified full filed ray tracing method of refraction-reflection type telescopic system |
| CN109307480A (en) * | 2018-09-30 | 2019-02-05 | 中国计量大学 | A kind of transmissive element multi-surface surface testing method |
Non-Patent Citations (10)
| Title |
|---|
| DAODANG WANG: "Transmitted wavefront testing with large dynamic range based on computer-aided deflectometry", 《JOURNAL OF OPTICS》 * |
| HAOSU XIAO: "Imaging quality evaluation of aerodynamically heated optical dome using ray tracing", 《APPLIED OPTICS》 * |
| N. ANANDA KRISHNA: "Discrete transfer method applied to radiative transfer in a variable refractive index semitransparent medium", 《JOURNAL OF QUANTITATIVE SPECTROSCOPY & RADIATIVE TRANSFER》 * |
| OSCAR E. CASTILLO: "Measurement of phase objects by the use of color phase-shifting technique", 《INFRARED REMOTE SENSING AND INSTRUMENTATION XXI》 * |
| TOMÁŠ PIKÁLEK: "Air refractive index measurement using low-coherence interferometry", 《APPLIED OPTICS》 * |
| 刘东: "条纹投影与相位偏折测量技术研究进展", 《红外与激光工程》 * |
| 张斌: "复杂流场光偏折层析的理论与实验研究", 《中国博士学位论文全文数据库 基础科学辑》 * |
| 徐平: "基于相位偏折术的大像差透射波前检测", 《仪器仪表学报》 * |
| 杨欢: "气动光学效应对红外成像影响的仿真分析", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
| 黄战华: "变折射率介质中光线追迹通用算法的研究", 《光学学报》 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116580076A (en) * | 2023-07-10 | 2023-08-11 | 中国空气动力研究与发展中心计算空气动力研究所 | Device, method, apparatus and medium for acquiring particle graph containing model boundary |
| CN116580076B (en) * | 2023-07-10 | 2023-09-12 | 中国空气动力研究与发展中心计算空气动力研究所 | Device, method, apparatus and medium for acquiring particle graph containing model boundary |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109870424B (en) | 2021-05-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Goldstein et al. | Optical systems for flow measurement: shadowgraph, schlieren, and interferometric techniques | |
| Meier | Computerized background-oriented schlieren | |
| Steinberg et al. | Measurements of turbulent premixed flame dynamics using cinema stereoscopic PIV | |
| Elliott et al. | Molecular filter based planar Doppler velocimetry | |
| Mantzaras | Three-dimensional visualization of premixed charge engine flames | |
| CN102323445B (en) | Method for measuring flowing velocity of recycle gas in multiple laser cavity | |
| CN102175178B (en) | System and method for measuring diffusion flame frontal surface three-dimensional structure of motion fire source | |
| Payri et al. | The effect of nozzle geometry over ignition delay and flame lift-off of reacting direct-injection sprays for three different fuels | |
| CN109883994A (en) | The measurement method of non-uniform dielectric field based on Hartmann's ray tracing | |
| CN106768890B (en) | Gray scale cosine distribution optical target simulator for modulation transfer function detection | |
| CN109870424A (en) | Hartmann's ray tracing method based on colored three steps transposition technology | |
| CN110261072B (en) | Dynamic target capable of simulating imaging in motion | |
| Trolinger | Flow visualization holography | |
| CN101694404B (en) | Lateral shearing interference image-forming device | |
| CN109883993A (en) | The measurement method of non-uniform dielectric field based on Hartmann's ray tracing | |
| CN109883995A (en) | The measuring system and method for non-uniform dielectric field based on Hartmann's ray tracing | |
| CN109883996A (en) | The measurement method of non-uniform dielectric field based on Hartmann's ray tracing | |
| Gomez et al. | Time-resolved volumetric (4D) laser induced fluorescence imaging of primary spray breakup | |
| CN105931292B (en) | A kind of multi-direction moire tomography method based on affine calibration | |
| CN201993091U (en) | Measuring system for three-dimensional structure of diffusion flame frontal area of movement fire source | |
| Amjad et al. | Three-dimensional density measurements of a heated jet using laser-speckle tomographic background-oriented schlieren | |
| Boxx et al. | On the mechanism responsible for extreme turbulence intensity generation in the Hi-pilot burner | |
| Desantes et al. | Characterization of local fuel concentrations and internal dynamics in DI Diesel sprays | |
| Raffel et al. | Applications: Flows at Different Temperatures | |
| Ristić | A-a view in the invisible |
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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20210511 |