CN106840039B - A kind of quick three-dimensional measurement method based on composite grating projection - Google Patents

A kind of quick three-dimensional measurement method based on composite grating projection Download PDF

Info

Publication number
CN106840039B
CN106840039B CN201710141583.3A CN201710141583A CN106840039B CN 106840039 B CN106840039 B CN 106840039B CN 201710141583 A CN201710141583 A CN 201710141583A CN 106840039 B CN106840039 B CN 106840039B
Authority
CN
China
Prior art keywords
phase
grating
echelon
wrapped
composite grating
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
Application number
CN201710141583.3A
Other languages
Chinese (zh)
Other versions
CN106840039A (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.)
Hefei University of Technology
Original Assignee
Hefei University of Technology
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 Hefei University of Technology filed Critical Hefei University of Technology
Priority to CN201710141583.3A priority Critical patent/CN106840039B/en
Publication of CN106840039A publication Critical patent/CN106840039A/en
Application granted granted Critical
Publication of CN106840039B publication Critical patent/CN106840039B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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/25—Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures by projecting a pattern, e.g. one or more lines, moiré fringes on the object
    • G01B11/254—Projection of a pattern, viewing through a pattern, e.g. moiré

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Length Measuring Devices By Optical Means (AREA)

Abstract

The present invention relates to the quick three-dimensional measurement methods projected based on composite grating, comprising: builds measuring system；The generation of composite grating is carried out, and is burned onto optical projection system；Composite grating projects measured object surface, and modulation generates the composite grating of new shape；New shape composite grating is acquired, and is transmitted in computer processing system；By five step phase shift method solution phases, cosine grating wrapped phase and echelon wrapped phase are obtained respectively；Echelon wrapped phase is subjected to rounding processing, and combines cosine grating wrapped phase that wrapped phase is launched into absolute phase；The reconstruction of object dimensional information is carried out by absolute phase non-linear relation corresponding with measured object depth.The present invention only needs to carry out the projection and acquisition of one group of five amplitude grating picture, and can be carried out the expansion of object phase by a five step phase shift methods；Method proposed by the present invention reduces the use of grating picture number, improves measuring speed with identical measurement accuracy.

Description

A kind of quick three-dimensional measurement method based on composite grating projection
Technical field
The present invention relates to quick three-dimensional field of measuring technique, especially a kind of quick three-dimensional based on composite grating projection is surveyed Amount method.
Background technique
In recent years, optical three-dimensional measurement technology because it has high-precision, with testee non-contacting advantage due to obtain it is fast The development of speed, with the development of computer technology, information technology and digitizing technique, by optical projection system, industrial camera and calculating The projecting structural optical three-dimension measuring system of machine processing system composition comes into being.
Projection measurement method based on Gray code is produced, based on the decoded throwing of neighborhood by many years research of numerous scholars Image measuring method, digital projection technique method based on phase etc..Wherein, based on the digital projection technique of phase because of its measurement speed Degree is fast, is able to carry out fine and close reconstruction, relative accuracy is higher, therefore has biggish development potentiality.Because of phase and depth There are non-linear relations, so most important step is that acquisition is continuous exhausted for the digital projection technique based on phase To phase.The different digital rasters with phase information is projeced on testee first in terms of realization, is then passed through Camera carries out data acquisition to the phase grating of projection, and the absolute phase of object, Jin Erke are obtained by solving phase processor accordingly It is a kind of to have the advantages that multiple method for three-dimensional measurement to obtain the depth information of object by non-linear relation.
By the research of lot of domestic and foreign scholar, the method for obtaining absolute phase now has two major classes method substantially, a kind of Method is to carry out phase acquisition using space neighborhood information, and such methods error is excessive, can not carry out good practical application.More More research is for another kind of method, and this kind of methods are to carry out phase unwrapping using time-domain information, in camera shooting speed Under the premise of limited and guarantee measurement accuracy, present relatively good projection three-dimensional measurement method at least needs to project six amplitude gratings Picture just can be carried out the expansion of absolute phase.Less grating picture has been used there are also some, but since grating is sparse etc. Reason, so that the precision for obtaining phase method is relatively low.
Summary of the invention
The purpose of the present invention is to provide one kind under the premise of keeping precision, effectively raises the speed of three-dimensional measurement Based on composite grating projection quick three-dimensional measurement method.
To achieve the above object, the invention adopts the following technical scheme: a kind of quick three based on composite grating projection Measurement method is tieed up, this method includes the steps that following order:
(a) measuring system being made of optical projection system, industrial camera and computer processing system is put up；
(b) generation of composite grating is carried out in computer processing system, and the composite grating is passed through into computer disposal System is burned onto optical projection system；
(c) for the composite grating by projection system projects to measured object surface, composite grating is high by measured object surface Degree modulation generates the composite grating of new shape；
(d) the new shape composite grating described with industrial camera acquisition, and the composite grating of the new shape of acquisition is transmitted Into computer processing system；
(e) in computer processing system, the composite grating of the new shape by industrial camera acquisition passes through five step phase shifts Method solution phase obtains cosine grating wrapped phase and echelon wrapped phase respectively；
(f) the echelon wrapped phase that will acquire carries out rounding processing, and wraps up phase in conjunction with the cosine grating Wrapped phase is launched into absolute phase by position；
(g) weight of object dimensional information is carried out by absolute phase non-linear relation corresponding with measured object depth It builds.
In the step (a), the measuring system has synchronous Trigger Function, can acquire the light projected away moment Grid are simultaneously sent to computer processing system.
In the step (b), the composite grating is the single cosine grating and single ladder light by specific period ratio What grid were combined under single Color Channel, the calculation formula of intensity profile is as follows:
Ik(x, y)=I'(x, y)+I " (x, y) cos [φr(x,y)-r×kπ/5]+I”'(x,y)cos[φs(x,y)-s× kπ/5] (1)
Wherein, (x, y) indicates that pixel coordinate, I ' (x, y) indicate composite grating background intensity, and I " (x, y) indicates cosine light The modulate intensity of grid, I " ' (x, y) indicate the modulate intensity of echelon, φrIndicate that cosine grating wrapped phase, r indicate cosine The scale factor of grating phase-shift phase, φsIndicate that echelon wrapped phase, s indicate the scale factor of echelon phase-shift phase, Ik (x, y) expression kth width composite grating, k=0,1,2,3,4.
In the step (e), the calculation formula of the five steps phase shift method is as follows:
Wherein, (x, y) indicates pixel coordinate, φr(x, y) indicates that cosine grating wrapped phase, r indicate cosine grating phase shift The scale factor of amount, φs(x, y) indicates that echelon wrapped phase, s indicate the scale factor of echelon phase-shift phase, Ik(x, Y) expression kth width composite grating, k=0,1,2,3,4.
In the step (f), the calculation formula of the rounding processing are as follows:
K (x, y)=round [5 φs(x,y)/2π] (3)
Wherein, φsAfter (x, y) indicates that echelon wrapped phase, round indicate that rounding function, K (x, y) indicate rounding Echelon wrapped phase；
The calculation formula of the absolute phase is as follows:
Φw(x, y)=φr(x,y)+K(x,y)·2π (4)
Wherein, φr(x, y) indicates that cosine grating wrapped phase, K (x, y) indicate the echelon wrapped phase after rounding, Φw(x, y) indicates absolute phase.
As shown from the above technical solution, the present invention has the advantages that first, the present invention only needs to carry out one group of five amplitude grating The projection and acquisition of picture, and can be carried out by a five step phase shift methods expansion of object phase；Second, the present invention proposes Method compared with traditional preferable time domain unpacking method, with identical measurement accuracy, The present invention reduces light The use of grid picture number, improves measuring speed.
Detailed description of the invention
Fig. 1 is the composite grating figure that the present invention uses；
Fig. 2 (a) is the whole composite grating figure of SNR=22 interference, and Fig. 2 (b) is no SNR composite grating interfered and SNR The local contrast figure of the composite grating sum of=22 interference；
Fig. 3 is cosine wrapped phase figure；
Fig. 4 is ladder wrapped phase figure；
Fig. 5 is the rounding ladder wrapped phase figure of emulation；
Fig. 6 is the absolute phase figure of emulation；
Fig. 7 is measuring system figure；
Fig. 8 is the pictorial diagram of acquisition；
Fig. 9 is the wrapped phase figure of measured object；
Figure 10 is the absolute phase figure of measured object；
Figure 11 is measured object three-dimensional information rendering figure.
Specific embodiment
Embodiment one
The present embodiment is simulated measurement, and all processes are handled in a computer, and here is specific simulation process:
Step 1: the composite grating that the present invention is designed passes through MATLAB Software Create, the composite grating intensity profile of design It is shown below, shape is as shown in Figure 1.
I1(x, y)=I'(x, y)+I " (x, y) cos [φr(x,y)-4π/5]+I”'(x,y)cos[φs(x,y)-8π/5]
I2(x, y)=I'(x, y)+I " (x, y) cos [φr(x,y)-2π/5]+I”'(x,y)cos[φs(x,y)-4π/5]
I3(x, y)=I'(x, y)+I " (x, y) cos [φr(x,y)]+I”'(x,y)cos[φs(x,y)]
I4(x, y)=I'(x, y)+I " (x, y) cos [φr(x,y)+2π/5]+I”'(x,y)cos[φs(x,y)+4π/5]
I5(x, y)=I'(x, y)+I " (x, y) cos [φr(x,y)+4π/5]+I”'(x,y)cos[φs(x,y)+8π/5]
Wherein, (x, y) indicates pixel coordinate, I'(x, y) indicate raster back color intensity, I " (x, y) indicates cosine grating Modulate intensity, I " ' (x, y) indicate the modulate intensity of echelon, φr(x, y) indicates cosine grating wrapped phase, φs(x,y) Indicate that echelon wrapped phase, the phase-shift phase of cosine grating are 2 π/5, the phase-shift phase of echelon is 4 π/5.
Step 2: the composite grating of generation being added to the interference of SNR=22, shown in image shape such as Fig. 2 (a)；SNR= The composite grating figure comparison that 22 composite grating figure is interfered with no SNR is as shown in Fig. 2 (b).
Step 3: the composite grating being interfered being handled with MATLAB software, is carried out at five step phase shift methods first Reason, treatment process are shown below:
Wherein, φr(x, y) indicates cosine grating wrapped phase, φs(x, y) indicates echelon wrapped phase, solution Cosine grating wrapped phase is as shown in figure 3, the echelon wrapped phase solved is as shown in Figure 4.
Step 4: echelon wrapped phase is subjected to rounding, shape such as Fig. 5 institute after rounding by process described in following formula Show.
K (x, y)=round [5 φs(x,y)/2π]
Wherein, round indicates rounding function, φsAfter (x, y) indicates that echelon wrapped phase, K (x, y) indicate rounding Echelon wrapped phase.
Step 5: by the echelon wrapped phase of rounding and cosine grating wrapped phase by process described in following formula into Row combines, and ultimately forms absolute phase, shape is as shown in Figure 6, it can be seen that the absolute phase solved in the case where SNR=22 Although there is a small amount of noise spot in position, overall, the absolute phase that solves or relatively good has extremely strong anti-dry Disturb ability.
Φw(x, y)=φr(x,y)+K(x,y)·2π
Wherein, Φw(x, y) indicates absolute phase, and K (x, y) indicates the ladder wrapped phase after rounding, φr(x, y) is indicated Cosine grating wrapped phase.
Embodiment two
The present embodiment is shooting measurement in kind, this measurement distance in kind is 500mm or so, detailed process is as follows institute Show:
Step 1: as shown in fig. 7, putting up the measurement system being made of optical projection system, industrial camera, computer processing system System, optical projection system, industrial camera are connected by data line and capture card with computer processing system.In this measuring system In, the LightCrafter 4500, industrial camera model IO of projector model TI (Texas Instrument) the company exploitation of selection The 2M360CL of company includes MATLAB programming software, operation projector and the software of industrial camera etc. in computer system.
Step 2: the composite grating of design is passed through into MATLAB Software Create, the composite grating gray scale of design such as following formula institute Show, shape is as shown in Figure 1.
I1(x, y)=I'(x, y)+I " (x, y) cos [φr(x,y)-4π/5]+I”'(x,y)cos[φs(x,y)-8π/5]
I2(x, y)=I'(x, y)+I " (x, y) cos [φr(x,y)-2π/5]+I”'(x,y)cos[φs(x,y)-4π/5]
I3(x, y)=I'(x, y)+I " (x, y) cos [φr(x,y)]+I”'(x,y)cos[φs(x,y)]
I4(x, y)=I'(x, y)+I " (x, y) cos [φr(x,y)+2π/5]+I”'(x,y)cos[φs(x,y)+4π/5]
I5(x, y)=I'(x, y)+I " (x, y) cos [φr(x,y)+4π/5]+I”'(x,y)cos[φs(x,y)+8π/5]
Wherein, (x, y) indicates pixel coordinate, I'(x, y) indicate raster back color intensity, I " (x, y) indicates cosine grating Modulate intensity, I " ' (x, y) indicate the modulate intensity of echelon, φr(x, y) indicates cosine grating wrapped phase, φs(x,y) Indicate that echelon wrapped phase, the phase-shift phase of cosine grating are 2 π/5, the phase-shift phase of echelon is 4 π/5.
Step 3: the composite grating of generation is thrown into optical projection system, and to object by computer system burning It penetrates.
Step 4: being reflected by the object composite grating back with industrial camera acquisition, and synchronous transfer is to computer disposal system In system, the collected composite grating of industrial camera is as shown in Figure 8.
Step 5: in computer processing system, collected composite grating being passed through into five step phase shift methods, treatment process It is shown below:
Wherein, φr(x, y) indicates cosine grating wrapped phase, φs(x, y) indicates echelon wrapped phase, solves Cosine grating wrapped phase and echelon wrapped phase are as shown in Figure 9.
Step 6: echelon wrapped phase is subjected to rounding by process described in following formula.
K (x, y)=round [5 φs(x,y)/2π]
Wherein, Round indicates rounding function, φsAfter (x, y) indicates that echelon wrapped phase, K (x, y) indicate rounding The echelon wrapped phase.
Step 7: by the echelon wrapped phase of rounding and cosine grating wrapped phase by process described in following formula into Row combines, and ultimately forms absolute phase, shape is as shown in Figure 10.
Φw(x, y)=φr(x,y)+K(x,y)·2π
Wherein, Φw(x, y) indicates object absolute phase, and K (x, y) indicates the ladder wrapped phase after rounding, φr(x,y) Indicate cosine grating wrapped phase.
Step 8: object three is carried out by the absolute phase and the one-to-one non-linear relation of depth and calibration process The reconstruction of information is tieed up, the object dimensional pattern of reconstruction is as shown in figure 11.
In conclusion the present invention only needs to carry out the projection and acquisition of one group of five amplitude grating picture, and pass through five steps Phase shift method can be carried out the expansion of object phase；Method proposed by the present invention compared with traditional preferable time domain unpacking method, With identical measurement accuracy, The present invention reduces the uses of grating picture number, improve measuring speed.

Claims (1)

1. a kind of quick three-dimensional measurement method based on composite grating projection, this method include the steps that following order:
(a) measuring system being made of optical projection system, industrial camera and computer processing system is put up；
(b) generation of composite grating is carried out in computer processing system, and the composite grating is passed through into computer processing system It is burned onto optical projection system；
(c) composite grating passes through measured object apparent height tune by projection system projects to measured object surface, composite grating System generates the composite grating of new shape；
(d) the new shape composite grating described with industrial camera acquisition, and the composite grating of the new shape of acquisition is transmitted to meter In calculation machine processing system；
(e) in computer processing system, the composite grating of the new shape by industrial camera acquisition passes through five step phase shift method solutions Phase obtains cosine grating wrapped phase and echelon wrapped phase respectively；
(f) the echelon wrapped phase that will acquire carries out rounding processing, and will in conjunction with the cosine grating wrapped phase Wrapped phase is launched into absolute phase；
(g) reconstruction of object dimensional information is carried out by absolute phase non-linear relation corresponding with measured object depth；? In the step (a), the measuring system has synchronous Trigger Function, can acquire the grating projected away moment and be sent to Computer processing system；
In the step (b), the composite grating is existed by the single cosine grating and single echelon of specific period ratio It is combined under single Color Channel, the calculation formula of intensity profile is as follows:
Ik(x, y)=I'(x, y)+I " (x, y) cos [φr(x,y)-r×kπ/5]+I”'(x,y)cos[φs(x,y)-s×kπ/ 5] (1)
Wherein, (x, y) indicates that pixel coordinate, I ' (x, y) indicate composite grating background intensity, and I " (x, y) indicates cosine grating Modulate intensity, I " ' (x, y) indicates the modulate intensity of echelon, φrIndicate that cosine grating wrapped phase, r indicate cosine grating The scale factor of phase-shift phase, φsIndicate that echelon wrapped phase, s indicate the scale factor of echelon phase-shift phase, Ik(x,y) Expression kth width composite grating, k=0,1,2,3,4；
In the step (e), the calculation formula of the five steps phase shift method is as follows:
Wherein, (x, y) indicates pixel coordinate, φr(x, y) indicates that cosine grating wrapped phase, r indicate cosine grating phase-shift phase Scale factor, φs(x, y) indicates that echelon wrapped phase, s indicate the scale factor of echelon phase-shift phase, Ik (x, y) table Show kth width composite grating, k=0,1,2,3,4；
In the step (f), the calculation formula of the rounding processing are as follows:
K (x, y)=round [5 φs(x,y)/2π] (3)
Wherein, φs(x, y) indicates that echelon wrapped phase, round indicate that rounding function, K (x, y) indicate the ladder after rounding Grating wrapped phase；
The calculation formula of the absolute phase is as follows:
Φw(x, y)=φr(x,y)+K(x,y)·2π (4)
Wherein, φr(x, y) indicates that cosine grating wrapped phase, K (x, y) indicate the echelon wrapped phase after rounding, Φw (x, y) indicates absolute phase.
CN201710141583.3A 2017-03-10 2017-03-10 A kind of quick three-dimensional measurement method based on composite grating projection Active CN106840039B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710141583.3A CN106840039B (en) 2017-03-10 2017-03-10 A kind of quick three-dimensional measurement method based on composite grating projection

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710141583.3A CN106840039B (en) 2017-03-10 2017-03-10 A kind of quick three-dimensional measurement method based on composite grating projection

Publications (2)

Publication Number Publication Date
CN106840039A CN106840039A (en) 2017-06-13
CN106840039B true CN106840039B (en) 2019-03-26

Family

ID=59143417

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710141583.3A Active CN106840039B (en) 2017-03-10 2017-03-10 A kind of quick three-dimensional measurement method based on composite grating projection

Country Status (1)

Country Link
CN (1) CN106840039B (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110375671B (en) * 2019-02-19 2021-04-02 西安理工大学 Self-correcting method for phase error
TWI699525B (en) * 2019-05-21 2020-07-21 中華學校財團法人中華科技大學 Three-dimensional phase shift defect detection method and system
CN110360950B (en) * 2019-06-04 2021-07-06 重庆大学 Five-step phase sequence black-white digital raster coding method

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR100910573B1 (en) * 2008-05-29 2009-08-04 선문대학교 산학협력단 3-d measuring method by using color-multi wavelength reference phase
CN101788274B (en) * 2010-03-15 2012-08-08 电子科技大学 Method for 3D shape measurement of colourful composite grating
CN101806587B (en) * 2010-04-29 2012-12-12 浙江师范大学 Optical three-dimensional measurement method with absolute phase measurement
CN101881605B (en) * 2010-06-02 2011-11-30 南京航空航天大学 Optical three-dimensional measuring method based on phase coding technology
CN103759673B (en) * 2014-01-21 2016-07-06 南京理工大学 De-packaging method time phase based on double frequency three GTG sinusoidal grating fringe projection
CN105180840B (en) * 2015-10-14 2018-03-09 苏州大学 The measuring method and its device of heavy caliber concave mirror surface face shape based on the reflection of colored Contour fringes
CN106441161B (en) * 2016-11-01 2018-12-11 哈尔滨工程大学 A kind of fast phase unwrapping method based on period coding

Also Published As

Publication number Publication date
CN106840039A (en) 2017-06-13

Similar Documents

Publication Publication Date Title
CN106931910B (en) A kind of efficient acquiring three-dimensional images method based on multi-modal composite coding and epipolar-line constraint
CN108955571B (en) The method for three-dimensional measurement that double frequency heterodyne is combined with phase-shift coding
CN103791856B (en) Phase solving and unwrapping method based on four structured light images
CN109489585B (en) Three-dimensional measurement method based on improved multi-frequency stripe structured light
CN100520287C (en) Method and apparatus for measuring shape of an object
WO2016188068A1 (en) Method and system for stereoscopic vision three-dimensional measurement taking computing laser speckles as texture
CN109253708A (en) A kind of fringe projection time phase method of deploying based on deep learning
CN106840039B (en) A kind of quick three-dimensional measurement method based on composite grating projection
CN102261896A (en) Method and system for measuring three-dimensional shape of object based on phase measurement
CN108195313A (en) A kind of high dynamic range method for three-dimensional measurement based on Intensity response function
CN102184566A (en) Micro projector mobile phone platform-based portable three-dimensional scanning system and method
CN109141291A (en) A kind of fast phase unwrapping algorithm
CN105066906A (en) Fast high dynamic range three-dimensional measurement method
CN102183214A (en) Light detection method for large-aperture aspherical mirror structure
CN104567730B (en) A kind of method that space-time binary coding produces sinusoidal light field
CN107167093A (en) The combined type measuring system and measuring method of a kind of laser line scanning and shadow Moire
CN110174079B (en) Three-dimensional reconstruction method based on four-step phase-shift coding type surface structured light
CN108759721A (en) A kind of 3 D measuring method and device for projecting and reflecting based on optical stripe
CN105303572B (en) Based on the main depth information acquisition method passively combined
CN110207620A (en) A kind of three-dimensional rebuilding method determining digital fringe projection structure optical orders by different frequency
CN104482877A (en) Motion compensation method and system in three-dimensional imaging of dynamic object
CN109523627A (en) A kind of three-dimensional rebuilding method based on Taylor's exponential expression moulding surface structure light
CN106840038B (en) A kind of double frequency phase shift and Gray's code character close method
CN103940370A (en) Target object three-dimensional information acquisition method based on periodical co-prime hybrid coding
CN102798354B (en) Binary stripe stack based sinusoidal grating generation method

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