CN107367244B - A kind of optimal striped sequence method of design twined based on time phase solution - Google Patents
A kind of optimal striped sequence method of design twined based on time phase solution Download PDFInfo
- Publication number
- CN107367244B CN107367244B CN201710600109.2A CN201710600109A CN107367244B CN 107367244 B CN107367244 B CN 107367244B CN 201710600109 A CN201710600109 A CN 201710600109A CN 107367244 B CN107367244 B CN 107367244B
- Authority
- CN
- China
- Prior art keywords
- frequency
- phase
- phase shift
- fringe
- noise
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 45
- 230000010363 phase shift Effects 0.000 claims abstract description 89
- 238000005259 measurement Methods 0.000 claims abstract description 18
- 238000005316 response function Methods 0.000 claims abstract description 18
- 230000008569 process Effects 0.000 claims description 27
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 4
- 230000009467 reduction Effects 0.000 claims description 3
- 230000008859 change Effects 0.000 claims description 2
- 239000002184 metal Substances 0.000 claims 1
- 238000001314 profilometry Methods 0.000 description 5
- 239000000919 ceramic Substances 0.000 description 3
- 238000010587 phase diagram Methods 0.000 description 3
- 230000002123 temporal effect Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 238000012795 verification Methods 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B11/00—Measuring arrangements characterised by the use of optical techniques
- G01B11/24—Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures
- G01B11/2433—Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures for measuring outlines by shadow casting
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Cosmetics (AREA)
- Length Measuring Devices By Optical Means (AREA)
Abstract
The invention discloses a kind of optimal striped sequence methods of design twined based on time phase solution, initially set up the noise model of the three-dimensionalreconstruction under the slight out-of-focus appearance of projecting apparatus system, and the frequency response function of projecting apparatus system is measured, curve matching, according to the noise model of foundation and fit come frequency response function determine coding optimal fringe frequency；Then the threshold condition during time phase solution twines is derived, and the frequency frame for deriving different phase shift step numbers according to the threshold condition is chosen optimal phase shift step number according to the formula and combined than formula；Threshold condition during finally being twined according to time phase solution derives that low frequency auxiliary fringe frequency with the relational expression between high frequency fringes frequency, according to the relational expression, determines the frequency of stripeds at different levels one by one.The present invention is detected by the frequency response function to projecting apparatus system, finds the optimal fringe frequency f for encodingopt, to improve the precision of three-dimensional measurement.
Description
Technical Field
The invention belongs to the field of fringe projection profilometry, and particularly relates to a method for designing an optimal fringe sequence based on time phase unwrapping.
Background
In the field of fringe projection profilometry, phase unwrapping is a very important process, and provides a powerful tool for the transition of wrapped phase to absolute phase, thereby providing guarantee for subsequent three-dimensional reconstruction. The phase unwrapping algorithm can be broadly divided into spatial phase unwrapping and temporal phase unwrapping. Spatial-domain phase unwrapping has been proposed to be developed to date, and has been in the past fifty years, and although it has been improved, the limitations of the algorithm itself cannot be avoided: in the treatment of complex and discontinuous objects with surfaces, it becomes difficult to adapt (Saldner HO, Huntley JM. "Temporal phase unwarping: application to surface profiling of discrete objects." Appl Opt,1997,36: 2770-. In order to break through the limitation of spatial domain phase unwrapping, a time domain phase unwrapping algorithm is gradually proposed to become a large hot spot of the phase unwrapping algorithm. In the short last two decades, there are tens of time domain phase unwrapping algorithms proposed. These algorithms can be broadly divided into three categories: one cycle, number theory and heterodyne (Chao Zuo. "Temporal phase unwrapting analysis profiler: A compatationview"). Compared with other two methods, the single period method has the strongest capacity of resisting time noise.
Although the time domain phase unwrapping algorithm can well process objects with complex or discontinuous surfaces, the time domain phase unwrapping additionally requires many auxiliary stripes (Yan Hu. "Real-time microscopic 3-D shape Measurement based on optimized pulse-width-modulation reconstruction" 2017 ") in order to obtain a high-precision three-dimensional reconstruction result, which undoubtedly seriously affects the efficiency of three-dimensional Measurement. Therefore, how to reasonably design the stripe sequence and effectively reduce the number of auxiliary stripes is a very important and deeply studied problem. In addition, in the model established in the past, the phenomenon that the projector has slight defocusing is almost ignored. And slight defocusing of the projector can obviously inhibit the modulation degree of the high-density sinusoidal stripes, thereby reducing the signal-to-noise ratio of the phase. Therefore, it is necessary and meaningful to find a coding stripe with an optimal frequency to achieve a three-dimensional measurement result with high precision. However, in the time domain phase unwrapping process, how to select the fringe frequency sequence to stabilize the time domain phase unwrapping and effectively avoid the transmission of time noise has not been proposed yet. The above three problems are three problems which are very important in fringe projection profilometry, are closely related to the accuracy, efficiency and system robustness of three-dimensional measurement, but are technically ignored and cannot be solved.
Disclosure of Invention
The invention aims to provide a method for designing an optimal fringe sequence based on fringe projection profilometry time phase unwrapping, so as to improve the precision, efficiency and robustness of three-dimensional reconstruction in fringe projection profilometry.
The technical solution for realizing the purpose of the invention is as follows: a method for designing an optimal fringe sequence based on time phase unwrapping comprises the steps of firstly establishing a three-dimensional reconstruction noise model of a projector system in a slight out-of-focus state, measuring and curve fitting a frequency response function of the projector system, and determining the encoded optimal fringe frequency according to the established noise model and the fitted frequency response function;
then deducing a threshold condition in the time phase unwrapping process, deducing a frequency-frame ratio formula of different phase shift steps according to the threshold condition, and selecting an optimal phase shift step combination according to the formula;
and finally, deducing a relational expression between the frequency of the low-frequency auxiliary stripes and the frequency of the high-frequency stripes according to a threshold condition in the time phase unwrapping process, and determining the frequency of each level of stripes one by one according to the relational expression.
Compared with the prior art, the invention has the following remarkable advantages: (1) finding the optimal fringe frequency f for encoding by detecting the frequency response function of the projector systemoptTo improve the accuracy of the three-dimensional measurement. (2) And deducing frequency-frame ratio formulas of different phase shift steps according to the threshold condition in the time phase unwrapping process, and selecting the optimal phase shift step according to the formula so as to improve the efficiency of three-dimensional measurement. (3) And deducing a mathematical relation between the frequency of the low-frequency auxiliary stripes and the frequency of the high-frequency stripes according to a threshold condition in the time phase unwrapping process, and determining the frequency of each level of stripes one by one according to the mathematical relation so as to improve the stability of three-dimensional measurement.
The present invention is described in further detail below with reference to the attached drawing figures.
Drawings
FIG. 1 is a schematic flow chart of the method for designing an optimal fringe sequence according to the present invention.
Fig. 2 is a diagram of the verification process of fig. 1 for phase shift combination 2-2-3.
FIG. 3 is a graph of phase results of the present invention measured on a standard ceramic plate.
Fig. 4 is a graph showing the results of the measurement of the standard ball according to the present invention.
Detailed Description
The invention relates to a method for designing an optimal fringe sequence based on time phase unwrapping, which comprises the steps of firstly establishing a three-dimensional reconstruction noise model of a projector system in a slight out-of-focus state, measuring and curve fitting a frequency response function of the projector system, and determining the encoded optimal fringe frequency according to the established noise model and the fitted frequency response function; then deducing a threshold condition in the time phase unwrapping process, deducing a frequency-frame ratio formula with different phase shift steps according to the threshold condition, and selecting the optimal phase shift step according to the formula; and finally, deducing a relational expression between the frequency of the low-frequency auxiliary stripes and the frequency of the high-frequency stripes according to a threshold condition in the time phase unwrapping process, and determining the frequency of each level of stripes one by one according to the relational expression.
With reference to fig. 1, the specific steps for implementing the method are as follows:
step one, establishing a three-dimensional reconstruction noise model of the projector system in a slight out-of-focus state.
①, establishing a wrapped phase noise model under the focusing state of the projector, the noise variance of the wrapped phase (i.e. the wrapped phase noise model) solved by the standard N-step phase shift algorithm and the 2+ 1-step phase shift algorithm is as follows:
wherein,is the variance of the noise, σ, wrapping the phase2Is the noise variance of the fringe pattern captured by the camera, N is the phase shift step number, and B is the modulation of the fringe pattern captured by the camera. For the 2+1 step phase shift algorithm, N is the number of phase shift steps of the fringes that provide the average intensity.
②, establishing a wrapped phase noise model of the projector in the out-of-focus state, the projected fringe modulation degree B (f) is inhibited in the out-of-focus state of the projector:
wherein, B0Is the modulation degree, sigma, of the sinusoidal fringes projected by the projector in perfect focusgIs a coefficient proportional to the projector defocus level, f is the frequency of the encoded fringes, and h (f) is the frequency response function of the projector system.
Therefore, the constant B in equation (1) should be replaced by B (f) in equation (2):
③, a noise model of the three-dimensional reconstruction of the projector in the out-of-focus state is established, the phase noise does not change in the process of unwrapping the phase into the absolute phase, but after the absolute phase is mapped into the three-dimensional coordinate, the noise variance of the three-dimensional coordinate can be greatly reduced, and the reduction degree is determined by the frequency of the coding stripes:
wherein,is the noise variance of the three-dimensional coordinates. In the formula (4), σ2Is the noise variance of the fringe pattern acquired by the camera, whose magnitude depends only on the measurement environment and the measurement equipment, and N is the number of phase shift steps. Although increasing the number of phase shift steps reduces noise, this reduction is very inefficient and increasing the number of phase shift steps can severely impact the efficiency of the measurement. Therefore, it is the best choice to focus on increasing the frequency of the fringes. However, the higher the fringe frequency, the better, due to the projector defocus. To reduce noise to the maximum, the maximum of the product of the fringe frequency f and the fringe modulation b (f) should be found.
And step two, measuring and curve fitting are carried out on the frequency response function of the projector system, and the optimal fringe frequency of the code is determined according to the established noise model and the fitted frequency response function.
① projecting sinusoidal stripes with different frequencies by the projector, projecting stripes with the frequency of 1 to the standard ceramic flat plate by the projector,where w is the lateral or longitudinal resolution of the projector and is acquired synchronously with the camera. The expression for these sinusoidal stripes is:
wherein, InIs the light intensity of the fringe pattern collected by the camera, a is the average light intensity, x is the horizontal or vertical pixel coordinate of the camera, n is the phase shift index, and n is 0, 1, 2.
② the sinusoidal stripes are acquired synchronously by the camera and the stripe modulation degree is calculated from the gray scale distribution of the stripe pattern:
selecting the average value of the modulation degrees corresponding to the 200 × 200 pixel positions in the central areaTo eliminate errors.
③ fitting the calculated fringe modulation degrees to the 13 calculatedPerforming third-order polynomial fitting to obtainAnalysis formula (II)Where a, b, c, d are the fitted coefficients, and then plottedFinding the maximum value of the curveAnd selectValue range ofStripe frequency range [ f ] corresponding to intervalmin，fmax]Inner frequency as the optimal code stripe frequency fopt。
And step three, deducing a threshold condition in the time phase unwrapping process, deducing a frequency-frame ratio formula with different phase shift steps according to the threshold condition, and selecting the optimal phase shift step according to the formula.
①, a threshold condition is found in the time phase unwrapping process, in order to ensure the correctness of the phase unwrapping process, the fringe order error must be guaranteed to be 0, namely:
wherein, the [ alpha ], [ beta ]]Is the operator of nearest rounding, Δ khIs the fringe order error, fhIs the frequency of the high frequency fringes, flIs the frequency of the low-frequency auxiliary fringes, Δ ΦlIs the noise of the low-frequency auxiliary phase,is noise that high frequency wraps the phase. To ensure that the stripe order error is 0, then:
equation (8) can be further converted to:
to ensure the reliability of the result, we can approximate the maximum of the noise with 4.5 times the standard deviation, so we can get the threshold condition in the time phase unwrapping process:
wherein,is the standard deviation of the noise of the low frequency auxiliary phase,is the standard deviation of the noise of the high frequency phase. The fringe sequence is designed to achieve the highest accuracy that can be achieved by the measurement system with the least number of fringes, simply the frequency foptThe wrapped phase corresponding to the stripe is accurately unwrapped using a phase shift method with a minimum number of stripe frames.
② definition of frequency-to-frame ratio for convenience of description, we will use (f, N) to denote a fringe pattern with frequency f, phase-shift step number N,shows in a striped pattern (f)lAnd N) the highest frequency of the fringe pattern corresponding to the wrapping phase that can be unwrapped when the fringe is an auxiliary fringe. To better analyze the problem, the concept of frequency-frame ratio (FFR) is proposed. When the phase shift is N and the frequency is flIf the high frequency f can be expanded when the stripes of (2) are auxiliary stripeshThe wrap phase corresponding to the stripe of which phase shift step number is 3 cannot be spread outhThe wrapped phase corresponding to +1 is the ratioIs the so-called frequency-to-frame ratio.
③ derive the frequency-frame ratio formula for different phase shift steps assuming 3 phase shift steps and f frequencylWhen the fringe of (2) is used as an auxiliary fringe, the number of phase shift steps of the fringe corresponding to the unwrapped wrapping phase is 3, and the frequency isThen, according to the formula (3) and the formula (10), it can be obtained:
similarly, when the phase shift step number is N, the frequency is flAs auxiliary stripesThe number of phase shift steps corresponding to the unwrapped wrapped phase is 3, and the frequency isThen, according to the formula (3) and the formula (10), it can be obtained:
from equation (12), we can obtain:
it can thus be obtained that the auxiliary stripe frequency is flThe frequency-frame ratio when the phase shift step number is N is:
wherein, the upper mark flIndicating the frequency of the assist fringe and the subscript N indicating the number of phase shift steps of the assist fringe. Assuming that a continuous phase obtained by 2+1 step phase shift is used as the auxiliary phase, the frequency f of the auxiliary fringe can be obtained according to the above analysislThe frequency of the unwrapping phase which can be unwrapped by the two phase shifting steps and the corresponding frequency-frame ratio are as follows:
as can be seen from equations (15) and (16), for the standard N-step phase shift, N is 3, but the frequency-frame ratio of the 2+1 step phase shift is larger than that of the three step phase shift.
④, selecting the optimal phase shift step number combination, according to the definition of the frequency frame ratio, it can be known that the larger the frequency frame ratio is, the larger the corresponding phase shift step number contributes to the time phase unwrapping, therefore, we should try to select the continuous phase obtained by the 2+1 phase shifting to unwrappe the high frequency wrapping phase, but at least one group of stripe pattern with 3 phase shift steps is needed to provide the average light intensity for the 2+1 phase shifting.
TABLE 1 partial optimal stripe combinations
And step four, deducing a relational expression between the low-frequency auxiliary stripe frequency and the high-frequency stripe frequency according to a threshold condition in the time phase unwrapping process and determining the frequency of each level of stripes one by one according to the relational expression.
① relationship between the frequency of the low frequency assist stripes and the frequency of the high frequency stripes, when the high frequency stripes correspond to a phase shift of 2+1 steps, and the low frequency stripes correspond to a phase shift of N steps,heelWill be maximized, thus:
as can be seen from the phase shift step combinations in Table 1, N is ≦ 4 for the optimal phase shift step combination. In addition to this, the present invention is,the size of the range is related to the defocusing degree of the projector system, and the specific value range can be determined according to the detected frequency response function. Therefore, we can consider that,the upper limit of (d) is a constant related to the defocus level of the projector system, called PDC, which can be expressed as:
thus, equation (10) can be:
in general, for a slightly out-of-focus projector system,and in order to improve the stability of the system, twice the noise variance of the detected low-frequency auxiliary phase is used as the noise variance in the actual measurement, the formula (19) can be further converted into:
the equation (20) can be used as a relational equation that must be satisfied to determine the frequencies of the low and high frequency fringes.
② the optimal fringe frequency sequence is determined step by step according to equation (20) firstDetermining the optimal phase shift step combination by using the total frame number of the fringe pattern (selecting an optimal phase shift step combination from the table 1 according to the sequence of the frame numbers from small to large); then, according to the selected phase shift step number combination, gradually determining the fringe frequency of each stage, wherein the determining process is as follows: detecting noise standard deviation of monocycle phaseWill f isl1 and detectedSubstituting into equation (20), f is calculatedh(ii) a Then detecting the frequency fhNoise standard deviation of phase ofWill f isl＝fhAnd detectedSubstituting into equation (20), f is calculatedh(ii) a Detecting and calculating one by one according to the sequence until finally calculating the frequency f corresponding to the last group of stripeshiUntil now.
③ comparison fhiAnd foptIf f is large or smallhi≥foptIf so, then the selected phase shift step combination does not meet the requirement, and if not, then the next phase shift step combination is selected according to the sequence in step ② to determine the fringe frequency sequence.
FIG. 1 is a schematic flow chart of the method for designing an optimal fringe sequence according to the present invention. Fig. 2 is a diagram of the verification process of fig. 1 for phase shift combination 2-2-3.
Fig. 3 is the phase results measured on a standard ceramic plate. FIGS. 3(a) - (c) are absolute phase diagrams of the frequencies of each stage of the fringe frequency series {1,8,180 }. As can be seen from the figure, 212 error points appear in the last absolute phase. FIGS. 3(d) - (f) are absolute phase diagrams of each stage frequency of the fringe frequency series {1,15,180}, which are designed according to the present invention, and it can be seen that no error point appears in the absolute phase of each stage frequency. FIGS. 3(g) - (i) are absolute phase diagrams of the stage frequencies of the fringe frequency series {1,25,180 }. As can be seen from the figure, 138 error points occurred in the last absolute phase. From this comparison, it can be seen that the striped sequence designed according to the present invention has a higher stability.
Fig. 4 is a measurement result for a standard ball. FIG. 4 shows a standard sphere in the fringe frequency series {1,15,90}, where,
{1,15,180} and {1,15,256}, the first row is the measurement of the standard sphere, the second row is the calculated radius profile, and the third row is the radius error map. As can be seen from the figure, the fringe sequence {1,15,180} designed according to the present invention has the highest measurement accuracy.
Claims (3)
1. A method for designing an optimal fringe sequence based on time phase unwrapping is characterized by firstly establishing a three-dimensional reconstruction noise model of a projector system in a slight out-of-focus state, measuring and curve fitting a frequency response function of the projector system, and determining the encoded optimal fringe frequency according to the established noise model and the fitted frequency response function;
then deducing a threshold condition in the time phase unwrapping process, deducing a frequency-frame ratio formula of different phase shift steps according to the threshold condition, and selecting an optimal phase shift step combination according to the formula;
finally, according to the threshold condition in the time phase unwrapping process, deriving a relational expression between the low-frequency auxiliary stripe frequency and the high-frequency stripe frequency, and according to the relational expression, determining the frequency of each level of stripes one by one;
the process of establishing the three-dimensional reconstruction noise model of the projector system in the slight out-of-focus state is as follows:
①, establishing a wrapped phase noise model under the projector focusing state, namely solving the noise variance of the wrapped phase through a standard N-step phase shift algorithm and a 2+ 1-step phase shift algorithm as follows:
wherein,is the variance of the noise, σ, wrapping the phase2The noise variance of the fringe pattern acquired by the camera, N is the phase shift step number, B is the modulation degree of the fringe pattern captured by the camera, and for the 2+1 step phase shift algorithm, N is the phase shift step number of the fringe providing the average light intensity;
②, establishing a wrapped phase noise model in the out-of-focus state of the projector, wherein the projected fringe modulation degree B (f) is inhibited in the out-of-focus state of the projector:
wherein, B0Is the modulation degree, sigma, of the sinusoidal fringes projected by the projector in perfect focusgIs a coefficient proportional to the projector defocus level, f is the frequency of the encoded fringes, h (f) is the frequency response function of the projector system; thus, the constant B in equation (1) is replaced by B (f) in equation (2):
③, a three-dimensional reconstruction noise model under the defocusing state of the projector is established, the phase noise does not change in the process of wrapping phase unwrapping into absolute phase, but after the absolute phase is mapped into a three-dimensional coordinate, the noise variance of the three-dimensional coordinate can be greatly reduced, and the reduction degree is determined by the frequency of the coding stripes:
wherein,is the noise variance of the three-dimensional coordinates, in equation (4), σ2The noise variance of the fringe pattern acquired by the camera is only related to the measuring environment and the measuring equipment, N is the phase shift step number, and the formula (4) is an established noise model of the three-dimensional reconstruction of the projector in the defocusing state;
deducing a frequency-frame ratio formula of different phase shift steps according to the threshold condition, and selecting an optimal phase shift step combination according to the formula, wherein the specific process comprises the following steps:
①, a threshold condition is determined during the time phase unwrapping process, in order to ensure the accuracy of the phase unwrapping process, the fringe order error must be guaranteed to be 0, namely:
wherein, the [ alpha ], [ beta ]]Is the operator of nearest rounding, Δ khIs the fringe order error, fhIs the frequency of the high frequency fringes, flIs the frequency of the low-frequency auxiliary fringes, Δ ΦlIs the noise of the low-frequency auxiliary phase,is noise of high frequency wrapped phase, in order to ensure that the fringe order error is 0, then:
equation (8) further translates to:
to ensure the reliability of the results, the maximum value of the noise is approximated with 4.5 times the standard deviation, thus obtaining the threshold condition in the time-phase unwrapping process:
wherein,is the standard deviation of the noise of the low frequency auxiliary phase,is the standard deviation of the noise of the high frequency phase;
② definition of frequency-to-frame ratio for convenience of description, a fringe pattern having a frequency f, a phase shift step number N is expressed by (f, N),shows in a striped pattern (f)lN) the highest frequency of the fringe pattern corresponding to the wrapping phase that can be unwrapped when the fringe is an auxiliary fringe;
when the phase shift is N and the frequency is flIf the high frequency f can be expanded when the stripes of (2) are auxiliary stripeshThe wrap phase corresponding to the stripe of which phase shift step number is 3 cannot be spread outhThe wrapped phase corresponding to +1 is the ratioIs the frequency to frame ratio;
③ deriving the frequency-frame ratio formula of different phase shift steps, assuming that the phase shift step number is 3 and the frequency is flWhen the fringe of (2) is used as an auxiliary fringe, the number of phase shift steps of the fringe corresponding to the unwrapped wrapping phase is 3, and the frequency isThen, according to the formula (3) and the formula (10), we obtain:
similarly, when the phase shift step number is N, the frequency is flWhen the fringe of (2) is used as an auxiliary fringe, the number of phase shift steps corresponding to the unwrapped wrapping phase is 3, and the frequency isThen, according to the formula (3) and the formula (10), we obtain:
from equation (12):
thereby obtaining an auxiliary fringe frequency of flThe frequency-frame ratio when the phase shift step number is N is:
wherein, the upper mark flDenotes the frequency of the auxiliary fringe, and the subscript N denotes the phase shift step number of the auxiliary fringe;
assuming that a continuous phase obtained by 2+1 step phase shift is used as an auxiliary phase, the analysis is performed according to the above method to obtain an auxiliary phaseStripe frequency of flThe frequency of the unwrapping phase which can be unwrapped by the two phase shifting steps and the corresponding frequency-frame ratio are as follows:
as can be seen from the formula (15) and the formula (16), for the standard N-step phase shift, there is the largest frequency-frame ratio when N is 3, but the frequency-frame ratio of the 2+1 step phase shift is larger than that of the three step phase shift;
④, selecting the optimal phase shift step number combination, according to the definition of the frequency frame ratio, when the frequency frame ratio is larger, the corresponding phase shift step number will make larger contribution to the time phase expansion, so the continuous phase obtained by 2+1 phase shift is selected as much as possible to expand the high frequency wrapping phase, but at least one group of stripe pattern with 3 phase shift steps is needed to provide the average light intensity for 2+1 phase shift, thereby obtaining the optimal phase shift step number, the phase shift step number of the highest frequency stripe is 3 or 4, and the phase shift step numbers of other frequencies are all 2.
2. The method of claim 1, wherein the determining and curve fitting are performed on a frequency response function of the projector system, and the encoded optimal fringe frequency is determined according to the established noise model and the fitted frequency response function, and the method comprises the following steps:
① projecting sinusoidal stripes of different frequencies by a projector, i.e. projecting stripes of a frequency ofThe expression of these sinusoidal stripes is:
where w is the transverse or longitudinal resolution of the projector, InIs the light intensity of the fringe pattern collected by the camera, a is the average light intensity, x is the horizontal or vertical pixel coordinate of the camera, n is the phase shift index, and n is 0, 1, 2;
② the sinusoidal stripes are acquired synchronously by a camera to calculate the stripe modulation degree, and the modulation degree of the sinusoidal stripes is calculated from the gray distribution of the stripe pattern:
selecting the average value of the modulation degrees corresponding to the 200 × 200 pixel positions in the central areaTo eliminate errors;
③ fitting the calculated fringe modulation degrees, i.e. 13 calculatedPerforming third-order polynomial fitting to obtainAnalysis formula (II)Where a, b, c, d are the fitted coefficients, and then plottedFinding the maximum value of the curveAnd selectValue range ofStripe frequency range [ f ] corresponding to intervalmin，fmax]Inner frequency as the optimal code stripe frequency fopt。
3. The method for designing an optimal fringe sequence based on time phase unwrapping as claimed in claim 1, wherein a relationship between a frequency of a low frequency auxiliary fringe and a frequency of a high frequency fringe is derived according to a threshold condition during the time phase unwrapping, and the frequencies of the fringes at each level are determined one by one according to the relationship, wherein the specific process is as follows:
① relationship between the frequency of the low frequency assist stripes and the frequency of the high frequency stripes, when the high frequency stripes correspond to a 2+1 step phase shift, and the low frequency stripes correspond to an N step phase shift,heelWill be maximized, thus:
for the optimal phase shift step number combination, N is less than or equal to 4; in addition to this, the present invention is,the value of the projection system is related to the defocusing degree of the projection system, and the specific value range is determined according to the detected frequency response function; therefore, the temperature of the molten metal is controlled,the upper limit of (A) is a constant related to the degree of defocus of the projector system, referred to as PDC, expressed as:
thus, equation (10) is:
in order to improve the stability of the system, twice the noise variance of the detected low-frequency auxiliary phase is used as the noise variance in the actual measurement, and the formula (19) is further converted into:
the formula (20) is used as a relational expression for judging the satisfaction between the frequencies of the low frequency stripes and the high frequency stripes;
② determining the optimal fringe frequency sequence step by step according to the formula (20), firstly determining the optimal phase shift step combination according to the total frame number of the fringe pattern, and then determining the fringe frequency of each stage step by step according to the selected phase shift step combination by detecting the noise standard deviation of the single period phaseWill f isl1 and detectedSubstituting into equation (20), f is calculatedh(ii) a Then detecting the frequency fhNoise standard deviation of phase ofWill f isl＝fhAnd detectedSubstituting into equation (20), f is calculatedh(ii) a Detecting and calculating one by one according to the sequence until finally calculating the frequency f corresponding to the last group of stripeshiUntil the end;
③ comparison fhiAnd foptIf f is large or smallhi≥foptIf the selected combination of phase shift steps meets the requirements, otherwise if the selected combination of phase shift steps fails to meet the requirements, the total number of frames of the fringe pattern should be increased, and the combination of phase shift steps determined again in accordance with step ② to determine the fringe frequency sequence.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710600109.2A CN107367244B (en) | 2017-07-21 | 2017-07-21 | A kind of optimal striped sequence method of design twined based on time phase solution |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710600109.2A CN107367244B (en) | 2017-07-21 | 2017-07-21 | A kind of optimal striped sequence method of design twined based on time phase solution |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107367244A CN107367244A (en) | 2017-11-21 |
| CN107367244B true CN107367244B (en) | 2019-06-25 |
Family
ID=60307136
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710600109.2A Active CN107367244B (en) | 2017-07-21 | 2017-07-21 | A kind of optimal striped sequence method of design twined based on time phase solution |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107367244B (en) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111207694B (en) * | 2020-01-13 | 2021-06-04 | 南昌航空大学 | Three-dimensional measurement method combining double-step phase shift method with phase coding |
| CN113959364B (en) * | 2021-10-22 | 2022-08-16 | 华中科技大学 | Three-frequency phase unwrapping method and measuring device |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103759673A (en) * | 2014-01-21 | 2014-04-30 | 南京理工大学 | Time phase unwrapping method based on double frequency three-gray-scale sinusoidal grating fringe projection |
| CN103791856A (en) * | 2014-01-21 | 2014-05-14 | 南京理工大学 | Phase solving and unwrapping method based on four-width optical grating stripe image |
| CN105093893A (en) * | 2015-05-20 | 2015-11-25 | 北京信息科技大学 | Space-time three-dimensional phase unwrapping method and device |
| CN105890540A (en) * | 2016-04-08 | 2016-08-24 | 山东师范大学 | Digital image correlation-based object out-of-plane deformation phase measurement method |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US9074993B2 (en) * | 2010-09-03 | 2015-07-07 | The Arizona Board Of Regents On Behalf Of The University Of Arizona | White light achromatic grating imaging polarimeter |
-
2017
- 2017-07-21 CN CN201710600109.2A patent/CN107367244B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103759673A (en) * | 2014-01-21 | 2014-04-30 | 南京理工大学 | Time phase unwrapping method based on double frequency three-gray-scale sinusoidal grating fringe projection |
| CN103791856A (en) * | 2014-01-21 | 2014-05-14 | 南京理工大学 | Phase solving and unwrapping method based on four-width optical grating stripe image |
| CN105093893A (en) * | 2015-05-20 | 2015-11-25 | 北京信息科技大学 | Space-time three-dimensional phase unwrapping method and device |
| CN105890540A (en) * | 2016-04-08 | 2016-08-24 | 山东师范大学 | Digital image correlation-based object out-of-plane deformation phase measurement method |
Non-Patent Citations (5)
| Title |
|---|
| Absolute fringe order calculation using optimised multi-frequency selection in full-field profilometry;C.E. Towers, D.P. Towers, J.D.C. Jones;《Optics and Lasers in Engineering》;20050731;全文 |
| Absolute phase recovery of three fringe patterns with selected spatial frequencies;Yi Ding,Jiangtao Xi,Yanguang Yu,Fuqin Deng;《Optics and Lasers in Engineering》;20150731;第70卷;全文 |
| Color fringe projection system based on optimum frequency selection;Zonghua Zhang, Catherine Towers, David P. Towers;《SPIE Optics + Photonics》;20060831;全文 |
| Frequency selection in absolute phase maps recovery with two frequency projection fringes;Yi Ding等;《Optics Express》;20120604;第20卷(第12期);全文 |
| 结构光三维测量中相位相关技术研究;吴勇辉;《中国优秀硕士学位论文全文数据库信息科技辑》;20131231(第S1期);全文 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107367244A (en) | 2017-11-21 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106840036A (en) | A kind of diadactic structure light optimization method suitable for fast three-dimensional appearance measuring | |
| KR101475382B1 (en) | Method for extracting self adaptive window fourie phase of optical three dimensionl measurement | |
| CN107367244B (en) | A kind of optimal striped sequence method of design twined based on time phase solution | |
| WO2018040017A1 (en) | Method and system for correcting distortion of projector lens based on adaptive fringes | |
| CN103292734B (en) | Gamma value scaling method in phase measuring system | |
| CN104048621B (en) | A kind of digital phase-shifting technique three-dimensional measurement technology of profiling phase error detecting method based on constant phase coding | |
| CN106840038B (en) | A kind of double frequency phase shift and Gray's code character close method | |
| CN112880589B (en) | Optical three-dimensional measurement method based on double-frequency phase coding | |
| CN113639644A (en) | Domain mapping simple gamma calculation method for fringe projection profile measurement | |
| CN117011478B (en) | Single image reconstruction method based on deep learning and stripe projection profilometry | |
| Zhang et al. | Improving the extracting precision of stripe center for structured light measurement | |
| CN112945089A (en) | Structured light coding method based on stripe width modulation | |
| Hu et al. | A combined binary defocusing technique with multi-frequency phase error compensation in 3D shape measurement | |
| Yu et al. | High sensitivity fringe projection profilometry combining optimal fringe frequency and optimal fringe direction | |
| CN114240916A (en) | Multi-polarized light point cloud data fusion method and device for appearance state of steel rail | |
| CN105588518B (en) | Three-dimensional appearance acquisition methods based on double angle multi-frequency fringe projections and device | |
| Zappa et al. | Fourier-transform profilometry calibration based on an exhaustive geometric model of the system | |
| CN110428459A (en) | A method of the Phase- un- wrapping based on numerical order coding | |
| Deng et al. | Optimal pulse width modulation technique combined with stair phase-coding method for absolute phase retrieval with projector defocusing | |
| CN109506590A (en) | A kind of boundary jump phase error method for rapidly positioning | |
| Savio et al. | Optical properties of spectacle lenses computed by surfaces differential quantities | |
| Duan et al. | Quasi-periodic phase coding for long-depth-range 3-D measurements of moving objects | |
| García-Isáis et al. | One shot profilometry using phase partitions | |
| CN104125873A (en) | Misalignment calculation system | |
| Takacs | Minimizing artifacts in analysis of surface statistics |
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 |