CN109870424B - Hartmann ray tracing method based on color three-step phase shifting technology - Google Patents

Hartmann ray tracing method based on color three-step phase shifting technology Download PDF

Info

Publication number
CN109870424B
CN109870424B CN201910163686.9A CN201910163686A CN109870424B CN 109870424 B CN109870424 B CN 109870424B CN 201910163686 A CN201910163686 A CN 201910163686A CN 109870424 B CN109870424 B CN 109870424B
Authority
CN
China
Prior art keywords
formula
light
projection screen
color
space
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.)
Expired - Fee Related
Application number
CN201910163686.9A
Other languages
Chinese (zh)
Other versions
CN109870424A (en
Inventor
孔明
王道档
单良
赵军
许新科
刘维
郭天太
徐平
龚志东
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Jiliang University
Original Assignee
China Jiliang University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China Jiliang University filed Critical China Jiliang University
Priority to CN201910163686.9A priority Critical patent/CN109870424B/en
Publication of CN109870424A publication Critical patent/CN109870424A/en
Application granted granted Critical
Publication of CN109870424B publication Critical patent/CN109870424B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Abstract

The invention discloses a Hartmann ray tracing method based on a color three-step phase shifting technology, which comprises the steps of calculating an image by adopting a color three-step phase shifting method and obtaining phase information; according to the invention, a color three-step phase shifting method is adopted, a three-dimensional reconstruction algorithm of an inverse integral curve of a three-dimensional inhomogeneous medium field is combined, and accurate tracking of the measured light of the inhomogeneous medium field and measurement of the transient refraction characteristic of the refractive index of a three-dimensional space are realized through the combined arrangement of a projection screen and a telecentric optical system, so that the measurement accuracy and efficiency are greatly improved; and the measurement system of the inhomogeneous medium field based on Hartmann ray tracing has precise overall design, high measurement precision, lower cost and wide application range, has important theoretical significance and engineering application value, and is suitable for popularization and application.

Description

Hartmann ray tracing method based on color three-step phase shifting technology
Technical Field
The invention relates to a measuring method, in particular to a Hartmann ray tracing method based on a color three-step phase shifting technology.
Background
The nonuniform medium field refers to a spatial field with nonuniform spatial physical characteristics formed by nonuniform spatial distribution of factors such as pressure, speed, temperature, density, substances and the like, and the reproduction of physical quantities such as temperature, pressure, air flow and the like can be realized through the measurement of the nonuniform medium field. Therefore, the method has very important scientific significance and wide application prospect in scientific research such as fluid mechanics, heat transfer science and the like and in various industries such as energy, environment, power, chemical industry and the like for three-dimensional measurement of the inhomogeneous medium field.
Chemical reaction flow, laminar flow, turbulent flow, rarefied flow, low-temperature flow, high-temperature gas dynamics, supersonic speed and hypersonic flow, wind tunnel experiment and water tunnel flow field are typical heterogeneous medium fields, detection of the heterogeneous medium fields is accurately achieved, effective parameters can be provided for chemical reaction processes such as combustion and the like, design and control of aerodynamic equipment and the like, optimal control strategies are achieved, and the method is widely applied to the fields of industrial boilers, power plant coal, gas turbines, aviation, automobile engine combustion, ship propellers, spraying operation and the like. For example, when the method is applied to monitoring a combustion temperature field of a boiler, unreasonable combustion conditions of a hearth are main factors causing operation accidents of the boiler, and the reasonableness and the unreasonable combustion states directly influence the heat efficiency of the boiler and the proportion of power supply and fuel loss. Therefore, in order to keep the combustion field of the hearth to operate uniformly and stably and improve the combustion efficiency, the on-line monitoring of the temperature distribution of the combustion field in the hearth is realized, and the on-line monitoring is a necessary measure for regulating and controlling the combustion state; the method can also provide basic measurement data for the design and improvement of aerodynamic equipment, for example, the measurement of a feather-like jet plume ejected by a rocket engine in a high and vacuum state, the plume flow pattern defect can cause the swinging, the thrust weakening, the unstable flight and the incapability of flying according to a preset orbit of the aircraft, and finally the aircraft is probably not capable of reaching a preset height and area in the flying process; and the data of a real-time mixed flow field, such as the proportion of oil, water and gas in the process of oil and gas extraction, can be provided for the extraction of oil and gas, and the strategy relates to the strategy of oil and gas extraction, and is the basis of the oil and gas extraction.
At present, the general measurement method for the inhomogeneous medium field is mainly realized by measuring the change of the refractive index of the inhomogeneous medium field. When light penetrates through the refractive index field of the inhomogeneous medium, two changes occur, namely light deviates from the original direction, and light generates relative phase change and phase shift. The refractive index method determines the refractive index based on these two characteristics. The refractive index method has the irreplaceable advantages of other methods, has no interference to an original field, belongs to non-insertion measurement, has high time resolution, can be used for researching a very quick transient process, and even can track the process of space inhomogeneous medium field change; the spatial resolution is high, and the method can be used for researching non-uniform medium fields of the space with extremely small temperature gradient, such as shock waves, small flames and the like; the device can simultaneously provide measurement of a spatially inhomogeneous medium field, but not point measurement, and has a large amount of information. Other major methods include schlieren, shadow, optical, holographic, moire, laser speckle photography, laser speckle shearing involving talbot interference, and the like. The shadow method and the schlieren method are measured according to the deflection effect of light; the optical interference method and the holographic interference method directly measure the phase change, and assume that the light rays are not deflected; moire deflection, laser speckle photography, laser speckle shearing interference and the like relate to the deflection angle of light rays directly measured by Taibo interference.
The schlieren method utilizes the phenomenon that deflection can be generated when light propagates in media with different internal densities, and the deflected light is shielded by using a knife edge device, so that the throughput of the light is changed, and the light and shade change can be observed on a receiving surface such as a light screen or a camera to realize refractive index measurement. In addition to the conventional schlieren method, the schlieren method commonly used at present includes a focus schlieren method, a background schlieren method, a rainbow schlieren method, and the like. The principle of the shadow method is the same as that of the schlieren method, but the measurement of the refractive index is realized by recording the displacement of a point on a projection screen caused by the deflection of light. The experimental principle and the system structure of the schlieren method and the shadow method are simple, but the precision of the shadow method is relatively low, the traditional schlieren method is easily influenced by the background brightness, the displacement of each point cannot be accurately determined by the background schlieren method, and the measurement precision is limited.
In general, a mach-zehnder interferometer, a double-mirror interferometer, a michelson interferometer, a raman interferometer, a schlieren interferometer, or the like is used as the optical interference method. The light emitted by the same light source is divided into two beams, one beam is coherent with the other beam of reference light after passing through the measuring area, and the measurement of the spatial refractive index and the non-uniform medium field information of the reaction space are realized by analyzing the phase change of the light penetrating through the measuring area. The mach-zehnder interferometer, the double-mirror interferometer, the michelson interferometer, the raman interferometer, and the like all employ an amplitude spectroscopy, and the schlieren interferometer employs a polarization spectroscopy. The holographic interference method is similar to common optical interference method in principle, but it adopts time light splitting method to record the same light beam on the same holographic plate at different time, then makes these wave fronts produce interference again, and utilizes the analysis of phase change of light passed through the measuring region to implement measurement of space refractive index and reflect the nonuniform medium field information of space. Interferometry is currently widely used in heat transfer, combustion and aerodynamic research, and is also the most basic and classical method of thermophysical measurements.
The moire deflection method realizes the measurement of refractive index in the inhomogeneous medium field through moire interference phenomenon, the light that the light source sent places the grating around two, form two grating fringe superimposed moire fringes, under general conditions, moire fringe is parallel or according to the regular orderly distribution of setting for, when measuring space is inhomogeneous medium field, light can be because the change of refractive index takes place the deflection, moire fringe can produce deformation or displacement, the displacement of this fringe of analysis can obtain the deflection angle of light, thereby calculate the distribution of space inhomogeneous medium field. The moire deflexion method reduces the requirements of the interference method on mechanical, optical components and environmental control, has wider application range, but the diffraction phenomenon in the moire deflexion method reduces the contrast ratio and the spatial resolution of the fringes.
The talbot interference is measured using the talbot effect. When parallel light perpendicularly irradiates a grating with a period d, an identical image of the grating appears on a certain plane behind the grating, which is called self-imaging, namely a Talbot effect, and an imaging plane is called a Talbot surface. And another identical grating is placed on the stage-surface of the former grating, the latter grating is parallel to and coaxial with the former grating, the two gratings are positioned in the direction of a vertical axis identically, and a zero light field is observed behind the latter grating at the moment. When a non-uniform medium field of a phase object is inserted between the two gratings, light is disturbed, the self-imaging of the first grating generates distortion, the distorted self-imaging interacts with the second grating, interference fringes can be observed on an observation screen behind the second grating, and the interference fringes can be examined to obtain information of a light deflection angle, namely the Bobby interference. The Talbot interference is mostly used for optical instruments and parameter measurement, and has been developed for measurement of inhomogeneous medium fields. However, the processing precision requirement of the manufacturing of the above interferometric measuring devices is very high, and the adjustment of the devices is time-consuming and labor-consuming, so that the devices are generally expensive.
In laser speckle photography, laser beams with good coherence are diffusely reflected on the surface of an object to generate countless tiny point light sources, coherent light emitted by the point light sources is spread in space and is coherent with each other, bright spots and dark spots with different distributions are generated on projection equipment at different positions, when the laser beams penetrate through a space non-uniform medium field, patterns generated on the projection equipment move according to a certain rule, the motion information of the patterns is tracked, the displacement of each scattered spot is obtained, and the change of the refractive index of the space non-uniform medium field can be calculated in an inversion mode. The technique is insensitive to vibrations and does not require prior knowledge of the gradient direction of the inhomogeneous medium field. The point-by-point analysis of speckle photography is time consuming and cumbersome and requires correction for the fringe width due to the halo effect.
The laser speckle shearing interference images speckles on a CCD of a camera through double holes, and a crystal oscillation plane of the CCD is slightly deviated from a conjugate image plane, so that the speckles passing through the double holes form a slightly defocused image on the CCD, namely shearing, and the speckles are coherent to form a stripe pattern, namely the speckle shearing interference. The laser speckle shearing interference method measures the deflection angle of light, and the formed coherent fringes are equal displacement gradient lines of the light and represent the distribution condition of spatial nonuniformity, so that the laser speckle shearing interference method can perform online visual measurement on a spatial nonuniform medium field. Laser speckle shear interferometry is commonly used for the measurement of strain, tilt, bending and other mechanical parameters, and is also being developed for the measurement of inhomogeneous media fields. The method reduces the calculation amount of point-by-point analysis in laser speckle photography, but still needs to correct the width of the halo effect stripe.
The spatial measurement methods for the inhomogeneous medium field have advantages and disadvantages, the interference method has high requirements on equipment processing and mounting precision, so that the cost is higher, and the measurement precision of the schlieren method and the shadow method is lower. Therefore, the invention provides a low-cost and high-precision three-dimensional refractive index measuring method, which is based on the Hartmann ray tracing principle and realizes accurate tracing of measuring rays through the combined arrangement of a projection screen and a telecentric optical system, thereby realizing the measurement of the refractive index in a three-dimensional space.
Disclosure of Invention
The technical problem to be solved by the invention is to overcome the defects of the prior art and provide a Hartmann ray tracing method based on a color three-step phase shifting technology, which realizes accurate tracing of measured rays and measurement of three-dimensional space refractive index through the combined arrangement of a projection screen and a telecentric optical system.
In order to solve the technical problem, the invention is solved by the following technical scheme:
the Hartmann ray tracing method based on the color three-step phase shifting technology comprises the following specific contents:
each light spot on the projection screen is taken as a point light source to emit spherical incident light to the measurement space, wherein one spherical incident light which forms a deflection angle theta with the optical axis penetrates through a non-uniform medium field in the measurement space to be refracted to form deflection light parallel to the optical axis; deducing and obtaining the angle corresponding relation between incident light and refracted light through a formula (A), and deducing and obtaining the relation between a deflection angle theta and image position deviation shown in a formula (B) according to the formula (A), so that the relation between the integral of the refractive index gradient of the non-uniform medium field and the image displacement quantity shown in the formula (C) is obtained through the formula (A) and the formula (B); calculating the image by adopting a color three-step phase shifting method to obtain phase information, and determining the optical refraction characteristics and the corresponding relation of the image generated after the light spot on the projection screen is refracted by the non-uniform medium field through the phase information so as to realize light ray tracing; wherein, formula (A), formula (B) and formula (C) are respectively as follows:
Figure GDA0002967864540000041
Figure GDA0002967864540000042
Figure GDA0002967864540000043
wherein z is the direction of the optical axis,
Figure GDA0002967864540000044
the amount of change in the X, Y, Z axis, n is the spatial index of refraction,
Figure GDA0002967864540000045
is the change of space refractive index, F is the focal length of the lens, delta is the displacement of the pixel point of the image on the imaging screen of the CCD camera, S1The distance between the image space focus of the object space telecentric diaphragm and the imaging screen of the CCD camera is L, and the distance between the projection screen and the emitting end of the measuring space is L.
Preferably, the inhomogeneous medium flowing inside the measuring space generates a temperature field and/or a pressure field, and measurement of inhomogeneous medium fields due to temperature, pressure and flow is realized. The non-uniform medium field comprises a temperature field and/or a pressure field generated by a liquid-liquid mixed medium or a gas-liquid mixed medium or an air medium flowing in the measuring space, the measuring space is provided with a temperature adjusting device for adjusting the temperature field existing in the measuring space, the non-uniform medium field formed by temperature, pressure and flow is measured, and different non-uniform medium fields formed by the liquid-liquid mixed medium, the gas-liquid mixed medium and the air medium are measured. The liquid-liquid mixed medium is an oil-water mixed liquid, the gas-liquid mixed medium is a mixed liquid of natural gas and oil, the gas medium is an air flow or other mixed gas flow, and the pressure field is generated by liquid-liquid mixed flow or gas flow; the temperature regulating device comprises a heating temperature control system and a constant temperature control system, and a constant temperature field is provided by the temperature regulating device.
Preferably, the color three-step phase shifting method is based on a three-step phase shifting method, and comprises the following specific contents:
information I to be modulated1、I2、I3、I4、I5、I6The color RGB light source is respectively modulated into R, G, B three different channels of color RGB through brightness subareas, so that R, G, B three different channels of color RGB on a projection screen and a CCD camera respectively output and collect three phase images in the transverse direction and the longitudinal direction and obtain three phase information, and the CCD camera transmits the three phase information to the projection screen; calculating RGB three-channel information M (RGB) acquired by CCD camera by formula (D)]', RGB three-channel information M and actual RGB three-channel information N ═ RGB]' transformation matrix A and background matrix M0=[R0G0B0]' conversion is achieved by equation (E); finally, measuring transient refraction characteristics according to the actual RGB three-channel information N obtained through calculation; wherein the projection screen displays modulation information I of three cosine waves with different phase shifts, and three cosine waves with different phase shifts in the transverse direction and the longitudinal direction respectively1、I2、I3、I4、I5、I6Respectively obtained by calculation through a formula (F) and a formula (G); calculating three phase images collected by a CCD camera through a formula (H) and a formula (I) respectively to obtain phase information of each pixel, and comparing the obtained phase information of each pixel with three cosine wave stripes with different phase shifts displayed on a projection screen, so that the optical refraction characteristic and the corresponding relation of the images generated after light spots on the projection screen are refracted through an uneven medium field are determined, and light ray tracing is realized; the formulas (D), (E), (F), (G), (H), (I) are shown below:
Figure GDA0002967864540000051
M=A·N+M0 (E),
Figure GDA0002967864540000052
Figure GDA0002967864540000061
wherein I' (x, y) is a background gray value and I "(x, y) is a modulation gray value;
Figure GDA0002967864540000062
Figure GDA0002967864540000063
wherein phi is1(x,y)、φ2And (x, y) are phase information.
In the invention, the refractive index change data of the inhomogeneous medium field existing in the measurement space is calculated by adopting an inverse integral curve three-dimensional reconstruction algorithm; the inverse integral curve three-dimensional reconstruction algorithm comprises the following specific contents: establishing a ill-conditioned matrix with the refractive index distribution of the non-uniform medium field as a variable, and then solving the ill-conditioned matrix, thereby converting the reconstruction problem into a solving problem of the ill-conditioned matrix; wherein, the ill-conditioned matrix shown in formula (J) is obtained by deforming and dispersing formula (C):
T=WDS2 (J),
wherein T is a displacement information matrix of image pixels, represents the displacement of the measuring light compared with the uniform medium field through the non-uniform medium field, and is directly obtained from the corresponding relation of the image; w is a model coefficient matrix, and is obtained by direct body reconstructionThe model coefficient calculation method is obtained by adopting a Monte Carlo method, and the Monte Carlo method is regarded as the method for solving the probability that any point in a voxel space falls in a light beam, so that a large amount of random sampling can be directly carried out in the voxel space, the probability that a sampling point in a sample space falls in the light beam is further calculated according to the spatial position relation of the sampling point and the light beam, and when the number of samples is large enough, the probability value is regarded as a model coefficient; d is a refractive index change vector matrix to be solved; s2Is a projection matrix, projecting the three-dimensional trajectory onto the offset of the distance; the calculation method of the formula (J) can specifically adopt an additive algebra reconstruction, a multiplicative algebra reconstruction or a filtering back projection reconstruction algorithm to solve, so as to realize the reconstruction of the three-dimensional inhomogeneous medium field.
The Hartmann light tracing method based on the color three-step phase shifting technology is realized by a measuring system of a non-uniform medium field based on Hartmann light tracing, the measuring system of the non-uniform medium field based on the Hartmann light tracing comprises a measuring space, a projection screen and an object-side telecentric optical system, wherein the measuring space is internally provided with the non-uniform medium field, the projection screen is used for emitting spherical incident light to the measuring space, the object-side telecentric optical system is used for refracting the spherical incident light emitted by a light spot on the projection screen through the non-uniform medium field in the measuring space to form deflection light parallel to an optical axis, screening and converging the deflection light into an image, and the measuring space is positioned between the projection screen and the object-side telecentric optical; calculating the image by adopting a color three-step phase shifting method, obtaining phase information, and determining the optical refraction characteristic and the corresponding relation of the image according to the phase information; calculating to obtain the change data of the refractive index of the inhomogeneous medium field existing in the measuring space by adopting an inverse integral curve three-dimensional reconstruction algorithm; the basic measurement principle of the measurement system of the inhomogeneous medium field based on Hartmann ray tracing is as follows: each light spot on the projection screen is taken as a point light source to emit spherical incident light rays to the measuring space, wherein one of the spherical incident light rays which form a deflection angle theta with the optical axis penetrates through a non-uniform medium field in the measuring space to be refracted to form deflection light rays parallel to the optical axis, the deflection light rays which are emitted from the measuring space and are parallel to the optical axis are emitted to the object-side telecentric optical system, and the object-side telecentric optical system screens the deflection light rays and converges the deflection light rays into an image; calculating the image by adopting a color three-step phase shifting method, obtaining phase information, and determining the optical refraction characteristic and the corresponding relation of the image according to the phase information so as to realize light ray tracing; calculating to obtain the change data of the refractive index of the inhomogeneous medium field flowing in the measuring space by adopting an inverse integral curve three-dimensional reconstruction algorithm, thereby realizing the real-time measurement of the three-dimensional space refractive index of the inhomogeneous medium field existing in the measuring space; according to the invention, a color three-step phase shifting method is adopted, a three-dimensional reconstruction algorithm of an inverse integral curve of a three-dimensional inhomogeneous medium field is combined, and accurate tracking of the measured light of the inhomogeneous medium field and measurement of the transient refraction characteristic of the refractive index of a three-dimensional space are realized through the combined arrangement of a projection screen and a telecentric optical system, so that the measurement accuracy and efficiency are greatly improved; the system has the advantages of precise overall design, high measurement precision, lower cost and wide application range, has important theoretical significance and engineering application value, and is suitable for popularization and application.
The object space telecentric optical system comprises a lens for refracting the deflection light rays emitted from the measuring space, a CCD camera for shooting the converged deflection light rays and obtaining an image, and an object space telecentric diaphragm for converging the deflection light rays emitted from the lens on the CCD camera, wherein the lens is positioned between the measuring space and the object space telecentric diaphragm, the object space telecentric diaphragm is positioned between the lens and the CCD camera, and the projection screen, the measuring space, the lens, the object space telecentric diaphragm and the CCD camera are arranged in a manner of being opposite to each other in sequence by taking an optical axis as a central axis. In the invention, each light spot on a projection screen is taken as a point light source to emit spherical incident light to a measurement space, wherein one spherical incident light which forms a deflection angle theta with an optical axis penetrates through a non-uniform medium field in the measurement space to be refracted to form deflection light parallel to the optical axis, the deflection light parallel to the optical axis emitted from the measurement space is emitted to a lens, the lens refracts the deflection light parallel to the optical axis again, the deflection light after the secondary refraction is screened by an object-side telecentric diaphragm and converged on a CCD camera, and the CCD camera shoots the converged deflection light and obtains an image, namely: the incident light rays emitted by the light spot on the projection screen pass through the measuring space in sequence and are refracted twice by the lens, the deflected light rays are screened by the object space telecentric stop and are converged on the CCD camera, so that the CCD camera shoots the converged deflected light rays and obtains an image.
The measuring system of the inhomogeneous medium field based on the Hartmann ray tracing also comprises a measuring ring, a projection screen and a lens are both arranged on the inner wall of the measuring ring, the measuring space is a circular area concentrically arranged with the measuring ring, the radius of the circular area is smaller than that of the measuring ring, an optical axis passes through the measuring ring and the circle center of the circular area, and the projection screen and the lens are symmetrically arranged by taking the circle center as the center. The measuring ring and the circular area which are concentrically arranged further improve the accuracy and the efficiency of measurement.
The invention relates to a projection screen, a lens, an object space telecentric diaphragm and a CCD camera, which are composed of 4-8 groups, wherein the 4-8 groups of projection screen, lens, object space telecentric diaphragm and CCD camera are uniformly distributed along the 360-degree annular direction of a measuring ring to form a plurality of groups of rear projection Hartmann systems and the cross parts of measuring light sections thereof, each group of rear projection Hartmann systems measures and obtains the linear displacement of light rays generated in the measuring direction according to the light ray tracing, and calculates and obtains the change data of the refractive index of the non-uniform medium field in the measured three-dimensional space by combining with a corresponding inverse integral curve three-dimensional reconstruction algorithm according to the light ray tracing result in a plurality of measuring directions, thereby further improving the measuring accuracy and efficiency.
In the invention, the angle corresponding relation between the incident light and the refracted light is obtained by derivation through the formula (A), and the relation between the deflection angle theta and the image position deviation shown in the formula (B) is obtained by derivation according to the formula (A), so that the relation between the integral of the gradient of the refractive index of the inhomogeneous medium field and the image displacement shown in the formula (C) is obtained through the formula (A) and the formula (B).
In the invention, the CCD camera is connected with the projection screen, the CCD camera transmits the shot and obtained image to the projection screen, and the image is displayed on the projection screen.
According to the invention, the inhomogeneous medium field flowing in the measuring space generates a temperature field and/or a pressure field, and the measuring space is provided with a temperature adjusting device for adjusting the temperature field existing in the measuring space, so that the inhomogeneous medium field formed by temperature, pressure and flow is measured. Wherein the inhomogeneous medium field is oil-water mixed liquid, and the pressure field is generated by the inhomogeneous medium field flowing; the temperature regulating device comprises a heating temperature control system and a constant temperature control system, and a constant temperature field is provided by the temperature regulating device.
Due to the adoption of the technical scheme, the invention has the remarkable technical effects that: by adopting a color three-step phase shifting method, combining a three-dimensional reconstruction algorithm of a reverse integral curve of a three-dimensional inhomogeneous medium field and combining a projection screen and a telecentric optical system, accurate tracking of measured light of the inhomogeneous medium field and measurement of transient refraction characteristics of refractive indexes of three-dimensional space are realized, and measurement accuracy and efficiency are greatly improved; the system has the advantages of precise overall design, high measurement precision, lower cost and wide application range, has important theoretical significance and engineering application value, and is suitable for popularization and application.
Drawings
Fig. 1 is a schematic structural diagram of a measurement system using a conventional hartmann method.
FIG. 2 is a schematic structural diagram of a measurement principle embodiment of the Hartmann ray tracing-based inhomogeneous medium field measurement system of the present invention.
FIG. 3 is a schematic structural diagram of an embodiment of a rear projection Hartmann measuring system composed of 4 groups of projection screens, lenses, an object space telecentric diaphragm and a CCD camera.
FIG. 4 is a view showing a projection screen of the present invention showing a-2 π/3 phase image in the lateral direction.
Fig. 5 is a 0-phase image displayed in the lateral direction of the projection screen of the present invention.
FIG. 6 is a 2 π/3 phase image displayed in the transverse direction of the projection screen of the present invention.
FIG. 7 is a view showing a-2 π/3 phase image displayed in the longitudinal direction of the projection screen of the present invention.
Fig. 8 is a 0-phase image displayed in the longitudinal direction of the projection screen of the present invention.
FIG. 9 is a 2 π/3 phase image displayed in the longitudinal direction of the projection screen of the present invention.
Fig. 10 shows RGB color phase images synthesized from fig. 4, 5, 6, 7, 8 and 9 displayed on the projection screen according to the present invention.
Detailed Description
The present invention will be described in further detail with reference to the accompanying drawings and examples.
The Hartmann ray tracing method based on the color three-step phase shifting technology, as shown in FIGS. 2 and 3, includes the following concrete contents:
each light spot on the projection screen 2 is taken as a point light source to emit spherical incident light to the measuring space 3, wherein one spherical incident light which forms a deflection angle theta with the optical axis 1 penetrates through a non-uniform medium field in the measuring space 3 to be refracted to form deflection light parallel to the optical axis 1, the deflection light which is emitted from the measuring space 3 and is parallel to the optical axis 1 is emitted to the lens 4, the lens 4 refracts the deflection light parallel to the optical axis 1 again, the deflection light after the secondary refraction is screened through the object-side telecentric diaphragm 5 and converged on the CCD camera 6, and the CCD camera 6 shoots the converged deflection light and obtains an image; deducing and obtaining the angle corresponding relation between incident light and refracted light through a formula (A), and deducing and obtaining the relation between a deflection angle theta and image position deviation shown in a formula (B) according to the formula (A), so that the relation between the integral of the refractive index gradient of the non-uniform medium field and the image displacement quantity shown in the formula (C) is obtained through the formula (A) and the formula (B); calculating the image by adopting a color three-step phase shifting method to obtain phase information, and determining the optical refraction characteristics and the corresponding relation of the image generated after the light spot on the projection screen 2 is refracted by the non-uniform medium field through the phase information so as to realize light ray tracing; wherein, formula (A), formula (B) and formula (C) are respectively as follows:
Figure GDA0002967864540000091
Figure GDA0002967864540000092
Figure GDA0002967864540000101
wherein z is the direction of the optical axis 1,
Figure GDA0002967864540000102
the amount of change in the X, Y, Z axis, n is the spatial index of refraction,
Figure GDA0002967864540000103
f is the focal length of the lens 4, and delta is the displacement of the pixel point of the image on the imaging screen of the CCD camera 6, S1Is the distance between the image-side focal point of the object-side telecentric diaphragm 5 and the imaging screen of the CCD camera 6, and L is the distance between the projection screen 2 and the emergent end of the measuring space 3.
In this embodiment, the inhomogeneous medium flowing inside the measurement space 3 generates a temperature field and/or a pressure field, thereby realizing measurement of the inhomogeneous medium field formed by temperature, pressure and flow. The inhomogeneous medium field comprises a temperature field and/or a pressure field generated by a liquid-liquid mixed medium or a gas-liquid mixed medium or an air medium flowing in the measuring space 3, the measuring space 3 is provided with a temperature adjusting device for adjusting the temperature field existing in the measuring space, the inhomogeneous medium field formed by temperature, pressure and flow is measured, and different inhomogeneous medium fields formed by the liquid-liquid mixed medium, the gas-liquid mixed medium and the air medium are measured. The liquid-liquid mixed medium is an oil-water mixed liquid, the gas-liquid mixed medium is a mixed liquid of natural gas and oil, the gas medium is an air flow or other mixed gas flow, and the pressure field is generated by liquid-liquid mixed flow or gas flow; the temperature regulating device comprises a heating temperature control system and a constant temperature control system, and a constant temperature field is provided by the temperature regulating device
In this embodiment, the color three-step phase shifting method includes the following specific contents:
information I to be modulated1、I2、I3、I4、I5、I6The images are respectively modulated into R, G, B three different channels of color RGB through brightness subareas, so that R, G, B three different channels of color RGB on the projection screen 2 and the CCD camera 6 respectively output and collect three phase images in the transverse direction and the longitudinal direction and obtain three phase information, and the CCD camera 6 transmits the three phase information to the projection screen 2; calculating RGB three-channel information M ═ RGB acquired by CCD camera 6 by formula (D)]', RGB three-channel information M and actual RGB three-channel information N ═ RGB]' transformation matrix A and background matrix M0=[R0G0B0]' conversion is achieved by equation (E); finally, measuring transient refraction characteristics according to the actual RGB three-channel information N obtained through calculation; wherein the projection screen 2 displays modulation information I of three different phase-shifted cosine waves in the transverse and longitudinal directions1、I2、I3、I4、I5、I6The three cosine waves with different phase shifts in the transverse direction and the longitudinal direction are respectively shown in fig. 4, 5, 6, 7, 8 and 9, and the RGB color phase image synthesized by the images in fig. 4, 5, 6, 7, 8 and 9 is a 2 pi/3 pitch phase image shown in fig. 10; calculating three phase images collected by the CCD camera 6 through a formula (H) and a formula (I) respectively to obtain phase information of each pixel, and comparing the obtained phase information of each pixel with three cosine wave stripes with different phase shifts displayed on the projection screen 2, so that the optical refraction characteristic and the corresponding relation of the images generated after the light spots on the projection screen 2 are refracted through the uneven medium field are determined, and the light ray tracing is realized; the formulas (D), (E), (F), (G), (H), (I) are shown below:
Figure GDA0002967864540000111
M=A·N+M0 (E),
Figure GDA0002967864540000112
Figure GDA0002967864540000113
wherein I' (x, y) is a background gray value and I "(x, y) is a modulation gray value;
Figure GDA0002967864540000114
Figure GDA0002967864540000115
wherein phi is1(x,y)、φ2And (x, y) are phase information.
In the embodiment, the refractive index change data of the inhomogeneous medium field existing in the measuring space 3 is obtained by calculating through an inverse integral curve three-dimensional reconstruction algorithm; the inverse integral curve three-dimensional reconstruction algorithm comprises the following specific contents: establishing a ill-conditioned matrix with the refractive index distribution of the non-uniform medium field as a variable, and then solving the ill-conditioned matrix, thereby converting the reconstruction problem into a solving problem of the ill-conditioned matrix; wherein, the ill-conditioned matrix shown in formula (J) is obtained by deforming and dispersing formula (C):
T=WDS2 (J),
wherein T is the displacement information matrix of the image point, W is the model coefficient matrix, D is the refractive index change vector matrix to be solved, S2Is a projection matrix.
The Hartmann ray tracing method based on the color three-step phase shifting technology of the embodiment is realized by a measuring system of a non-uniform medium field based on Hartmann ray tracing, and the measuring system of the non-uniform medium field based on the Hartmann ray tracing comprises a measuring space 3 with the non-uniform medium field inside, a projection screen 2 for emitting spherical incident rays to the measuring space 3, and an object-side telecentric optical system for refracting the spherical incident rays emitted by a light spot on the projection screen 2 through the non-uniform medium field inside the measuring space 3 to form deflected rays parallel to an optical axis 1 for screening and converging into an image, wherein the measuring space 3 is positioned between the projection screen 2 and the object-side telecentric optical system; calculating the image by adopting a color three-step phase shifting method, obtaining phase information, and determining the optical refraction characteristic and the corresponding relation of the image according to the phase information; and calculating to obtain the change data of the refractive index of the inhomogeneous medium field existing in the measuring space 3 by adopting an inverse integral curve three-dimensional reconstruction algorithm.
In this embodiment, the object-side telecentric optical system includes a lens 4 for refracting the deflected light rays emitted from the measurement space 3, a CCD camera 6 for shooting the converged deflected light rays and obtaining an image, and an object-side telecentric stop 5 for screening the deflected light rays emitted from the lens 4 and converging on the CCD camera 6, wherein the lens 4 is located between the measurement space 3 and the object-side telecentric stop 5, the object-side telecentric stop 5 is located between the lens 4 and the CCD camera 6, the projection screen 2, the measurement space 3, the lens 4, the object-side telecentric stop 5, and the CCD camera 6 all use the optical axis 1 as a central axis and are arranged in sequence. In the invention, each light spot on a projection screen is taken as a point light source to emit spherical incident light to a measurement space, wherein one spherical incident light which forms a deflection angle theta with an optical axis penetrates through a non-uniform medium field in the measurement space to be refracted to form deflection light parallel to the optical axis, the deflection light parallel to the optical axis emitted from the measurement space is emitted to a lens, the lens refracts the deflection light parallel to the optical axis again, the deflection light after the secondary refraction is screened by an object-side telecentric diaphragm and converged on a CCD camera, and the CCD camera shoots the converged deflection light and obtains an image, namely: compared with the traditional Hartmann method shown in figure 1, the Hartmann screen 8 is replaced by the optical design of the combination of the projection screen, the lens, the object space telecentric diaphragm and the CCD camera, the track of imaging light is limited, and the quick tracing of light is realized.
In this embodiment, the system for measuring a non-uniform medium field based on hartmann ray tracing further includes a measuring ring 7, the projection screen 2 and the lens 4 are both disposed on the inner wall of the measuring ring 7, the measuring space 3 is a circular region concentrically disposed with the measuring ring 7, the radius of the circular region is smaller than that of the measuring ring 7, the optical axis 1 passes through the centers of the measuring ring 7 and the circular region, and the projection screen 2 and the lens 4 are symmetrically disposed with the center of the circle as the center.
In this embodiment, the projection screen 2, the lens 4, the object space telecentric stop 5, and the CCD camera 6 are composed of 4 groups, and the 4 groups of the projection screen 2, the lens 4, the object space telecentric stop 5, and the CCD camera 6 are uniformly distributed along the 360-degree annular direction of the measurement ring 7.
In this embodiment, the corresponding relationship between the incident light and the refracted light is obtained by derivation according to the formula (a), and the relationship between the deflection angle θ and the image position deviation shown in the formula (B) is obtained by derivation according to the formula (a), so that the relationship between the integral of the refractive index gradient of the inhomogeneous medium field and the image displacement shown in the formula (C) is obtained by the formula (a) and the formula (B).
In this embodiment, the CCD camera 6 is connected to the projection screen 2, and the CCD camera 6 transmits the image captured and obtained to the projection screen 2, and the image is displayed on the projection screen 2.
In this embodiment, the inhomogeneous medium field flowing inside the measurement space generates a temperature field and/or a pressure field, and the measurement space is provided with a temperature adjusting device for adjusting the temperature field existing inside the measurement space, thereby realizing measurement of the inhomogeneous medium field formed by temperature, pressure and flow. Wherein the inhomogeneous medium field is oil-water mixed liquid, and the pressure field is generated by the inhomogeneous medium field flowing; the temperature regulating device comprises a heating temperature control system and a constant temperature control system, and a constant temperature field is provided by the temperature regulating device.
In summary, the above-mentioned embodiments are only preferred embodiments of the present invention, and all equivalent changes and modifications made in the claims of the present invention should be covered by the claims of the present invention.

Claims (2)

1. The Hartmann ray tracing method based on the color three-step phase shifting technology is characterized in that: the method comprises the following specific contents:
each light spot on the projection screen (2) is taken as a point light source to emit spherical incident light rays to the measuring space (3), wherein one of the spherical incident light rays which form a deflection angle theta with the optical axis (1) penetrates through a non-uniform medium field in the measuring space (3) to be refracted to form a deflection light ray parallel to the optical axis (1), the deflection light ray parallel to the optical axis (1) emitted from the measuring space (3) is emitted to the lens (4), the lens (4) refracts the deflection light ray parallel to the optical axis (1) again, the deflection light ray after being refracted again is screened through the object space telecentric diaphragm (5) and converged on the CCD camera (6), and the CCD camera (6) shoots the converged deflection light ray and obtains an image; deducing and obtaining the angle corresponding relation between incident light and refracted light through a formula (A), and deducing and obtaining the relation between a deflection angle theta and image position deviation shown in a formula (B) according to the formula (A), so that the relation between the integral of the refractive index gradient of the inhomogeneous medium field and the image displacement shown in the formula (C) is obtained through the formula (A) and the formula (B); calculating the image by adopting a color three-step phase shifting method to obtain phase information, and determining the optical refraction characteristics and the corresponding relation of the image generated after the light spot on the projection screen (2) is refracted by the uneven medium field through the phase information so as to realize light ray tracing; wherein, formula (A), formula (B) and formula (C) are respectively as follows:
Figure FDA0002967864530000011
Figure FDA0002967864530000012
Figure FDA0002967864530000013
wherein z is the direction of the optical axis (1),
Figure FDA0002967864530000014
the amount of change in the X, Y, Z axis, n is the spatial index of refraction,
Figure FDA0002967864530000015
f is the focal length of the lens (4), delta is the displacement of the pixel point of the image on the imaging screen of the CCD camera (6), S is the change of the space refractive index1The distance between the image space focus of the object space telecentric diaphragm (5) and the imaging screen of the CCD camera (6) is L, and the distance between the projection screen (2) and the ejection end of the measurement space (3) is L;
the color three-step phase shifting method comprises the following specific contents:
information I to be modulated1、I2、I3、I4、I5、I6The images are respectively modulated into R, G, B three different channels of color RGB through brightness subareas, so that R, G, B three different channels of color RGB on a projection screen (2) and a CCD camera (6) respectively output and collect three phase images in the transverse direction and the longitudinal direction and obtain three phase information, and the CCD camera (6) transmits the three phase information to the projection screen (2); RGB three-channel information M (RGB) acquired by CCD camera (6) is calculated through formula (D)]', RGB three-channel information M and actual RGB three-channel information N ═ RGB]' transformation matrix A and background matrix M0=[R0G0B0]' conversion is achieved by equation (E); finally, measuring transient refraction characteristics according to the actual RGB three-channel information N obtained through calculation; wherein the projection screen (2) displays modulation information I of three differently phase shifted cosine waves in both the transverse and longitudinal directions1、I2、I3、I4、I5、I6Respectively obtained by calculation through a formula (F) and a formula (G); the CCD camera (6) is collected by a formula (H) and a formula (I) respectivelyThe three phase images are calculated to obtain the phase information of each pixel, and the obtained phase information of each pixel is compared with three cosine wave stripes with different phase shifts displayed on the projection screen (2), so that the optical refraction characteristic and the corresponding relation of the image generated after the light spot on the projection screen (2) is refracted through the uneven medium field are determined, and the light ray tracing is realized; the formulas (D), (E), (F), (G), (H), (I) are shown below:
Figure FDA0002967864530000021
M=A·N+M0 (E),
Figure FDA0002967864530000022
Figure FDA0002967864530000023
wherein I' (x, y) is a background gray value and I "(x, y) is a modulation gray value;
Figure FDA0002967864530000024
Figure FDA0002967864530000025
wherein phi is1(x,y)、φ2And (x, y) are phase information.
2. The Hartmann ray tracing method based on the color three-step phase shifting technology according to claim 1, characterized in that: the inhomogeneous medium flowing inside the measuring space (3) generates a temperature field and/or a pressure field.
CN201910163686.9A 2019-03-05 2019-03-05 Hartmann ray tracing method based on color three-step phase shifting technology Expired - Fee Related CN109870424B (en)

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 CN109870424A (en) 2019-06-11
CN109870424B true 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)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116580076B (en) * 2023-07-10 2023-09-12 中国空气动力研究与发展中心计算空气动力研究所 Device, method, apparatus and medium for acquiring particle graph containing model boundary

Family Cites Families (26)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE3917388C1 (en) * 1989-05-29 1990-11-29 Rainer 8000 Muenchen De Thiessen
JPH10289479A (en) * 1997-04-10 1998-10-27 Tdk Corp Optical recording medium
DE19949029C2 (en) * 1999-10-11 2002-11-21 Innovatis Gmbh Method and device for characterizing a culture fluid
DE10140399C2 (en) * 2001-08-17 2003-08-07 Fraunhofer Ges Forschung Optical analysis method and optical detector device
JP2004028843A (en) * 2002-06-27 2004-01-29 Nikon Corp Ray tracing method
JP4641981B2 (en) * 2006-07-21 2011-03-02 富士フイルム株式会社 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
CN100494984C (en) * 2007-07-06 2009-06-03 云南大学 Method for accurately measuring trace quantity liquid refractivity
CN101285855B (en) * 2008-01-31 2010-06-16 西北工业大学 Self-focusing type optical current mutual inductor
JP5206974B2 (en) * 2009-02-12 2013-06-12 信越化学工業株式会社 Pattern formation method
JP4770945B2 (en) * 2009-02-27 2011-09-14 ソニー株式会社 Variable focal length lens system and imaging apparatus
JP4968966B2 (en) * 2009-12-07 2012-07-04 キヤノン株式会社 Refractive index distribution measuring method and measuring apparatus
CN101769821A (en) * 2010-02-04 2010-07-07 北京理工大学 Lens refractive index and thickness measuring method and device based on differential confocal technology
CN102183307B (en) * 2011-03-17 2013-05-15 北京航空航天大学 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
CN102975383B (en) * 2012-10-18 2015-07-01 西安交通大学 Preparation method for microwave device with gradient refractive indexes based on mixed liquid medium
CN105138759B (en) * 2015-08-17 2018-02-06 哈尔滨工业大学 A kind of transmission method of light in variable refractivity flat board side window
CN106501214B (en) * 2016-10-21 2019-05-03 厦门大学嘉庚学院 The complex refractivity index measurement method of absorbing medium based on real refraction horn cupping
CN107401977B (en) * 2017-08-15 2019-04-30 合肥工业大学 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
CN108507495B (en) * 2018-03-19 2019-12-31 中国计量大学 Free-form surface detection method based on reverse Hartmann detection
CN108444877B (en) * 2018-06-11 2024-02-23 浙江大学 Phase particle interference imaging method and device for measuring liquid drops
CN109029932B (en) * 2018-08-03 2020-06-09 长春理工大学 Full-aperture full-view-field ray tracing method of catadioptric telescopic system
CN109307480B (en) * 2018-09-30 2020-08-14 中国计量大学 Method for detecting multi-surface shape of transmission element

Also Published As

Publication number Publication date
CN109870424A (en) 2019-06-11

Similar Documents

Publication Publication Date Title
CN102323445B (en) Method for measuring flowing velocity of recycle gas in multiple laser cavity
CN105043946A (en) Dual-wavelength based scattering angle self-calibration full-field rainbow measurement method and device
Tian et al. Aero-optical wavefront measurement technique based on BOS and its applications
CN110763159A (en) Optical deflection microscopic surface measuring device and method
CN109870424B (en) Hartmann ray tracing method based on color three-step phase shifting technology
CN109883994B (en) Hartmann ray tracing-based non-uniform medium field measuring method
CN114689170A (en) Device and method for measuring single-mode perfect vortex rotation with large topological charge value
CN109883996B (en) Hartmann ray tracing-based non-uniform medium field measuring method
CN1187588C (en) Space three-dimensional flow field real-time measurement system
CN109883993B (en) Hartmann ray tracing-based non-uniform medium field measuring method
CN109883995B (en) Hartmann ray tracing-based non-uniform medium field measurement system and method
CN106225730B (en) The big working distance autocollimation of portable combined zeroing high-precision laser and method
CN104792269A (en) Calculation method for optical fiber end face height insensitive to linear phase-shift errors
CN106323197B (en) The big working distance autocollimation of portable array zeroing high-precision laser and method
CN106017362B (en) A kind of big working distance autocollimation of portable high dynamic precision and method
Nakamura et al. Speckle beam-oriented schlieren technique
CN106225725B (en) The big working distance autocollimation of portable array zeroing laser and method
CN108413872A (en) Three-dimensional dimension precision measurement method based on Fabry-Perot multiple-beam interference
CN114518181A (en) Flame temperature distribution measuring system and method based on wavelength modulation shearing interference
CN112504164A (en) Measuring device and method capable of dynamically measuring surface shape of planar optical element
CN106091990B (en) The big working distance autocollimation of portable array zeroing high dynamic precision and method
CN106225726B (en) The big working distance autocollimation of array zeroing high-precision laser and method
CN106017363B (en) A kind of big working distance autocollimation of high dynamic precision and method
CN208012761U (en) A kind of wide spectrum high-resolution spectra interferometric measuring means
Grunwald et al. Orbital Angular Momentum Encoded Beam Tracking and Wavefront Sensing

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

Granted publication date: 20210511

CF01 Termination of patent right due to non-payment of annual fee