Machine vision methods for identifying collinear sets of points from an image

Information

  • Patent Grant
  • 6075881
  • Patent Number
    6,075,881
  • Date Filed
    Tuesday, March 18, 1997
    27 years ago
  • Date Issued
    Tuesday, June 13, 2000
    24 years ago
Abstract
The invention provides a two-phase method for identifying sets of collinear points in data points, e.g., derived from an image. During a first phase, the "angle-finding" phase, the angular orientation of the lines defined by the sets of points is estimated. During a second phase, the "line-finding" phase, sets of collinear points aligned at that angle are identified. The invention has application, for example, in finding parallel and perpendicular point sets in an image.
Description

BACKGROUND OF THE INVENTION
The invention pertains to machine vision and, more particularly, to methods for identifying sets of collinear points from an image.
Machine vision refers to the automated analysis of an image to determine characteristics of objects and other features shown therein. It is often employed in automated manufacturing lines, where images of components are analyzed to determine placement and alignment prior to assembly. Machine vision is also used for quality assurance. For example, in the pharmaceutical and food packing industries, images of packages are analyzed to insure that product labels, lot numbers, "freshness" dates, and the like, are properly positioned and legible.
In many machine vision applications, it is essential to identify sets of points that define lines and, more particularly, parallel or perpendicular lines. This is typically critical in determining the position and orientation of an imaged object. In electronic circuit board assembly, for example, integrated circuit chips and sockets must be precisely positioned before they can be soldered into place. Metallic solder pads feature predominantly in images of these components and, as such, are often used in determining chip location and orientation.
Objects whose image is to be analyzed by machine vision typically include calibration targets. Often a cross-shaped symbol or parallel array of dots, such a target facilitates determining the position and orientation of the object in the image. A calibration target can also be used to correlate the image reference frame with that of the motion stage or conveyor belt on which the imaged object is placed.
Regardless of whether they are formed by solder pads, calibration targets, or otherwise, there are several known approaches to the problem of finding parallel and perpendicular lines from an image. One of the best known techniques is the Hough transform. According to that technique, every point in the set "votes" into an accumulator array in parameter space. Parameter pairs with the sufficiently many votes are identified as lines. Though oft used, the Hough transform can require significant computation resources and time to carry out these tasks.
For small data sets, the art also suggests exhaustively partitioning the set of points into all possible combinations of candidate lines, and evaluating the "goodness" of each one, e.g., with a least squares test. This procedure is impractical for large data sets because the computational complexity grows rapidly with the number of points.
An object of this invention is to provide an improved machine vision methods and, particularly, improved methods for identifying sets of collinear points from an image.
A related object is to provide methods for identifying parallel and perpendicular lines from an image.
A further object is to provide such methods as can be readily implemented on conventional digital data processors or other conventional machine vision analysis equipment.
Yet still another object of the invention is to provide such methods that can rapidly analyze images without undue consumption of resources.
SUMMARY OF THE INVENTION
The foregoing objects are among those attained by the invention, which provides a two-phase method for identifying sets of collinear points in data points, e.g., extracted from an image. During a first, "angle-finding" phase, the angular orientation of the predominant line defined by the sets of points is estimated. During a second, "line-finding" phase, sets of collinear points aligned at that angle (or at angles offset, e.g., by 90.degree. therefrom) are identified.
More particularly, one aspect of the invention provides a machine vision method for identifying one or more sets of collinear points in data points, e.g., extracted from an image. The method includes finding an estimated angle of the sets by taking projections of the data points at multiple angles, e.g., at 10.degree. increments between 0.degree. and 170.degree., and estimating the angle of the predominant sets from the projection(s) with the best-formed, or highest density, peak(s).
The method further includes identifying the sets of collinear points by taking a projection of the data points at the estimated angle and determining the locations of the sets from the locations of the peaks in that projection. Once those sets have been identified, characteristics of the lines they define can be output, e.g., for use by other machine vision tools or apparatus. Those characteristics can include, for example, the numbers of lines, the angles and locations of the lines, the points that make up the lines, etc.
In a related aspect, the invention provides a method as described above in which the estimated angle determined in the angle-finding phase is a function of the two (or more) candidate angles that resulted in projections with the highest density peaks. Thus, for example, where projections at 40.degree. and 50.degree. produce the highest density peaks, the estimated angle is interpolated between those angles, e.g., in accord with the relative density or height of those peaks.
In another related aspect of the invention, a method as described above utilizes a filter to determine the density of the peaks in the projections and, more particularly, to determine which of the projections has the highest density peaks. The filter can be, for example, a "boxcar" filter sized in accord with the estimated width of the peaks that the sets of collinear points are expected to produce in the projections.
Further aspects of the invention provide for operator determination of the candidate angles utilized in determining the estimated angle. Alternatively, those angles can be determined in accord with the expected length and spacing of the lines formed by the sets of collinear points.
According to another aspect of the invention, the estimated angle can be determined in successive approximations. A first determination of the estimated angle can be made by taking projections of the data points at a first set of candidate angles covering a first angular range at a first angular resolution. Second and successive determinations can be made over narrower angular ranges covering successively smaller ranges, e.g., bracketing the estimate from each prior pass, with finer angular resolutions. For example, taking projections at 10.degree. increments between 0.degree. and 170.degree. can provide a first determination of the estimated angle. Based on that determination, a second determination can be made by taking projections at 1.degree. increments between 45.degree. and 55.degree., and so forth.
As in the angle-finding phase, a further aspect of the invention provides that the line-finding phase can use a filter, such as a boxcar filter, to find the highest density peaks in the projection taken at the estimated angle. That filter, too, can be sized according to the estimated width of the peaks that the sets of collinear points are expected to produce in the projection.
These and other aspects of the invention are evident in the drawings and in the detailed description that follows.





BRIEF DESCRIPTION OF THE DRAWINGS
A better understanding of the invention may be attained by reference to the drawings in which:
FIG. 1 depicts a machine vision system for practice of the invention;
FIG. 2 is a flow chart summarizing operation of a method according to the invention for identifying sets of collinear points in data points, e.g., extracted from an image;
FIG. 3 is a flow chart detailing operation of an angle-finding phase of a method according to the invention;
FIG. 4 is a flow chart detailing operation of an line-finding phase of a method according to the invention;
FIGS. 5A-5B depict projections of collinear point sets at different candidate angles during an angle-finding phase of a method according to the invention;
FIGS. 6A-6G depict the boxcar filtering of a projection array during an angle-finding phase of a method according to the invention; and
FIGS. 7A-7C depict identification of lines corresponding to collinear sets of points during a line-finding phase of a method according to the invention.





DETAILED DESCRIPTION OF THE ILLUSTRATED EMBODIMENT
FIG. 1 illustrates a machine system 10 of the type in which the invention is practiced. The system 10 includes a capturing device 16, such as a conventional video camera or scanner, that generates an image of an object 12 that bears one or more sets of collinear points--in this case, two sets of collinear points 14 arranged perpendicularly to one another. Image data (or pixels) generated by the capturing device 16 represent, in the conventional manner, the image intensity (e.g., color or brightness) of each point in the scene at the resolution of the capturing device 16.
The digital image data is transmitted from capturing device 16 via a communications path 18 to an image analysis system 20. This can be a conventional digital data processor, or a vision processing system of the type commercially available, for example, from the assignee hereof, Cognex Corporation, programmed in accord with the teachings hereof to identify sets of collinear points from the image data.
The image analysis system 20 may have one or more central processing units 22, main memory 24, input-output system 26, and disc drive (or other mass storage device) 28, all of the conventional type.
The system 20 and, more particularly, central processing unit 22, is configured by programming instructions according to the teachings hereof for identifying sets of collinear points, as described in further detail below. Those skilled in the art will appreciate that, in addition to implementation on a programmable digital data processor, the methods and apparatus taught herein can be implemented in special purpose hardware.
FIG. 2 is a flow chart summarizing operation of a method according to the invention for identifying sets of collinear points from an image. In step 28, the method inputs data points to be analyzed. These may be extracted from an image (e.g., by a vision processing tool, such as by connectivity analysis, normalized correlation, etc.), or they may be artificially derived, supplied by the operator, or otherwise obtained. The data points contain sets of collinear points, as well as non-collinear points and noise. The sets of collinear points may define a single line or sets of parallel and/or perpendicular lines. The sets may also define sets of lines at angles other than 0.degree. (parallel) and 90.degree. (perpendicular).
In step 30, data structures used in the angle-finding and line-finding phases of the method are allocated and initialized. This results in creation of pointers that are passed to functions used in these phases. For example, pointers to the image data points may be cached in the private data structure, so that those points themselves are not be destroyed until the private structure itself is destroyed.
As those skilled in the art will appreciate on reading the description below, different data points may be specified for use in the angle-finding and line-finding phases. For example, points taken from blob centroids that are close to image edges may give skewed results in the angle-finding phase, even though those points may be used in the line-finding phase.
The projection bins used in the angle-finding and line-finding phases are preferably allocated and initialized during the initialization phase. Adequate numbers of bins are allocated for the widest range of points possible, and for the smallest bin size possible.
In step 32, the method performs an angle-finding step to estimate the angular orientation of predominant line defined by sets of collinear points in the image. Details of this step are shown in FIG. 3 and described below.
In the case of an image with a single set of collinear points, the estimated angle is the angle of the line defined by those points. In the case of an image with multiple sets of parallel collinear points, the estimated angle is the angle of the predominant line defined by those sets. In the case of an image with one or more perpendicular sets of collinear points, the estimated angle is the angle of the predominant line defining one axis of the sets, e.g., the predominant angle of the "horizontal" lines or of the "vertical" lines.
In step 34, the method performs a line-finding step to identify sets of collinear points in the data points aligned with the estimated angle. Details of this step are shown in FIG. 4 and described below.
In step 36, the method outputs a characteristic of the sets identified in step 34. This can be to another machine vision tool or to a human operator. Those characteristics include the numbers of lines (or, more properly, line segments), the angles of the lines, the coordinates of the data points defining those lines, the locations of those lines, e.g., distance and direction from the image origin, etc.
Although the discussion herein is principally directed to finding parallel or perpendicular lines, those skilled in the art will appreciate that the teachings herein can be applied in finding collinear sets of points, regardless of their angular orientation to one another. Thus, by way of non-limiting example, it will be appreciated that the invention can be applied to finding points defining edges of similarly triangular or octagonal objects from an image.
Furthermore, although the discussion herein is directed to finding sets of collinear points from an image, it will appreciated that the teachings are equally applicable to finding those sets in portions of an image, or windowed image.
FIG. 3 is a flow chart depicting the methodology of a preferred angle-finding phase in a system according to the invention. In this phase, the method finds the angle of the predominant or strongest line in a set of data points. To do so, the entire set of points is projected at multiple angles. At each projection angle, a count of how many points project at a particular distance is accumulated in an array of projection bins. This one-dimensional accumulator space is then filtered to find the maximum density of projected points. The projection angle that produces a peak with the highest density is the best projection angle: the strongest line is aligned with that angle.
In step 40, the method establishes a bin size and filter size for projections taken during the angle-finding phase. In a preferred embodiment, this information is provided by the operator, e.g., via a keyboard. Alternatively, it can be pre-stored in data or in any other known manner.
Both the angle-finding phase and the line-finding phase (discussed below) operate by projecting the data points into bins, then filtering the bins with an unweighted filter to find the maximum density of projected points. It is assumed that the maximum density of projected points under the filter is caused by sets of collinear points, or lines. Therefore, a filter is chosen that is neither too large (so points from other lines or noise points are not included under the filter) nor too small (so all points from a complete line are included under the filter).
The number of bins under the filter is selected so the filter only covers a few bins. More bins under the filter will slow performance without adding extra accuracy to the answer. Ideally, the peaks resulting of projections of each set of collinear points in the data points should fall in a handful, say, three to five, bins. This avoids having to use very large filters.
The filter size is preferably large enough to capture any peaks resulting from projection of a line in the data points. For example, if the angle-finding phase takes projections every x degrees, the worst difference between the true angle of the sets of collinear points and the closest sample projection is x/2. The filter size is, accordingly, set to (line length) * sin(x/2). In this regard, the line length is assumed to be the shortest expected length of a complete line in the image. If the number chosen for this length is too small, chance alignments of points might be mistaken as a line. If it is too large, lines with only a few points might be missed.
In step 42, the method determines an angular range and resolution for projections taken of the data points in order to determine the estimated angle. In one embodiment, this information is input by the operator, e.g., via a keyboard. Alternatively, it can be pre-stored in data or source code, or in any other known manner.
In order to search the entire coordinate space, the range can be selected to be 0.degree. to 170.degree., and the angle increment can be selected to be 10.degree.. The angular resolution or increment is be chosen small enough so that the "real" maximum density in the projection is greater than the expected "chance" maximum density resulting from the projection of overlapping sets of collinear points.
Generally speaking, it is more efficient to run the angle finding phase at several search resolutions, iteratively narrowing the range to the correct angle. For example, if searching all possible angles to an accuracy of 0.25.degree., one could search between 0.degree. and 170.degree. at 0.5.degree. angle increments, causing 340 projections of the input points. It is more efficient, however, to make three searches. An initial search is made between 0.degree. and 170.degree., at 10.degree. increments, to find the angle to the closest 5 degrees. A second search is then made in a +/-10.degree. range bracketing the found angle at 3.degree. increments. A third search is made in a +/-1.5.degree. range bracketing the second found angle at 0.5.degree. increments. By taking three passes, each input point is projected only 32 times.
In order to determine the minimum number of samples that must be evaluated to get 0.5.degree. accuracy, the following analysis may be used.
To begin, S is defined as the number of sample projections used for a given angular range. For example, at an increment of 10 degrees, S=5 samples are needed to cover the search range from 10 degrees to 50 degrees. It will be noted that, although there are five samples, the range is divided into only four 10.degree. subranges.
L is defined as the number of search levels, i.e., the number of different angular ranges that will be used. R is defined as the number of retained increments, i.e., the number of subranges that are re-searched at each successive search level. In a preferred embodiment, R is three. F is defined as the fraction of the search range that will be re-searched at the next successive search level. In a preferred embodiment, where R is three, f=3/(S-1). More generally, f=R/(S-1).
In order to reduce the initial search range by a factor of 100 from, say, 50 degrees to 0.5 degrees), f.sup.L =1/100, or (R/(S-1)).sup.L =1/100. In order to reduce the search range using the fewest total sample projections, the total number of sample projections, S*L, must be minimized. These conditions are solved in the following tables, which show obvious minima for S*L
______________________________________Using 3.0 retained incrementsEvaluation table for final ratio = 100Number of Samples Levels Needed Total Evaluations______________________________________ 5 16.0078 80.04 6 9.0152 54.09 7 6.6439 46.51 8 5.4351 43.48 9 4.6952 42.2610 4.1918 41.9211 3.8250 42.0712 3.5444 42.5313 3.3219 43.1914 3.1406 43.9715 2.9895 44.84______________________________________
In step 44, the method takes a projection of the data points at the first of the candidate angles in the angular range established in step 42. The projection, which is a well known function in the machine vision arts, can be performed using any conventional projection vision tool. One technique for taking a projection is described in Horn, Robot Vision (M.I.T. 1986), Chapter 3, Section 3.3, pages 53-57, the teachings of which are incorporated herein by reference.
In yet another preferred embodiment, the projection of a point (x, y) into a bin at angle "theta" is determined in accord with the relation:
bin#=(x*cos+y*sin)/bs,
where
bs=bin.sub.-- size,
sin=sin(theta) and
cos=cos(theta).
This expression can be rearranged as follows:
bin#=x*(cos/bs)+y*(sin/bs);
=(x/k)*(cos*k/bs)+(y/k)*(sin*k/bs);
={(x/k)*[(2 32)*cos*k/bs]+(y/k)*[(2 32)*sin*k/bs]}*(2 -32);
This latter expression is advantageous in that the quantities in square brackets need to be computed only once for a given projection angle. (x/k) and (y/k) can be computed once per field of view (FOV). A more important reason is that, given an appropriate choice of `k`, the following four quantities can always be expressed as 32-bit integers with no overflow and little or no loss of precision:
(x/k), (y/k), [(2 32)*cos*k/bs], [(2 32)*sin*k/bs]
This makes the math easy to do using only two integer multiplies and a couple of integer additions.
To pick an appropriate value for k, one must consider what values make sense for the quantities (x/k) and (y/k), as well as the quantities [(2 32)*cos*k/bs] and [(2 32)*sin*k/bs]. In a preferred embodiment for use in machine vision inspection of surface mounted devices, k is chosen as one-fiftieth of a pixel.
FIGS. 5A-5B illustrate projections taken at two candidate angle (80.degree., 90.degree.) of three sets of collinear points 70, each comprising five points. At each angle, the points 70 project into different combinations of bins of projection arrays 72A, 72B. Those different combinations result in peaks of different densities.
Thus, as shown in FIG. 5A, the five points in each of the three sets 70 project into respective groups of bins in array 72A. This is represented by the three groups 74A of projected-points, which are represented in the drawing by groups squares, each of which designates the projection of a single respective point 70 into the array 72A. Each of the groups, which is referred to herein as a peak, has a maximum density of two, as described below. Those skilled in the art will appreciate that these "peaks" represent portions of the array 72A which exhibit higher density of projected-points relative to the ambient level thereof (e.g., resulting from noise and non-collinear sets of points). For example, groups 74A and 74B represents peaks, as that term is used herein, insofar as they exhibit higher density of projected-points relative to the ambient level.
Points 70 project into different peaks when projected at different angles. Thus, as shown in FIG. 5B, the points project into peaks 74B of density three when projected at 90.degree..
Referring to FIG. 3, in step 46, the method determines the density of the peak resulting from the projection of step 44. As evident in the discussion above, the density of peaks will be the highest for in those projections that are taken at angles most closely aligned with the angle of the predominant set of collinear points in the image.
Referring to FIGS. 5A-5B, the density of the peaks 74A-74B is determined by passing a filter over projection arrays 72A-72B. This is illustrated in FIGS. 6A-6G, which show movement of a boxcar filter 76 over one of the peaks 74B in projection array 72B.
As the filter 76 moves, the method counts how many projected-points (i.e., squares) it covers. For example, in FIG. 6A, the filter 76 is not over the peak 74B; hence, it appears to have a count and density of zero. As the filter 76 is moved, however, it becomes apparent that the peak has a definite density. Particularly, in FIG. 6D, the filter 76 is centered over the peak 74B, revealing its true count and density of three.
Those skilled in the art will appreciate that techniques other than those shown in FIGS. 6A-6G and discussed above can be used to determine which of the projections is best aligned with the sets of collinear points 70 of FIGS. 5A-5B. For example, filters other than boxcar filters can be passed over the projection arrays 72A-72B. In other embodiments of the invention, filters are not used in determining the quality of the peaks. Instead, only the heights of the peaks are tracked.
Referring back to step 46, of FIG. 3, the method stores an indication of the maximum density of the projection taken in step 44. That information is later used in determining which of the projections was taken along an angle best aligned with the sets of collinear points in the image.
In step 48, the method determines whether further candidate angles are required to cover the angular range established in step 42. If so, the method loops back to step 44. From there, steps 44 and 46 are repeated using the next angle at the designated angular resolution. If no further candidate angles are required, the method proceeds to step 50.
In step 50, the method determines which of the projections was taken at an angle best aligned with the sets of collinear points in the image. This is accomplished by comparing the peak densities determined in step 46 for each of the projections. In one embodiment of the invention, the estimated angle of the collinear points is assumed to be equal to the candidate angle of the single projection with the highest density peak. In a preferred embodiment, the estimated angle is interpolated from the candidate angles associated with the two highest density peaks, as follows: ##EQU1## where, .THETA..sub.MAX is candidate angle associated with the highest density peak;
.THETA..sub.INCR is the angular difference between the candidate angles associated with the two highest density peaks;
dMAX is the density of the highest density peak; and
dOTHER is the density of the second highest density peak;
In step 52, the method determines whether any further passes of the angle finding technique will be made to improve the resolution of the estimated angle. If so, the method loops back to step 42, where an angular range and resolution are set as described above to bracket the previously determined estimated angle. If not further passes are required, the method utilizes the estimated angle in the line finding phase.
FIG. 4 is a flow chart depicting the methodology of a preferred line-finding phase in a system according to the invention. In this phase, the method identifies the sets of collinear points in the data points aligned with the estimated angle (i.e., the angle of the sets of predominant collinear points, as determined in the angle-finding phase), as well as lines at 90.degree. to the estimated angle.
It does this by projecting the points at two angles: the input angle and the input angle plus 90 degrees. Lines are located at peaks in the projection. However, this phase is interested in finding all lines at the angle, not just the strongest line. It also gathers additional information about lines: the identities (e.g., coordinates) of the points that make them up, the number of points that make them up, and/or their distances along the line.
In step 54, the method determines whether to identify sets of collinear points aligned with the estimated angle or whether to identify those at 90.degree. to that angle. This determination can be made based on user input, preprogrammed instructions, etc. In a preferred embodiment, the method finds both types of lines: first it finds those aligned with the estimated angle, then it finds those at 90.degree. to that angle.
Depending on the flow effected by step 54, in step 56 the method takes a projection of the data points at the estimated angle, while in step 58 the method takes a projection at the estimated angle plus 90.degree.. This projection, which is conducted as described above, is illustrated in FIG. 7A. There, the sets of collinear points 78 are projected at the estimated angle (90.degree.) into projection array 80. Peaks 82 arising from the projection are shown as groupings of squares in that array 80.
Regardless of the angle at which the projection is taken, in step 60 of FIG. 4 the method applies a filter to identify all peaks in the projection and, thereby, to identify all sets of points aligned with the angle of the projection. The filter applied in this step can be identical to that discussed above in connection with the angle-finding phase; however, it need not be. This is illustrated in FIG. 7B, where boxcar filters 86 of the type described above are shown centered above the peaks 82. Though illustrated otherwise for the sake of convenience, in a preferred embodiment, a single boxcar filter 86 is passed over the entire projection array 80.
To avoid "false positives," the filtering step 60 of FIG. 4 preferably imposes a minimum threshold on the density and/or height required for peak to be considered as being associated with a line. Thus, for example, if a peak has a density of four or less (i.e., indicating that there are four or less points in the corresponding set that are aligned with the projection angle), it may be ignored for purposes of the line-finding phase. Likewise, a peak may be ignored if its height falls below a minimum threshold.
The filtering step 60 preferably stores identifying information for the sets of collinear points that contributed to the peaks that exceed the aforementioned thresholds. In a preferred embodiment, this is effected by storing the coordinates of each set of data points that contributed to those peaks.
In step 62, the method merges sets of collinear points that are sufficiently close to constitute a single line. This can be accomplished as a final phase of the filtering step, i.e., by comparing the distances between the peaks identified by the filter. It is preferably accomplished by comparing the distances between the points that make up the respective sets of collinear points. It can, alternatively, be accomplished by identifying the lines defined by those points and, in turn, comparing the distances between the lines.
The effect of merging step 62 is shown in FIG. 7C. There a single set of the collinear points 78 is identified with line segment 88, while two closely-spaced sets are associated with line segment 90.
In step 64, the method determines whether further sets of collinear points are to be identified in the data points. Thus, for example, if a first pass of the line-finding phase was dedicated to finding sets of collinear points aligned with the estimated angle, step 64 can route flow back to step 54 to find sets aligned with the estimated angle plus 90.degree..
At the conclusion of the line-finding phase, the method returns, e.g., to the operator or to another machine vision tool, the identity of the line segments 88 and 90 identified during that phase. This can include the slopes of the line segments, the coordinates of the points that make them up, the length of the lines segments, etc.
Described herein is a method for identifying sets of collinear points from an image meeting the objects set forth above. It will be appreciated that the embodiments shown in the drawings and described above are merely illustrative and that other embodiments incorporating modifications thereto fall within the scope of the invention. Thus, for example, it will be appreciated that the invention can be used to find not only parallel and perpendicular collinear point sets, but also collinear points sets at other angles. In this regard, for example, step 58 of FIG. 4 can be utilized to take projections at the estimated angle plus 45.degree., 90.degree., 135.degree., or any other angle.
Claims
  • 1. A machine vision method for identifying one or more sets of collinear points from data points, the method comprising the steps of:
  • A. searching among projections taken at a plurality of candidate angles determine an estimated angle of the sets of collinear points;
  • B. identifying the sets of collinear points from a projection of the data points at an angle based on the estimated angle; and
  • C. outputting a characteristic of a line formed by at least one of the sets of collinear points.
  • 2. A method according to claim 1, wherein step (A) includes taking a projection of the data points at each of the candidate angles, and identifying one or more of the projections that produce highest density peaks.
  • 3. A method according to claim 2, wherein step (A) includes the step of determining the estimated angle of the sets of collinear points from the candidate angles of plural projections that produce the highest density peaks.
  • 4. A method according to claim 3, wherein step (A) includes the step of determining the estimated angle of the sets of collinear points by interpolating between the candidate angles of two projections that produce the two highest density peaks.
  • 5. A method according to claim 2, wherein step (A) includes the step of filtering the projections to find those that produce the highest density peaks.
  • 6. A method according to claim 5, wherein step (A) includes the step of using a filter sized in accord with an estimated width of a peak produced by a projection of the sets of collinear points at the candidate angle.
  • 7. A method according to claim 1, wherein step (A) includes the step of defining the candidate angles in accord with any of (i) operator input, and (ii) an expected length and spacing of lines defined by the collinear sets of points.
  • 8. A method according to claim 1, wherein step (A) includes the steps of
  • making a first determination of the estimated angle of the sets of collinear points by taking projections of the data points at a first set of candidate angles, that first set of candidate angles covering a first angular range at a first angular resolution, and
  • making successive determinations of the estimated angle of the sets of collinear points by taking projections of the data points at successive sets of candidate angles, each successive set of candidate angles covering a successively narrow angular range and being at a successively finer angular resolution than a prior set of candidate angles.
  • 9. A method according to claim 1, wherein step (B) includes the steps of taking a projection of the data points at the estimated angle and finding one or more peaks in that projection corresponding to respective sets of collinear points.
  • 10. A method according to claim 9, wherein step (B) includes the step of filtering the projection to find those peaks.
  • 11. A method according to claim 10, wherein step (B) includes the step of using a filter sized in accord with an estimated width of a peak resulting from taking a projection of the set of collinear points at the estimated angle.
  • 12. A method according to claim 9, wherein step (B) includes the step of merging into one set of collinear points two or more sets of collinear points that correspond to adjacent peaks in the projection.
  • 13. A method according to claim 1, wherein step (C) includes the step of outputting, as the characteristic of a line formed by a set of collinear points, any of an angle of that line, a position of that line, and identities or locations of points making up that line.
  • 14. A machine vision method for identifying one or more sets of collinear points from data points, the method comprising the steps of:
  • A. determining an estimated angle of the sets of collinear points by taking projections of the data points at a plurality of candidate angles, and by searching among the projection to find one or more projections that produce highest density peaks,
  • B. identifying one or more sets of collinear points by taking a projection of the data points at an angle based on the estimated angle, and finding one or more peaks in that projection corresponding to respective sets of collinear points, and
  • C. outputting a characteristic of a line formed by at least one of the sets of collinear points.
  • 15. A machine vision method for identifying a set of collinear points from an data points, the method comprising the steps of:
  • A. making a first determination of an estimated angle of the sets of collinear points by taking projections of the data points at a first set of candidate angles covering a first angular range at a first angular resolution, and finding one or more projections producing highest density peaks,
  • B. making one or more successive determinations of the estimated angle by taking successive projections of the data points at successive sets of candidate angles, and finding one or more projections producing highest density peaks for each such set, each set of candidate angles covering a successively narrow angular range and being at a successively finer angular resolution than a prior set of candidate angles,
  • C. identifying one or more sets of collinear points by taking a projection of the data points at an angle based on the estimated angle, and finding one or more peaks in the projection corresponding to respective sets of collinear points, and
  • D. outputting a characteristic of a line formed by at least one of the sets of collinear points.
  • 16. A method according to claim 15, wherein any of steps (A)-(C) include the step of filtering the projection to find one or more peaks therein.
  • 17. A method according to claim 16, wherein including the step of using a filter sized in accord with an estimated width of a peak resulting from taking a projection of the corresponding set of collinear points.
US Referenced Citations (208)
Number Name Date Kind
3816722 Sakoe et al. Jun 1974
3936800 Ejiri et al. Feb 1976
3967100 Shimomura Jun 1976
3968475 McMahon Jul 1976
3978326 Shimomura Aug 1976
4011403 Epstein et al. Mar 1977
4115702 Nopper Sep 1978
4115762 Akiyama et al. Sep 1978
4183013 Agrawala et al. Jan 1980
4200861 Hubach et al. Apr 1980
4254400 Yoda et al. Mar 1981
4286293 Jablonowski Aug 1981
4300164 Sacks Nov 1981
4385322 Hubach et al. May 1983
4435837 Abernathy Mar 1984
4441124 Heebner et al. Apr 1984
4441206 Kuniyoshi et al. Apr 1984
4519041 Fant et al. May 1985
4534813 Williamson et al. Aug 1985
4541116 Lougheed Sep 1985
4570180 Baier et al. Feb 1986
4577344 Warren et al. Mar 1986
4581762 Lapidus et al. Apr 1986
4606065 Beg et al. Aug 1986
4617619 Gehly Oct 1986
4630306 West et al. Dec 1986
4641349 Flom et al. Feb 1987
4688088 Hamazaki et al. Aug 1987
4706168 Weisner Nov 1987
4707647 Coldren et al. Nov 1987
4728195 Silver Mar 1988
4730260 Mori et al. Mar 1988
4731858 Grasmueller et al. Mar 1988
4736437 Sacks et al. Apr 1988
4742551 Deering May 1988
4752898 Koening Jun 1988
4758782 Kobayashi Jul 1988
4764870 Haskin Aug 1988
4771469 Wittenburg Sep 1988
4776027 Hisano et al. Oct 1988
4782238 Radl et al. Nov 1988
4783826 Koso Nov 1988
4783828 Sadjadi Nov 1988
4783829 Miyakawa et al. Nov 1988
4809077 Norita et al. Feb 1989
4831580 Yamada May 1989
4860374 Murakami et al. Aug 1989
4860375 McCubbrey et al. Aug 1989
4876457 Bose Oct 1989
4876728 Roth Oct 1989
4891767 Rzasa et al. Jan 1990
4903218 Longo et al. Feb 1990
4907169 Lovoi Mar 1990
4912559 Ariyoshi et al. Mar 1990
4914553 Hamada et al. Apr 1990
4922543 Aknhlbom et al. May 1990
4926492 Tanakan et al. May 1990
4932065 Feldgajer Jun 1990
4953224 Ichinose et al. Aug 1990
4955062 Terui Sep 1990
4959898 Landman et al. Oct 1990
4962423 Yamada et al. Oct 1990
4972359 Silver et al. Nov 1990
4982438 Usami et al. Jan 1991
5012402 Akiyama Apr 1991
5012524 LeBeau Apr 1991
5024419 Davis Jun 1991
5046190 Daniel et al. Sep 1991
5054096 Beizer Oct 1991
5060276 Morris et al. Oct 1991
5063608 Soege Nov 1991
5073958 Imme Dec 1991
5081656 Baker et al. Jan 1992
5081689 Meyer et al. Jan 1992
5086478 Kelly-Mahaffey et al. Feb 1992
5090576 Menten Feb 1992
5091861 Geller et al. Feb 1992
5091968 Higgins et al. Feb 1992
5093867 Hori et al. Mar 1992
5113565 Cipolla et al. May 1992
5115309 Hang May 1992
5119435 Berkin Jun 1992
5124622 Kawamura et al. Jun 1992
5133022 Weideman Jul 1992
5134575 Takagi Jul 1992
5143436 Baylor et al. Sep 1992
5145432 Midland et al. Sep 1992
5151951 Ueda et al. Sep 1992
5153925 Tanioka et al. Oct 1992
5155775 Brown Oct 1992
5159281 Hedstrom et al. Oct 1992
5159645 Kumagai Oct 1992
5164994 Bushroe Nov 1992
5168269 Harlan Dec 1992
5179419 Palmquist et al. Jan 1993
5185810 Freischlad Feb 1993
5185855 Kato et al. Feb 1993
5189712 Kajiwara et al. Feb 1993
5206820 Ammann et al. Apr 1993
5216503 Paik Jun 1993
5225940 Ishii et al. Jul 1993
5230027 Kikuchi Jul 1993
5243607 Masson et al. Sep 1993
5253306 Nishio Oct 1993
5253308 Johnson Oct 1993
5265173 Griffin et al. Nov 1993
5271068 Ueda et al. Dec 1993
5287449 Kojima Feb 1994
5297238 Wang et al. Mar 1994
5297256 Wolstenholme et al. Mar 1994
5299269 Gaborski et al. Mar 1994
5307419 Tsujino et al. Apr 1994
5311598 Bose et al. May 1994
5315388 Shen et al. May 1994
5319157 Nakahashi et al . Jun 1994
5327156 Masukane et al. Jul 1994
5329469 Watanabe Jul 1994
5337262 Luthi et al. Aug 1994
5337267 Colavin Aug 1994
5363507 Nakayama et al. Nov 1994
5367439 Mayer et al. Nov 1994
5367667 Wahlquist et al. Nov 1994
5371690 Engel et al. Dec 1994
5388197 Rayner Feb 1995
5388252 Dreste et al. Feb 1995
5398292 Aoyama Mar 1995
5418892 Aghajan et al. May 1995
5432525 Maruo et al. Jul 1995
5440699 Farrand et al. Aug 1995
5455870 Sepai et al. Oct 1995
5455933 Schieve et al. Oct 1995
5471312 Watanabe et al. Nov 1995
5475766 Tsuchiya et al. Dec 1995
5477138 Efjavic et al. Dec 1995
5481712 Silver et al. Jan 1996
5485570 Bushboom et al. Jan 1996
5491780 Fyles et al. Feb 1996
5495424 Tokura Feb 1996
5495537 Bedrosian et al. Feb 1996
5500906 Picare et al. Mar 1996
5506617 Parulski et al. Apr 1996
5506682 Pryor Apr 1996
5511015 Flockencier Apr 1996
5519840 Maties et al. May 1996
5526050 King et al. Jun 1996
5528703 Lee Jun 1996
5532739 Garakani et al. Jul 1996
5539409 Mathews et al. Jul 1996
5548326 Michael Aug 1996
5550763 Michael Aug 1996
5553859 Kelly et al. Sep 1996
5557410 Huber et al. Sep 1996
5557690 O'Gorman et al. Sep 1996
5566877 McCormack Oct 1996
5568563 Tanaka et al. Oct 1996
5574668 Beaty Nov 1996
5574801 Collet-Beillon Nov 1996
5581632 Koljonen et al. Dec 1996
5583949 Smith et al. Dec 1996
5583954 Garakani Dec 1996
5586058 Aloni et al. Dec 1996
5592562 Rooks Jan 1997
5594859 Palmer et al. Jan 1997
5602937 Bedrosian et al. Feb 1997
5608490 Ogawa Mar 1997
5608872 Schwartz et al. Mar 1997
5640199 Garakani et al. Jun 1997
5640200 Michael Jun 1997
5642158 Petry, III et al. Jun 1997
5647009 Aoki et al. Jul 1997
5657403 Wolff et al. Aug 1997
5672334 Nichani et al. Sep 1997
5676302 Petry, III Oct 1997
5696848 Patti et al. Dec 1997
5715369 Spoltman et al. Feb 1998
5717785 Silver Feb 1998
5724439 Mizuoka et al. Mar 1998
5742037 Scola et al. Apr 1998
5751853 Michael May 1998
5754679 Koljonen et al. May 1998
5757956 Koljonen et al. May 1998
5761326 Brady et al. Jun 1998
5768443 Michael e tal. Jun 1998
5793899 Wolff et al. Aug 1998
5796386 Lioscomb et al. Aug 1998
5796868 Dutta-Choudhury Aug 1998
5801966 Ohashi Sep 1998
5805722 Cullen et al. Sep 1998
5809658 Jackson et al. Sep 1998
5818443 Schott Oct 1998
5825483 Michael et al. Oct 1998
5825913 Rostami et al. Oct 1998
5835622 Koljonen et al. Nov 1998
5845007 Ohashi et al. Dec 1998
5848189 Pearson et al. Dec 1998
5850466 Schott Dec 1998
5859923 Petry, III et al. Jan 1999
5861909 Garakani et al. Jan 1999
5872870 Michael Feb 1999
5878152 Sussman Mar 1999
5900975 Sussman May 1999
5901241 Koljonen et al. May 1999
5909504 Whitman Jun 1999
5912768 Sissom et al. Jun 1999
5912984 Michael et al. Jun 1999
5918196 Jacobson Jun 1999
5933523 Drisko et al. Aug 1999
5943441 Michael Aug 1999
Foreign Referenced Citations (11)
Number Date Country
0527 632 A2 Feb 1993 EPX
0 777 381 A2 Nov 1996 EPX
WO 9522137 Aug 1995 WOX
WO 9521376 Aug 1995 WOX
WO 9722858 Jun 1997 WOX
WO 9721189 Jun 1997 WOX
WO 9724693 Jul 1997 WOX
WO 9724692 Jul 1997 WOX
WO 9852349 Nov 1998 WOX
WO 9915864 Apr 1999 WOX
WO 9859490 Dec 1999 WOX
Non-Patent Literature Citations (20)
Entry
Grimson, W. Eric L. and Huttenlocher, Daniel P., "On the Sensitivity of the Hough Transform for Object Recognition", May 1990, IEEE Transactions on Pattern Analysis and machine Intelligence, vol. 12, No. 3.
Medina-Mora et al. (1981) An Incremental Programming Environment, IEEE Transactions on Software Eng. SE-7:472-482.
Teitelbaum et al. (19810 The Cornell Program Synthesizer: A Syntax-Directed Programming Environment, Communications of the ACM 24:563-573.
Newletter from Acquity Imaging, Inc., "Remote Vision Support Package--The Phones Are Ringing!," 1 page.
PictureTel Corporation Product Brochure "Picturetel Live PCS 100(tm) Personal Visual Communications System," 3 pp. (1993).
PictureTel Corporation Procuct Brochure "Picturetel System 1000: Completer Videoconferencing for cost Sensitive Applications," 4 pp. (1993).
PictureTel Corporation Product Brochure, "Picturetel System 4000(tm) A Family of Models to Fit Your Application From Offices to Boardrooms, Classrooms, and Auditoriums," 4 pp. (1993).
Symantec Corporation, "The Norton pcAnywhere User's Guide," Table of Contents 8 pp; Introduction of pcAnywhere Technology pp i-vii; Chapter 7--Sessions; pp. 191-240 (1991).
Bursky,Dave, "CMOS Four-Chip Set Process Images at 20-MHz Data Rates," Electronic Design May 28, 1987, pp. 39-44.
Plessey Semiconductors, Preliminary Information, May 1986, Publication No. PS2067, May 1986, pp. 1-5.
NEC Electronics Inc., PD7281 Image Pipelined Processor, Product Information Brochure, pp. 2-169-2-211.
Horn, Berthold Klaus Paul. "Robot Vision", The Massachusetts Institute for Technology, 1986.
Racca Roberto G., Stephenson Owen, and Clements Reginald, M. High-speed video analysis system using multiple shuttered charge-coupled device imagers and digital storage. Optical Engineering (Jun. 1992) 31:6.
Gevorkian David Z., Astola Jaakko T., and Atourian Samvel M. "Improving gil-Werman Algorithm for Running lMin and Max Filters" IEEE Transactions on Pattern Analysis and Machie Intelligence, vol. 19, No. 5, May 1997, pp. 526-529.
Gil, Joseph and Werman Michael. "Computing 2-D Min, Median, and Max Filters" IEEE transactions on Pattern Analysis and Machine Intelligence, vol. 15, No. 5, may 1993, pp. 504-507.
Chapetr 3: "Guidelines for Developing MMX Code," Intel.
Chapter 4: "MMX Code Development Strategy," Intel.
Chapter 5: "MMX Coding Techniques," Intel.
Chapter 3: "Optimization Techniques for Integer Blended Code," Intel.
Rosenfeld, Azriel. "Computer Vision: Basic Principles," Proceedings of the IEEE. vol. 76, No. 8, Aug. 1988, pp. 863-868.