US 6,477,272 B1Grant

Object recognition with co-occurrence histograms and false alarm probability analysis for choosing optimal object recognition process parameters

Issue Date:2002-11-05

•27 Claims

•8 Drawing Sheets

Abstract

This invention is directed toward an object recognition system and process that identifies the location of a modeled object in a search image. This involves first capturing model images of the object whose location is to be identified in the search image. A co-occurrence histogram (CH) is then computed for each model images. A model image CH is computed by generating counts of every pair of pixels whose pixels exhibit colors that fall within the same combination of a series of pixel color ranges and which are separated by a distance falling within the same one of a series of distance ranges. Next, a series of search windows, of a prescribed size, are generated from overlapping portions of the search image. A CH is also computed for each of these search windows using the pixel color and distance ranges established for the model image CHs. A comparison between each model image CH and each search window CH is conducted to assess their similarity. A search window that is associated with a search window CH having a degree of similarity to one of the model image CHs which exceeds a prescribed search threshold is designated as potentially containing the object being sought. This designation can be presumed final, or further refined. This system and process requires that the size of the search window, color ranges and distance ranges be chosen ahead of time. The choice of these parameters can be optimized via a false alarm analysis.

Metadata

Assignee

Microsoft Corporation

Inventors

John Krumm
Peng Chang

Application Information

Application Number:US 09/336,274

Filing Date:1999-06-18

Priority Date:1999-06-18

Art Unit:7

Classifications

IPC:

G06K 900

Field of Search:

382165382168382170382194382217382218382224382305

Patent Drawings (8 sheets)

Description

Background

[0002] 1. Technical Field

[0003] The invention is related to a computer-implemented object recognition system and process, and more particularly, to such a system and process employing co-occurrence histograms (CH) for finding an object in a search image.

[0004] 2. Background Art

[0005] Object recognition in images is typically based on a model of the object at some level of abstraction. This model is matched to an input image which has been abstracted to the same level as the model. At the lowest level of abstraction (no abstraction at all), an object can be modeled as a whole image and compared, pixel by pixel, against a raw input image. However, more often unimportant details are abstracted away, such as by using sub-templates (ignoring background and image position), normalized correlation (ignoring illumination brightness), or edge features (ignoring low spatial frequencies). The abstraction itself is embodied in both the representation of the object and in the way it is matched to the abstracted image. For instance, Huttenlocher et al. [1] represent objects as simple edge points and then match with the Hausdorff distance. While the edge points form a completely rigid representation, the matching allows the points to move nonrigidly.

[0006] One interesting dimension of the aforementioned abstraction is rigidity. Near one end of this dimension are the several object recognition algorithms that abstract objects into a rigid or semi-rigid geometric juxtaposition of image features. These include Hausdorff distance [1], geometric hashing [2], active blobs [3], and eigenimages [4, 5]. In contrast, some histogram-based approaches abstract away (nearly) all geometric relationships between pixels. In pure histogram matching, e.g. Swain & Ballard [6], there is no preservation of geometry, just an accounting of the number of pixels of given colors. The technique of Funt & Finlayson [7] uses a histogram of the ratios of neighboring pixels, which introduces a slight amount of geometry into the representation.

[0007] Abstracting away rigidity is attractive, because it allows the algorithm to work on non-rigid objects and because it reduces the number of model images necessary to account for appearance changes due to scaling and viewpoint change. One can start with a geometrically rigid approach and abstract away some rigidity by using geometric invariants [8], loosening the matching criteria [1], or explicitly introducing flexibility into the model [3]. On the other hand, one can start with a method like Swain & Ballard's color indexing [6], which ignores all geometry, and add some geometric constraints. For example, some histogram-based approaches, most of which are used to find images in a database rather than to find an object in an image, have employed attempts to add spatial information to a regular color histogram. Included among these are Huang et al. [9] where a “color correlogram” is used to search a database for similar images, or Pass and Zabih [10] where “color coherence vectors” are employed that represent which image colors are part of relatively large regions of similar color.

[0008] The dilemma comes in deciding how much to abstract away. The goal is to ignore just enough details of the object's appearance to match all anticipated images of the object, but not so many details that the algorithm generates false matches. The present invention addresses this issue.

[0009] It is noted that in the preceding paragraphs, as well as in the remainder of this specification, the description refers to various individual publications identified by a numeric designator contained within a pair of brackets. For example, such a reference may be identified by reciting, “reference [1]” or simply “[1]”. Multiple references will be identified by a pair of brackets containing more than one designator, for example, [4, 5]. A listing of the publications corresponding to each designator can be found at the end of the Detailed Description section.

Summary

[0010] This invention is directed toward an object recognition system. and process that identifies the location of a modeled object in an image. Essentially, this involves first creating model images of the object whose location is to be identified in the search image. As the object may be depicted in the search image in any orientation, it is preferred that a number of model images be captured, of which each shows the object from a different viewpoint. Ideally, these model images would be taken from viewpoints spaced at roughly equal angles from each other around the object. In addition, multiple images could be taken at each angular viewpoint where each is captured at a different distance away from the object being modeled. This latter method would better model an object whose distance from the camera capturing the search image is unknown. One way of determining how far apart each angular viewpoint should be to ensure a good degree of match between the object in the search image and one of the model images is to require a high degree of matching between each model image. Similarly, it is desirable that a good degree of match exist between adjacent model images taken at the same angular viewpoint, but at different distances from the object, for the same reasons.

[0011] Each model image is processed to compute a color co-occurrence histogram (CH). This unique histogram keeps track of the number of pairs of certain “colored” pixels that occur at particular separation distances in the model image. In this way, geometric information is added to the histogram to make it more selective of the object in the search image, than would a simple color histogram. In generating the CH there are two parameters that are selected ahead of time. These parameters are the color ranges and the distance ranges. Choosing the color ranges involves dividing the total possible color values associated with the pixels of the model images into a series of preferably equal sized color ranges. Choosing the distance rangesiinvolves selecting a set of distance ranges, for example (1-2 pixels), (2-3 pixels), (3-4 pixels) . . . , up to a prescribed maximum separation distance between pairs of pixels that will be checked. The size of each color and distance range (& so ultimately the total number of different ranges) is preferably selected :to optimize the search process as will be discussed later. Essentially, these ranges should be made small enough to ensure a high selectivity in finding the object in the search image and to minimize false matches. On the other hand, the ranges should be large enough to reduce the amount of processing and to allow enough flexibility in the match that small changes in the shape of the object in the search image (e.g. due to flexure, a slight change in viewpoint, or a different zoom level or scale) do not prevent an object from being recognized.

[0012] Once the color and distance ranges are established, each pixel in the model image is quantized in regards to its color by associating it to the color range which includes the actual color value of the pixel. The model images are quantized in regards to distance by associating each possible unique, non-ordered, pair of pixels in a model image to one of the distance ranges based on the actual distance separating the pixel pair. The CH is then generated by establishing a count of the number of pixel pairs in: a model image which exhibit the same “mix” of color ranges and the same distance range.

[0013] The search image (i.e., the image to be searched for the modeled object) is first cordoned into a series of preferably equal sized sub-images or search windows. These search windows preferably overlap both side-to-side and up-and-down. In the tested embodiment of the present invention, the overlap was set to one-half the width and height of the search window. The size of the search window is preferably as large as possible so as to minimize the search process. However, it is also desired that the search window not be so large that false matches or alarms occur. It should be noted that the size of the search window, as well as the size of the aforementioned color and distance ranges, can be optimized via a unique false alarm analysis. This analysis will be discussed in greater detail below.

[0014] A CH is calculated for each search window in the same manner that they were generated in connection with the model images. In particular, the same color & distance ranges are employed in the quantization steps of the process. Each of these search window CHs is compared to, each of the model image CHs. Essentially, this comparison assesses the similarity between the histograms. The results of this comparison can be handled in two ways. In a first version, a similarity value threshold is established and any comparison that exceeds the threshold would cause the associated search window to be declared as potentially containing the object. This version would allow identifying each location of the object in a search image which depicts more than one of the objects. In a second version, the threshold would still be used, but only the greatest of the similarity values exceeding the threshold would be chosen to indicate a search window potentially contains the object. This version would be especially useful where it is known only one of the objects exists in the image. It could also be useful in that the threshold could be made lower and the object still located, for example in a noisy image.

[0015] It is noted the preferred approach for assessing similarity between a search window CH and a model image CH is via an intersection analysis. In the context of the CHs, an intersection analysis essentially entails comparing the count in each corresponding color mix/distance category between a search window CH and a model image CH, and then identifying the smaller count. The identified smaller counts from each color mix/distance category are added together to produce a similarity value. It is noted that two matching images of an object will have a larger similarity value than non-matching images because the smallest count from each category will be nearly as large as the larger count, whereas the smaller count in non-matching images are likely to be significantly smaller than the larger value. Thus, the sum of the smaller counts from matching images should be larger than the sum of counts from nonmatching images.

[0016] Once the search window (or windows) potentially containing the object has been identified, the location can be fined tuned. This is accomplished by repeating the comparison process using the particular model image CH involved in the match and search window CHs generated from search windows respectively created by moving the identified search window up by one pixel row, down by one pixel row, to the left by one pixel column, and to the right by one pixel column, respectively. The search window associated with the largest similarity value of the five search windows (i.e., the original, up, down, right, & left) is selected. This refinement process is repeated until a newly selected window is associated with the maximum similarity value, unless the original search window already exhibits the maximum similarity value. This search window selected at the end of the refinement process is declared as containing the object.

[0017] The foregoing object recognition process employing color co-occurence histograms required that some of the search parameters be chosen ahead of time. Specifically, the size of the search window, color ranges and distance ranges must be selected. With enough trial-and-error experimentation, it is possible to discover reasonable values for these parameters. Such experimentation is not preferred because it can be tedious, and gives only a vague sense of the sensitivity of the search process to changes in these parameters. However, the choice of the parameters can be optimized via the unique false alarm analysis of the present invention.

[0018] The analysis employs a numerical algorithm that approximates a mathematical model for estimating the false alarm probability of the object recognition process. In the context of the CH object recognition system, the false alarm probability is the probability that the intersection of a search window CH and a model image CH (I_pm) will exceed the search threshold (αI_mm) when the object is not actually in the search window. Intuitively, the false alarm rate would be expected to increase as the search window grows in size, because there will be more opportunities for accidentally accounting for all the entries in the model CH. Also, the false alarm probability would be expected to decrease with an increase in the number of colors (n_c) and the number of distances (n_d), because these increases lead to more specific models that are less likely to be accidentally accounted for by a random background. Pulling these parameters in the opposite direction is the desire for a speedy algorithm and for one that does not required too many model viewpoints of each object. The algorithm would run faster if the search window were larger, because the whole image could then be searched faster. The intersection computation would be faster with fewer colors and distances. Fewer distances would also tend to generalize a model of a particular viewpoint to match those of more disparate viewpoints, and it would match better with objects that flex. The false alarm probability helps arbitrate between these desires, making it possible to set the parameters for the fastest possible execution and most tolerant matching without undue risk of false alarm.

[0019] The process of computing the false alarm probability begins by computing the probability of occurrence of a given CH on a random background region. To this end, the number n_kof ways that a distance interval can occur in the search window must be ascertained. For instance, a distance interval of 5 pixels can only occur in a limited number of ways in a given size image. This number can be computed using a simple Matlab program. For a fixed distance interval there are n_c(n_c+1)/2 possible unique, nonordered color pairs with corresponding bins in the CH, each pair occurring with probability p_ij. Each bin contains n_ij^kcounts. The probability of a partition of the n_kcolor pairs into n_c(n_c+1)/2 bins with n_ij^kin each bin is given by the multinomial distribution:

f (n_{11}^{k}, n_{12}^{k}, K, n_{ij}^{k}, K, n_{n_{c}, n_{c}}^{k}) = (\begin{matrix} n_{k} \\ n_{11}^{k}, n_{12}^{k}, K, n_{ij}^{k}, K, n_{n_{c}, n_{c}}^{k} \end{matrix}) * p_{11}^{n_{11}^{k}} p_{12}^{n_{12}^{k}} K p_{ij}^{n_{ij}^{k}} K p_{n_{c} n_{c}}^{n_{n_{c}, n_{c}}^{k}}

[0020] where

(\begin{matrix} n_{k} \\ n_{11}^{k}, n_{12}^{k}, K, n_{ij}^{k}, K, n_{n_{c}, n_{c}}^{k} \end{matrix}) = \frac{n_{k}!}{n_{11}^{k}! n_{12}^{k}! K n_{ij}^{k}! K n_{n_{c}, n_{c}}^{k}!}

[0021] The function f(n₁₁^k,n₁₂^k,K,n_ij^k,K,n_n^k,ⁿ) will be abbreviated as f(n_c,n_k)for convenience. This function f(n_c,n_k) is the probability of observing a given set of co-occurrences in a distance range (k−1,k). As it is assumed (albeit not totally true) that these probabilities are independent with respect to k, the probability of seeing a particular CH in the search window would be

P (CH) = \prod_{k = 1}^{n_{d}} f (n_{c}, n_{k})

[0022] If the intersection of the model CH and the image region CH exceeds a prescribed threshold without the object being there, then this is a false alarm. Ideally, to compute the probability of a false alarm, a list of all the CHs whose intersection with the model CH exceeds the threshold would be made. The probability of each CH in the list is computed as described above and summed to get the false alarm probability. However, the list is much too long. Instead the aforementioned multinomial distribution is approximated as a multivariate Gaussian distribution and integrated. This simplifies the problem from summing values in an enormous list to integrating a multidimensional Gaussian. Specifically, f(n_c,n_k) can be approximated by the ordinate of a multidimensional Gaussian:

f (n_{c}, n_{k}) \approx g_{k} (n) = \frac{1}{{(2 π)}^{m / 2} {| \sum |}^{1 / 2}} \exp [- \frac{1}{2} {(n - μ)}^{T} \sum^{- 1} (n - μ)]

[0023] where

[0024] m=n_c(n_c+1)/2−1

[0025] n=(n₁₁^k,n₁₂^k,K,n_ij^k,K,n_n_−1,n^k)^Ti≦j

[0026] μ=n_k(p₁₁,p₁₂,K,p_ij,K,p_n_−1,n)^Ti≦j

[0027] and where the inverse covariance is

\sum^{- 1} = {a_{rs}} = {\begin{matrix} \frac{1}{n_{k}} (\frac{1}{q_{r}} + \frac{1}{q^{*}}) & if r = s \\ \frac{1}{n_{k} q^{*}} & if r \neq s \end{matrix} (r, s) \in [1, 2, K, m]

[0028] The Gaussian g_k(n) can be integrated to give probabilities of sets of CHs occurring in the image background. The integration limits are given in terms of the number of co-occurrences of particular color pairs. However, the Gaussian only applies to a single specified distance range (k−1,k). Thus, it is still necessary to list all the CHs that could cause a false alarm. These would ideally be represented in terms of integration limits for the Gaussian approximation, but this list is too complex to characterize in this way. However, by simplifying the definition of a false alarm from “any background CH whose intersection with the model exceeds the threshold T=αI_mm” to “any background CH, each of whose entries n_ij^kexceeds a threshold T_ij^k=αm_ij^k”, integrating the Gaussian gives the probability of false alarm for all color pairs at one distance range (k−1,k). Since it can be is presumed from

P (CH) = \prod_{k = 1}^{n_{d}} f (n_{c}, n_{k})

[0029] that the co-occurrences at different distances are independent, the probability of a false alarm can be approximated as

P (false alarm) \approx \prod_{k = 1}^{n_{d}} \int_{Ω_{k}} g_{k} (n)

[0030] where Ω_kis the integration region for CH distance range (k−1,k). The integration region Ω_kis unbounded above and bounded below as n₁₁^k≧αm₁₁^k,n₁₂^k≧αm₁₂^k,K,n_ij^k≧αn_ij^k,K,n_n_−1,n^k. The actual integration is preferably performed using a Monte Carlo technique.

[0031] The foregoing approximation of the probability of a false alarm occurring can then be used to plot probability curves for one of the parameters (i.e., either the colors ranges, the distance ranges, or the search window size) at different values thereof, while holding the other parameters constant. The goal of this procedure is to find the optimum value. In the case of the color and distance ranges this optimum will occur at a point where the ranges are as large as possible while at the same time exhibiting an acceptable probability of false alarms. In the case of the search window size, the optimum is the largest window size possible with the same acceptable prom!bability of false alarms. It is noted that the plotting process can be repeated for various combinations of constant parameter values to find the curve that gives the largest color or distance range, or the largest window size, at an acceptably low false alarm probability. Preferably, the new values of the constant parameters would be chosen by performing an identical probability analysis using one of the constant parameters as the varying parameter and holding the others constant. Once initial values for each parameter have been computed the process is repeated iteratively until a set of mutually optimal values is established. It is also noted that the results are specific to the particular model image employed in the calculation, and so could be repeated for each model image. However, it was found in a tested embodiment that the computed optimal parameter values from one model image worked reasonably well for all the other models as well.

[0032] In view of the foregoing description, it is evident that the color.CH is an effective way to represent objects for recognition in images. By keeping track of pairs of pixels, it allows a variable amount of geometry to be added to the regular color histogram. This in turn allows the object recognition process to work in spite of confusing background clutter and moderate amounts of occlusion and object flexing. Specifically, by adjusting the distances over which the co-occurrences are checked, the sensitivity of the process can be adjusted to account for geometric changes in the object's appearance such as caused by viewpoint change or object flexing. The CH is also robust to partial occlusions, because it is not required that the image account for all the co-occurrences of the model.

[0033] In addition to the just described benefits, other advantages of the present invention will become apparent from the detailed description which follows hereinafter when taken in conjunction with the drawing figures which accompany it.

Description of the Drawings

[0034] The specific features, aspects, and advantages of the present invention will become better understood with regard to the following description, appended claims, and accompanying drawings where:

[0035] The file of the patent contains at least one drawing executed in color. Copies of the patent with color drawing(s) will be provided by the U.S. Patent and Trademark Office upon request and payment of the necessary fee.

[0036] FIG. 1 is a diagram depicting a general purpose computing device constituting an exemplary system for implementing the present invention.

[0037] FIG. 2 is a block diagram of an overall object recognition process for identifying the location of a modeled object in a search image according to the present invention.

[0038] FIG. 3 is a diagram graphically depicting two pixels separated by a vector in an image plane.

[0039] FIG. 4 is a block diagram of a process for accomplishing the rmodel image co-occurrence histogram creation module of the overall process of FIG. 3.

[0040] FIG. 5 is a block diagram of a process for accomplishing the co-occurrence histogram comparison module of the overall process of FIG. 3.

[0041] FIG. 6 is a block diagram of a process for refining the location of a search window containing the sought after object according to the present invention.

[0042] FIG. 7 is a diagram graphically depicting some of the ways that a distance of five pixels can occur in a 6×6 search window.

[0043] FIGS. 8(a) through 8(l) are images depicting various model images of a figurine.

[0044] FIG. 9 is a graph plotting false alarm probability as a function of search window size. The left, middle and right curves represent the use of 6, 8 and 10 color ranges respectively. All the curves are based on 12 distance ranges.

[0045] FIG. 10 is a graph plotting false alarm probability as a function of the number of color ranges. The left, middle and right curves represent the use of 14, 12 and 10 distance ranges respectively. All the curves are based on a search window scaling factor of 3.

[0046] FIGS. 11 and 12 are images showing the results of employing the object recognition process according to the present invention and the model images shown in FIGS. 8(a) through 8(l). FIG. 11 shows that the location of the figurine was found in spite of background clutter and partial occlusion. FIG. 12 shows that the location of the figurine was found in spite of being flexed such that it was in a standing position with one arm outstretch instead of sitting with arms at the side as in the model images.

Detailed Description of the Preferred Embodiments

[0047] In the following description of the preferred embodiments of the present invention, reference is made to the accompanying, drawings which form a part hereof, and in which is shown by way of illustration specific embodiments in which the invention may be practiced. It is understood that other embodiments may be utilized and structural changes may be made without departing from the scope of the present invention.

[0048] FIG. 1 and the following discussion are intended to provide a brief, general description of a suitable computing environment in :which the invention may be implemented. Although not required, the invention will be described in the general context of computer-executable instructions, such as program modules, being executed by a personal computer. Generally, program modules include routines, programs, objects, components, data structures, etc. that perform particular tasks or implement particular abstract data types. Moreover, those skilled in the art will appreciate that the invention may be practiced with other computer system configurations, including hand-held devices, multiprocessor systems, microprocessor-based or programmable consumer electronics, network PCs, minicomputers, mainframe computers, and the like. The invention may also be practiced in distributed computing environments where tasks are performed by remote processing devices that are linked through a communications network. In a distributed computing environment, program modules may be located in both local and remote memory storage devices.

[0049] With reference to FIG. 1, an exemplary system for implementing the invention includes a general purpose computing device in the form of a conventional personal computer 20, including a processing unit 21, a system memory 22, and a system bus 23 that couples various system components including the system memory to the processing unit 21. The system bus 23 may be any of several types of bus structures including a memory bus or memory controller, a peripheral bus, and a local bus using any of a variety of bus architectures. The system memory includes read only memory (ROM) 24 and random access memory (RAM) 25. A basic input/output system 26 (BIOS), containing the basic routine that helps to transfer information between elements within the personal computer 20, such as during start-up, is stored in ROM 24. The personal computer 20 further includes a hard disk drive 27 for reading from and writing to a hard disk, not shown, a magnetic disk drive 28 for reading from or writing to a removable magnetic disk 29, and an optical disk drive 30 for reading from or writing to a removable optical disk 31 such as a CD ROM or other optical media. The hard disk drive 27, magnetic disk drive 28, and optical disk drive 30 are connected to the system bus 23 by a hard disk drive interface 32, a magnetic disk drive interface 33, and an optical drive interface 34, respectively. The drives and their associated computer-readable media provide nonvolatile storage of computer readable instructions, data structures, program modules and other data for the personal computer 20. Although the exemplary environment described herein employs a hard disk, a removable magnetic disk 29 and a removable optical disk 31, it should be appreciated by those skilled in the art that other types of computer readable media which can store data that is accessible by a computer, such as magnetic cassettes, flash memory cards, digital video disks, Bernoulli cartridges, random access memories (RAMs), read only memories (ROMs), and the like, may also be used in the exemplary operating environment.

[0050] A number of program modules may be stored on the hard disk, magnetic disk 29, optical disk 31, ROM 24 or RAM 25, including an operating system 35, one or more application programs 36, other program modules 37, and program data 38. A user may enter commands and information into the personal computer 20 through input devices such as a keyboard 40 and pointing device 42. Of particular significance to the present invention, a camera 55 (such as a digital/electronic still or video camera, or film/photographic scanner) capable of capturing a sequence of images 56 can also be included as an input device to the personal computer 20. The images 56 are input into the computer 20 via an appropriate camera interface 57. This interface 57 is connected to. the system bus 23, thereby allowing the images to be routed to and stored in the RAM 25, or one of the other data storage devices associated with the computer 20. However, it is noted that image data can be input into the computer 20 from any of the aforementioned computer-readable media as well, without requiring the use of the camera 55. Other input devices (not shown) may include a microphone, joystick, game pad, satellite dish, scanner, or the like. These and other input devices are often connected to the processing unit 21 through a serial port interface 46 that is coupled to the system bus, but may be connected by other interfaces, such as a parallel port, game port or a universal serial bus (USB). A monitor 47 or other type of display device is also connected to the system bus 23 via an interface, such as a video adapter 48. In addition to the monitor, personal computers typically include other peripheral output devices (not shown), such as speakers and printers.

[0051] The personal computer 20 may operate in a networked environment using logical connections to one or more remote computers, such as a remote computer 49. The remote computer 49 may be an:other personal computer, a server, a router, a network PC, a peer device or other common network node, and typically includes many or all of the elements described above relative to the personal computer 20, although only a memory storage device 50 has been illustrated in FIG. 1. The logical connections depicted in FIG. 1 include a local area network (LAN) 51 and a wide area network (WAN) 52. Such networking environments are commonplace in offices, enterprise-wide computer networks, intranets and the Internet.

[0052] When used in a LAN networking environment, the personal computer 20 is connected to the local network 51 through a network interface or adapter 53. When used in a WAN networking environment, the personal computer 20 typically includes a modem 54 or other means for establishing communications over the wide area network 52, such as the Internet. The modem 54, which may be internal or external, is connected to the system bus 23 via the serial port interface 46. In a networked environment, program modules depicted relative to the personal computer 20, or portions thereof, may be stored in the remote memory storage device. It will be appreciated that the network connections shown are exemplary and other means of establishing a communications link between the computers may be used.

[0053] The exemplary operating environment having now been discussed, the remaining part of this description section will be devoted to a description of the program modules embodying an object recognition process that identifies the location of a modeled object in an image according to the present invention.

[0054] 1.0 The Obiect Recognition Process.

[0055] The general process will be explained first. Then, each step in the process will be discussed in detail in the sections that follow. Referring to FIG. 2, the first step 200 in the general object recognition process involves capturing model images of the object that will be searched for in an image (hereinafter referred to as the search image). These model images are captured from a plurality of “viewpoints” around the object. For the purposes of the present description the term “viewpoint” refers to a viewing direction toward the object, rather than a specific point in space. Next, in step 202, a co-occurrence histogram (CH) is computed for each of the model images. Then, in step 204, a series of search windows, of a prescribed size, are generated from overlapping portions of the search image. A co-occurrence histogram is computed for each of the search windows in step 206. A comparison between the CH of each of the model images and the CH of each of the search windows is then conducted to assess the similarity between each pair of CHs (step 208). Finally, in step 210, any search window that is associated with a search window CH having a degree of similarity to one of the model image CHs which exceeds a prescribed search threshold is designated as potentially containing the object being sought. This designation can be presumed final, or refined as will be discussed later.

[0056] 1.1 Creating Model Images.

[0057] As described above, the object recognition process first involves creating model images of the object whose location is to be identified in the search image. As the object may be depicted in the search image in any orientation, it is preferred that a number of model images be captured, of which each shows the object from a different viewpoint. Ideally, these model images would be taken from viewpoints spaced at equal angles from each other around the object.

[0058] The angle separating each model image viewpoint should be selected so that a high degree of matching (e.g., by exceeding the aforementioned search threshold) exists between one of the model images and a portion of the search image actually containing the modeled object. On the other hand, the aforementioned comparison process becomes more onerous as the number of model images increases. Therefore, it is preferred that the angle be large enough to minimize the number of model images required, while still providing the desired degree of match. One way of determining this optimal separation angle is to require a good degree of matching between each adjacent model image viewpoint, as measured by an appropriate threshold. Specifically, this can be accomplished by selecting a series of candidate separation angles, and for each candidate, assessing the degree of similarity:between at least one pair of model images having adjacent viewpoints separated by the candidate separation angle. The largest candidate separation angle that has an assessed degree of similarity exceeding a prescribed model image separation angle threshold is designated as the separation angle to be employed in capturing the model images of the object. The aforementioned model image separation angle threshold is chosen so that the number of model images required is minimized, while still ensuring the model images cover every possible orientation of the sought after object in the search image.

[0059] Another aspect of the model image creation process is the distance away from the object that the images are captured. If it is generally known how far away from the camera the object is in the search image, then the model images could be captured at approximately the same distance away from the object being modeled. However, if this distance in not known, then it would be desirable to capture more than one model image at each viewpoint angle, where each image is captured at a different distance away from the object. In this way, the model images are more likely to accurately model the object of interest in the search image. It should be noted that the aforementioned similarity test used to ensure a good degree of matching between model, images taken at adjacent angular viewpoints could also be employed to ensure a good degree of matching between model images taken at the same angular viewpoint, but at different distances away from the object. Specifically, this can be accomplished by selecting a series of candidate separation distances, and for each candidate, assessing the degree of similarity between at least, one pair of adjacent model images separated by the candidate separation distance. The largest candidate separation distance that has an assessed degree of similarity exceeding a prescribed model image separation distance threshold is designated at the aforementioned separation distance to be employed in capturing the model images of the object. Similar to the aforementioned angle threshold, the model image separation distance threshold is chosen so that the number of model images required is minimize, while still ensuring the model images cover every possible distance away for the sought after object that might be encountered in the search image. In this way, it is more likely that -one of the model images associated with a particular angular viewpoint will match (i.e., exceed the aforementioned search threshold) an object in the search image having a similar angular viewpoint.

[0060] 1.2. Computing Co-occurrence Histograms (CHs) For Each Model Image.

[0061] As discussed in the description of the general object recognition process, a co-occurrence histogram (CH) is computed for each model image. It is noted that for the following description pixel “color” will be employed as the basis for computing the CH. The pixel “color” can be defined using any conventional characterization of the color space. Thus, each model image will be represented as a color CH. However, it is not intended that the present invention be limited to CHs based only on pixel color. Rather, it is envisioned that various other pixel characterizations could also be employed. For example, pixel gray levels could be used for gray scale images, or pixel depth values might be used.

[0062] Essentially, a color CH is a count of the number of occurrences of pairs of color pixels c₁=(R₁,G₁,B₂)and c₂=(R₂,G₂,B₁) separated by a vector in the image plane (Δx,Δy), as depicted in FIG. 3. It is noted that for (Δx,Δy)=(0,0), the color CH reduces to a regular color histogram (where c₁=c₂). For convenience of notation, the color co-occurrence histogram will be written symbolically as CH(c₁,c₂,Δx,Δy).

[0063] In order to make the CHs invariant to rotation in the image plane, the direction of (Δx,Δy) is ignored and only the magnitude d={square root over ((Δx)²+(Δy)²)} is tracked. In addition, the pixel colors-are quantized into a set of n_crepresentative colors C={c₁,c₂,K,c_n}, and the distances are quantized into a set of n_ddistance ranges D={(0,1),(1,2),K,(n_d−1,n_d)}. As to the quantization of the pixel colors, in general the process involves separately determining the colors of each model image. This is preferably accomplished using a standard nearest neighbor k-means algorithm, where the number of means is n_c, and where the Euclidian distance in (R,G,B) is employed. The whole gamut of determined colors is then divided into a set of color ranges to produce a series of quantized colors. Each pixel is then assigned to a quantized color by mapping its actual color from the gamut to its most similar (nearest) quantized color in the set C={c₁,c₂,K,c_n}. As for the distance ranges, it is noted that they can be of equal size as indicated above, or can also be adjusted in size so that each range is smaller or larger, as desired, depending on what works best for the application. The quantization in color and distance means the CH is represented by CH(i,j,k), where i and j index the two colors from set C, and k indexes the distance range from set D. The question of how to choose the number of color ranges and distance ranges will be discussed in connection with a false alarm rate analysis later in this description.

[0064] In practice, a model image CH is computed via the steps illustrated in FIG. 4. In the first step 400, the colors exhibited by the pixels of the model image are determined, as indicated above. The overall gamut of possible pixel “colors” exhibited by the pixels is then divided into a series of discrete color ranges (step 401). Preferably these color ranges are equal in size, although they need not be made so. The model image is quantized into the n_ccolors. This is done as shown in step 402 by comparing the raw color exhibited by a pixel to the set of quantized colors and assigning the pixel to the nearest quantized color.

[0065] The distance ranges are also established, either at this point, or more preferably, prior to beginning the process. This involves choosing a maximum distance of interest and dividing it into a series of distance ranges (step 404). For instance, in a tested embodiment a maximum distance range of 12 pixels was chosen and divided into equal sized distance ranges each being one pixel wide (i.e., (1-2 pixels), (2-3 pixels), (3-4 pixels) . . . , (11-12 pixels)). Every possible unique, non-ordered pixel pair is then identified in step 406. For each of the identified pixel pairs, the next step 408 is to quantize the distance between the pixels making up the pair by assigning the pair to the appropriate one of the previously selected distance ranges.

[0066] Now that each pixel pair has been quanitzed as to the colors of the pixel making it up and the distance separating the pixels, the next step 410 in the process is to categorize each pixel pair by the assigned color range of each constituent pixel, as well as its assigned distance range. A count is then established for the number of pixel pairs exhibiting, the same mix of assigned quantized colors and having the same assigned distance range (step 412).

[0067] 1.3 Generating Search Windows and Search Window CHs.

[0068] In searching an image for an object, the image is scanned for a window that gives a CH similar to one of the model CHs. This is repeated for each of the model images. The image windows from which the model and image CHs are computed need not be the same size. In fact it is desirable to make the image search window as large as reasonably possible to increase the speed of scanning the whole image. However, the search window size should not be made so large that a substantial risk of a false match or alarm is introduced, Thus, it is desirable to choose a search window size which is as large as possible, while at the same time not being so large as to risk false alarms. The issue of choosing such a search window size is also addressed later using the aforementioned false alarm analysis.

[0069] The search windows are preferably generated by cordoning off the search image into a series of equally-sized regions having the chosen search window size and shape configuration. These regions preferably overlap both side-to-side and up-and-down, where possible (i.e. the edges of the image will preclude overlapping in some directions for regions adjacent the edge). In the tested embodiment of the present invention, the overlap was set to one-half the width and height of the search window.

[0070] It is noted that while rectangular search windows are envisioned to be the most practical shape, other shapes may be feasible depending on the search image being searched. For example, triangular, polygonal, circular or even elliptical search windows might be employed. For the purposes of explaining the present invention a rectangular search window will be assumed.

[0071] Once a search window has been generated, its CH is computed. This can occur as each search window is identified, or after all the windows have been generated. Regardless, a CH is calculated for each search window in the same manner that they were generated in connection with the model images. In particular, the same quantized colors & distance ranges are employed in the quantization steps of the process. Specifically, for those pixels of the search window having a color falling within one of the previously computed pixel color ranges, the pixel's color is quantized by assigning the pixel to the pixel color range into which its color falls. Then, as with the model image CHs, every possible unique, non-ordered pair of pixels in the search window is identified. Next, the distance separating the two pixels of those pixel pairs having separation distances that fall within one of the previously computed distance ranges is quantized by assigning the pixel pair to the distance range into which its separation distance falls. Finally, each pixel pair is categorized as was done with the model image CHs, and a count of the number of categorized pixel pairs having the same pixel color mix/distance category is established.

[0072] 1.4 Comparing Model Images And Search Windows.

[0073] Image and model CHs, CH_p(i,j,k) and CH_m(i,j,k) respectively, are preferably compared by computing their intersection, i.e.,

\begin{matrix} I_{pm} = \sum_{i = 1}^{n_{c}} \sum_{j = 1}^{n_{c}} \sum_{k = 1}^{n_{d}} \min [({CH}_{p} (i, j, k), {CH}_{m} (i, j, k)] & (1) \end{matrix}

[0074] This intersection indicates how well the image CH accounts for the model CH. For example, if the image accounts for all the entries in the model CH, then the intersection will be equal to the sum of all the entries in the model CH, or I_mm.

[0075] Specifically, as illustrated in FIG. 5, the count in each corresponding color mix/distance category is compared between a search window CH and a model image CH (step 500). The purpose is of this comparison is to identify the smaller of the two counts in each category, as indicated in step 502. The identified smaller counts from each color mix/distance category are added together to produce a similarity value (step 504). The larger this similarity value, the greater the similarity between the associated model image and search window.

[0076] It is noted that two matching images of an object will have a larger similarity value than non-matching images because the smallest count from each category will be nearly as large as the larger count, whereas the smaller count in non-matching images is likely to be significantly smaller than the larger value. Thus, the sum of the smaller counts from matching images should be larger than the sum for nonmatching images.

[0077] The intersection method is an attractive measure of similarity, because if the image CH contains extra entries from the background, they will not reduce the quality of the match. However, if the image CH has large values for all its entries, say due to examining a large region in the image, it could accidentally account for the entries in the model CH, causing a false alarm. As indicated earlier, the effect of the search window size, as well as other parameters, is examined based on the false alarm probability in the next section.

[0078] 1.5 Findinq a Match Between A Model Image And A Search Window.

[0079] An object is potentially considered to be contained in a particular portion of the search image associated with a search window if the previously described comparison between a model image and the search window yields a similarity value that exceeds a search threshold T. For instance, if it is desired that the image account for at least a fraction α of the co-occurrences in the model, then T=αI_mm. The exact value of α will vary depending on the degree of accuracy required for the application for which the present invention is being employed. However, in a tested embodiment it was determined that an α of 0.9 produced satisfactory results.

[0080] Once it is determined that the intersection of a model and image search window CH exceeds the prescribed search threshold, it can be handled in two ways. In a first version, any intersection that exceeds the threshold would cause the associated search window to be declared as potentially containing the object. This version would allow identifying each location of the object in a search image which depicts more than one of the objects. In a second version, the threshold would still be used, but only the greatest of the similarity values exceeding the threshold would be chosen to indicate a search window potentially contains the object. This version would be especially useful where it is known only one of the objects exists in the image. It could also be useful in that the threshold could be made lower and the object still located, for example in a noisy image.

[0081] 1.6 Refining the Obeect Recognition Process.

[0082] The search window associated with a similarity value exceeding the prescribed search threshold (either any search window or the search window having the greatest similarity value exceeding the threshold, as the case may be) can be considered to contain the object being sought without further processing. However, depending on the application, it might be, desirable to further refine the exact location of the object within the search image. This refinement can be achieved using a two-part process. This two-part process begins as described above to find one or more search windows potentially containing the object being sought.

[0083] Once such a search window (or windows if looking for multiple objects in the search image) has been identified, a second, higher resolution, search is performed to find the exact location of the object in the search image. The high resolution search is essentially a pixel-by-pixel hill-climbing search. Referring to FIG. 6, the high resolution search is accomplished in connection with each of the identified search windows by first designating the identified search window as the primary search window (step 600). Then, in step 602, additional search windows are generated from portions of the search image respectively defined by moving the primary search window up by one pixel row, down by one pixel row, to the left by one pixel column, and to the right by one pixel column. This produces a total of five search windows in most cases, i.e. the primary and four translated windows. Of course, if the primary search window is located at the edge of the search image, it will not be possible to generate all four of the translated search windows and the refinement process would be conducted with fewer search windows (i.e., three or four).

[0084] The refinement process continues with the computation of a color CH for each of the additional search windows (step 604) and assessing the similarity between the CHs of each of the additional search windows and the CH of the model image associated with the original search window (step 606). These steps are performed in the same manner as described in connection with the process used to identify the original search window. The result of the foregoing steps is to produce a similarity value for each of the translated search windows. It is next determined if any of the similarity values associated with the translated search windows exceeds the similarity value previously computed for the primary search window, as indicated by step 608. If none of the similarity values associated with the additional search windows exceeds the similarity value of the primary search window, then in accordance with step 610, the primary search window is declared as actually containing the object being sought. If, however, one of the similarity values computed for the additional search windows does exceed that of the primary search window, then that particular additional search window is designated the current primary search window (step 612) and the process of steps 600 through 610 or 612 is repeated using the newly designated primary search window in place of the previously designated primary search window. Thus, the refinement process is repeated, unless the originally designated primary search window is associated with the maximum similarity value, and continues until a newly designated primary search window is associated with the maximum similarity value. This fine tunes the location of a search window containing the sought after object. Either way, a currently designated primary search window having the maximum similarity value in comparison to the additional search windows is declared to contain the object being sought (see step 610).

[0085] 1.7 Differences In Scale.

[0086] There is a concern that scale or zoom differences between the model images and the search windows could skew the results of the object recognition process. For example, if the scale of the model images is different from that of the search windows, the distance ranges used to categorize pixel pairs in the respective images would not be the same, even if the object was contained in the search window. However, it is noted that the distance ranges would be proportional to each other. Thus, the scale problem could be avoided by normalizing the distance ranges between the images. This can be accomplished by designating the model image size to be associated with a scale of one. Then, the relative scale of the search windows taken from the image being searched would be used as a scaling factor to equalize the distance ranges. For example, suppose the search windows taken from the image being searched have a scale twice that of the model images. The model images would be quantized for the distance ranges as described previously. However, these ranges would be doubled in size when quantizing the search windows. In this way the number of ranges stays the same. When comparing the images, the number of pixel pairs in the unaltered distance ranges associated with the model image CH are compared to the number of pixel pairs in the corresponding doubled distance ranges in the search image CH. The result of this scale compensating process is that approximately the same number of pixel pairs will be associated with the “corresponding” distance ranges in the compared model image and search window when these images match.

[0087] 2.0 Using False Alarm Probability To Choose Obiect Recoqnition Process Parameters.

[0088] With enough trial-and-error experimentation, it is possible to discover reasonable values for the number of colors n_c, the number of distances n_d, and the search window size. However, this experimentation is tedious, and it gives only a vague sense of the sensitivity of the algorithm to changes in these parameters. To remedy this, a mathematical model can be created to approximate the false alarm probability and a numerical algorithm derived to approximate the mathematical model. This provides a principled way of choosing the process's adjustable parameters.

[0089] The false alarm probability is the probability that the intersection I_pmwill exceed the match threshold αI_mmwhen the object is not actually in the search window. Intuitively, it is expected that the false alarm rate will increase as the search window grows, because there will be more opportunities for accidentally accounting for all the entries in the model CH. Also, it is expected that the false alarm probability will decrease as the number of colors n_cand the number of distances n_dis increased, because these increases lead to more specific models that are less likely to be accidentally accounted for by a random background. Pulling these parameters in the opposite direction is desired for a speedy algorithm and for one that does not required too many model viewpoints of each object. The algorithm would run faster if the search window were larger, because the whole image could then be searched faster. The intersection computation would be faster with fewer colors and,distances. Fewer distances would tend to generalize a model of a particular viewpoint to match those of more disparate viewpoints, and it would match better with objects that flex. The false alarm probability helps arbitrate between these desires, letting the parameters to be set for the fastest possible execution and most tolerant matching without undue risk of false alarm.

[0090] The next subsection 2.1 shows the derivation of the probability of a given CH occurring on a random background. This is then used as discussed in subsection 2.2 to compute the probability that a random background will match a model CH. And finally, in subsection 2.3 the analysis is used to choose parameters for a real object.

[0091] 2.1 Probability Of A Given CH.

[0092] The derivation begins by computing the probability of occurrence of a given CH on a random background region of size h×w. The image is quantized into n_ccolors, with color i occurring with probability p_i. The co-occurrence histogram looks at pairs of colors, so it is observed that the probability of seeing colors i and j at (x_i,y_i) and (x_j,y_j) respectively is p_ip_j. The CH would not distinguish the case where the two colors exchanged positions, so it is noted that the probability of colors i and j being at (x_i,y_i) and (x_j,y_j) respectively or (x_j,y_j) and (x_i,y_i) respectively is

\begin{matrix} p_{ij} = {\begin{matrix} 2 p_{i} p_{j} & if i \neq j \\ p_{i} p_{j} & if i = j \end{matrix} i, j \in {1, 2, K, n_{c}} . & (2) \end{matrix}

[0093] For lack of any reason to believe otherwise at this stage, it will be assumed that the n_ccolors are uniformly distributed in the image background, giving p_i=1/n_c. The probabilities of the unique pairs sum to one as expected:

\begin{matrix} \sum_{i = 1}^{n_{c}} \sum_{j = 1}^{i} p_{ij} = 1 & (3) \end{matrix}

[0094] Computing the probability of a given CH is complex, because there is no obvious way to break the probability into a series of independent probabilities that can be multiplied together. An approximation will be made that the probability of a juxtaposition of colors in one distance interval is independent of the probability in any other distance interval. In reality there is some correlation between distances. For instance, if there are many pairs of red pixels at one separation distance, it would be expected to see many of the same color pairs at nearby separation distances.

[0095] For a given distance interval (k−1,k), there are many different pairs of pixels separated by that amount in an h×w search window. For instance, the diagram shown in FIG. 7 depicts some of the ways that a distance of five pixels can occur in a search window of size 6×6. The variable nk defines the number of ways that a distance interval (k−1,k) can occur in the search window. This number n_kis computed using a standard Matlab program.

[0096] For a fixed distance interval (k−1,k) there are n_c(n_c+1)/2 possible unique, nonordered color pairs with corresponding bins in the CH, each pair occurring with probability p_ij. Each bin contains n_ij^kcounts (n_ij^k=CH(i,j,k), i≦j), with

n_{ij}^{k} \geq 0, n_{k} = \sum_{i = 1}^{n_{c}} \sum_{j = 1}^{i} n_{ij}^{k},

[0097] and i,jε{1,2,K,n_c}. The probability of a partition of particular color pairs into bins with n_ij^kin each bin is p₁₁ⁿp₁₂ⁿKp_ijⁿKp_n_nⁿ. This partition can happen in

(\begin{matrix} n_{k} \\ n_{11}^{k}, n_{12}^{k}, K, n_{ij}^{k}, K, n_{n_{c}, n_{c}}^{k} \end{matrix})

[0098] ways, each of which produces an identical CH. Thus, the probability of a partition of the n_kcolor pairs into n_c(n_c+1)/2 bins with n_ij^kin each bin is given by the multinomial distribution:

\begin{matrix} f (n_{11}^{k}, n_{12}^{k}, K, n_{ij}^{k}, K, n_{n_{c}, n_{c}}^{k}) = (\begin{matrix} n_{k} \\ n_{11}^{k}, n_{12}^{k}, K, n_{ij}^{k}, K, n_{n_{c}, n_{c}}^{k} \end{matrix}) * p_{11}^{n_{11}^{k}} p_{12}^{n_{12}^{k}} K p_{ij}^{n_{ij}^{k}} K p_{n_{c} n_{c}}^{n_{n_{c}, n_{c}}^{k}} & (4) \end{matrix}

[0099] where

(\begin{matrix} n_{k} \\ n_{11}^{k}, n_{12}^{k}, K, n_{ij}^{k}, K, n_{n_{c}, n_{c}}^{k} \end{matrix}) = \frac{n_{k}!}{n_{11}^{k}! n_{12}^{k}! K n_{ij}^{k}! K n_{n_{c}, n_{c}}^{k}!}

[0100] The function f(n₁₁^k,n₁₂^k,K,n_ij^k,n_n_,n^k) will be abbreviated as f(n_c,n_k) in the following description.

[0101] The function f(n_c,n_k) is the probability of observing a given set of co-occurrences in a distance range (k−1,k). The aforementioned assumption says that these probabilities are independent with respect to k, so the probability of seeing a particular CH in the search window is

\begin{matrix} P (CH) = \prod_{k = 1}^{n_{d}} f (n_{c}, n_{k}) & (5) \end{matrix}

[0102] 2.2 False Alarm Probability Approximation.

[0103] If the intersection of the model CH and the image region CH exceeds the threshold T without the object being there, then this is a false alarm. Ideally, to compute the probability of a false alarm, a list of all the CHs whose intersection with the model CH exceeds the threshold would be made. Equation (5) would then be used to compute the probability of each CH in the list and sum the probabilities to get the false alarm probability. However, this list is much too long.

[0104] Instead, the multinomial distribution is approximated as a multivariate Gaussian distribution and integrated. This simplifies the problem from summing values in an enormous list to integrating a multidimensional Gaussian.

[0105] From [12] it is known that Equation (4) can be approximated by the ordinate of a multidimensional Gaussian:

\begin{matrix} f (n_{c}, n_{k}) \approx g_{k} (n) = \frac{1}{{(2 π)}^{m / 2} {| Σ |}^{1 / 2}} \exp [- \frac{1}{2} {(n - μ)}^{T} Σ^{- 1} (n - μ)] & (6) \end{matrix}

[0106] where

[0107] m=n_c(n_c+1)/2−1

[0108] n=(n₁₁^k,n₁₂^k,K,n_ij^k,K,n_n_−1,n^k)^Ti≦j

[0109] μ=n_k(p₁₁,p₁₂,K,p_ij,K,p_n_−1,n)^Ti≦j

[0110] It is noted that the dimensionality of these vectors is one less than the number of color pairs n_c(n^c+1)/2. This is because once all but the last color co-occurrence bin has been filled, the last bin must contain exactly enough counts to make the total n_k.

[0111] To specify the inverse covariance Σ⁻¹the co-occurrence probabilities are renamed: q₁=p₁₁, q₂=p₁₂, . . . , q_n_(n_+1)/2−1=p_n_−1,n, and q*=p_n_,n. Then the inverse covariance is

Σ^{- 1} = {a_{rs}} = {\begin{matrix} \frac{1}{n_{k}} (\frac{1}{q_{r}} + \frac{1}{q^{*}}) if r = s \\ \frac{1}{n_{k} q^{*}} if r \neq s \end{matrix} (r, s) \in [1, 2, K, m]

[0112] The Gaussian g_k(n) can be integrated to give probabilities of sets of CHs occurring in the image background. The integration limits are given in terms of the number of co-occurrences of particular color pairs. The Gaussian only applies to a single specified distance range (k−1,k).

[0113] The problem of listing all the CHs that could cause a false alarm must still be addressed. These would ideally be represented in terms of integration limits for the Gaussian approximation. However, this list is too complex to characterize in this way. Thus, the definition of a false alarm is simplified from “any background CH whose intersection with the model! exceeds the threshold T=αI_mm” to “any background CH, each of whose entries n_ij^kexceeds a threshold T_ij^k=αm_ij^k.” Here m_ij^k=I_m(i,j,k) is the model CH, n_ij^k=CH_p(i,j,k) is the search window CH, and 0<α≦1 is the fraction of the histogram counts that we require to declare the model found in the image. This tends to underestimate the probability of a false alarm, but experiments show that it is still approximately correct.

[0114] Integrating the Gaussian gives the probability of false alarm for all color pairs at one distance range (k−1,k). From Equation (5), the co-occurrences at different distances are considered to be independent. Thus, the probability of a false alarm can be approximated as:

\begin{matrix} P (false alarm) \approx \prod_{k = 1}^{n_{d}} \int_{Ω_{k}} g_{k} (n) & (7) \end{matrix}

[0115] where Ω_kis the integration region for CH distance range (k−1,k). The region Ω_kis unbounded above and bounded below as n₁₁^k≧αm₁₁^k,n₁₂^k≧αm₁₂_k,K,n_ij^k≧αn_ij^k,K,n_n_−1,n^k≧αm_n_−1,n. The actual integration is preferably performed using a standard Monte Carlo technique.

[0116] The next subsection shows plots of the false alarm probability to demonstrate how the analysis can be used to pick some of the algorithm's operating parameters.

[0117] 2.3. Using False Alarm Probabilitv To Choose Parameters.

[0118] The present object recognition system and process requires the specification of three parameters: the number of colors n_c, the number of distances n_d, and the size of the search window h×w. This section shows how a false alarm analysis can be used to guide the choice of these parameters for a model image—namely the one shown in FIG. 8(h).

[0119] The analysis begins by deciding the size of the search window. The graph in FIG. 9 illustrates the effect on the false alarm rate of scaling the search window size from the original size of the model image (i.e., 100×60) up to 5 times this size in each dimension, while fixing the other two parameters of interest—namely setting n_c=6,8, and 12 respectively, and n_d=12. The false alarm rate increases as expected, staying essential zero if the scaling factor is less than three. As the goal in the selection of a search window size is to use the largest window possible without risking a significant possibility of false alarms, a search window size of 300×180 (i.e., 3 times the model image size) is selected.

[0120] Next, the optimal number of colors n, is chosen for the aforementioned k-means color quantization. Using a search window size of 3 times the model window (as chosen above) and n_d=14,12, and 10 respective distance ranges, the false alarm probability as a function of n_cis shown in FIG. 10. The false alarm probability decreases with more colors, because more colors make the model more specific and less prone to matching a random background. Based on this plot, n_c=8 is chosen because that is where the false alarm rate first goes to nearly zero.

[0121] Finally, the optimal number of distances n_dis chosen using the search window size of 3 times the model window and n_c+8 colors. The false alarm probability as a function of n_ddecreases with more distances because like more colors, more distance ranges makes the model more specific and less prone to matching a random background. In a tested embodiment, it was found that n_d=12 should be chosen because that is where th,e false alarm rate first goes to nearly zero.

[0122] It is noted that the plots of the false alarm probability, shown in FIGS. 9 and 10, illustrate how the probability varies as two of the three parameters of interest are fixed, while the third variable was varied. It is also noted that in the process of choosing the optimum search window size, n_cand n_dwere fixed at the values later found to be nearly optimal for these parameters. The same was true for n_din the process used to select the optimal number of colors n_c. This was no accident. The process of selecting the optimal parameters actually required a series of iterations through the parameters and associated plots to find a mutually good-working set of values. Only the final iteration is depicted in FIGS. 9 and 10.

[0123] An alternate analysis approach can also be employed in selecting the optimal, search parameters. Generally, this alternate approach involves searching through variations of all parameters simultaneously rather then fixing all but one of the parameters as described above. This simultaneous search is accomplished by exhaustively computing the false alarm rate for a substantial number of combinations of candidate search parameters or by using a standard numerical search routine to find the best combination.

[0124] The foregoing analysis approaches will result in a reasonable set of parameters to use for the first model image. Ideally a similar analysis would be performed for the model images of other viewpoints. However, it was found that using the same parameter values for all the model images produced satisfactory results. Thus, it may not be necessary to perform the analysis on all the model images if some small deviation from the optimal parameter values is acceptable. It is interesting to note that for the autocorrelograms used in [9], the authors used 64 colors and 4 distances, compared to the 8 colors and 12 distances selected in the foregoing analysis for the present system and process. It is believed that using the parameter values from [9] would likely make the present process miss many good matches due to unavoidable color shifts induced by inter-reflections with surrounding objects.

[0125] 3.0 Experiments And Results.

[0126] Referring to FIGS. 8(a)-(l), a series of model images are shown of a figurine object taken as it was rotating on a motorized turntable. A long focal length lens (75 mm on a C-mount camera) was chosen so that the camera could be placed far away from the objects (about eight meters). This reduced the effects of distance-induced scale change as the objects move around in the scene. The camera's automatic gain control and automatic white balance were also disabled to help make the images consistent. The model images with captured with a black background to make segmentation easier.

[0127] As discussed previously, an important number to choose for the present system and process is the number of image models to use for a given object. For the object depicted in FIGS. 8(a)-(l), we used 12 images to cover 360° of rotation. This choice was based on an analysis of how well the model images match themselves. We began by taking model images of the object every 5° for a total of 72 images. Adjacent images were then compared to ensure a good degree of matching existed. It was found, however, that as few as 12 model images would still provide a good degree of matching between the adjacent model images. Thus, it was decided that 12 model images were adequate for modeling the full set of viewpoints around 360° for the tested embodiment.

[0128] It is noted that while only one angle of the object's point of view was varied for testing, it would be possible to generalize to an arbitrary point of view using more model images. Thus, instead of using just 12 model images taken around an essentially horizontal viewing circle, more images taken about other non-horizontal viewing circles would also be captured. For example, to ensure the object depicted in FIGS. 8(a)-(l) could be recognized from any angle, about 12²=144 model images around the viewing sphere would be necessary, i.e., 12 equally spaced images around each of 12 equally spaced viewing circles.

[0129] For test images, a bright, multi-colored background was used to make clutter. It was found that inter-reflections from the colored background caused a significant color shift on the objects, which the present system and process overcame by using a rough quantization of colors (i.e., n_c=8).

[0130] Some typical results are shown in FIGS. 11 and 12. FIG. 11 shows that the present system and process finds the figurine modeled with the images of FIGS. 8(a)-(l) despite significant occlusion. In FIG. 12, the system and process was able to find the aforementioned figurine standing with one arm outstretched, even though the model images (of FIGS. 8(a)-(l)) show it sitting with arms at the sides. This demonstrates that the system and process can tolerate some flexing in an object.

References

[0131] [1] D. P. Huttenlocher, G. A. Klanderman, and W. J. Ricklidge, “Comparing Images Using the Hausdorff Distance,” IEEE Transactions on Pattern Analysis and Machine Intelligence, vol. 15, pp. 850-863, 1993.

[0132] [2] Y. Lamdan and H. J. Wolfson, “Geometric Hashing: A General and Efficient Model-Based Recognition Scheme;,” presented at Second International Confer ene on Computer Vision, Tampa, Fla., 1988.

[0133] [3] S. Sclaroff and J. Isidoro, “Active Blobs,” presented at Sixth International Conference on Computer Vision, Bombay, India, 1998.

[0134] [4] M. Turk and A. Pentland, “Eigenfaces for Recognition,” Journal of Cognitive Neuroscience, vol. 3, pp. 71-86, 1991.

[0135] [5] H. Murase and S. K. Nayar, “Visual Learning and Recognition of 3-D Objects from Appearance,” International Journal of Computer Vision, vol. 14, pp. 5-24, 1995.

[0136] [6] M. J. Swain and D. H. Ballard, “Color Indexing,” International Journal of Computer Vision, vol. 7, pp. 11-32, 1991.

[0137] [7] B. V. Funt and G. D. Finlayson, “Color Constant Color Indexing,” IEEE Transactions on Pattern Analysis and Machine Intelligence, vol. 17, pp. 522-529, 1995.

[0138] [8] D. Forsyth, J. L. Mundy, A. Zisserman, C. Coelho, A. Heller, and C. Rothwell, “Invariant Descriptors for 3-D Object Recognition and Pose,” IEEE Transactions on Pattern Analysis and Machine Intelligence, vol. 13, pp. 971-991, 1991.

[0139] [9] J. Huang, S. R. Kumar, M. Mitra, W.-J. Zhu, and R. Zabih, “Image Indexing Using Color Correlograms,” presented at IEEE Conference on Computer Vision and Pattern Recognition, San Juan, Puerto Rico, 1997.

[0140] [10] G. Pass and R. Zabih, “Histogram Refinement for Content-Based Image Retrieval,” presented at IEEE Workshop on Applications of Computer Vision, Sarasota, Fla., 1996.

[0141] [11] R. M. Haralick, K. Shanmugam, and 1. Dinstein, “Textural Features for Image Classification,” IEEE Transactions on, Systems, Man, and Cybernetics, vol. SMC-3, pp. 610-621, 1973.

[0142] [12] O. Kempthorne and L. Folks, in Probability, Statistics, and Data Analysis. Ames, Iowa: Iowa State Univeristy Press, 1971, pp. 108-121.

Claims

Post-issuance certificates present

The text shown below may not reflect later post-issuance changes.

What is claimed is:

1. A computer-implemented process for determining optimal search parameters for use in an object recognition procedure designed to find an object in a search image, said process comprising using a computer to perform the following acts:

capturing model images of the object from a plurality of viewpoints around the object;

computing a co-occurrence histogram (CH) for each model image, wherein a model image CH is computed by generating counts of every pair of pixels whose pixels exhibit a prescribed pixel characteristic that fall within the same combination of a series of pixel characteristic ranges and which are separated by a distance falling within the same one of a series of distance ranges;

generating a plurality of search windows each comprising a portion of the search image;

computing a CH for each search window, wherein a search window CH is computed by generating counts of every pair of pixels whose pixels exhibit a prescribed pixel characteristic that fall within the same combination of said series of pixel characteristic ranges and which are separated by a distance falling within the same one of said series of distance ranges;

assessing a degree of similarity between each model image CH and each of the search window CHs; and

designating a search window associated with a search window CH having a degree of similarity to one of the model image CHs which exceeds a prescribed search threshold as potentially containing the object being sought; and wherein,

said series of pixel characteristic ranges, said series of distance ranges, and the size of the search windows constitute search parameters which are optimized by ensuring these parameters are large enough to minimize the processing required to compute and assess the similarity of the model image and search window CHs, while at the same time producing an acceptable risk of a false designation that a search window potentially contains the object being sought.

2. The process of claim 1, wherein the act of optimizing the search parameters comprises the acts of:

modeling the probability that a random background region of the search image will be characterized by a particular model image CH;

deriving a numerical algorithm from said probability model that approximates the probability that a search window CH associated with the background region will match the model image CH, thereby constituting a false match, whenever a particular set of said search parameters are employed;

computing the probability of a false match for a plurality of different values of a first of said search parameters using the numerical algorithm, while keeping the values of a second and a third of said search parameters constant at a prescribed value;

determining the particular value of the first search parameter that will result in the most efficient processing of the object recognition procedure while at the same time having a false match probability that is less than a prescribed false alarm threshold; and

designating said particular value of the first search parameter as a potential optimal first search parameter.

3. The process of claim 2, wherein the act of optimizing the search parameters further comprises the acts of:

computing the probability of a false match for a plurality of different values of the second of said search parameters using the numerical algorithm, while using the particular value of the first search parameter designated as being the potential optimal first search parameter and keeping the value of the third of said search parameters constant at a prescribed value;

determining the particular value of the second search parameter that will result in the most efficient processing of the object recognition procedure while at the same time having a false match probability that is less than the prescribed false alarm threshold; and

designating said particular value of the second search parameter as a potential optimal second search parameter.

4. The process of claim 3, wherein the act of optimizing the search parameters further comprises the acts of:

computing the probability of a false match for a plurality of different values of the third of said search parameters using the numerical algorithm, while employing the particular values of the first and second search parameters designated as being the potential optimal first search parameter and potential optimal second search parameter, respectively;

determining the particular value of the third search parameter that will result in the most efficient processing of the object recognition procedure while at the same time having a false match probability that is less than the prescribed false alarm threshold; and

designating said particular value of the third search parameter as a potential optimal third search parameter.

5. The process of claim 4, further comprising the act of designating the potential optimal first, second and third search parameters as actual optimal parameters.

6. The process of claim 4, further comprising the acts of:

(a) re-computing the probability of a false match for a plurality of different values of the first of said search parameters using the numerical algorithm, while employing the particular values of the second and third search parameters designated as being the potential optimal second and third search parameters, respectively;

(b) re-determining the particular value of the first search parameter that will result in the most efficient processing of the object recognition procedure while at the same time having a false match probability that is less than the prescribed false alarm threshold;

(c) designating said re-determined value of the first search parameter as the current potential optimal third search parameter;

(d) re-computing the probability of a false match for a plurality of different values of the second of said search parameters using the numerical algorithm, while employing current values of the first and third search parameters designated as being the potential optimal first and third search parameters, respectively;

(e) re-determining the particular value of the second search parameter that will result in the most efficient processing of the object recognition procedure while at the same time having a false match probability that is less than the prescribed false alarm threshold;

(f) designating said re-determined value of the second search parameter as the current potential optimal third search parameter;

(g) re-computing the probability of a false match for a plurality of different values of the third of said search parameters using the numerical algorithm, while employing current values of the first and second search parameters designated as being the potential optimal first and second search parameters, respectively;

(h) re-determining the particular value of the third search parameter that will result in the most efficient processing of the object recognition procedure while at the same time having a false match probability that is less than the prescribed false alarm threshold; and

(i) designating said re-determined value of the third search parameter as the current potential optimal third search parameter;

(j) repeating acts (a) through (i) until a change in the designated current potential optimal first, second and third search parameter values has not varied substantially from the immediately preceding designated values for these search parameters; and

(k) declaring the last designated values of the first, second and third search parameters as actual optimal first, second and third search parameters, respectively.

7. The process of claim 1, wherein the act of optimizing the search parameters comprises the acts of:

modeling the probability that a random background region of the search image will be characterized by a particular model image CH;

computing the probability of a false rmatch for a substantial number of combinations of candidate search parameters using the numerical algorithm;

determining the particular combination of search parameters that will result in the most efficient processing of the object recognition procedure while at the same time having a false match probability that is less than a prescribed false alarm threshold; and

designating said particular combination of search parameters as an optimal set of search parameters.

8. The process of claim 1, wherein the prescribed pixel characteristic is the pixel color.

9. The process of claim 1, wherein the prescribed false alarm threshold is chosen to be as large as possible, while still providing an acceptable degree of risk of a false match.

10. A system for determining optimal search parameters for use in an object recognition procedure designed to find an object in a search image, comprising:

a general purpose computing device;

a computer program comprising program modules executable by the computing device, wherein the computing device is directed by the program modules of the computer program to,

capture model images of the object from a plurality of viewpoints around the object;

compute a co-occurrence histogram (CH) for each model image, wherein a model image CH is computed by generating counts of every pair of pixels whose pixels exhibit a prescribed pixel characteristic that fall within the same combination of a series of pixel characteristic ranges and which are separated by a distance falling within the same one of a series of distance ranges;

generate a plurality of search windows each comprising a portion of the search image;

compute a CH for each search window, wherein a search window CH is computed by generating counts of every pair of pixels whose pixels exhibit a prescribed pixel characteristic that fall within the same combination of said series of pixel characteristic ranges and which are separated by a distance falling within the same one of said series of distance ranges;

assess a degree of similarity between each model image CH and each of the search window CHs; and

designate a search window associated with a search window CH having a degree of similarity to one of the model image CHs which exceeds a prescribed search threshold as potentially containing the object being sought; and wherein,

said series of pixel characteristic ranges, said series of distance ranges, and the size of the search windows constitute search parameters which are optimized by ensuring these parameters are large enough to minimize the processing required to compute arid assess the similarity of the model image and search window CHs, while at the same time producing an acceptable risk of a false designation that a search window potentially contains the object being sought.

11. The system of claim 10, wherein thea program module for optimizing the search parameters comprises sub-modules for:

modeling the probability that a random background region of the search image will be characterized by a particular model image CH;

designating said particular value of the first search parameter as a potential optimal first search parameter.

12. The system of claim 11, wherein the program module for optimizing the search parameters further comprises sub-modules for:

designating said particular value of the second search parameter as a potential optimal second search parameter.

13. The system of claim 12, wherein the program module for optimizing the search parameters further comprises sub-modules for:

designating said particular value of the third search parameter as a potential optimal third search parameter.

14. The system of claim 13, further comprising a program module for designating the potential optimal first, second and third search parameters as actual optimal parameters.

15. The system of claim 13, further comprising program modules for:

(a) re-computing the probability of a false match for a plurality of different values of the first of said search parameters using the numerical algorithm, while employing the particular values of:the second and third search parameters designated as being the potential optimal second and third search parameters, respectively;

(c) designating said re-determined value of the first search parameter as the current potential optimal third search parameter;

(d) re-computing the probability of; a false match for a plurality of different values of the second of said search parameters using the numerical algorithm, while employing current values of the first and third search parameters designated as being the potential optimal first and third search parameters, respectively;

(f) designating said re-determined value of the second search parameter as the current potential optimal third search parameter;

(i) designating said re-determined value of the third search parameter as the current potential optimal third search parameter;

(k) declaring the last designated values of the first, second and third search parameters as actual optimal first, second and third search parameters, respectively.

16. The system of claim 10, wherein the program module for optimizing the search parameters comprises sub-modules for:

modeling the probability that a random background region of the search image will be characterized by a particular model image CH;

computing the probability of a false match for a substantial number of combinations of candidate search parameters using the numerical algorithm;

designating said particular combination of search parameters as an optimal set of search parameters.

17. The system of claim 10, wherein the prescribed pixel characteristic is the pixel color.

18. The system of claim 10, wherein the prescribed false alarm threshold is chosen to be as large as possible, while still providing an acceptable degree of risk of a false match.

19. A computer-readable memory for determining optimal search parameters for use in an object recognition procedure designed to find an object in a search image, comprising:

a general purpose computing device;

a computer program comprising program modules executable by the computing device, wherein the computing device is directed by the program modules of the computer program to,

capture model images of the object from a plurality of viewpoints around the object;

generate a plurality of search windows each comprising a portion of the search image;

assess a degree of similarity between each model image CH and each of the search window CHs; and

20. The computer-readable memory of claim 19, wherein the program module for optimizing the search parameters comprises sub-modules for:

modeling the probability that a random background region of the search image will be characterized by a particular model image CH;

designating said particular value of the first search parameter as a potential optimal first search parameter.

21. The computer-readable memory of claim 20, wherein the program module for optimizing the search parameters further comprises sub-modules for: computing the probability of a false match for a plurality of different values of the second of said search parameters using the numerical algorithm, while using the particular value of the first search parameter designated as being the potential optimal first search parameter and keeping the value of the third of said search parameters constant at a prescribed value;

designating said particular value of the second search parameter as a potential optimal second search parameter.

22. The computer-readable memory of claim 21, wherein the program module for optimizing the search parameters further comprises sub-modules for:

designating said particular value of the third search parameter as a potential optimal third search parameter.

23. The computer-readable memory of claim 22, further comprising a program module for designating the potential optimal first, second and third search parameters as actual optimal parameters.

24. The computer-readable memory of claim 22, further comprising program modules for:

(c) designating said re-determined value of the first search parameter as the current potential optimal third search parameter;

(d) re-computing the probability of a false match for a plurality of different values of the second of said search pararmeters using the numerical algorithm, while employing current values of the first and third search parameters designated as being the potential optimal first and third search parameters, respectively;

(f) designating said re-determined value of the second search parameter as the current potential optimal third search parameter;

(i) designating said re-determined value of the third search parameter as the current potential optimal third search parameter;

(k) declaring the last designated values of the first, second and third search parameters as actual optimal first, second and third search parameters, respectively.

25. The computer-readable memory of claim 19, wherein the program module for optimizing the search parameters comprises sub-modules for:

modeling the probability that a random background region of the search image will be characterized by a particular model image CH;

computing the probability of a false match for a substantial number of combinations of candidate search parameters using the numerical algorithm;

designating said particular combination of search parameters as an optimal set of search parameters.

26. The computer-readable memory of claim 19, wherein the prescribed pixel characteristic is the pixel color.

27. The computer-readable memory of claim 19, wherein the prescribed false alarm threshold is chosen to be as large as possible, while still providing an acceptable degree of risk of a false match.

Patent Citations (3)

Patent	Date	Inventor	Cited By
US5612928(A)	1997-03-01	Haley et al.	Examiner
US5671290(A)	1997-09-01	Vaidyanathan	Examiner
US5982944(A)	1999-11-01	Vaidyanathan et al.	Examiner

Non-Patent Literature (14)

D. P. Huttenlocher, G. A. Klanderman, and W. J. Ricklidge, “Comparing Images Using the Hausdorff Distance,” IEEE Transactions on Pattern Analysis and Machine Intelligence, vol. 15, pp. 850-863, 1993.Applicant
Y. Lamdan and H. J. Wolfson, “Geometric Hashing: A General and Efficient Model-Based Recognition Scheme,” presented at Second International Conference on Computer Vision, Tampa, Florida, 1988.Applicant
S. Sclaroff and J. Isidoro, “Active Blobs,” presented at Sixth International Conference on Computer Vision, Bombay, India, 1998.Applicant
M. Turk and A. Pentland, “Eigenfaces for Recognition,” Journal of Cognitive Neuroscience, vol. 3, pp. 71-86, 1991.Applicant
H. Murase and S. K. Nayar, “Visual Learning and Recognition of 3-D Objects from Appearance,” International Journal of Computer Vision, vol. 14, pp. 5-24, 1995.Applicant
M. J. Swain and D. H. Ballard, “Color Indexing,” International Journal of Computer Vision, vol. 7, pp. 11-32, 1991.Applicant
B. V. Funt and G. D. Finlayson, “Color Constant Color Indexing,” IEEE Transactions on Pattern Analysis and Machine Intelligence, vol. 17, pp. 522-529, 1995.Applicant
D. Forsyth, J. L. Mundy, A. Zisserman, C. Coelho, A Heller, and C. Rothwell, “Invariant Descriptors for 3-D Object Recognition and Pose,” IEEE Transactions on Pattern Analysis and Machine Intelligence, vol. 13, pp. 971-991, 1991.Applicant
J. Huang, S. R. Kumar, M. Mitra, W.-J. Zhu. and R. Zabih, “Image Indexing Using Color Correlograms,” presented at IEEE Conference on Computer Vision and Pattern Recognition, San Juan, Puerto Rico, 1997.Applicant
G. Pass and R. Zabih, “Histogram Refinement for Content-Based Image Retrieval,” presented at IEEE Workshop on Applications of Computer Vision, Sarasota, Florida, 1996.Applicant
R. M. Haralick, K. Shanmugam and I. Dinstein, “Textural Features for Image Classification,” IEEE Transactions on Systems, Man. and Cybernetics, vol. SMC-3, pp. 610-621, 1973.Applicant
O. Kempthorne and L. Folks, in Probability, Statistics, and Data Analysis, Ames, Iowa: Iowa State University Press, 1971, pp. 108-121.Applicant
Ettelt, E.; Schmidt, G., Optimized template trees for appearance based object recognition, Oct. 11-14, 1998, pp. 4536-4541 vol. 5, ISBN: 0-7803-4778-1, IEEE Catalog No.: 98CH36218.Applicant
Ennesser, F; Medioni, G., Finding Waldo, or focus of attention using local color information, pp. 805-809, vol. 17, Issue: 8, ISSN: 0162-8828.Examiner

Abstract

Description

Background

[0002] 1. Technical Field

[0004] 2. Background Art

Summary

f (n_{11}^{k}, n_{12}^{k}, K, n_{ij}^{k}, K, n_{n_{c}, n_{c}}^{k}) = (\begin{matrix} n_{k} \\ n_{11}^{k}, n_{12}^{k}, K, n_{ij}^{k}, K, n_{n_{c}, n_{c}}^{k} \end{matrix}) * p_{11}^{n_{11}^{k}} p_{12}^{n_{12}^{k}} K p_{ij}^{n_{ij}^{k}} K p_{n_{c} n_{c}}^{n_{n_{c}, n_{c}}^{k}}

[0020] where

(\begin{matrix} n_{k} \\ n_{11}^{k}, n_{12}^{k}, K, n_{ij}^{k}, K, n_{n_{c}, n_{c}}^{k} \end{matrix}) = \frac{n_{k}!}{n_{11}^{k}! n_{12}^{k}! K n_{ij}^{k}! K n_{n_{c}, n_{c}}^{k}!}

P (CH) = \prod_{k = 1}^{n_{d}} f (n_{c}, n_{k})

f (n_{c}, n_{k}) \approx g_{k} (n) = \frac{1}{{(2 π)}^{m / 2} {| \sum |}^{1 / 2}} \exp [- \frac{1}{2} {(n - μ)}^{T} \sum^{- 1} (n - μ)]

[0023] where

[0024] m=n_c(n_c+1)/2−1

[0025] n=(n₁₁^k,n₁₂^k,K,n_ij^k,K,n_n_−1,n^k)^Ti≦j

[0026] μ=n_k(p₁₁,p₁₂,K,p_ij,K,p_n_−1,n)^Ti≦j

[0027] and where the inverse covariance is

\sum^{- 1} = {a_{rs}} = {\begin{matrix} \frac{1}{n_{k}} (\frac{1}{q_{r}} + \frac{1}{q^{*}}) & if r = s \\ \frac{1}{n_{k} q^{*}} & if r \neq s \end{matrix} (r, s) \in [1, 2, K, m]

P (CH) = \prod_{k = 1}^{n_{d}} f (n_{c}, n_{k})

[0029] that the co-occurrences at different distances are independent, the probability of a false alarm can be approximated as

P (false alarm) \approx \prod_{k = 1}^{n_{d}} \int_{Ω_{k}} g_{k} (n)