Source: http://www.docstoc.com/docs/55655654/Method-And-System-For-Image-Processing-For-Structured-Light-Profiling-Of-A-Part---Patent-7302109
Timestamp: 2013-05-19 17:24:25
Document Index: 420902808

Matched Legal Cases: ['art. 14', 'art. 19', 'art 16', 'art 16', 'art 16', 'art 16']

Method And System For Image Processing For Structured Light Profiling Of A Part - Patent 7302109 Abstract
Inventors:Qingying Hu; Kevin George Harding; Joseph Benjamin Ross; Peter William Lorraine
Primary Examiner:Kassa; Yosef
Attorney, Agent of firm:Clarke; Penny A.Patnode; Patrick K.
Filed Date: 8/28/2003
US Patent Application: 10/652,366
Patent number: 7302109
Primary Class:Image Analysis (382/260)
Other US Classes:Facsimile And Static Presentation Processing (358/520)Facsimile And Static Presentation Processing (358/538)Image Analysis (382/274)Image Analysis (382/282)Image Analysis (382/284)
International Classes:G06K 9/40(20060101)
Categories:Legal > Patents > Image Analysis Tags:Method and system for image processin..., Hu, et al., Qingying Hu, Kevin George Harding, Joseph Benjamin Ross, Peter William Lorraine, Application number 10 652-366, Image Analysis, Facsimile And Static Presentation Pro..., image processing, image data, patent search, distribution function, present invention relates, storage medium, digital image, Patent Community, Page Range, English Español
4875777 Off-axis high accuracy structured light profilerHarding10/1/1989
5430810 Real time implementation of the hough transformSaeki7/1/1995
5528339Color image reproduction of scenes with color enhancement and preferential tone mappingBuhr et al.6/1/1996
5606390 Visual-line detecting device and optical apparatus having the sameArai et al.2/1/1997
5644390 Intensity distribution simulating methodYasuzato7/1/1997
5999840 System and method of registration of three-dimensional data setsGrimson et al.12/1/1999
6005984 Process and apparatus for extracting and recognizing figure elements using division into receptive fields, polar transformation, application of one-dimensional filter, and correlation between plurality of imagesKawakami et al.12/1/1999
6483641 Apparatus and methods relating to spatially light modulated microscopyMacAulay11/1/2002
6633683 Apparatus and method for adaptively reducing noise in a noisy input image signalDinh et al.10/1/2003
6639597 Visibility splatting and image reconstruction for surface elementsZwicker et al.10/1/2003
6782137 Digital image display improvement system and methodAvinash8/1/2004
7916898Method and system for identifying an edge of a cropAnderson3/29/2011
7969583System and method to determine an object distance from a reference point to a point on the object surfaceBendall, et al.6/28/2011
Claims What is claimed is: 1. An image processing method for structured light profiling, said method comprising: sampling an image of a structured light pattern to obtain an intensity distribution; selecting a plurality of sets of sampled points from the intensity distribution, wherein each of the respective sets comprises a plurality of sampled points; fitting each of the sets of sampled points to a respective Gaussian distribution function; andfiltering the Gaussian distribution functions to select a representative distribution function for the intensity distribution. 2. The image processing method of claim 1, further comprising extracting a center of the representative distribution function. 3. The image processing method of claim 2, wherein said extracting the center comprises using zero-crossing. 4. The image processing method of claim 1, wherein each of the sets of sampled points comprises at least three sampled points. 5. The image processing method of claim 1, wherein each of the distribution functions comprises a Gaussian distribution. 6. The image processing method of claim 1, further comprising extracting a respective center for each of the distribution functions, wherein said filtering includes using the centers. 7. The image processing method of claim 6, wherein said extracting comprises using zero-crossing. 8. The image processing method of claim 6, wherein said filtering includes rejecting a plurality of outliers. 9. The image processing method of claim 6, wherein said filtering includes selecting a median of the centers, wherein the representative distribution function corresponds to the median. 10. The image processing method of claim 9, wherein said filtering further includes using a histogram to determine the median. 11. The image processing method of claim 9, wherein said filtering further includes sorting the centers to determine the median. 12. The image processing method of claim 1, wherein said structured light pattern comprises at least one laser stripe, and wherein said intensity distribution is an intensity profile across one of the laser stripes. 13. The image processing method of claim 1, further comprising acquiring the image of the structured light pattern of a part. 14. The image processing method of claim 13, further comprising reconstructing a three dimensional profile of the part, using the representative distribution function. 15. An image processing method for structured light profiling, said method comprising: sampling an image of a structured light pattern to obtain an intensity distribution; selecting a plurality of sets of sampled points from the intensitydistribution, wherein each of the sets comprises at least three sampled points; fitting each of the sets of sampled points to respective Gaussian distribution functions; extracting a center for each of the Gaussian distribution functions; filteringthe Gaussian distribution functions by using the centers to select a representative distribution function for the intensity distribution. 16. The image processing method of claim 15, wherein said extracting comprises using zero-crossing. 17. The image processing method of claim 15, wherein said filtering includes selecting a median of the centers, wherein the representative distribution function corresponds to the median. 18. An image processing method for reconstructing a three dimensional profile of a part, said method comprising: projecting at least one laser beam on the part; acquiring an image of a structured light pattern of the part; sampling the imageto obtain an intensity distribution; selecting a plurality of sets of sampled points from the intensity distribution, wherein each of the sets comprises at least three sampled points; fitting each of the sets of sampled points to a respective Gaussiandistribution function; extracting a center for each of the Gaussian distribution functions; filtering the Gaussian distribution functions by using the centers to select a representative distribution function for the intensity distribution; and usingthe representative distribution function for the intensity distribution to generate the three dimensional profile of the part. 19. The image processing method of claim 18, wherein said filtering includes selecting a median of the centers, wherein the representative distribution function corresponds to the median. 20. A system for obtaining a three dimensional profile of a part using a structured light pattern, said system comprising: a source of structured light positioned at a predetermined distance from a part, wherein the source projects a beam ofstructured light to illuminate the part; at least one imaging device configured to acquire an image of a structured light pattern of the part, wherein the at least one imaging device is positioned such that an angle of view of the imaging device isdifferent from an angle of illumination of the source; and a processor coupled to the at least one imaging device, wherein the processor is configured for sampling the image of a structured light pattern to obtain an intensity distribution; selecting aplurality of sets of sampled points from the intensity distribution, and wherein each of the sets comprises at least three sampled points; fitting each of the sets of sampled points to respective Gaussian distribution functions; extracting a center foreach of the Gaussian distribution functions; filtering the Gaussian distribution functions by using the centers to select a representative distribution function for the intensity distribution; and reconstructing a three dimensional profile of the part,using the representative distribution function. 21. The system of claim 20, wherein said processor is further configured for selecting a median of the respective centers, wherein the representative distribution function corresponds to the median. 22. A computer readable medium for storing and/or transmitting instructions that, when executed by a computer, perform a method for image processing for structured light profiling, said method comprising: sampling an image of a structured lightpattern to obtain an intensity distribution; selecting a plurality of sets of sampled points from the intensity distribution, wherein each of the sets comprises a plurality of sampled points; fitting each of the sets of sampled points to a respectivedistribution function; and filtering the distribution functions to select a representative distribution function for the intensity distribution. 23. The computer readable medium of claim 22, wherein each of the sets of sampled points comprises at least three sampled points. 24. The computer readable medium of claim 22, wherein each of the distribution functions comprises a Gaussian distribution. 25. The computer readable medium of claim 22, wherein the method further includes extracting a respective center for each of the distribution functions, wherein said filtering includes using the centers. 26. The computer readable medium of claim 25, wherein the extracting comprises using zero-crossing. 27. The computer readable medium of claim 25, wherein the filtering includes selecting a median of the centers, wherein the representative distribution function corresponds to the median. Description BACKGROUND OF THE INVENTIONThe present invention relates generally to structured light profiling of a part and more specifically to methods and systems for image processing to obtain a three dimensional profile of a part using a structured light source.In structured light applications, a three dimensional profile of a part is obtained by using a laser source and at least one video camera. A laser source emanates a laser beam and spreads into multiple (one or more) laser stripes which strikes asurface of a part. These stripes are viewed from one or more video cameras from an angle other than the illumination angle. Typically, applications for structured light profiling include obtaining shape profile information about the part or generatinga 3D contour map of the part.The limitations of using laser stripes to obtain accurate profile information are mainly attributed to sampling error and the noise associated with the laser because the center of a laser stripe may not be imaged at the center of the pixel of thecamera and may not be the detected intensity peak. Sampling error occurs while locating the center of the laser stripe on the image. There are image processing techniques such as maximum intensity, intensity center, Gaussian fitting and zero-crossingwhich attempt to extract the relevant information from the laser stripe. The associated problem with several of these techniques is that it gives the location of the highest peak, which is not the true center of the stripe. Current correctiontechniques for this error include neighborhood averaging over neighboring pixels and doing a weighted average or using fitting methods, but these techniques also fail to adequately address the sampling error.The noise associated with the laser primarily takes the form of laser speckle, which is the oscillation of the intensity profile for a laser when it is reflected from the surface of the part and is caused by coherency of the laser. One way toreduce the speckle noise is by choosing an appropriate viewing system. By changing the size of the aperture, the size of the speckle changes, the larger the aperture, the smaller the size of the speckle. However, in this case there's a depth of fieldtradeoff.Therefore there is need for an improved image processing technique for reducing the speckle noise and also the sampling error in structured light applications.BRIEF DESCRIPTION OF THE INVENTIONBriefly, in accordance with one aspect of the present invention, an image processing method for structured light profiling includes sampling an image of a structured light pattern to obtain an intensity distribution. A number of sets of sampledpoints are then selected from the intensity distribution, where each of the respective sets includes a number of sampled points. Each of the sets of sampled points is then fitted to a respective distribution function. Finally, the distributionfunctions are filtered to select a representative distribution function for the intensity distribution.In accordance with another aspect of the present invention, a system for obtaining a three dimensional profile of a part using a structured light pattern includes a source of structured light positioned at a predetermined distance from a part,where the source projects a beam of structured light to illuminate the part. The system also includes at least one imaging device configured to acquire an image of a structured light pattern of the part, where the at least one imaging device ispositioned such that an angle of view of the imaging device is different from an angle of illumination of the source. A processor is coupled to the at least one imaging device, and the processor is configured for a variety of tasks including samplingthe image of a structured light pattern to obtain an intensity distribution; selecting a number of sets of sampled points from the intensity distribution, where each of the sets includes at least three sampled points; fitting each of the sets of sampledpoints to respective Gaussian distribution functions; extracting a center for each of the Gaussian distribution functions; filtering the Gaussian distribution functions by using the centers to select a representative distribution function for theintensity distribution; and reconstructing a three dimensional profile of the part, using the representative distribution function. BRIEF DESCRIPTION OF THE DRAWINGSThese and other features, aspects, and advantages of the present invention will become better understood when the following detailed description is read with reference to the accompanying drawings in which like characters represent like partsthroughout the drawings, wherein:FIG. 1 illustrates a system for obtaining a three dimensional profile of a part using a structured light source;FIG. 2 illustrates a set of graphs for three different orientations of a curved surface being viewed by two cameras;FIG. 3 is a graphical representation of a number of sampling points at three different locations of an intensity distribution;FIG. 4 is a graphical representation of sampled intensity profiles for the three locations illustrated in FIG. 3;FIG. 5 is a graphical representation of a sampled intensity profile with speckle noise and sampling error;FIG. 6 is a flowchart that illustrates one aspect of an image processing method for structured light profiling for use in the system of FIG. 1;FIG. 7 illustrates a flowchart for determining the position of laser stripe by using a center of the representative distribution function according to an aspect of the image processing technique; andFIG. 8 illustrates a flowchart for determining the position of laser stripe by using centers of distribution functions according to an aspect of the image processing method.DETAILED DESCRIPTION OF THE INVENTIONFIG. 1 illustrates a system 10 for obtaining a three dimensional profile of a part 16 using a structured light pattern 18. The system includes a source 12 of structured light 14 positioned at a predetermined distance from a part. The source 12projects a beam of structured light 14 to illuminate the part 16. In the illustrated embodiment, a laser is used as the structured light source 14. However other mediums like white light may also be used. At least one imaging device 20 is configuredto acquire an image 28 of a structured light pattern of the part 16 by viewing along the direction 22. The imaging device 20 is positioned such that an angle of view (.alpha.) of the imaging device 20 is different from an angle of illumination (.beta.)of the source 12. A processor 26 is coupled through 24 to the imaging device 20 and is configured for reconstructing the three dimensional profile of the part using the methods described hereinbelow.As would be appreciated by those skilled in the art, the intensity profile for a flat surface in the case of illumination with Gaussian profile by source 12 and any viewing angle of the imaging device 20, for any tilt angle (.theta.) of thesurface of the part 16, remains a Gaussian distribution. For such a symmetric intensity profile on a flat surface, it will be well appreciated by those skilled in the art that sampling of the symmetric point and applying zero-crossing which is a secondderivative on the profile, will yield the center corresponding to the maximum intensity.However, in the case of a complex part, where the surface is not flat, the intensity profile is no longer symmetric. In case of a curved surface, the curvature of the surface may be approximated by a cylinder to analyze the intensity profilesfor a Gaussian distribution. It is found that the degree of asymmetry .delta. (the maximum distance difference between two points with the same intensity to the maximum intensity point) is small and hence the effect on image due to this asymmetry issmall. Therefore the intensity profile near a cylindrical neighborhood in a curved surface can also be assumed to be a Gaussian distribution. In this case, considering that the maximum intensity point is sampled, the image bias after applyingzero-crossing is also small, of the order of 10.sup.-2. This is shown by a set of graphs in FIG. 2, for three different orientations of a curved surface at angles -40 degrees, 0 degrees and 40 degrees being viewed by two cameras, with the views beingdepicted as a left view and a right view respectively. Also in FIG. 2, {right arrow over (V)} is the viewing vector unit determined by viewing angle .alpha., {right arrow over (N)} is the surface normal unit determined by surface tilt angle .theta. and{right arrow over (B)} is the vector unit of the beam center line determined by beam divergent angle .beta.. The degree of asymmetry .delta. is shown in FIG. 2 for the left and right views for different tilt angles of the curved surface.However, if the maximum intensity point is not sampled, i.e. there is sampling error, zero-crossing does not give accurate results and leads to a larger image bias as shown by the set of graphical representations in FIG. 3 and FIG. 4. FIG. 3depicts a number of sampling points 43 at three different locations 44, 46, 48 of an intensity distribution 36. FIG. 4 depicts the sampled intensity profiles 50, 52 and 54 for the locations 44, 46 and 48 respectively. These profiles 52 and 54 (as isclear from FIG. 4) are not a true representation of the intensity profile 36 due to the sampling error. Further, speckle also leads to enhanced image bias error, as it causes the intensity profile to change as shown in FIG. 5. Therefore, if there isspeckle noise, the intensity profile 36 of FIG. 4 will appear as the sampled intensity profile 58 with sampling error and added speckle noise at 56.In order to reduce sampling error and speckle noise, an imaging processing method is provided. FIG. 6 illustrates one aspect of an image processing method for structured light profiling for use in the system 10 of FIG. 1. The method includesacquiring the image of the structured light pattern of a part. The structured light pattern includes at least one laser stripe, and the intensity distribution is an intensity profile across one of the laser stripes. Referring to the flowchart of FIG.6, the method starts at 60 with an image 28, which is sampled to obtain an intensity distribution of the structured light pattern at step 62. Next, at step 64, a number of sets of sampled points are selected from the intensity distribution, where eachof the respective sets includes a number of sampled points. In one example each of the sets of sampled points includes at least three sampled points. In Step 66 each of the sets of sampled points is fit to a respective distribution function. In oneexample each of the distribution functions is a Gaussian distribution. It is found that fitting the intensity profile for a curved surface with a Gaussian distribution leads to reducing sampling error. In Step 68 the distribution functions are filteredto select a representative distribution function with reduced effect of speckle noise for the intensity distribution. The method concludes at 69 with an output of the representative distribution function, which is further processed to reconstruct thethree dimensional profile of the part.FIG. 7 illustrates additional aspects of the image processing method. For this embodiment, the image processing method also includes extracting a center of the representative distribution function, which is shown as step 72 in FIG. 7. The input70 to the process illustrated in FIG. 7 is the representative distribution function, and the output 74 is the center of the representative distribution function depicting the position of the laser stripe on the complex part. The center of therepresentative distribution function may be extracted, for example, using zero-crossing. The exemplary fitting method, according to an aspect of the present technique, fits sets of three sampled points to respective distribution functions, in order tocalculate the center of the representative distribution function, as mentioned above. As would be appreciated by one skilled in the art, more than three sampled points can be used as well.FIG. 8 illustrates another embodiment of the image processing method. For this embodiment, a center is extracted for each of the distribution functions as shown in step 82 of FIG. 8. The centers may be extracted using zero-crossing. The input80 to this process will be the distribution functions, which are generated at step 66 in FIG. 6. The filtering step includes using these centers in the subsequent step 84. For example, multiple centers may be calculated using different sets of threepoints from the intensity profile. If the multiple centers are close, the average or median can be used as a fitted center of the Gaussian distribution. Namely, for this example, the filtering step includes selecting a median of the centers, and therepresentative distribution function corresponds to the median of the centers. A histogram may be used to determine the median. Alternatively, the centers may be sorted to determine the median. As will be appreciated by one skilled in the art, thesefiltering techniques select the representative center for the distribution functions that depicts the position of laser stripe center on the complex part.In case there is a speckle which can be observed as an outlier on the intensity profile, the corresponding center using the outlier can be easily discarded and the speckle noise can be addressed. This is indicated in FIG. 8, in which step 84 offiltering includes rejecting a number of distribution functions that involve outliers. A filtering example is discussed with reference to FIG. 5. As shown in FIG. 5, there is an added speckle at sampling point 4 (total 9 points) and the calculatedseven centers are, 0, -6.06, 0.16, -0.32, 0, 0, 0. For this filtering example, a closest group amongst the calculated centers is selected and zero is selected. The statistical analysis methods can further use histograms to screen speckle noise usingthe median value. Beneficially, this fitting requires only three sampled points at a time to calculate multiple values for the centers which can be used for filtering without compromising computation speed.As noted above, the flowcharts illustrated in FIG. 6, FIG. 7, and FIG. 8 describe the aspects of the method discussed hereinabove. The foregoing flow charts also show the functionality and operation of the method and the system forreconstructing a three dimensional profile of a part. In this regard, each block/component represents a module, segment, or portion of code, which comprises one or more executable instructions for implementing the specified logical functions. It shouldalso be noted that in some alternative implementations, the functions noted in the blocks may occur out of the order noted in the figures or, for example, may in fact be executed substantially concurrently or in the reverse order, depending upon thefunctionality involved. Also, one of ordinary skill in the art will recognize that additional blocks may be added. Furthermore, the functions can be implemented in programming languages such as C++ or JAVA; however, other languages can be used.The various embodiments and aspects of the invention described above comprise an ordered listing of executable instructions for implementing logical functions. The ordered listing can be embodied in any computer-readable medium for use by or inconnection with a computer-based system that can retrieve the instructions and execute them. In the context of this application, the computer-readable medium can be any means that can contain, store, communicate, propagate, transmit or transport theinstructions. The computer readable medium can be an electronic, a magnetic, an optical, an electromagnetic, or an infrared system, apparatus, or device. An illustrative, but non-exhaustive list of computer-readable mediums can include an electricalconnection (electronic) having one or more wires, a portable computer diskette (magnetic), a random access memory (RAM) (magnetic), a read-only memory (ROM) (magnetic), an erasable programmable read-only memory (EPROM or Flash memory) (magnetic), anoptical fiber (optical), and a portable compact disc read-only memory (CDROM) (optical).Note that the computer readable medium may comprise paper or another suitable medium upon which the instructions are printed. For instance, the instructions can be electronically captured via optical scanning of the paper or other medium, thencompiled, interpreted or otherwise processed in a suitable manner if necessary, and then stored in a computer memory.The various aspects of the technique described hereinabove have utility in industrial as well as medical environments. The methods can be used for non-contact measurement of complex parts e.g. aircraft parts for inspection, in the extrusionprocess in the steel industry and other high temperature manufacturing environments where contact measurement of a part is difficult. These methods are also useful in medical fields for surgery planning, where these may be used for profiling thedifferent parts of a human body to have precision in surgery.Although only certain features of the invention have been illustrated and described herein, many modifications and changes will occur to those skilled in the art. It is, therefore, to be understood that the appended claims are intended to coverall such modifications and changes as fall within the true spirit of the invention.* * * * *
Register and get access to “Method And System For Image Processing For Structured Light Profiling Of A Part - Patent 7302109”
Login to continue downloading: “Method And System For Image Processing For Structured Light Profiling Of A Part - Patent 7302109”
7+ Documents Related to “Method And System For Image Processing For Structured Light Profiling Of A Part - Patent 7302109”