CN102798354A - Binary stripe stack based sinusoidal grating generation method - Google Patents

Binary stripe stack based sinusoidal grating generation method Download PDF

Info

Publication number
CN102798354A
CN102798354A CN2012103107143A CN201210310714A CN102798354A CN 102798354 A CN102798354 A CN 102798354A CN 2012103107143 A CN2012103107143 A CN 2012103107143A CN 201210310714 A CN201210310714 A CN 201210310714A CN 102798354 A CN102798354 A CN 102798354A
Authority
CN
China
Prior art keywords
scale
picture
pictures
sinusoidal grating
binary
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
CN2012103107143A
Other languages
Chinese (zh)
Other versions
CN102798354B (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.)
Northwestern Polytechnical University
Original Assignee
Northwestern Polytechnical University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Northwestern Polytechnical University filed Critical Northwestern Polytechnical University
Priority to CN201210310714.3A priority Critical patent/CN102798354B/en
Publication of CN102798354A publication Critical patent/CN102798354A/en
Application granted granted Critical
Publication of CN102798354B publication Critical patent/CN102798354B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Length Measuring Devices By Optical Means (AREA)

Abstract

The invention discloses a binary stripe stack based sinusoidal grating generation method. The method provided by the invention comprises the following steps of: firstly generating a sinusoidal grating picture, and decomposing one sinusoidal grating picture into 8 binary pictures through a binary method, wherein each point of each binary picture represents a bit of each point of the original sinusoidal grating after binary, and the value only can be 1 or 1; and through projecting the 8 binary pictures, finally combining the 8 binary pictures into one sinusoidal grating picture. Through the method provided by the invention, a composite plane sinusoidal grating protecting picture E can be obtained, the range of the pixel value of each point of the projecting picture E is from 0 to 255, and is in sinusoidal distribution, and thus errors caused by the nonlinear response of a projector and a camera can be prevented.

Description

Sinusoidal grating generation method based on the stack of scale-of-two striped
Technical field
The invention belongs to projector and camera nonlinear response association area.Relate to the synthetic problem of sine streak grating.Be specially a kind of sinusoidal grating generation method based on the stack of scale-of-two striped.
Background technology
Optical three-dimensional measurement has widespread use in various fields such as industry detection automatically, production quality control, reverse engineer, biomedicine, virtual reality, the reproduction of the cultural relics, anthropological measuring.This huge application demand has been impelled the fast development of multiple optical measuring technique.Along with computer technology, digital picture are obtained the development of equipment and optical device, and a lot of three dimensional optical measuring technology have got into the stage of ripeness of commercial application, are also continuing to bring out with stylish method for three-dimensional measurement.Have noncontact, high precision, characteristics such as high-speed based on the optical 3-dimensional surface shape measurement method of fringe projection, be widely used in fields such as machine vision, the detection of face shape, profiling in kind, production quality control, biomedicine.
The measuring error of structured light measurement system has two main sources: phase displacement error and nonsinusoidal waveform error.Phase displacement error is because due to phase shift step pitch inaccurate, it usually is inevitably, but can the real-time phase shift correction technology of employing reduces in accurate phase changer and the measuring process through adopting.Along with the development of digital display technique, business-like digital projector (DLP) is used widely in structured light measurement system.Digital projector projection grating can be eliminated phase displacement error; But because design reasons; Projector and camera all have the character of gray scale nonlinear distortion, make that taking the sinusoidal grating projected fringe that obtains does not have good sine property, can bring phase error to measurement thereupon.The nonlinear distortion that projector has is called gradation distortion, and the input and output that are projector present the exponential function response relation.Gradation distortion will become the principal element that influences the measuring system measuring accuracy.
The existing method that reduces the non-linear measuring error that causes of gamma comprises: it is non-linear to utilize mathematical method directly to handle gamma; Increase the phase shift step number and reduce error; Adopt look-up table to preestablish error; Based on the spline curve fitting method.These methods can improve measuring accuracy greatly, but also have some defectives: the spline-fitting method need be through iteration repeatedly when phase calculation, and operand is bigger; Increase the increase that the phase shift step number means projection image quantity; Look-up table needs the plenty of time in the process of labeling projection appearance gamma value.The present invention proposes a kind ofly to this problem, solve nonlinear problem based on binary pattern.
Summary of the invention
The technical matters that solves
For solving the problem that prior art exists, the present invention proposes a kind of sinusoidal grating generation method based on the stack of scale-of-two striped.Can overcome projector and camera nonlinear response and theoretical scale-of-two picture stack generated sinusoidal grating through binary law.
Technical scheme
The present invention at first generates a sinusoidal grating picture, again through binary method, 1 sinusoidal grating picture is decomposed into 8 scale-of-two pictures.According to binary decomposition principle, each point of the former sinusoidal grating of each some representative of every binary picture carries out a position behind the scale-of-two, and value can only be 0 or 1.Through 8 scale-of-two pictures of projection, 8 scale-of-two pictures synthesize sinusoidal grating pictures the most at last.
Technical scheme of the present invention is:
Said a kind of sinusoidal grating generation method based on the stack of scale-of-two striped is characterized in that: may further comprise the steps:
Step 1: the sinusoidal grating striped picture A that generates at least two cycles;
Step 2: according to the scale-of-two reconstruction formula:
I ( x , y ) = Σ i = 1 8 2 ( 8 - i ) M i ( x , y )
Sinusoidal grating striped picture A reconstruct is decomposed into 8 scale-of-two picture { M i(i=1,2 ..., 8), wherein (x, the y) coordinate figure of pixel among the expression picture A, I (x, the y) gray-scale value of expression sinusoidal grating picture A each point, M iThe scale-of-two picture that expression reconstruct obtains, 8 scale-of-two pictures are formed set of pictures B; At scale-of-two picture M iIn, work as M i(x, y)=0 o'clock, (x, gray-scale value y) is 0 to this point, works as M i(x, y)=1 o'clock, (x, gray-scale value y) they are 255 to this point;
Step 3: every scale-of-two picture among the set of pictures B is projected on the testee respectively, and take the projection image on the testee respectively through the CCD camera, obtain 8 projection images, 8 projection images are formed set of pictures C; Utilize big Tianjin threshold method that every projection image among the set of pictures C is done binary conversion treatment, 8 scale-of-two pictures that obtain after the processing are formed set of pictures D;
Step 4: 8 scale-of-two pictures among the set of pictures D are carried out reverse overlap-add procedure according to the scale-of-two reconstruction formula in the step 2, thereby obtain the sinusoidal grating image E of generation.
Said a kind of sinusoidal grating generation method based on the stack of scale-of-two striped, it is characterized in that: the angle of the shooting visual angle of CCD camera and projector optical axis is 20 ° ~ 40 ° in the step 3.
Beneficial effect
Can obtain the plane sinusoidal grating projection image E that opening and closing become through the present invention, the scope of the pixel value of projection image E each point is 0 to 255, and becomes Sine distribution, thereby has avoided the nonlinear response of projector and camera and the error that causes.
Description of drawings
Fig. 1 is sinusoidal grating figure.
Fig. 2 (a) obtains the sine streak raster pattern for through the projection camera plane.
Fig. 2 (b) is an image middle line intensity profile situation map.
The picture of Fig. 2 (c) for obtaining after the Fast Fourier Transform (FFT).
Fig. 3 is 8 binary pictures and corresponding sinusoidal grating figure.
Fig. 4 (a) is after taking 8 binary pictures, synthetic sine streak raster pattern.
Fig. 4 (b) is an image middle line intensity profile situation map.
The figure of Fig. 4 (c) for obtaining after the Fast Fourier Transform (FFT).
Embodiment
Below in conjunction with specific embodiment the present invention is described:
Based on the sinusoidal grating generation method of scale-of-two striped stack, may further comprise the steps in the present embodiment:
Step 1: the sinusoidal grating striped picture A that generates at least two cycles; Shown in accompanying drawing 1, the generation one-period is 3 sinusoidal grating striped picture A in the present embodiment.
In order to compare; (model C P270 BenQ) projects to standard sine grating picture on the testee, takes the projection image on the testee through CCD camera (model MVD-500SM) to use resolution to be set to the projector of 1024*768; Obtaining resolution is the 1024*768 projection image; Shown in accompanying drawing 2 (a), wherein the angle of the shooting visual angle of camera and projector optical axis is 20 ° ~ 40 °, gets 30 ° in the present embodiment.In order to contrast, get the intensity profile of middle row, shown in accompanying drawing 2 (b), and picture carried out Fast Fourier Transform (FFT), the result is shown in accompanying drawing 2 (c).
Step 2: according to the scale-of-two reconstruction formula:
I ( x , y ) = Σ i = 1 8 2 ( 8 - i ) M i ( x , y )
Sinusoidal grating striped picture A reconstruct is decomposed into 8 scale-of-two picture { M i(i=1,2 ..., 8), wherein (x, the y) coordinate figure of pixel among the expression picture A, I (x, the y) gray-scale value of expression sinusoidal grating picture A each point, M iThe scale-of-two picture that expression reconstruct obtains, M iIn i represent binary position, M i(x, value y) is 0 or 1; 8 scale-of-two pictures are formed set of pictures B; At scale-of-two picture M iIn, work as M i(x, y)=0 o'clock, (x, gray-scale value y) is 0 to this point, works as M i(x, y)=1 o'clock, (x, gray-scale value y) they are 255 to this point; Sinusoidal grating striped picture A resolves into 8 binary images shown in accompanying drawing (3);
Step 3: use resolution to be set to projector (the model C P270 of 1024*768; BenQ) project to every scale-of-two picture among the set of pictures B on the testee respectively; And take the projection image on the testee respectively through CCD camera (model MVD-500SM), and obtaining resolution is 8 projection images of 1024*768,8 projection images are formed set of pictures C; Wherein the angle of the shooting visual angle of camera and projector optical axis is 20 ° ~ 40 °, gets 30 ° in the present embodiment.Utilize big Tianjin threshold method that every projection image among the set of pictures C is done binary conversion treatment, 8 scale-of-two pictures that obtain after the processing are formed set of pictures D;
Step 4: 8 scale-of-two pictures among the set of pictures D are carried out reverse overlap-add procedure according to the scale-of-two reconstruction formula in the step 2, thereby obtain the sinusoidal grating image E of generation, shown in accompanying drawing 4 (a).In order to compare, get the intensity profile of middle row, shown in accompanying drawing 4 (b), and picture carried out Fast Fourier Transform (FFT), the result is shown in accompanying drawing 4 (c).
The image of the sinusoidal grating that the sinusoidal grating image through contrast standard and 8 scale-of-two picture stacks generate can be found clearly that the pixel value that adopts the scale-of-two picture to handle each point in the sinusoidal grating image that generates is 0 to 255, and become Sine distribution.And directly the pixel minimum of the sinusoidal grating image each point of projection is not 0, and does not become Sine distribution.Institute of the present invention elaboration method has overcome the non-linear corresponding problem of projector and camera.

Claims (2)

1. sinusoidal grating generation method based on scale-of-two striped stack is characterized in that: may further comprise the steps:
Step 1: the sinusoidal grating striped picture A that generates at least two cycles;
Step 2: according to the scale-of-two reconstruction formula:
I ( x , y ) = Σ i = 1 8 2 ( 8 - i ) M i ( x , y )
Sinusoidal grating striped picture A reconstruct is decomposed into 8 scale-of-two picture { M i(i=1,2 ..., 8), wherein (x, the y) coordinate figure of pixel among the expression picture A, I (x, the y) gray-scale value of expression sinusoidal grating picture A each point, M iThe scale-of-two picture that expression reconstruct obtains, 8 scale-of-two pictures are formed set of pictures B; At scale-of-two picture M iIn, work as M i(x, y)=0 o'clock, (x, gray-scale value y) is 0 to this point, works as M i(x, y)=1 o'clock, (x, gray-scale value y) they are 255 to this point;
Step 3: every scale-of-two picture among the set of pictures B is projected on the testee respectively, and take the projection image on the testee respectively through the CCD camera, obtain 8 projection images, 8 projection images are formed set of pictures C; Utilize big Tianjin threshold method that every projection image among the set of pictures C is done binary conversion treatment, 8 scale-of-two pictures that obtain after the processing are formed set of pictures D;
Step 4: 8 scale-of-two pictures among the set of pictures D are carried out reverse overlap-add procedure according to the scale-of-two reconstruction formula in the step 2, thereby obtain the sinusoidal grating image E of generation.
2. according to the said a kind of sinusoidal grating generation method based on the stack of scale-of-two striped of claim 1, it is characterized in that: the angle of the shooting visual angle of CCD camera and projector optical axis is 20 ° ~ 40 ° in the step 3.
CN201210310714.3A 2012-08-28 2012-08-28 Binary stripe stack based sinusoidal grating generation method Expired - Fee Related CN102798354B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210310714.3A CN102798354B (en) 2012-08-28 2012-08-28 Binary stripe stack based sinusoidal grating generation method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210310714.3A CN102798354B (en) 2012-08-28 2012-08-28 Binary stripe stack based sinusoidal grating generation method

Publications (2)

Publication Number Publication Date
CN102798354A true CN102798354A (en) 2012-11-28
CN102798354B CN102798354B (en) 2014-11-05

Family

ID=47197558

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210310714.3A Expired - Fee Related CN102798354B (en) 2012-08-28 2012-08-28 Binary stripe stack based sinusoidal grating generation method

Country Status (1)

Country Link
CN (1) CN102798354B (en)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104315996A (en) * 2014-10-20 2015-01-28 四川大学 Method for realizing fourier transform profilometry by using binary encoding strategy
CN104567730A (en) * 2015-01-15 2015-04-29 四川大学 Method for generating light field with sine structure by means of space-time binary encoding
CN108180836A (en) * 2014-06-27 2018-06-19 佳能株式会社 position detecting device
CN109029294A (en) * 2018-08-21 2018-12-18 合肥工业大学 Based on the Fast gray striped synthetic method for focusing two-value pattern
CN110375673A (en) * 2019-07-01 2019-10-25 武汉斌果科技有限公司 A kind of big depth of field two-value defocus method for three-dimensional measurement based on multifocal optical projection system
CN114941999A (en) * 2022-07-22 2022-08-26 南京信息工程大学 Binary coding stripe design method for structured light projection

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1928491A (en) * 2005-09-07 2007-03-14 北京航空航天大学 Grating projection and recognition method
CN101089551A (en) * 2007-07-11 2007-12-19 华中科技大学 3D measuring method based on two-dimensional code
KR20090116560A (en) * 2008-05-07 2009-11-11 선문대학교 산학협력단 Grating moving apparatus and method for measuring three-dimensional using moire
CN101813462A (en) * 2010-04-16 2010-08-25 天津理工大学 Three-dimensional feature optical measuring system controlled by uniprocessor and measuring method

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1928491A (en) * 2005-09-07 2007-03-14 北京航空航天大学 Grating projection and recognition method
CN101089551A (en) * 2007-07-11 2007-12-19 华中科技大学 3D measuring method based on two-dimensional code
KR20090116560A (en) * 2008-05-07 2009-11-11 선문대학교 산학협력단 Grating moving apparatus and method for measuring three-dimensional using moire
CN101813462A (en) * 2010-04-16 2010-08-25 天津理工大学 Three-dimensional feature optical measuring system controlled by uniprocessor and measuring method

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
马少鹏: "应用时域相位解包方法的三维形貌测量系统", 《光学技术》, vol. 28, no. 5, 30 September 2002 (2002-09-30), pages 395 - 400 *

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108180836A (en) * 2014-06-27 2018-06-19 佳能株式会社 position detecting device
CN104315996A (en) * 2014-10-20 2015-01-28 四川大学 Method for realizing fourier transform profilometry by using binary encoding strategy
CN104315996B (en) * 2014-10-20 2018-04-13 四川大学 The method that Fourier transform profilometry is realized with binary coding strategy
CN104567730A (en) * 2015-01-15 2015-04-29 四川大学 Method for generating light field with sine structure by means of space-time binary encoding
CN109029294A (en) * 2018-08-21 2018-12-18 合肥工业大学 Based on the Fast gray striped synthetic method for focusing two-value pattern
CN110375673A (en) * 2019-07-01 2019-10-25 武汉斌果科技有限公司 A kind of big depth of field two-value defocus method for three-dimensional measurement based on multifocal optical projection system
CN110375673B (en) * 2019-07-01 2021-04-16 武汉斌果科技有限公司 Large-depth-of-field binary out-of-focus three-dimensional measurement method based on multi-focus projection system
CN114941999A (en) * 2022-07-22 2022-08-26 南京信息工程大学 Binary coding stripe design method for structured light projection

Also Published As

Publication number Publication date
CN102798354B (en) 2014-11-05

Similar Documents

Publication Publication Date Title
CN107607060B (en) A kind of phase error compensation method applied in the measurement of grating tripleplane
CN102798354B (en) Binary stripe stack based sinusoidal grating generation method
CN101881605B (en) Optical three-dimensional measuring method based on phase coding technology
CN104897083B (en) Three-dimensional rapid measurement method for raster projection based on defocusing phase-unwrapping of projector
Falcao et al. Plane-based calibration of a projector-camera system
CN202074952U (en) Three dimensional profiles and colored textures acquisition system based on single camera-single overhead projector
CN105806259A (en) Three-dimensional measuring method based on binary grating out-of-focus projection
CN104567730B (en) A kind of method that space-time binary coding produces sinusoidal light field
CN104596439A (en) Speckle matching and three-dimensional measuring method based on phase information aiding
CN103400366A (en) Method for acquiring dynamic scene depth based on fringe structure light
CN104677308A (en) Three-dimensional scanning method for high-frequency two-value strip
CN109141291A (en) A kind of fast phase unwrapping algorithm
CN101236066B (en) Projection grating self-correction method
CN105528770A (en) Projector lens distortion correcting method
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
CN103791856A (en) Phase solving and unwrapping method based on four-width optical grating stripe image
CN105606038A (en) Gamma non-linear correction method of phase measurement profilometry and system thereof
CN110207620A (en) A kind of three-dimensional rebuilding method determining digital fringe projection structure optical orders by different frequency
CN109945802A (en) A kind of structural light three-dimensional measurement method
CN106500626A (en) A kind of mobile phone stereoscopic imaging method and three-dimensional imaging mobile phone
CN106840038B (en) A kind of double frequency phase shift and Gray's code character close method
CN105043301A (en) Grating strip phase solving method used for three-dimensional measurement
CN106091988A (en) A kind of compensation method of high-temperature forging three-dimensional dimension certainty of measurement
CN103942802A (en) Method for obtaining depth of structured light dynamic scene on basis of random templates

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: 20141105

Termination date: 20170828

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