CN1302999A - Method for reconstructing 3D contour of digital projection based on phase-shifting method - Google Patents
Method for reconstructing 3D contour of digital projection based on phase-shifting method Download PDFInfo
- Publication number
- CN1302999A CN1302999A CN01104163A CN01104163A CN1302999A CN 1302999 A CN1302999 A CN 1302999A CN 01104163 A CN01104163 A CN 01104163A CN 01104163 A CN01104163 A CN 01104163A CN 1302999 A CN1302999 A CN 1302999A
- Authority
- CN
- China
- Prior art keywords
- phase
- projection
- projection pattern
- projector
- reconstructing
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 36
- 239000004744 fabric Substances 0.000 claims description 18
- 230000008569 process Effects 0.000 claims description 8
- 238000013500 data storage Methods 0.000 claims description 3
- 230000010363 phase shift Effects 0.000 abstract description 8
- 238000012545 processing Methods 0.000 abstract description 4
- 230000008901 benefit Effects 0.000 abstract description 2
- 230000000694 effects Effects 0.000 abstract description 2
- 230000008859 change Effects 0.000 description 7
- 238000003384 imaging method Methods 0.000 description 6
- 244000309464 bull Species 0.000 description 5
- 238000005259 measurement Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 230000003287 optical effect Effects 0.000 description 4
- 238000013461 design Methods 0.000 description 3
- 238000006073 displacement reaction Methods 0.000 description 3
- 238000001914 filtration Methods 0.000 description 3
- 230000000007 visual effect Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 210000001747 pupil Anatomy 0.000 description 2
- 230000001105 regulatory effect Effects 0.000 description 2
- 241001270131 Agaricus moelleri Species 0.000 description 1
- 238000011960 computer-aided design Methods 0.000 description 1
- 235000009508 confectionery Nutrition 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000012467 final product Substances 0.000 description 1
- 238000005286 illumination Methods 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 238000003754 machining Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000003908 quality control method Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
Images
Landscapes
- Length Measuring Devices By Optical Means (AREA)
- Image Processing (AREA)
Abstract
A method based on phase shifting for reconstructing 3D contour of digital projection includes such steps as pre-compiling projection pattern generation program, projecting the patterns onto the surface of object to be tested at time intervals by projector controlled by computer, taking the deformed image of each pattern by camera, then collecting, storing and processing the image data, and reconstructing the 3D contour based on the phase information of each point on the object and the geometric relation between projector, camera and object. Its advantages are the simplified projector and phase-shift unit, and high effect.
Description
The invention belongs to three-dimensional sensing and metering field, particularly the computer reconstruction technology of three-dimensional object surface pattern.
The computer reconstruction technology of three-dimensional object surface profile is one of focus of at present domestic and international three-dimensional sensing and quantitative study, and it is widely used in fields such as Aero-Space, machine-building, online detection and quality control, medical diagnosis, computer aided design/computer aided machine (CAD/CAM), robotic vision system.Phase-shifting method is the most ripe, the most reliable a kind of method in all three-dimensional appearance reconstructing method, it have calculated amount less, can be with reach very advantage such as high sensitivity and precision and anti-static noise than coarse grating, be the preferred option under the situation of conditions permit.
The principle of work of phase-shifting method is described below:
Shown in Fig. 1 a and Fig. 1 b is N step phase-shifting method.The sinusoidal grating striped that it will project to body surface moves N time, utilizes discrete phase shifting method to obtain the modulated image of the N width of cloth (N>3), goes out the position according to the N-PHASE algorithm computation and distributes mutually, utilizes geometric relationship to try to achieve the height profile of body surface at last.
As shown in Figure 1a, sinusoidal light field 13 is with incident angle θ
iIncide on the reference planes 11, ccd image sensor 12 is at θ
vReceive on the direction by the light field after reference planes 11 modulation, wherein z is the light field cycle that receives.The cycle of supposing the sinusoidal light field of projection is x, and the plane moved Δ h, and shown in Fig. 1 b, at this moment z just has a variation delta z: Δ z=Δ h[tan (θ
i)+tan (θ
v)] cos (θ
v) (1) this can cause a phase change Δ φ in the light field in modulation:
As can be seen from the above equation, the height change h in somewhere on the reference planes can be tried to achieve by the phase shift Δ φ of modulation light field.When ccd image sensor during perpendicular to reference planes, θ
vBe zero, so height change is:
As long as we get access to Δ φ, just can reconstruct the three-D profile of object.
When a light intensity of being regulated out by sinusoidal grating projected on the diffuse reflection body surface by two-dimensional projection's pattern of Sine distribution, the modulated image that is detected by imaging system can be described as: I
i(x, y)=a (x, y) [(l+m (and x, y) cos (φ (and x, y)+φ
i)] (4)
Wherein, (x y) represents bias light intensity to a, and (x y) represents fringe contrast to m, and (x y) has represented the feature of candy strip to phase function φ, thereby has comprised the information of profiling object surface, φ
iBe starting phase angle, N is the phase shift number of times.By regulating grating, can produce the sinusoidal projection pattern of different initial phases.Mobile grating N time (N 〉=3), the distance of each mobile P/N (P is the cycle of sinusoidal grating), then producing phase change is the N width of cloth projection pattern of 2 π/N.Note image after the modulation of each width of cloth pattern with ccd image sensor, have according to the N-PHASE algorithm:
So Δ φ can ask.
Traditional phase-shifting method measuring system is utilized sinusoidal grating to produce the projection light field and is shone on the testee, drive grating by stepper motor and move (N 〉=3) N time, the distance of each mobile P/N (P is the cycle of sinusoidal grating), produce the sinusoidal projection pattern of the different initial phases of the N width of cloth, note the deformation pattern that each width of cloth projection pattern is formed by the body surface modulation successively with video camera then; Computing machine is gathered video image by image pick-up card, and image data storage is got up; Handle these view data with N step phase-shifting method algorithm (N-PHASE algorithm), obtain the phase information of each point on the object; Reconstruct the three-D profile of object according to the geometric relationship between the phase information of each point on the object and projector, gamma camera, the testee.Because a variety of causes such as step-out of machining and stepper motor, the grating displacement of expectation usually can not accurately obtain, and this has brought very mistake to measurement result.Traditional phase changer as shown in Figure 2, the distance that mostly adopts the next real-time monitoring grating of accurate micro-displacement sensor to move, the distance that reality is moved feeds back to computing machine, and then does corresponding compensation according to deviation by servoamplifier drive stepping motor drive worktable.Though this method can be controlled phase shift accurately, the micro-displacement sensor of working stability costs an arm and a leg on the market, and volume is bigger, makes the design and installation of this cover feedback control system very inconvenient.
The objective of the invention is to for overcoming traditional phase-shifting method measuring system phase changer cost height, design, the complicated weak point of control, a kind of method of digital projection has been proposed, utilize computing machine to realize projection light field and phase shift from software, simplify projection, phase shift section in the traditional measurement system greatly, obtained result preferably.
The method for reconstructing 3 D contour of digital projection that the present invention proposes is characterized in that: work out the projection pattern generator program in advance, with generating the sinusoidal grating of fixed skew projection pattern is arranged more than three width of cloth; Move described program, the digital projector that computerizeds control projects projection pattern on the testee surface at interval successively by certain hour; Video camera is noted the deformation pattern that each width of cloth projection pattern is formed by the body surface modulation successively; Computing machine is gathered video image by image pick-up card, and image data storage is got up; With these view data of N step phase-shifting method algorithm process, obtain the phase information of each point on the object; Reconstruct the three-D profile of object according to the geometric relationship between the phase information of each point on the object and projector, gamma camera, the testee.
The generation of above-mentioned projection pattern can realize by various programming languages.Since on the projection pattern that generates from left to right light intensity change by the sinusoidal rule cycle, so can use computing machine programming Control projector, from left to right draw a series of straight lines by sinusoidal grey scale change successively and realize sinusoidal projection light field pattern.
The processing of steps such as the collection of projected image of the present invention, processing all can be adopted traditional disposal route.
The camera acquisition image is for different image pick-up cards, and its implementation also may be different.For example,, can programme so to be allowed to condition at and gather a few width of cloth images during every width of cloth pattern throws, select for use the best width of cloth of effect to enter down the Flame Image Process in step as the modulation image of this width of cloth projection pattern correspondence if used image pick-up card is programmable.If image pick-up card can not be programmed, can adopt the method for artificial control computer to carry out image acquisition.
The present invention is by adopting the method for digital projection, and it is higher effectively to have overcome traditional phase-shifting method measuring system phase changer cost, and design, the complicated deficiency of control have obtained good result.
Brief Description Of Drawings:
Fig. 1 is the phase-shifting method principle schematic.
Fig. 2 moves control chart for the grating of traditional phase changer.
Fig. 3 is an embodiment of the invention measurement mechanism synoptic diagram.
Fig. 4 is an embodiment of the invention measurement mechanism geometric relationship synoptic diagram.
Fig. 5 is the projection pattern automatically generating program process flow diagram of the embodiment of the invention.
Fig. 6 is the Flame Image Process process flow diagram of the embodiment of the invention.
A kind of computer control projector that the present invention proposes realizes that the embodiment of the method for reconstructing 3 D contour of projection light field and phase shift is described in detail as follows in conjunction with the accompanying drawings:
The present invention is embodiment with 4 width of cloth projection pattern, the measuring system of present embodiment as shown in Figure 3, projector 32 adopts EPSON7100 type LCD projector, resolution 1024 * 768 pixels, illumination 650 lumens; Imageing sensor 35 is quick logical TMC-56 type colored CCD gamma cameras, 640 * 480 pixels, 256 grades of gray scales (B/W output); Computing machine 31 performances are Pentium 233MHz, have programmable video input card, and 33,34,36 is respectively projection pattern, testee and the video image that collects.
The projection that above-mentioned measuring system adopts and the geometric relationship of imaging system as shown in Figure 4, P
2And I
2Emergent pupil and the entrance pupil of representing optical projection system and imaging system respectively, O are the intersection points of optical projection system and imaging system optical axis, and reference planes are vertical with the optical axis of imaging system.Testee should be placed in about on the point-symmetric position of O as far as possible, and the visual field of projection, imaging is evenly distributed.Under this situation, in the computing formula (3)
, remaining measures Δ φ with phase-shifting method exactly.
The specific implementation step of the method for reconstructing 3 D contour of digital projection of present embodiment comprises: at first, establishment projection pattern generator program, its step as shown in Figure 5, defining variable at first, initialize; Define drawing area then, finish initialization; Adopt Visual Basic language to recirculate initial phase according to each width of cloth image by the setting-out of sinusoidal gray scale rule again by three.The projection pattern of definition control earlier number, every width of cloth projection pattern projection number of cycles, interior setting-out number of each projection cycle and their control threshold value variablees separately are respectively n, i, j, c, b, a, compose zero for n, i, j.Define a time-delay variable, assignment is 5 seconds in the present embodiment.Define 0.7 times of full frame width, full frame altitude range is a drawing area, and drawing area is provided with between two parties, finishes the drawing initialization.Present embodiment is drawn 128 lines by each in grey scale change cycle in drawing area, 7 cycles of picture draw altogether, so only needs be taken as 127 to a among Fig. 6, and b is taken as 6, and c is taken as 3, realizes that with Visual Basic Programming with Pascal Language this three recirculates and get final product.
Secondly, the program that operation is finished above, the LCD projector that computerizeds control projects the projection pattern that generates on the body surface successively by 5 second time interval setting.Read the Programming Notes of programmable vision capture card, work out corresponding program, during this period of time at interval in, allow the CCD gamma camera note the deformation pattern that each width of cloth projection pattern is formed by the body surface modulation successively, and they be stored as the BMP file.
At last, carry out Flame Image Process, its process as shown in Figure 6.At first will read in the image file of acquisition, wherein gradation data is saved as matrix, present embodiment uses the imread () function under the Matlab directly to read in the BMP file, and gradation data wherein is saved as matrix; Next the data that obtain are carried out filtering, remove noise wherein, present embodiment has utilized under the Matlab powerful wavelet filtering and Fourier transform kit to finish this work; When reading in successively, after Filtering Processing crosses each width of cloth view data, obtaining the phase information of each point on the object, obtain the height of each point by the geometric relationship of system by formula (3), thereby reconstruct the three-D profile of object according to N-PHASE algorithmic formula (6).
Experiment shows that method of the present invention has been simplified projection, the phase shift section in the traditional measurement system greatly, has obtained gratifying result.
Claims (1)
1, a kind of method for reconstructing 3 D contour of digital projection based on phase-shifting method is characterized in that may further comprise the steps:
(1) works out the projection pattern generator program in advance, the sinusoidal grating of fixed skew projection pattern is arranged more than three width of cloth with generating;
(2) move described program, the digital projector that computerizeds control projects projection pattern on the testee surface at interval successively by certain hour;
(3) video camera is noted the deformation pattern that each width of cloth projection pattern is formed by the body surface modulation successively;
(4) computing machine is gathered video image by image pick-up card, and image data storage is got up;
(5), obtain the phase information of each point on the object with these view data of N step phase-shifting method algorithm process;
(6) reconstruct the three-D profile of object according to the geometric relationship between the phase information of each point on the object and projector, gamma camera, the testee.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB011041633A CN1149388C (en) | 2001-02-23 | 2001-02-23 | Method for reconstructing 3D contour of digital projection based on phase-shifting method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB011041633A CN1149388C (en) | 2001-02-23 | 2001-02-23 | Method for reconstructing 3D contour of digital projection based on phase-shifting method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1302999A true CN1302999A (en) | 2001-07-11 |
CN1149388C CN1149388C (en) | 2004-05-12 |
Family
ID=4653721
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB011041633A Expired - Fee Related CN1149388C (en) | 2001-02-23 | 2001-02-23 | Method for reconstructing 3D contour of digital projection based on phase-shifting method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN1149388C (en) |
Cited By (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1302257C (en) * | 2003-10-09 | 2007-02-28 | 四川大学 | Method for realizing Fourier transfrom contour by generating pi phase shift by two-colour grating template |
CN1312461C (en) * | 2004-06-23 | 2007-04-25 | 香港理工大学 | Reconstruction system and method for sheet three-dimensional surface of flexible body |
CN100408972C (en) * | 2006-07-27 | 2008-08-06 | 西安交通大学 | Three dimension object contour phase measuring method based on double frequency color strip projection |
CN100416466C (en) * | 2006-05-31 | 2008-09-03 | 北京航空航天大学 | Single-eye vision semi-matter simulating system and method |
CN100417231C (en) * | 2006-05-31 | 2008-09-03 | 北京航空航天大学 | Three-dimensional vision semi-matter simulating system and method |
CN101299227B (en) * | 2008-06-27 | 2010-06-09 | 北京中星微电子有限公司 | Multi-person game system and method based on three-dimensional reconstruction |
CN101238488B (en) * | 2005-08-04 | 2010-08-04 | 皇家飞利浦电子股份有限公司 | Method and system for producing or reconstructing 3-D image |
CN101256672B (en) * | 2008-03-21 | 2011-10-12 | 北京中星微电子有限公司 | Object image depth restruction apparatus based on video camera apparatus as well as projecting apparatus thereof |
CN102261896A (en) * | 2011-04-19 | 2011-11-30 | 长春东瑞科技发展有限公司 | Method and system for measuring three-dimensional shape of object based on phase measurement |
CN102589481A (en) * | 2012-03-13 | 2012-07-18 | 周波 | Design of sinusoidal grating under modulation of optimal width |
CN102589477A (en) * | 2012-02-14 | 2012-07-18 | 乔付 | Design of sinusoidal grating in optic three-dimensional outline measurement |
CN102589480A (en) * | 2012-03-08 | 2012-07-18 | 乔付 | Design of sinusoidal grating in structured light three-dimensional vision measurement |
US8331688B2 (en) | 2009-01-07 | 2012-12-11 | International Business Machines Corporation | Focus-based edge detection |
CN104567721A (en) * | 2015-01-23 | 2015-04-29 | 清华大学 | Continuous shear interference measurement method |
CN105547540A (en) * | 2015-12-09 | 2016-05-04 | 清华大学 | Coherent gradient sensitive interference method for real-time space phase shift |
CN107810384A (en) * | 2015-06-17 | 2018-03-16 | 惠普德国股份有限公司 | Fringe projection method, fringe projector apparatus and computer program product |
CN108784836A (en) * | 2018-06-20 | 2018-11-13 | 安徽医科大学第附属医院 | Based on image processing system in the calm management of optimization and regional block orthopaedics anesthesia art |
-
2001
- 2001-02-23 CN CNB011041633A patent/CN1149388C/en not_active Expired - Fee Related
Cited By (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1302257C (en) * | 2003-10-09 | 2007-02-28 | 四川大学 | Method for realizing Fourier transfrom contour by generating pi phase shift by two-colour grating template |
CN1312461C (en) * | 2004-06-23 | 2007-04-25 | 香港理工大学 | Reconstruction system and method for sheet three-dimensional surface of flexible body |
CN101238488B (en) * | 2005-08-04 | 2010-08-04 | 皇家飞利浦电子股份有限公司 | Method and system for producing or reconstructing 3-D image |
CN100416466C (en) * | 2006-05-31 | 2008-09-03 | 北京航空航天大学 | Single-eye vision semi-matter simulating system and method |
CN100417231C (en) * | 2006-05-31 | 2008-09-03 | 北京航空航天大学 | Three-dimensional vision semi-matter simulating system and method |
CN100408972C (en) * | 2006-07-27 | 2008-08-06 | 西安交通大学 | Three dimension object contour phase measuring method based on double frequency color strip projection |
CN101256672B (en) * | 2008-03-21 | 2011-10-12 | 北京中星微电子有限公司 | Object image depth restruction apparatus based on video camera apparatus as well as projecting apparatus thereof |
CN101299227B (en) * | 2008-06-27 | 2010-06-09 | 北京中星微电子有限公司 | Multi-person game system and method based on three-dimensional reconstruction |
US8331688B2 (en) | 2009-01-07 | 2012-12-11 | International Business Machines Corporation | Focus-based edge detection |
US8509562B2 (en) | 2009-01-07 | 2013-08-13 | International Business Machines Corporation | Focus-based edge detection |
CN102261896A (en) * | 2011-04-19 | 2011-11-30 | 长春东瑞科技发展有限公司 | Method and system for measuring three-dimensional shape of object based on phase measurement |
CN102589477A (en) * | 2012-02-14 | 2012-07-18 | 乔付 | Design of sinusoidal grating in optic three-dimensional outline measurement |
CN102589480A (en) * | 2012-03-08 | 2012-07-18 | 乔付 | Design of sinusoidal grating in structured light three-dimensional vision measurement |
CN102589480B (en) * | 2012-03-08 | 2013-12-11 | 黑龙江科技学院 | Design of sinusoidal grating in structured light three-dimensional vision measurement |
CN102589481A (en) * | 2012-03-13 | 2012-07-18 | 周波 | Design of sinusoidal grating under modulation of optimal width |
CN102589481B (en) * | 2012-03-13 | 2013-12-11 | 黑龙江科技学院 | Design of sinusoidal grating under modulation of optimal width |
CN104567721A (en) * | 2015-01-23 | 2015-04-29 | 清华大学 | Continuous shear interference measurement method |
CN107810384A (en) * | 2015-06-17 | 2018-03-16 | 惠普德国股份有限公司 | Fringe projection method, fringe projector apparatus and computer program product |
US10801834B2 (en) | 2015-06-17 | 2020-10-13 | Hewlett-Packard Development Company, L.P. | Fringe projection for determining topography of a body |
CN107810384B (en) * | 2015-06-17 | 2020-11-20 | 惠普德国股份有限公司 | Stripe projection method, stripe projection apparatus, and computer program product |
CN105547540A (en) * | 2015-12-09 | 2016-05-04 | 清华大学 | Coherent gradient sensitive interference method for real-time space phase shift |
CN108784836A (en) * | 2018-06-20 | 2018-11-13 | 安徽医科大学第附属医院 | Based on image processing system in the calm management of optimization and regional block orthopaedics anesthesia art |
Also Published As
Publication number | Publication date |
---|---|
CN1149388C (en) | 2004-05-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN1149388C (en) | Method for reconstructing 3D contour of digital projection based on phase-shifting method | |
CN102188290B (en) | Device and method for acquiring 3-D surface profile image data of tooth | |
US7286246B2 (en) | Method and apparatus for non-contact three-dimensional surface measurement | |
EP1994503B1 (en) | Depth-varying light fields for three dimensional sensing | |
Lindner et al. | Lateral and depth calibration of PMD-distance sensors | |
US8050461B2 (en) | Depth-varying light fields for three dimensional sensing | |
JP3525964B2 (en) | 3D shape measurement method for objects | |
US10706570B2 (en) | System and method to acquire the three-dimensional shape of an object using a moving patterned substrate | |
JP7319913B2 (en) | METHOD FOR MEASURING Z-HEIGHT VALUE OF WORKPIECE SURFACE USING MACHINE VISION INSPECTION SYSTEM | |
EP3064895A1 (en) | Linear illumination device | |
JPH0338524B2 (en) | ||
CN111536904A (en) | Three-dimensional measurement method and system based on structural illumination and storage medium | |
Jia et al. | Two-step triangular-pattern phase-shifting method for three-dimensional object-shape measurement | |
Forest Collado | New methods for triangulation-based shape acquisition using laser scanners | |
EP3839599A1 (en) | Metrology system with transparent workpiece surface mode | |
Traxler et al. | Experimental comparison of optical inline 3D measurement and inspection systems | |
WO2004097335A1 (en) | Shape and deformation measurements of large objects by fringe projection | |
CN111982026A (en) | Grating projection three-dimensional measuring device and measuring method for high-reflection object curved surface | |
JP2006170744A (en) | Three-dimensional distance measuring instrument | |
JPH07260451A (en) | Three dimensional shape measuring system | |
US11928851B2 (en) | Computer-implemented method, computer program, data processing system and device for determining the reflectance behavior of a surface of an object, and storage medium with instructions stored thereon for determining the reflectance behavior of a surface of an object | |
Uyanik et al. | A method for determining 3D surface points of objects by a single camera and rotary stage | |
Kammel | Topography reconstruction of specular surfaces from a series of gray-scale images | |
Gray et al. | Simulation of a tire inspection system | |
CN215114433U (en) | Structured light three-dimensional shape measuring device based on optical phase element |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C06 | Publication | ||
PB01 | Publication | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
C19 | Lapse of patent right due to non-payment of the annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |