A template is analyzed to determine a step size for searching within a search area. first, the template is padded with pixels to increase its size. Cross-correlation between the padded template and the original template leads to identification of an effective step size. step sizes for each of horizontal and vertical axes are derived. Windows of the search area, selected based on the step size, then are tested in a fast search stage by correlating the template to the window. Any tested window which has a correlation coefficient exceeding a specific value is a local match. A full search of the vicinity of the local match then is performed for all potential windows within an area bounded by one step to either side of the local match along either axis. The location(s) corresponding to the highest correlation(s) exceeding the threshold value are matches.
|
4. A method for locating a match between a first template of data and a search area of data, the first template including a first plurality of data points, the search area including a second plurality of data points exceeding in number the first plurality of data points, the method comprising the steps of:
padding the first template with additional data to enlarge the template into a padded template, wherein a first window area within the padded template is formed by the first template and includes the first plurality of data points; correlating the first template to window areas within the padded template along a first axis to derive a first step size; stepping through the search area using the first step size for selecting window areas within the search area; and evaluating the selected window areas for locating a match between the first template of data and the search area of the data.
3. An apparatus for estimating a location of a first template of data within a search area of data, the first template including a first plurality of data points, the search area including a second plurality of data points exceeding in number the first plurality of data points, the apparatus comprising:
a processor which pads the first template with additional data to enlarge the template into a padded template, wherein a first window area within the padded template is formed by the first template and includes the first plurality of data points; a correlator which correlates the first template to window areas within the padded template to derive a first step size; a selector which steps though the search area using the first step size to select window areas within the search area; and an evaluator which evaluates the selected window areas for locating a match between the first template of data and the search area of the data.
1. A method for determining a fixed step size which is used for stepping through a search area to select a subset of window areas within a search area to be evaluated for locating a match between a subject template of data and the search area of data, the subject template including a first plurality of data points, the search area including a second plurality of data points exceeding in number the first plurality of data points, the method comprising:
creating a created template from the subject template without use of the search area data; and correlating the subject template to window areas within the created template to derive the fixed step size for stepping through the search area; wherein said creating comprises padding the subject template with additional data to enlarge the subject template into the created template, wherein a first window area within the created template is formed by the subject template and includes the first plurality of data points.
2. A device for determining a fixed step size which is used for stepping through a search area to select a subset of window areas within a search area to be evaluated for locating a match between a subject template of data and the search area of data, the subject template including a first plurality of data points, the search area including a second plurality of data points exceeding in number the first plurality of data points, the device comprising:
means for creating a created template from the subject template without use of the search area data; and means for correlating the subject template to window areas within the created template to derive the fixed step size for stepping through the search area; wherein the creating means comprises means for padding the subject template with additional data to enlarge the subject template into the created template, wherein a first window area within the created template is formed by the subject template and includes the first plurality of data points.
5. The method of
correlating the first template to window areas within the padded template along a first axis about the first window area to derive a second step size.
6. The method of
stepping through the search area using either one or both of the first step size and the second step size for selecting window areas within the search area; and correlating each one of the window areas among the selected subset of window areas to the first template.
7. The method of
correlating each one of the window areas among the selected subset of window areas to the first template.
8. The method of
identifying as a potential template match each of the correlated window areas for which the correlation to the first template results in a correlation coefficient exceeding a predetermined value.
9. The method of
respectively correlating the first template to one window area of a subset of potential window areas within the search area for each one window area within the subset of potential window areas, wherein the subset of potential window areas is selected from all potential window areas within the search area using either one or both of the first step size and the second step size, and wherein each respective correlation performed in the step of respectively correlating the first template to said one window area of the subset results in a correlation coefficient; for each of the respective correlations between the first template and said one window area within the subset of potential window areas, comparing the resulting correlation coefficient to a predetermined value.
10. The method of
identifying as a potential template match each said one window area for which the corresponding correlation results in a correlation coefficient which exceeds a predetermined value.
11. The method of
12. The method of
deriving a padding value from the first plurality of data points; using the padding value as a value for each data point added to the first template to form the padded template.
13. The method of
14. The method of
padding the first template with additional pixels to enlarge the template into a padded template.
15. The method of
selecting a second subset of window areas in the search area in the vicinity of a given potential template match; correlating each one of the window areas among the selected second subset of window areas to the first template.
16. The method of
identifying as a template match one or more of the correlated window areas among the second subset and the window area corresponding to the local template match for which the correlation to the first template results in a correlation coefficient which exceeds a threshold value.
17. The method of
18. The method of
19. The apparatus of
means identifying as a potential template match each of the correlated window areas for which the correlation to the first template results in a correlation coefficient exceeding a predetermined value; means for selecting a subset of window areas in the search area in the vicinity of a given potential template match; and means for correlating each one of the window areas among the selected subset of window areas to the first template.
20. The apparatus of
means for identifying as a template match one or more of the correlated window areas among the subset and the window area corresponding to the local template match for which the correlation to the first template results in a correlation coefficient which exceeds a threshold value.
|
This application is a continuation of U.S. patent application Ser. No. 09/216,692 filed on Dec. 18, 1998 (now U.S. Pat. No. 6,301,387).
This invention relates to template matching within a data domain, and more particularly to a method for locating a given data template within a data domain.
Template matching in the context of an image search is a process of locating the position of a subimage within an image of the same, or more typically, a larger size. The subimage is referred to as the template and the larger image is referred to as the search area. The template matching process involves shifting the template over the search area and computing a similarity between the template and the window of the search area over which the template lies. Another step involves determining a single or a set of matched positions in which there is a good similarity measure between the template and the search area window.
A common technique for measuring similarity in template matching and image registration is cross-correlation. A correlation measure is determined between the template and respective windows of the search area to find the template position which has maximum correlation. For a two-dimensional search area the correlation function generally is computed for all translations of the template within the search area. A statistical correlation measure is a common approach in which window areas are spatially convolved with the template using spatial filter functions. Because this approach is extremely expensive in terms of computation time, a more common computer implementation is to use a sum of absolute differences.
Rosenfeld et al., in "Coarse-Fine Template Matching," IEEE Transactions on Systems, Man and Cybernetics (February 1977, pp. 104-107) describe an approach where a `reduced-resolution` template is used during a first, coarse evaluation stage. The template is divided into blocks of equal size (e.g., `m` pixels per block). The average of each block is computed. For each pixel of the search area an average also is calculated over a neighborhood of the same size as the reduced-resolution template (e.g., m pixels). The average absolute difference between each template block average and the picture neighborhood average then is computed for each pixel of the search area. If the average absolute difference for any pixel of the search area is below a threshold value, then a possible match has been identified. Next, the full resolution template is compared to a window of the search area about each pixel point where the average absolute difference in the prior coarse evaluation step was below the threshold value. This fine evaluation step identifies if there actually is a good correlation.
Goshtasby et al. in "A Two-Stage Correlation Approach to Template Matching," IEEE Transaction on Pattern Analysis and Machine Intelligence, (Vol. PAMI-6, No. 3, May 1984), note the need for an accurate threshold value for the first stage evaluation. They describe a method for deriving the threshold value based upon sub-template size and false dismissal probability.
The coarse-fine or two stage method subsample the template to match with the image. The task of subsampling the template is not a trivial task and contributes significant processing cost. In addition, the false alarms result in wasted, or an ineffective use of, processing time. Accordingly, there is a need for a more efficient method of template matching.
In the area of motion estimation for digital video and multimedia communications a three stage correlation strategy is used. In a first step, a search step size of 4 is used. Once a maximum point is found, the step size is reduced to 2 to evaluate the neighborhood of the previously determined point to choose the next search point. The third step is to search all neighboring points to find the best match. This approach speeds up the search process, but also has a high probability of mismatches or suboptimal matches. It also has difficulty handling cases in which multiple match points occur. Thus, there is a need for a more reliable, fast search method for correlating a template to windows of a search area.
According to the invention, a correlation auto-predictive search method is used to compare a template to windows of a search area. The location(s) where the template has the highest correlation coefficient with the underlying window is selected as a match for the template. Local maximum criteria or other criteria then are used to select one or more match points within the search area. The principle of a correlation auto-predictive search as conceived by the inventors is (1) to extract statistical information from the template itself to determine the search step size, and (2) to perform fast searching based on this extracted information.
According to one aspect of the invention, during a first analytical step, autocorrelation is performed on the template to generate desired statistics. To use autocorrelation, the original template is padded with additional pixels to increase the template size. In one approach, where the search area is assumed to be periodic, circular padding is used. In such approach the padded template is an array of copies of the original template. This increases the template size to the search area (image) size.
In another approach linear padding is used in which pixels are added around the original template to increase the size of the template to the search area (image) size. According to an aspect of this invention, a mean pixel value of the original template is used as a padding constant (i.e., pixel value for the added pixels). Alternatively, a value of zero or another fixed value may be used as the padding constant for the padded pixels.
After generating the padded template, cross-correlation is performed between the padded template and the original template. The autocorrelation is highest at the center of the padded template as this area is formed by the original template. This corresponds to a peak in a graph of the autocorrelation of the padded template to original template. The width of the peak, either along a horizontal direction of the padded template, or along a vertical direction of the padded template, may be measured. The height of the maximum peak is 1∅ The horizontal width is taken as the distance along the horizontal axis between autocorrelation values of 0.5 to each side of the maximum peak. Similarly, the vertical width is taken as the distance along the vertical axis between autocorrelation values of 0.5 to each side of the maximum peak. Such value, 0.5, is referred to herein as the cut value. The cut value may differ.
According to another aspect of the invention, the autocorrelation between the padded template and the original template is not calculated for every point of the padded template. At the center of the padded template, the correlation is known to be 1.0 because the original template is located at such center of the padded template. The correlation then is derived about the center of the padded template along both horizontal and vertical axes. As the correlations are derived during this stepping along the axes, there comes a point where the correlation decreases to the cut value. Along the horizontal axis, there is a cut value reached to either direction of center. The horizontal distance between these two locations where the correlation has decreased to the cut value is the horizontal width. Further correlations along such axis need not be derived. Along the vertical axis, there also is a cut value reached to either direction of center. The vertical distance between these two locations where the correlation has decreased to the cut value is the vertical width. Further correlations along such vertical axis need not be derived. Thus, correlation coefficients are derived only for the steps along the axes away from center, and only to the step where the cut value is reached.
Next, horizontal step size and vertical step size are derived from the horizontal width and vertical width, respectively. In one embodiment the horizontal step size is 0.5 times the horizontal width. Similarly, a vertical step size is 0.5 times the vertical width. These step sizes are the correlative auto-predictive search (CAPS) step sizes. No additional correlation values need be derived between the padded template and the original template. The CAPS step sizes, then are used for template matching between the original template and the search area.
According to another aspect of this invention, a fast search then is performed between the template and the search area using the derived step sizes. Then, for correlations having a correlation coefficient exceeding a specific value, a full search is performed locally in each area where the fast search resulted a correlation coefficient exceeding the select value.
According to another aspect of the invention, the fast search is performed as a set of correlation between the original template and the search area. Specifically, a correlation is performed between the template and a window area within the search area. The set of correlations is selected by choosing window areas based upon the step size. For example, one window is the center of the search area. A positive or negative step then is taken along an axis using the corresponding horizontal width or vertical width to derive a correlation for another window. Any of the correlations which result in a correlation coefficient exceeding a specific value is considered a local match point.
According to another aspect of this invention, the specific value used to identify a local match during the fast search is the cut value times a threshold value. The cut value is the same cut value used during the first analytical step described above, to derive statistics from the template. The threshold value is assigned based upon image characteristics. Typical threshold values are between 0.8 and 0.9.
One or more locations are identified as local match points based upon the whether the correlation coefficient between the template and that location exceed the specific value (e.g., cut value times threshold value).
According to another aspect of the invention, a full search then is performed in the vicinity of any location which is a local match. A full search of such vicinity encompasses performing a correlation between the template and every potential search area window between the local match location window and the windows at the prior and next step in each of the horizontal and vertical axes. For example, if the horizontal step size is 3 pixels and the vertical step size is 4 pixels, then correlations are performed for windows ±1 pixel and ±2 pixels along the horizontal axis and ±1 pixel, ±2 pixels and ±3 pixels along the vertical axis. In addition correlations are performed for windows off the axes within the area delineated by the step sizes. Thus, the full search of the vicinity of the local match for this example includes 34 correlations between the template and the search area. Any locations among the local match locations and the locations tested during the full search of the vicinity which exceed the threshold value are considered template matches. In some embodiments, the only the location having the highest correlation is considered a match. In other embodiments there may be multiple matches. Thus, the top matches or all matches above the threshold are selected as resultant matches.
One advantage of the invention is that template matches are found more quickly and with greater reliability than prior correlation search methods. In particular, this search methodology is more tolerant of noise and offsets of the template as demonstrated empirically by forming a search area from copies of templates altered by low pass filtering or Gaussian noise. These and other aspects and advantages of the invention will be better understood by reference to the following detailed description taken in conjunction with the accompanying drawings.
Exemplary Host Computer System
The functions of the present invention preferably are performed by programmed digital computers of the type which are well known in the art, an example of which is shown in
Overview
Referring to
To reduce the number of windows 42 that the template 38 is compared with, an effective step size is derived from the template. According to a 2-dimensional implementation embodiment, a step size along a first axis 44 is derived and a step size along a second axis 46 is derived. Rather then compare the template to every possible window of the search area 40, the template 38 is moved along either or both of the first axis 44 and second axis 46 by the corresponding first axis step size or second axis step size.
Once the desired step sizes are derived, then the template 38 is compared to the various windows 42 of the search area 40 at the step size increments during a fast search process. In one embodiment the comparison is a correlation function of the template 38 and the window 42 and results in a correlation coefficient. Any window 42 in which the correlation coefficient with the template 38 is found to exceed a specific value is a local match for the template. In a preferred embodiment the specific value is the cut value times a threshold value.
Next, a full search then is performed in the vicinity of any location which is a local match. A full search of such vicinity encompasses performing a correlation between the template and every potential search area window between the local match location window and the windows at the prior and next step in each of the horizontal and vertical axes. For example, if the horizontal step size is 3 pixels and the vertical step size is 4 pixels, then correlations are performed for windows ±1 pixel and ±2 pixels along the horizontal axis and ±1 pixel, ±2 pixels and ±3 pixels along the vertical axis. In addition correlations are performed for windows off the axes within the area delineated by the step sizes. Thus, the full search of the vicinity of the local match for this example includes (2*2+1) * (2*3+1)-1=34 correlations between the template and the search area. Any locations among the local match locations and the locations tested during the full search of the vicinity which exceed the threshold value are considered template matches. In some embodiments, the only the location having the highest correlation is considered a match. In other embodiments there may be multiple matches. Thus, the top matches or all matches above the threshold are selected as resultant matches.
Determining Step Size
To determine effective step sizes, the template 38 itself is analyzed. Referring to
Referring again to
Referring to
Rather than perform a correlation for each potential window along the first axis 44, correlations are performed for windows along the axis 44 away from the center window in each direction 56, 58 until a window is identified in such direction where the correlation coefficient intersects the cut-off value. For two dimensional analysis, there is a cut-off point found in each direction from the center window 54c along the first axis 44. The distance between those two windows in data points is the width along the first axis.
Referring to
In steps 66 and 68 the correlations are repeated along the second axis 46 in directions 70, 72 to find a width along the second axis 46. Referring again to
Fast Search
Once the CAPS step sizes have been derived, a fast search is performed comparing the template 38 to the search area 40. It is a fast search in the sense that not every potential window of the search area is compared to the template. Referring to
Referring to
At step 80 the shifting along the first axis 44 and testing of windows is performed for a template center point repositioned over every y-th row of data points. Specifically, once the initial row of the search area has been tested, the template 38 is moved along the second axis 46 to another row that is y data points away, where y is the second axis step size. This next row then is tested by shifting along the first axis 44 using the first axis step size. A correlation is performed at each iteration. Then another row is tested which is y data points away along the second axis 46. In this manner the template is shifted by the second step size along the second axis 46 and by the first step size along the first axis 44 to select windows to be tested during the fast search. For example, in a search area which is 400 pixels by 400 pixels, and where the first axis step size is four and the second axis step size is four, there are 100 * 100=10,000 windows tested during the fast search.
Of the tested windows, at step 82 the window location for any correlation which resulted in a correlation coefficient which is greater than or equal to the product of the cut value times a predetermined threshold value is considered a local match. In a preferred embodiment the cut value is the same for each axis. Where the cut value used along one axis differs from the cut value used along the other axis, either cut value may be used. Alternatively, an average of the cut values may be used. The threshold value is a predetermined value and signifies the minimum correlation coefficient acceptable to designate a window as being a match for the template. Typical values are 0.8 and 0.9. The specific value may vary based upon the search area or type of date. The specific value may be determined empirically for different types of data or search area characteristics.
Local Full Search
Once the fast search is complete (or during the course of the fast search), a local full search is performed about each of the local matches. Referring to
Referring to
A correlation is performed between the template 38 and each window in the vicinity of the local match. For the vicinity shown in
Implementing the Correlation Function
The correlation coefficient for a correlation between two data sets `a` and `b` is defined below. The data set `a` is the template 38. The data set `b` is a window of the padded template 52 (or of a rotational offset of the padded template) for the process of finding the CAPS step sizes. The data set `b` is a window of the search area 40 (or of a rotational offset of the search area) for the process of identifying candidate locations, potential template matches or template matches. Each of data sets `a` and `b` may be a matrix, image or another set of data points. The correlation coefficient, corr is:
which may be simplified to
where E(x)=expected value of data set (x)
sd(x)=standard deviation of data set (x)
and corr is between -1.0 and +1∅
Meritorious and Advantageous Effects
One advantage of the invention is that template matches are found more quickly and with greater reliability than prior correlation search methods. In particular, this search methodology is more tolerant of noise and offsets of the template as demonstrated empirically by forming a search area from copies of templates altered by lowpass filtering or Gaussian noise.
Although a preferred embodiment of the invention has been illustrated and described, various alternatives, modifications and equivalents may be used. For example, although the fast search is described as performing respective correlations between the template 38 and a subset of window areas, an estimate of the template 38 may be used instead followed by additional correlation using the full template for correlations resulting in a coefficient exceeding a prescribed value.
Further, although the full local search is described to include performing a correlation for each potential window in the vicinity of the local match (e.g., step size equal 1), an intermediate level search may be performed instead using a step size less than the first step size and second step size, but greater than 1. The step size for the intermediate level search may be determined in the same manner as the for the fast search--based on template characteristics. For example, two cut values are used during the process of identifying step sizes. One cut value for a given axis is used to determine the step size for the fast search. The other cut value for the same axis is used to determine the step size for the intermediate level search. The step size for the fast search is to be larger than the step size for the intermediary level search.
In operation, the intermediate level search is performed in the vicinity of all local matches resulting from the fast search, where the vicinity is bounded by the first step size and second step size of the fast search. Any of the windows tested in the intermediate level search which has a correlation to the first template exceeding a prescribed value is an intermediate level search local match. The prescribed value is selected either independent of or in relation to the cut values and threshold value previously described. Next, a full search is performed in the vicinity of the intermediate level search local matches, where the vicinity is bounded by the step sizes used in the intermediate level search. Any one or more of the local matches and full search correlations which have a correlation to the first template exceeding the threshold value are template matches.
As used herein the term vicinity of a window refers to an area bounded by lines one step size away or less and parallel to the search area axes.
Therefore, the foregoing description should not be taken as limiting the scope of the inventions which are defined by the appended claims.
Kim, Yongmin, Sun, Shijun, Park, HyunWook
Patent | Priority | Assignee | Title |
7016535, | Jul 19 2001 | Fujitsu Limited | Pattern identification apparatus, pattern identification method, and pattern identification program |
7054491, | Nov 16 2001 | STMicroelectronics, Inc.; STMicroelectronics, Inc | Scalable architecture for corresponding multiple video streams at frame rate |
7116824, | Oct 23 2001 | STMICROELECTRONICS FRANCE | Detection of a pattern in a digital image |
7206778, | Dec 17 2001 | AVOLIN, LLC | Text search ordered along one or more dimensions |
8228990, | Jan 16 2008 | Sony Corporation; Sony Electronics Inc. | Template matching scheme using multiple predictors as candidates for intra-prediction |
8238691, | Sep 09 2005 | GRASS VALLEY LIMITED | Method of and apparatus for image analysis |
8485559, | Jul 24 2006 | THALES DIS FRANCE SAS | Document authentication using template matching with fast masked normalized cross-correlation |
Patent | Priority | Assignee | Title |
5359513, | Nov 25 1992 | Arch Development Corporation | Method and system for detection of interval change in temporally sequential chest images |
5495537, | Jun 01 1994 | Cognex Corporation | Methods and apparatus for machine vision template matching of images predominantly having generally diagonal and elongate features |
5784108, | Dec 03 1996 | ZAPEX TECHNOLOGIES INC | Apparatus for and method of reducing the memory bandwidth requirements of a systolic array |
5911001, | Oct 25 1994 | Fuji Machine Mfg. Co., Ltd. | Image processing apparatus |
5943442, | Jun 12 1996 | Nippon Telegraph and Telephone Corporation | Method of image processing using parametric template matching |
6014181, | Oct 13 1997 | Sharp Laboratories of America, Inc | Adaptive step-size motion estimation based on statistical sum of absolute differences |
6075557, | Apr 17 1997 | Sharp Kabushiki Kaisha | Image tracking system and method and observer tracking autostereoscopic display |
6122320, | Mar 14 1997 | Cselt-Centro Studi e Laboratori Telecomunicazioni S.p.A. | Circuit for motion estimation in digitized video sequence encoders |
Executed on | Assignor | Assignee | Conveyance | Frame | Reel | Doc |
Aug 23 2001 | University of Washington | (assignment on the face of the patent) | / |
Date | Maintenance Fee Events |
Nov 29 2006 | M1551: Payment of Maintenance Fee, 4th Year, Large Entity. |
Dec 16 2010 | M1552: Payment of Maintenance Fee, 8th Year, Large Entity. |
Nov 24 2014 | M1553: Payment of Maintenance Fee, 12th Year, Large Entity. |
Date | Maintenance Schedule |
Jun 24 2006 | 4 years fee payment window open |
Dec 24 2006 | 6 months grace period start (w surcharge) |
Jun 24 2007 | patent expiry (for year 4) |
Jun 24 2009 | 2 years to revive unintentionally abandoned end. (for year 4) |
Jun 24 2010 | 8 years fee payment window open |
Dec 24 2010 | 6 months grace period start (w surcharge) |
Jun 24 2011 | patent expiry (for year 8) |
Jun 24 2013 | 2 years to revive unintentionally abandoned end. (for year 8) |
Jun 24 2014 | 12 years fee payment window open |
Dec 24 2014 | 6 months grace period start (w surcharge) |
Jun 24 2015 | patent expiry (for year 12) |
Jun 24 2017 | 2 years to revive unintentionally abandoned end. (for year 12) |