CN104111039A - Calibrating method for randomly placing fringe projection three-dimensional measuring system - Google Patents

Calibrating method for randomly placing fringe projection three-dimensional measuring system Download PDF

Info

Publication number
CN104111039A
CN104111039A CN201410389633.6A CN201410389633A CN104111039A CN 104111039 A CN104111039 A CN 104111039A CN 201410389633 A CN201410389633 A CN 201410389633A CN 104111039 A CN104111039 A CN 104111039A
Authority
CN
China
Prior art keywords
measuring system
scaling
scaling board
height
coordinate
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201410389633.6A
Other languages
Chinese (zh)
Other versions
CN104111039B (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.)
University of Electronic Science and Technology of China
Original Assignee
University of Electronic Science and Technology of China
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 University of Electronic Science and Technology of China filed Critical University of Electronic Science and Technology of China
Priority to CN201410389633.6A priority Critical patent/CN104111039B/en
Publication of CN104111039A publication Critical patent/CN104111039A/en
Application granted granted Critical
Publication of CN104111039B publication Critical patent/CN104111039B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Length Measuring Devices By Optical Means (AREA)

Abstract

The invention provides a calibrating method for randomly placing a fringe projection three-dimensional measuring system. The calibrating method comprises the following steps of: 1, designing a calibrating plate, wherein the calibrating plate is a chessboard pattern with equal intervals; 2, placing the calibrating plate into a measuring system view field for CCD (Charge Coupled Device) camera calibration; 3, placing the chessboard pattern on any different positions in a measurement field, projecting a phase shifting sine fringe pattern to the surface of the chessboard pattern on each placement position, and then collecting images and calculating the distribution phi of an unwrapped phase of each pixel point positioned on the surface of the chessboard pattern and the distance h from the unwrapped phase to a reference surface; 4 establishing three polynomials to express the relationship between the unwrapped phase phi and a height h, between the height h and a horizontal coordinate X and between the height h and a horizontal coordinate Y; 5 converting the horizontal coordinates (X, Y) which are from the unwrapped phase to actual height data and correspond to the height data by utilizing calibrated polynomial coefficient to complete the calibration of the three-dimensional measuring system.

Description

For putting arbitrarily the scaling method of fringe projection three-dimension measuring system
Technical field
The present invention relates to the scaling method of the optical three-dimensional measurement system based on fringe projection, be specially and adopt fringe projection commercial measurement body surface phase place, set up the relation of launching between phase place and height and height and horizontal ordinate, for the technical field of object dimensional surface shape measurement.
Background technology
Along with the fast development in the fields such as precision optical machinery processing, automobile making, precision optics processing, semicon industry, exact instrument manufacture, Design of Dies, reverse-engineering, people wish accurately to measure precision parts, workpiece object day by day urgently, traditional 2-D data measuring equipment can not meet the production requirement of modernization increasingly, robotization, and three-dimensional measurement technology is arisen at the historic moment.But high accuracy three-dimensional measuring technique be unable to do without system calibrating technology.
The internal and external parameter of current existing three-dimensional measuring systems calibration method or dependence calibration for cameras unit and projecting cell, calculate the relative position of the two and the geometric parameter of whole measuring system, replace manual measurement, this kind of method execution difficulty and the parameter degree of accuracy obtaining are not high yet, thereby finally cannot obtain high accuracy three-dimensional data; Or the three-dimensional parameter that relies on accurate standard component carries out matching, set up the relation of phase place and height, this kind method manufacturer's standard part is very difficult and cost is expensive; Or rely on accurately mobile and positioning and demarcating plate, set up the relation of phase place and height, simultaneously high to measuring system structural requirement, must ensure that camera optical axis is vertical with reference surface, thus nominal time length and be confined to complete in laboratory environment.And the three-dimensional information of pixel and altitude information formation, the true three-dimension information that cannot provide horizontal ordinate and altitude information to form all can only be provided above-mentioned three class scaling methods.In applicant's range of search, the typical scaling method of three-dimension measuring system can be seen following pertinent literature information:
" demarcation of 3-d shape measurement system " (Calibration of a three-dimensional shape measurement system that 1.Q.Y.Hu etc. deliver, Opt.Eng.42.487-493 (2003)), propose one projector has been equivalent to CCD camera, demarcated the internal and external parameter of CCD camera and projector simultaneously.But due to the image of CCD camera and digital projector lens distortion, this procedure is complicated and cannot provide system external parameter accurately, thereby the three dimension system precision of demarcating is not high.
" for the least square scaling method of fringe projection technology of profiling " (Least-squares calibration method for fringe projection profilometry that 2.H.W.Guo etc. deliver, Opt.Eng.44,033603 (2005)), study and utilized accurate mobile horizontal shift platform movement High Accuracy Flat to obtain nominal data, then data are carried out to least square fitting, obtain phase height relation.The accurate application restric-tion that moves horizontally platform this kind of method can only in laboratory, complete, be not suitable for practical application.
3.Z.H.Zhang Deng " a kind of simply, flexibly the phase calculation scaling method based on the 3-D imaging system " (Simple delivering, flexible calibration of phase calculation-based three-dimensional imaging system, Opt.Lett.36,1257-1259 (2011)), studied a kind of simple calibrating method that utilizes chessboard grid and demarcate blank, scaling board can at random, at random be placed on any position in visual field.But this method requires camera optical axis vertical with reference surface, has in actual applications limitation.
" a kind of flexible calibration technology based on fringe projection technology of profiling 3-D imaging system " (Flexible calibration technique for fringe-projection-based three-dimensional imaging that 4.Minh Vo etc. delivers, Opt.Lett.35,3192-3194 (2010)), study a kind of three-dimensional scaling method for any system architecture.The method requires low to the structure setting of measuring system, but also can only provide phase height data relationship, can not provide the scaling method of horizontal ordinate.
Can be seen and be had calibration technique or process complexity, be difficult to realize by above-mentioned document, or equipment cost costliness, the practicality of not having, or system structure design be required harsh, or to provide be that three-dimensional data is comprehensive not.Therefore, the full coordinate of fringe projection three-dimension measuring system is demarcated, and particularly demarcation in actual applications, is a difficult problem for a not yet fine solution.
Summary of the invention
For the deficiencies in the prior art, the technical matters that quasi-solution of the present invention is determined is, provide a kind of for putting arbitrarily the scaling method of fringe projection three-dimension measuring system, this demarcation is based on fringe projection technology and CCD camera calibration technology, can determine fast, accurately, simply the full coordinate relation between phase place and height, height and horizontal ordinate, particularly be applicable to the actual demarcation of engineering and use, and with low cost.
The present invention to achieve these goals, by the following technical solutions:
For putting arbitrarily a scaling method for fringe projection three-dimension measuring system, the method comprises:
Step 1, design are demarcated dull and stereotyped, and scaling board is equidistant alternate chessboard grid, and has diffuse reflection surface;
Step 2, place described scaling board in any diverse location of measurement field, gather the scaling board image of each placement location, then to several phase shift sine streaks figure of scaling board surface projection of each placement location, and gather several sine streak images of the scaling board reflection of each placement location, wherein scaling board image is for the inner parameter matrix A of calibration for cameras and the external parameter matrix [R of this position, T], sine streak image is used to phase shift algorithm and phase-unwrapping algorithm and solves the expansion PHASE DISTRIBUTION Φ of this placement location;
Step 3, the volume coordinate (X of each pixel (m, n) in camera coordinates system that calculates the some correspondence on the scaling board of position according to the inside and outside portion parameter matrix obtaining in step 2 c, Y c, Z c);
Step 4, in the diverse location described in step 2, choose arbitrarily a position as with reference to plane, and according to the volume coordinate (X of this position obtaining in step 3 c, Y c, Z c) matching three dimensions plane equation;
Step 5, obtain the volume coordinate of other position scaling boards according to step 3, utilize point to calculate the height value h of each point to reference surface to the space length formula of plane, and the horizontal ordinate (X, Y) of the vertical point of the vertical line equation computer memory point that utilizes space plane in reference planes;
Step 6, set up a polynomial expression, for expressing phase place Φ and the relation between h highly launched:
h ( m , n ) = 1 + a 1 ( m + n ) · Φ ( m , n ) a 2 ( m , n ) + a 3 ( m , n ) · Φ ( m , n )
In formula, h (m, n) is the height that different spaces that pixel (m, n) is corresponding is put reference planes, the expansion PHASE DISTRIBUTION that Φ (m, n) is this pixel, a 1(m, n), a 2(m, n), a 3(m, n) is multinomial coefficient to be solved, and launches gained height value h in phase place and step 5 utilize least square method to solve multinomial coefficient according to step 2 gained;
Step 7, set up a polynomial expression respectively, for expressing the relation between height h (m, n) and horizontal ordinate X, Y:
X = Σ k = 0 M b k ( m . n ) · h k ( m , n ) Y = Σ k = 0 M c k ( m , n ) · h k ( m , n )
In formula, (X, Y) is intersection point horizontal ordinate, and M is the maximal value of polynomial expression level time k, and is positive integer, b k(m, n), c k(m, n) is multinomial coefficient to be solved, and utilizes least square method to solve multinomial coefficient according to gained height value h in step 5 and horizontal ordinate (X, Y).
Step 8, utilize the multinomial coefficient of step 6 and step 7 gained, conversion launches phase place to actual three-dimensional data, completes the demarcation of three-dimension measuring system.
In technique scheme, the scaling board in described step 1 is equidistant alternate square chessboard grid, and has diffuse reflection surface.
In technique scheme, camera calibration in step 2, obtains inner parameter matrix A and external parameter matrix [R, T], and detailed process is as follows:
1), first scaling board is placed on to the multiple optional positions in viewing field of camera, and gathers the image of each position; 2), by the coordinate X of the angle point of each position scaling board medium square w, Y wwith corresponding pixel coordinate (m, n) substitution internal and external parameter compute matrix, calculate the external parameter matrix [R, T] of inner parameter matrix A and each position;
m n 1 = A × X C Y C Z C = A × [ R , T ] = A × R 11 R 12 T 1 R 21 R 22 T 2 R 31 R 32 T 3 × X W Y W 1
(X in formula c, Y c, Z c) be the volume coordinate in camera coordinate system, s is scale factor, m, n is pixel coordinate, R11, R12, R21, R22, R31, R32 is the element of R in external parameter matrix, T1, T2, T3 are the element of T in external parameter matrix.
In technique scheme, in step 3, calculate the volume coordinate (X of each pixel (m, n) in camera coordinates system c, Y c, Z c), detailed process is as follows:
1), first define matrix G, as follows
G = g 1 g 2 g 3 g 4 g 5 g 6 g 7 g 8 g 9 = [ R , T ] - 1 × A - 1 = R 11 R 12 T 1 R 21 R 22 T 2 R 31 R 32 T 3 - 1 × A - 1 ,
2), by matrix below pixel (m, n) substitution, obtain volume coordinate (X c, Y c, Z c).
X C Y C Z C = 1 mg 7 + ng 8 + g 9 × A - 1 × m n 1
In technique scheme, in described step 4, adopt the three dimensions plane equation of least square fitting reference planes: a X C r + bY C r + Z C r + c = 0 .
Wherein for the volume coordinate of each point in reference planes in camera coordinate system, (a, b, c) is the equation coefficient for the treatment of matching, and (a, b, 1) also claims the normal vector of reference surface.
In technique scheme, the bar graph gathering in described step 2 can be expressed as:
I(X,Y)=A(X,Y)+B(X,Y)cos[Φ(X,Y)]
Wherein I (X, Y) light distribution of expression collected by camera, A (X, Y) be background light intensity, B (X, Y) is that degree of modulation distributes, Φ (X, Y) be and the PHASE DISTRIBUTION of body surface height correlation, (X, Y) is the two-dimentional horizontal ordinate of deforming stripe image.
The present invention possesses following beneficial effect:
What the present invention was the most outstanding is propose first a kind of for putting arbitrarily the scaling method of fringe projection three-dimension measuring system.The method has advantages of simply, flexible, cost is low, can accurately provide the relation between PHASE DISTRIBUTION and object height data and altitude information and horizontal ordinate (X, Y) simultaneously.Disclosed by the invention for putting arbitrarily the scaling method experiments of measuring of fringe projection three-dimension measuring system, success, obtained efficiently and accurately the three-dimensional data information of object, the method has the advantages such as easy to operate, low cost, spirit and property are strong with respect to additive method.
Brief description of the drawings
Fig. 1 is that system arranges any fringe projection three-dimension measuring system figure;
Fig. 2 is that this measuring system is demarcated gridiron pattern pattern used;
Fig. 3 is the bar graph on the scaling board of digital camera record;
Fig. 4 is the deforming stripe figure of the face mask surface of digital camera record;
Fig. 5 is the expansion PHASE DISTRIBUTION of the face mask that demodulates;
Fig. 6 is the horizontal 2-D data figure of face mask;
Fig. 7 is the three-dimensional data figure of face mask;
Fig. 8 is the measurement result to 4mm position;
Fig. 9 is the measurement result to 16mm position;
Figure 10 is workflow diagram.
Embodiment
The invention discloses a kind of for putting arbitrarily the scaling method of fringe projection three-dimension measuring system, the method by fringe projection technology of profiling (FPP) and in conjunction with CCD camera calibration technology for fringe projection 3 d shape measure.Below in conjunction with the drawings and specific embodiments, working of an invention scheme is illustrated.
Step 1: experimental measurement system is set: system mainly comprises camera, projector, scaling board, object under test, control computer and support.The distribution of projector, camera, scaling board is random, and relative position is any, as shown in Figure 1.Adjustment System makes digital camera can photograph the candy strip on scaling board and object under test surface.
Step 2: scaling board is positioned over to the optional position in measuring system visual field.In each position, produce and project the sinusoidal phase shift candy strip of several standards by computer control projector, digital camera is taken the bar graph by demarcating surface reflection.Fig. 2 and Fig. 3 have shown the scaling board and the bar graph that in the method, use.
Step 3: to the bar graph of each position acquisition, the continuous phase that utilizes phase shift algorithm and space phase method of deploying to demodulate this scaling board surface, position distributes.The full coordinate that completes measuring system according to camera calibration technology and volume coordinate, distance relation is demarcated, and solves multinomial coefficient a 1(m, n), a 2(m, n), a 3(m, n) and b k(m, n), c kthe value of (m, n), creates the coefficient storage unit of corresponding each location of pixels, sets up the phase height of each pixel position, highly-horizontal ordinate relation, completion system staking-out work in pixel index mode.
Step 4: produced and projected the sinusoidal phase shift candy strip of several standards by computer control projector, digital camera is taken by the deforming stripe figure of determinand surface reflection.Utilize phase shift algorithm and space phase method of deploying to demodulate the continuous phase distributed intelligence on determinand surface.
Step 5: the phase height relation providing according to full coordinate scaling method obtains object under test surface elevation information, then provide height-horizontal ordinate relation according to full coordinate scaling method and obtain the horizontal ordinate that each pixel is corresponding.Finally complete the measurement of three-dimensional data.Fig. 5 is the expansion PHASE DISTRIBUTION demodulating, and this phase place and object height information are closely related.Fig. 6 represents that Fig. 7 represents the object Shape ' after reconstruction according to the three-dimensional data of demarcating gained phase height relation and height-horizontal ordinate relation and obtain body surface.Can find out from Fig. 6, Fig. 7, Fig. 8, Fig. 9, a kind of full coordinate scaling method for any fringe projection three-dimension measuring system that this patent proposes has higher accuracy and validity, and simple to operate.
Embodiment 1
For a full coordinate scaling method for any fringe projection three-dimension measuring system, it is characterized in that comprising following step:
A., experimental measurement system is set: system mainly comprises camera, projector, scaling board, object under test, control computer and support.The distribution of projector, camera, scaling board is random, and relative position is any.Adjustment System makes digital camera can photograph the candy strip on scaling board and object under test surface.
B. scaling board is positioned over to the multiple optional positions (in this example, scaling board being positioned over arbitrarily to 20 diverse location places) in measuring system visual field.In each position, produce and project the sinusoidal phase shift candy strip of 10 width standards by computer control projector, digital camera is taken the bar graph by demarcating surface reflection, as shown in Figure 3.This example adopts the white equidistant alternate square grid of 18 × 27 indigo plants, and the length of side of grid is 10mm, as shown in Figure 2.In this example, CCD camera is Prosilica GT1660C, and camera uses the tight shot (Computar M5018-MP2) of 50mm.
C. the bar graph to each position acquisition, utilizes 10 step phase shift algorithm and space phase method of deploying to demodulate the continuous phase distributed intelligence on this scaling board surface, position.The full coordinate that completes measuring system according to camera calibration technology and volume coordinate, distance relation is demarcated, and solves multinomial coefficient a 1(m, n), a 2(m, n), a 3(m, n) and b k(m, n), c kthe value of (m, n), creates the coefficient storage unit of corresponding each location of pixels, sets up the phase height of each pixel position, highly-horizontal ordinate relation, completion system staking-out work in pixel index mode.
D. produced and projected the sinusoidal phase shift candy strip of 10 width standards by computer control projector, digital camera is taken by the deforming stripe figure of determinand surface reflection.In this example, adopt face mask as object under test, as shown in Figure 4.The deforming stripe light intensity gathering can be expressed as:
I(X,Y)=A(X,Y)+B(X,Y)cos[Φ(X,Y)]
Wherein I (X, Y) light distribution of expression collected by camera, A (X, Y) be background light intensity, B (X, Y) is that degree of modulation distributes, Φ (X, Y) be and the PHASE DISTRIBUTION of body surface height correlation, (X, Y) is the two-dimentional horizontal ordinate of deforming stripe image.
E. the stripe pattern to gained, utilizes 10 step phase shift algorithm and space phase method of deploying to demodulate the continuous phase distributed intelligence Φ (X, Y) on determinand surface, as shown in Figure 5.
F. obtain, after continuous Φ (X, Y), utilizing the phase height relation of the demarcation in step C, can obtain the three-dimensional altitude information on object under test (face mask) surface.Step C solves coefficient a 1(m, n), a 2(m, n), a 3after (m, n), phase place height relationships is as follows:
h ( m , n ) = 1 + a 1 ( m + n ) · Φ ( m , n ) a 2 ( m , n ) + a 3 ( m , n ) · Φ ( m , n )
Here Φ (m, n)=Φ (X, Y), because in measuring system, the corresponding different lateral coordinates of each pixel.
G. calculating after altitude information h (m, n), utilizing height-horizontal ordinate relation of demarcating in step C, can obtain the horizontal 2-D data of object under test (face mask), as shown in Figure 6.Step C solves coefficient b k(m, n), c k(m, n) is rear, highly-horizontal ordinate relation is as follows:
X = Σ k = 0 M b k ( m . n ) · h k ( m , n ) Y = Σ k = 0 M c k ( m , n ) · h k ( m , n )
Here M is the maximal value of polynomial expression level time k, and we allow M=3.
The horizontal ordinate (X, Y) of the relevant position that the altitude information h (m, n) that H. comprehensive step F obtains and step G obtain, completes the three-dimensional measurement to object under test (face mask), as shown in Figure 7.
Embodiment 2
For the phase height relation of checking the inventive method demarcation and accuracy and the validity of height-horizontal ordinate relation, in measurement field, relative reference face position, accurately positioning and demarcating plate is in 4mm and two positions of 16mm.Its measurement result as shown in Figure 8 and Figure 9.Measurement result is got to average, be respectively 4.053mm and 16.057mm, its corresponding standard deviation is 0.028 and 0.033.To the measurement of horizontal ordinate, we adopt the spacing of grid angle point on scaling board as assessment.At 4mm place, the average of directions X is 9.968mm, and the average of Y-direction is 9.977mm, and its corresponding standard deviation is 0.062 and 0.075.At 16mm place, the average of directions X is 10.061mm, and the average of Y-direction is 10.070mm, and its corresponding standard deviation is 0.064 and 0.073.
To sum up, the present invention has higher degree of accuracy and validity, the visible last process flow diagram of its staking-out work process.

Claims (5)

1. for putting arbitrarily a scaling method for fringe projection three-dimension measuring system, the method comprises:
Step 1, design are demarcated dull and stereotyped, and scaling board is equidistant alternate chessboard grid, and has diffuse reflection surface;
Step 2, place described scaling board in any diverse location of measurement field, gather the scaling board image of each placement location, then to several phase shift sine streaks figure of scaling board surface projection of each placement location, and gather several sine streak images of the scaling board reflection of each placement location, wherein scaling board image is for the inner parameter matrix A of calibration for cameras and the external parameter matrix [R of this position, T], sine streak image is used to phase shift algorithm and phase-unwrapping algorithm and solves the expansion PHASE DISTRIBUTION Φ of this placement location;
Step 3, the volume coordinate (X of each pixel (m, n) in camera coordinates system that calculates the some correspondence on the scaling board of position according to the inside and outside portion parameter matrix obtaining in step 2 c, Y c, Z c);
Step 4, in the diverse location described in step 2, choose arbitrarily a position as with reference to plane, and according to the volume coordinate (X of this position obtaining in step 3 c, Y c, Z c) matching three dimensions plane equation;
Step 5, obtain the volume coordinate of other position scaling boards according to step 3, utilize point to calculate the height value h of each point to reference surface to the space length formula of plane, and the horizontal ordinate (X, Y) of the vertical point of the vertical line equation computer memory point that utilizes space plane in reference planes;
Step 6, set up a polynomial expression, for expressing phase place Φ and the relation between h highly launched:
h ( m , n ) = 1 + a 1 ( m + n ) · Φ ( m , n ) a 2 ( m , n ) + a 3 ( m , n ) · Φ ( m , n )
In formula, h (m, n) is the height that different spaces that pixel (m, n) is corresponding is put reference planes, the expansion PHASE DISTRIBUTION that Φ (m, n) is this pixel, a 1(m, n), a 2(m, n), a 3(m, n) is multinomial coefficient to be solved, and launches gained height value h in phase place and step 5 utilize least square method to solve multinomial coefficient according to step 2 gained;
Step 7, set up a polynomial expression respectively, for expressing the relation between height h (m, n) and horizontal ordinate X, Y:
X = Σ k = 0 M b k ( m . n ) · h k ( m , n ) Y = Σ k = 0 M c k ( m , n ) · h k ( m , n )
In formula, (X, Y) is intersection point horizontal ordinate, and M is the maximal value of polynomial expression level time k, and is positive integer, b k(m, n), c k(m, n) is multinomial coefficient to be solved, and utilizes least square method to solve multinomial coefficient according to gained height value h in step 5 and horizontal ordinate (X, Y).
Step 8, utilize the multinomial coefficient of step 6 and step 7 gained, conversion launches phase place to actual three-dimensional data, completes the demarcation of three-dimension measuring system.
2. according to claim 1 a kind of for putting arbitrarily the scaling method of fringe projection three-dimension measuring system, it is characterized in that: the scaling board in described step 1 is equidistant alternate square chessboard grid, and has diffuse reflection surface.
3. according to claim 1 a kind of for putting arbitrarily the scaling method of fringe projection three-dimension measuring system, it is characterized in that: camera calibration in step 2, obtain inner parameter matrix A and external parameter matrix [R, T], detailed process is as follows:
1), first scaling board is placed on to the multiple optional positions in viewing field of camera, and gathers the image of each position; 2), by the coordinate X of the angle point of each position scaling board medium square w, Y wwith corresponding pixel coordinate (m, n) substitution internal and external parameter compute matrix, calculate the external parameter matrix [R, T] of inner parameter matrix A and each position;
m n 1 = A × X C Y C Z C = A × [ R , T ] = A × R 11 R 12 T 1 R 21 R 22 T 2 R 31 R 32 T 3 × X W Y W 1
(X in formula c, Y c, Z c) be the volume coordinate in camera coordinate system, s is scale factor, m, n is pixel coordinate, R11, R12, R21, R22, R31, R32 is the element of R in external parameter matrix, T1, T2, T3 are the element of T in external parameter matrix.
4. according to claim 1 a kind of for putting arbitrarily the scaling method of fringe projection three-dimension measuring system, it is characterized in that: in step 3, calculate the volume coordinate (X of each pixel (m, n) in camera coordinates system c, Y c, Z c), detailed process is as follows:
1), first define matrix G, as follows
G = g 1 g 2 g 3 g 4 g 5 g 6 g 7 g 8 g 9 = [ R , T ] - 1 × A - 1 = R 11 R 12 T 1 R 21 R 22 T 2 R 31 R 32 T 3 - 1 × A - 1 ,
2), by matrix below pixel (m, n) substitution, obtain volume coordinate (X c, Y c, Z c).
X C Y C Z C = 1 mg 7 + ng 8 + g 9 × A - 1 × m n 1
5. according to claim 1 a kind of for putting arbitrarily the scaling method of fringe projection three-dimension measuring system, it is characterized in that: the three dimensions plane equation that adopts least square fitting reference planes in described step 4: a X C r + bY C r + Z C r + c = 0 .
Wherein for the volume coordinate of each point in reference planes in camera coordinate system, (a, b, c) is the equation coefficient for the treatment of matching, and (a, b, 1) also claims the normal vector of reference surface.
CN201410389633.6A 2014-08-08 2014-08-08 For arbitrarily putting the scaling method of fringe projection three-dimension measuring system Expired - Fee Related CN104111039B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410389633.6A CN104111039B (en) 2014-08-08 2014-08-08 For arbitrarily putting the scaling method of fringe projection three-dimension measuring system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410389633.6A CN104111039B (en) 2014-08-08 2014-08-08 For arbitrarily putting the scaling method of fringe projection three-dimension measuring system

Publications (2)

Publication Number Publication Date
CN104111039A true CN104111039A (en) 2014-10-22
CN104111039B CN104111039B (en) 2016-08-24

Family

ID=51707927

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410389633.6A Expired - Fee Related CN104111039B (en) 2014-08-08 2014-08-08 For arbitrarily putting the scaling method of fringe projection three-dimension measuring system

Country Status (1)

Country Link
CN (1) CN104111039B (en)

Cited By (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106767533A (en) * 2016-12-28 2017-05-31 深圳大学 Efficient phase three-dimensional mapping method and system based on fringe projection technology of profiling
CN107170010A (en) * 2017-05-11 2017-09-15 四川大学 System calibration method, device and three-dimensional reconstruction system
CN107610183A (en) * 2017-09-25 2018-01-19 天津大学 New striped projected phase height conversion mapping model and its scaling method
CN108362205A (en) * 2017-11-14 2018-08-03 沈阳工业大学 Space ranging method based on fringe projection
CN109405765A (en) * 2018-10-23 2019-03-01 北京的卢深视科技有限公司 A kind of high accuracy depth calculation method and system based on pattern light
CN109443214A (en) * 2018-12-19 2019-03-08 广东工业大学 A kind of scaling method of structured light three-dimensional vision, device and measurement method, device
WO2020016812A1 (en) * 2018-07-18 2020-01-23 Universidad De Cartagena 3d reconstruction method and device for measuring papules on the skin
CN111023995A (en) * 2019-11-18 2020-04-17 西安电子科技大学 Three-dimensional measurement method based on random two-frame phase shift fringe pattern
CN111043989A (en) * 2019-12-16 2020-04-21 电子科技大学 Sinusoidal fringe field projection module based on liquid crystal negative
CN112050751A (en) * 2020-07-17 2020-12-08 深圳大学 Projector calibration method, intelligent terminal and storage medium
CN113465548A (en) * 2021-08-23 2021-10-01 江南大学 Calibration and precision evaluation method of fringe projection three-dimensional measurement system
CN113566709A (en) * 2021-08-26 2021-10-29 苏州小优智能科技有限公司 Calibration method and device of structured light measurement system and electronic equipment
CN114018175A (en) * 2021-10-22 2022-02-08 湖南长步道光学科技有限公司 Real-time three-dimensional shape measurement method and system
CN114234847A (en) * 2021-12-08 2022-03-25 苏州恒视智能科技有限公司 Grating projection system and automatic correction and compensation method for grating phase shift height measurement
CN114459383A (en) * 2022-02-28 2022-05-10 嘉兴市像景智能装备有限公司 Calibration method based on sine stripe phase shift profilometry and implementation device
CN114485404A (en) * 2022-01-30 2022-05-13 嘉兴市像景智能装备有限公司 Calibration compensation method of height mapping based on path

Family Cites Families (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0617920B2 (en) * 1989-03-18 1994-03-09 株式会社島津製作所 Measurement data processor
JPH10281737A (en) * 1997-04-04 1998-10-23 Nikon Corp Method for synthesizing wave front
CN1155799C (en) * 2001-12-10 2004-06-30 西安交通大学 Laser linear scanning three-dimensional measurement double liquid knife virtual grid mapping calibrating method and equipment
CN101334270B (en) * 2008-07-25 2010-06-23 西安交通大学 Laser line scanning feeler geometric transformation calibration and curved face interpolation correcting method and apparatus
CN101718860B (en) * 2009-11-27 2012-01-11 电子科技大学 Non-cooperative target radio positioning method with unknown signal wave type
CN101788274B (en) * 2010-03-15 2012-08-08 电子科技大学 Method for 3D shape measurement of colourful composite grating
CN101949693B (en) * 2010-08-03 2012-01-11 河北工业大学 Method for calibrating three-dimensional imaging system
CN103868460B (en) * 2014-03-13 2016-10-05 桂林电子科技大学 Binocular stereo vision method for automatic measurement based on parallax optimized algorithm
CN103940371A (en) * 2014-05-12 2014-07-23 电子科技大学 High-precision three-dimensional shape measurement method for jump object

Cited By (24)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106767533A (en) * 2016-12-28 2017-05-31 深圳大学 Efficient phase three-dimensional mapping method and system based on fringe projection technology of profiling
CN106767533B (en) * 2016-12-28 2019-07-05 深圳大学 Efficient phase-three-dimensional mapping method and system based on fringe projection technology of profiling
CN107170010A (en) * 2017-05-11 2017-09-15 四川大学 System calibration method, device and three-dimensional reconstruction system
CN107610183B (en) * 2017-09-25 2020-12-15 天津大学 Calibration method of fringe projection phase height conversion mapping model
CN107610183A (en) * 2017-09-25 2018-01-19 天津大学 New striped projected phase height conversion mapping model and its scaling method
CN108362205A (en) * 2017-11-14 2018-08-03 沈阳工业大学 Space ranging method based on fringe projection
WO2020016812A1 (en) * 2018-07-18 2020-01-23 Universidad De Cartagena 3d reconstruction method and device for measuring papules on the skin
CN109405765A (en) * 2018-10-23 2019-03-01 北京的卢深视科技有限公司 A kind of high accuracy depth calculation method and system based on pattern light
CN109443214B (en) * 2018-12-19 2021-03-16 广东工业大学 Calibration method and device, measurement method and device for structured light three-dimensional vision
CN109443214A (en) * 2018-12-19 2019-03-08 广东工业大学 A kind of scaling method of structured light three-dimensional vision, device and measurement method, device
CN111023995A (en) * 2019-11-18 2020-04-17 西安电子科技大学 Three-dimensional measurement method based on random two-frame phase shift fringe pattern
CN111023995B (en) * 2019-11-18 2021-08-06 西安电子科技大学 Three-dimensional measurement method based on random two-frame phase shift fringe pattern
CN111043989A (en) * 2019-12-16 2020-04-21 电子科技大学 Sinusoidal fringe field projection module based on liquid crystal negative
CN112050751A (en) * 2020-07-17 2020-12-08 深圳大学 Projector calibration method, intelligent terminal and storage medium
CN112050751B (en) * 2020-07-17 2022-07-22 深圳大学 Projector calibration method, intelligent terminal and storage medium
CN113465548B (en) * 2021-08-23 2022-06-07 广东维正科技有限公司 Calibration and precision evaluation method of fringe projection three-dimensional measurement system
CN113465548A (en) * 2021-08-23 2021-10-01 江南大学 Calibration and precision evaluation method of fringe projection three-dimensional measurement system
CN113566709A (en) * 2021-08-26 2021-10-29 苏州小优智能科技有限公司 Calibration method and device of structured light measurement system and electronic equipment
CN114018175A (en) * 2021-10-22 2022-02-08 湖南长步道光学科技有限公司 Real-time three-dimensional shape measurement method and system
CN114234847A (en) * 2021-12-08 2022-03-25 苏州恒视智能科技有限公司 Grating projection system and automatic correction and compensation method for grating phase shift height measurement
CN114234847B (en) * 2021-12-08 2024-01-30 苏州恒视智能科技有限公司 Grating projection system and grating phase shift height measurement automatic correction compensation method
CN114485404A (en) * 2022-01-30 2022-05-13 嘉兴市像景智能装备有限公司 Calibration compensation method of height mapping based on path
CN114459383A (en) * 2022-02-28 2022-05-10 嘉兴市像景智能装备有限公司 Calibration method based on sine stripe phase shift profilometry and implementation device
CN114459383B (en) * 2022-02-28 2023-12-15 嘉兴市像景智能装备有限公司 Calibration method based on sine stripe phase shift profilometry

Also Published As

Publication number Publication date
CN104111039B (en) 2016-08-24

Similar Documents

Publication Publication Date Title
CN104111039A (en) Calibrating method for randomly placing fringe projection three-dimensional measuring system
CN103398660B (en) For obtaining the structured light vision sensor parameter calibration method of weld bead height information
CN100489446C (en) Method for measuring three-dimensional contour based on phase method
CN103292710B (en) A kind of distance measurement method applying binocular vision vision range finding principle
CN103267491B (en) The method and system of automatic acquisition complete three-dimensional data of object surface
CN103528543B (en) System calibrating method in a kind of optical grating projection three-dimensional measurement
CN104376558B (en) Cuboid-based intrinsic parameter calibration method for Kinect depth camera
CN101706263B (en) Three-dimensional surface measurement method and measurement system
CN101949693B (en) Method for calibrating three-dimensional imaging system
CN102032878B (en) Accurate on-line measurement method based on binocular stereo vision measurement system
CN103559735B (en) A kind of three-dimensional rebuilding method and system
Zhao et al. Calibration for stereo vision system based on phase matching and bundle adjustment algorithm
CN103530880A (en) Camera calibration method based on projected Gaussian grid pattern
CN104596439A (en) Speckle matching and three-dimensional measuring method based on phase information aiding
CN103606147B (en) Multiple stage is not total to visual field and measures the coordinate system conversion scaling method of camera
CN103559707B (en) Based on the industrial fixed-focus camera parameter calibration method of motion side's target earnest
CN104061879A (en) Continuous-scanning structured light three-dimensional surface shape perpendicular measuring method
TWI573984B (en) System and method of adjusting matching image
CN102564348A (en) Systematic geometric demarcation method for reflection three-dimensional measurement of stripe
CN105069789B (en) Structure light dynamic scene depth acquisition methods based on coding grid template
CN104266608A (en) Field calibration device for visual sensor and calibration method
Yang et al. Flexible and accurate implementation of a binocular structured light system
CN103942830B (en) Directly utilize and there is the method that the phase place of nonlinearity erron realizes scene three-dimensional reconstruction
CN102261896A (en) Method and system for measuring three-dimensional shape of object based on phase measurement
CN105403173A (en) Calibration mirror of light field deflection measuring system and application method of calibration mirror

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20160824