Enhanced background model employing object classification for improved background-foreground segmentation

[0001] The present invention relates to background-foreground segmentation performed by computer systems, and more particularly, to the generation and maintenance of improved background models for use in background-foreground segmentation.

[0002] Background-foreground segmentation is a well-known computer vision based technique for detecting objects in the field of view of a stationary camera. Initially, a system learns a scene during a training phase while no objects are present. A background model of the scene is built during the training phase using a sequence of images captured from the scene. Thereafter, during normal operation, new images are compared with the background model. Pixel positions with significant deviation from the background model are classified as foreground pixels, while the remaining pixels are labeled as background pixels. The output of the algorithm is generally a binary image depicting the silhouette of the foreground objects found in the scene.

[0003] Conventional background-foreground segmentation techniques perform well for segmenting and tracking people and other objects in open outdoor areas, such as a parking lot, or enclosed, spacious facilities, such as warehouses, office spaces, or subway platforms. These scenes, however, are quite different from those of a typical home. For example, a residential environment typically contains many objects in a small area. In addition, many objects in a residential environment are non-rigid, such as garments and curtains, or deformable, such as furniture and blinds (or both), and people tend to frequently vary their pose in a residential environment, such as between a standing, sitting and laying down position.

[0004] Most existing background-foreground segmentation techniques do not perform well in the presence of occlusions of lower body parts in cluttered environments, non-upright body poses, and spontaneous movement of large background objects such as doors, chairs, and tables. A need therefore exists for a method and apparatus for generating and maintaining improved background models for use in background-foreground segmentation.

[0005] Generally, a method and apparatus are disclosed for generating and maintaining improved background models for use in background-foreground segmentation. According to one aspect of the invention, enhanced background models are employed to improve the segmentation of the foreground and background. The background model is modified in accordance with the present invention to contain an indication of objects that are typically stationary. Thereafter, if an object moves and has been previously identified as an object that is typically stationary, the object will not unnecessarily be identified as part of the foreground during background-foreground segmentation.

[0006] In an exemplary implementation, moving objects are classified into two sets. A first set includes objects that typically move independently, such as humans and animals, and a second set includes objects that are typically stationary, such as furniture. Generally, once an object is assigned to the second (stationary object) set, the object will remain in the background, even if the object is moved (normally, movement of the object would cause the object to become part of the foreground).

[0007] The enhanced background models of the present invention can be employed to improve background-foreground segmentation. During a training procedure, the parameters of the background model is estimated. When movement of an object is detected, the object is classified into one of the two sets. During normal operation, the likelihood probability of captured images is estimated using the background models. Pixels with low likelihood probability of belonging to the background scene that have not previously been flagged as stationary objects are classified as foreground, while the rest are labeled as background.

[0008] A more complete understanding of the present invention, as well as further features and advantages of the present invention, will be obtained by reference to the following detailed description and drawings.

[0013] Referring now to FIG. 1, a video processing system 120 is shown that generates and maintains a background model in accordance with preferred embodiments of the present invention for background-foreground segmentation. Video processing system 120 is shown interoperating with one or more cameras, such as camera 105, through video feed 107, a Digital Versatile Disk (DVD) 110 and a network 115. Video processing system 120 comprises a processor 130, a medium interface 135, a network interface 140, and a memory 145. Memory 145 comprises image grabber 150, an input image 155, a background model generation/maintenance process 200/300, a background model object classification process 400, probability tables 165, a global threshold 180, a pixel threshold 195, and a segmented image 190.

[0014] In the exemplary embodiment, the background model is embodied as a set of probability tables 165, discussed below, that comprise a plurality of probability tables 170-11 through 170-HW. One probability table 170-11 is shown in FIG. 1 comprising entries 175-11 through 175-NM. A background model is generally built for each camera to facilitate the fast segmentation of foreground and background. In order to prevent the need for background models to detect every object (for treatment as a foreground object), the background model is modified in accordance with the present invention to contain an indication of objects that are typically stationary. Thereafter, if an object moves and has been previously identified as an object that is typically stationary, the object will not unnecessarily become part of the foreground during background-foreground segmentation.

[0015] In one implementation, moving objects are classified into one of two sets. A first set includes objects that typically move independently, such as humans and animals. When an object from the first set moves in an image, the object should become part of the foreground. The second set includes objects that are typically stationary, such as furniture. Generally, once an object is assigned to the second set, the object will remain in the background, even if the object is moved (normally, movement of the object would cause the object to become part of the foreground). In this manner, background objects do not unnecessarily become part of the foreground when they are moved.

[0016] As shown in FIG. 1, video processing system 120 couples video feed 107 from camera 105 to image grabber 150. Image grabber 150 “grabs” a single image from the video feed 107 and creates input image 155, which will generally be a number of pixels. Illustratively, input image 155 comprises H pixels in height and W pixels in width, each pixel having 8 bits for each of red, green, and blue (RGB) information, for a total of 24 bits of RGB pixel data. Other systems may be used to represent an image, but RGB is commonly used.

[0017] The background-foreground segmentation process 200, 300 generates the background model and performs background-foreground segmentation. Background-foreground segmentation process 200 is used during normal operation of video processing system 120, while background-foreground segmentation process 300 is used during training to generate the background model. It is expected that one single process will perform both processes 200 and 300, and that the single process will simply be configured into either normal operation mode or training mode. However, separate processes may be used, if desired.

[0018] During normal operation of video processing system 120, the background-foreground segmentation process 200 uses the probability tables 165 to determine likelihood probabilities for each of the H×W pixels in input image 155. Each of the likelihood probabilities is compared with the pixel threshold 195. If the likelihood probability is below pixel threshold 195, then the pixel is assumed to belong to the background. It is also possible to modify probability models used by the background-foreground segmentation process 200 to allow video processing system 120 to assume that a pixel belongs to the background if the likelihood probability for the pixel is greater than the pixel threshold 195. It is even possible for the video processing system 120 to assign a pixel to the background if the likelihood probability for the pixel is within a range of pixel thresholds. However, it will be assumed herein, for simplicity, that a pixel is assumed to belong to the background if the likelihood probability is below a pixel threshold 195.

[0019] During normal operation, the background-foreground segmentation process 200 determines the segmented image 190 from the input image by using the probability tables 165 and the pixel threshold 195. Additionally, probability models (not shown) are used by the background-foreground segmentation process 200 to determine the likelihood probability for each pixel. Preferred probability models are discussed below in detail. These probability models are “built into” the background-foreground segmentation process 200 (and 300) in the sense that the background-foreground segmentation process 200 performs a series of steps in accordance with the models. In other words, the background-foreground segmentation process 200 has its steps defined, at least in part, by a probability model or models. For the sake of simplicity, the probability model used to perform the background-foreground segmentation and the background-foreground segmentation process will be considered interchangeable. However, it should be noted that the background-foreground segmentation process, while performing the steps necessary to determine probabilities according to a model, may have additional steps not related to determining probabilities according to a model. For instance, retrieving data from input image 155 and storing such data in a data structure is one potential step that is not performed according to a probability model.

[0020] During training, the background-foreground segmentation process 300 defines and refines the probability tables 170-11 through 170-HW (collectively, “probability tables 170” herein). Preferably, there is one probability table for each pixel of input image 155. Each probability table will have an M×N matrix, illustrated for probability table 170-11 as entries 175-11 through 175-NM (collectively, “entries 175” herein). There will be M global states and N Gaussian modes for each pixel. Generally, each probability table 170 will start with one global state and one Gaussian mode and, after training, contain the M×N entries 175.

[0021] During training, global threshold 180 is used by the background-foreground segmentation process 300 to determine whether a state should be added or parameters of a selected state modified. The pixel threshold 195 is used, during training, to determine whether another Gaussian mode should be added or whether parameters of a selected Gaussian mode should be adjusted.

[0022] It should be noted that the exemplary background-foreground segmentation process 300 allows training to be incremental. The exemplary background-foreground segmentation process 300 allows parameters of the model to be adjusted every time an image is passed to the model or after a predetermined number of images have been passed to the model. The former is preferred although the latter is possible.

[0023] As is known in the art, the methods and apparatus discussed herein may be distributed as an article of manufacture that itself comprises a computer-readable medium having computer-readable code means embodied thereon. The computer-readable program code means is operable, in conjunction with a computer system such as video processing system 120, to carry out all or some of the steps to perform the methods or create the apparatuses discussed herein. The computer-readable medium may be a recordable medium (e.g., floppy disks, hard drives, compact disks such as DVD 110 accessed through medium interface 135, or memory cards) or may be a transmission medium (e.g., a network 115 comprising fiber-optics, the world-wide web, cables, or a wireless channel using time-division multiple access, code-division multiple access, or other radio-frequency channel). Any medium known or developed that can store information suitable for use with a computer system may be used. The computer-readable code means is any mechanism for allowing a computer to read instructions and data, such as magnetic variations on a magnetic medium or height variations on the surface of a compact disk, such as DVD 110.

[0024] Memory 145 will configure the processor 130 to implement the methods, steps, and functions disclosed herein. The memory 145 could be distributed or local and the processor 130 could be distributed or singular. The memory 145 could be implemented as an electrical, magnetic or optical memory, or any combination of these or other types of storage devices. The term “memory” should be construed broadly enough to encompass any information able to be read from or written to an address in the addressable space accessed by processor 130. With this definition, information on a network, such as network 115 accessed through network interface 140, is still within memory 145 of the video processing system 120 because the processor 130 can retrieve the information from the network. It should also be noted that all or portions of video processing system 120 may be made into an integrated circuit or other similar device, such as a programmable logic circuit.

[0025] Now that a system has been discussed, probability models will be discussed that can provide global and local pixel dependencies and incremental training.

[0026] Probability Models

[0027] In a preferred probabilistic framework, images (i.e., two-dimensional array of pixel appearances) are interpreted as samples drawn from a high-dimensional random process. In this process, the number of pixels of the image defines the number of dimensions. More formally, let I={I_x,yεΘ^WH} represent an image of W×H pixels with values in the observation space Θ (i.e., RGB values at 24 bits per pixel).

[0028] The probability distributions associated with that random process, P(I|Ω), would capture the underlying image-generating process associated with both the scene and the imaging system. This includes the colors and textures present in the scene as well as the various sources of image variations such as motion in the scene, light changes, auto-gain control of the camera, and other image variations.

[0029] Most conventional algorithms model the images of a scene assuming each of the pixels as independent from each other. In practice, the image-formation processes and the physical characteristics of typical scenes impose a number of constraints that make the pixels very much inter-dependant in both the global sense (i.e., the whole image or a series of images) as well as in the local sense (i.e., regions within the image).

[0030] The proposed exemplary model exploits the dependency among the pixels within the images of a scene by introducing a hidden process ξ that captures the global state of the observation of the scene. For example, in the case of a scene with several possible illumination settings, a discrete variable ξ could represent a pointer to a finite number of possible illumination states.

[0031] A basic idea behind the proposed model is to separate the model term that captures the dependency among the pixels in the image from the one that captures the appearances of each of the pixels so that the problem becomes more tractable. That is, it is beneficial to compute the likelihood probability of the image from:

[0032] $\begin{array}{cc}P\left(I|\Omega \right)=\sum _{\forall \xi }^{}P\left(I|\xi ,\Omega \right)P\left(\xi |\Omega \right),& \left\{1\right\}\end{array}$
where P(ξ|Ω) represents the probability of the global state of the scene, and P(I|ξ, Ω) represents the likelihood probability of the pixel appearances conditioned to the global state of the scene ξ. Note that as the dependency among the pixels is captured by the first term, it is reasonable to assume that, conditioned to the global state of the scene ξ, the pixels of the image I are independent from each other. Therefore, Equation {1} can be re-written as:

[0033] $\begin{array}{cc}P\left(I|\Omega \right)=\sum _{\forall \xi }^{}P\left(\xi |\Omega \right)\prod _{\forall \left(x,y\right)}^{}P\left(I_{x,y}|\xi ,\Omega \right),& \left\{2\right\}\end{array}$
where P(I_x,y|ξ, Ω) represents the probability used to model the (x, y) pixel of the image I.

[0034] Depending upon the complexity of the model used to capture the global state of the observation of a scene, namely P(ξ|Ω), the implemented process would be able to handle different types of imaging variations present in the various application scenarios. For example, it is feasible to implement a background-foreground segmentation process robust to the changes due to the auto-gain control of a camera, if a parameterized representation of the gain function is used in the representation of ξ.

[0035] In the interest of simplicity, each of the pixel values conditioned to a global state ξ, P(I_x,y|ξ, Ω), is modeled using a mixture-of-Gaussian distribution with full covariance matrix in the three-dimensional RGB-color space. More formally, one can use the following equation:

[0036] $P\left(I_{x,y}|\xi ,\Omega \right)=\sum _{\forall \alpha }^{}P\left({\alpha }_{x,y}\right)N\left(I;{\stackrel{\_}{I}}_{\alpha ,x,y},{\sum }_{\alpha ,x,y}\right),$
where Ī_α,x,y and Σ_α,x,y are the mean value and the covariance matrix of the α-th mixture-of-Gaussian mode for the (x, y) pixel. These parameters are a subset of the symbolic parameter variable Ω used to represent to whole image model.

[0037] Note that previous research has shown that other color spaces are preferable to deal with issues such as shadows, and this research may be used herein if desired. However, the present description will emphasize modeling the global state of the scene.

[0038] The global state of the observation of a scene is preferably modeled using a discrete variable ξ={1, 2, . . . , M} that captures global and local changes in the scene, so that Equation {2} becomes the following:

[0039] $\begin{array}{cc}P\left(I|\Omega \right)=\sum _{\forall \xi }^{}P\left({\alpha }_{x,y}\right)\prod _{\forall \left(x,y\right)}^{}\sum _{\forall \alpha }^{}P\left({\alpha }_{x,y}\right)N\left(I;{\stackrel{\_}{I}}_{\alpha ,x,y},{\sum }_{\alpha ,x,y}\right).& \left\{3\right\}\end{array}$

[0040] Note the difference between the described model and the traditional mixture of Gaussians. The model of the present invention uses a collection of Gaussian distributions to model each pixel in connection to a global state, as opposed to a mixture-of-Gaussian distribution that models each of the pixels independently.

[0041] Equation 3 can be re-written as the following:

[0042] $\begin{array}{cc}P\left(I|\Omega \right)=\sum _{\forall \xi }^{}\prod _{\forall \left(x,y\right)}^{}\sum _{\forall \alpha }^{}G\left(\xi ,{\alpha }_{x,y}\right)N\left(I;{\stackrel{\_}{I}}_{\alpha ,x,y},{\sum }_{\alpha ,x,y}\right),& \left\{4\right\}\end{array}$
where the term

[0043] $G\left(\xi ,{\alpha }_{x,y}\right)={P\left(\xi |\Omega \right)}^{\frac{1}{\mathrm{WH}}}P\left({\alpha }_{x,y}\right)$
can be simply treated as M×N matrixes associated with each of the pixel positions of the image model. In this example, M is the number of global states, and N is the number of Gaussian modes. In the example of FIG. 1, the M×N matrixes are stored in probability tables 165, where there is one M×N matrix 170 for each pixel.

[0044] Segmentation Procedure

[0045] Assuming that one of the proposed models, discussed above, has been successfully trained from a set of image observations from a scene, the segmentation procedure of a newly observed image is simply based on maximum likelihood classification. Training is discussed in the next section.

[0046] An exemplary segmentation procedure is shown as method 200 of FIG. 2. Method 200 is used by a system during normal operation to perform background-foreground segmentation. As noted above, training has already been performed.

[0047] Method 200 begins in step 210 when an image is retrieved. Generally, each image is stored with 24 bits for each pixel of the image, the 24 bits corresponding to RGB values. As described above, other systems may be used, but exemplary method 200 assumes RGB values are being used.

[0048] Given the test image, I^t, the segmentation algorithm determines (step 220) the global state ξ* that maximizes the likelihood probability of the image given the following model:

[0049] $\begin{array}{cc}{\xi }^{*}=\arg \max P\left(\xi |\Omega \right)\prod _{\forall \left(x,y\right)}^{}P\left(I_{x,y}^t|\xi ,\Omega \right).& \left\{5\right\}\end{array}$

[0050] Then, the background-foreground segmentation is performed on each pixel independently using the individual likelihood probability, but only considering the most likely global state ξ*. To perform this step, a pixel is selected in step 230. The individual likelihood probability for each pixel is determined for the most likely global state (step 240), and the likelihood probability is used in the following equation to determine whether each pixel should be assigned to the background or foreground (step 250):

[0051] $\begin{array}{cc}{s}_{x,y}=\{\begin{array}{cc}1& P\left({\xi }^{*}|\Omega \right)P\left(I_{x,y}^t|{\xi }^{*},\Omega \right)<P_{\mathrm{th}}\\ 0& \mathrm{otherwise}\end{array}\forall \left(x,y\right),& \left\{6\right\}\end{array}$
where s={s_x,y∀(x, y)} represents a binary image of the background-foreground segmentation, where non-zero pixels indicate foreground objects. Basically, Equation {6} states that, if the likelihood probability for a pixel is less than a pixel threshold (step 250=YES), then a further test is performed during step 255 to determine if the pixel has been previously flagged as part of the background by the background model object classification process 400 in accordance with the present invention, otherwise (step 250=NO) the pixel is assigned to the background (step 270).

[0052] If it is determined during step 255 that the pixel has been previously flagged as part of the background by the background model object classification process 400, then the pixel is assigned to the background during step 270. If, however, it is determined during step 255 that the pixel has not been previously flagged as part of the background by the background model object classification process 400, then the pixel is assigned to the foreground (step 260).

[0053] Equation {6} is performed for each pixel of interest (generally, all the pixels in an image). Thus, in step 280, if all pixels in the image have been assigned to the background or foreground (step 280=NO), then the method 200 ends, otherwise (step 280=YES) the method continues in step 230 and Equation 6 is performed for a newly selected pixel.

[0054] Note how it is possible for process 200 to successfully classify a pixel as foreground even in the case that its color value is also modeled as part of the background under a different global state. For example, if a person wearing a red shirt walks by in the back of the scene during the training procedure, the red color would be captured by one of the mixture-of-Gaussian modes in all the pixels hit by that person's shirt. Later during testing, if that person walks again in the back of the scene (of course, roughly following the same path) he or she will not be detected as foreground. However, if that person comes close to the camera, effectively changing the global state of the scene, his or her red shirt would be properly segmented even in the image regions in which red has been associated with the background.

[0055] As an additional example, consider the case in which a part of the background looks (i) black under dark illumination in the scene, and (ii) dark green when the scene is properly illuminated. The models of the present invention, which exploit the overall dependency among pixels, would be able to detect black objects of the background when the scene is illuminated, as well as green foreground objects when the scene is dark. In traditional models, both black and green would have been taken as background colors, so that those objects would have been missed completely.

[0056] Training Procedure

[0057] Offline training of the proposed models with a given set of image samples (e.g., a video segment) is straightforward using the Expectation-Maximization (EM) algorithm. For example, the parameters of the individual pixel models,

[0058] $P\left(I_{x,y}^t|{\xi }^{*},\Omega \right),$
could be initialized randomly around the mean of the observed training data, while the probability of the individual states could be initialized uniformly. Then, using EM cycles, all the parameters of the model would be updated to a local-maximum solution, which typically is a good solution. The EM algorithm is a well known algorithm and is described, for instance, in A. Dempster, N. Laird, and D. Rubin, “Maximum Likelihood From Incomplete Data via the EM Algorithm,” J. Roy. Statist. Soc. B 39:1–38 (1977), the disclosure of which is hereby incorporated by reference.

[0059] The exemplary training procedure described in FIG. 3 incrementally trains the models and automatically determines the appropriate number of global states. Incremental training of the models allows the processes to run continuously over long periods of time, in order to capture a complete set of training samples that include all the various image variations of the modeled scene. The automatic determination of the number of global states minimizes the size of the model, which, in turn, reduces the memory requirements of the process and speeds up the background-foreground segmentation procedure.

[0060] An exemplary training process is shown in FIG. 3. This exemplary training process comprises an incremental procedure in which an unlimited number of training samples can be passed to the model. Every time a new sample image is passed to the model (i.e., a new image I^t passed to the model in step 305), the process 300 first executes an expectation step (E-step, from the EM algorithm) determining the most likely global state ξ* (step 310) and the most likely mixture-of-Gaussian mode, α_x,y, of each pixel of the image (step 315). Note that these steps are similar to steps in the segmentation procedure process 200.

[0061] In step 320, the likelihood probability of the same image for the selected state is determined. Then, depending upon the value of the likelihood probability of the sample image for the selected global state (step 325), the process 300 selects between adjusting the parameters of the selected state (step 335) or adding a new one (step 330). If the likelihood probability of the sample image for the selected state is greater than a global threshold (step 325=YES), then the parameters of the selected state are adjusted (step 335). If the likelihood probability of the sample image for the selected state is less than or equal to a global threshold (step 325=NO), then a new state is added (step 330).

[0062] In step 340, the individual likelihood probabilities of the selected mixture-of-Gaussian modes for each pixel position are determined. Then, depending upon the individual likelihood probabilities of the selected mixture-of-Gaussian modes for each pixel position, the algorithm selects between adjusting the selected modes or adding new ones. To do this, in step 345, a pixel is selected. If the individual likelihood probability of the selected mixture-of-Gaussian modes for this pixel position is greater than a pixel threshold (step 350=YES), then the selected mode is adjusted (step 360), otherwise (step 350=NO) a new mode is added (step 355). If there are more pixels (step 365=YES), the method 300 continues in step 345, otherwise (step 365=NO), the method continues in step 370. If there are more sample images to process (step 370=YES), the method 300 continues in step 305, otherwise (step 370=NO) the method ends.

[0063] Note that two thresholds are used in the exemplary training method 300: one for the decision at each pixel position, and the other for the decision about the global state of the image.

[0064] Each mixture-of-Gaussian mode of every pixel position preferably keeps track of the total number of samples used to compute its parameters, so that when a new sample is added, the re-estimation of the parameters is carried out incrementally. For example, means and covariances of the mixture-of-Gaussian modes are simply updated using:

[0065] ${\stackrel{\_}{I}}_{\alpha ,x,y}=\frac{1}{\left(1+K_{\alpha ,x,y}\right)}\left\{I_{x,y}^t+K_{\alpha ,x,y}{\stackrel{\_}{I}}_{\alpha ,x,y}\right\},{\sum }_{\alpha ,x,y}=\frac{1}{K_{\alpha ,x,y}}\left[{\left(I_{x,y}^t-{\stackrel{\_}{I}}_{\alpha ,x,y}\right)}^{\prime }\left(I_{x,y}^t-{\stackrel{\_}{I}}_{\alpha ,x,y}\right)+\left(1-K_{\alpha ,x,y}\right){\sum }_{\alpha ,x,y}\right],$
where k_α,x,y is the number of samples already used for training that mixture-of-Gaussian mode.

[0066] Similarly, each global state keeps track of the total number of samples used for training, so that when a sample is added, the probability tables G(ξ, α_x,y) are updated using the usage statistics of the individual states and mixture-of-Gaussian modes, considering the addition of the new sample.

[0067] Beneficially, the overall model is initialized with only one state and one mixture-of-Gaussian mode for each pixel position. Also, a minimum of 10 samples should be required before a global state and/or a mixture-of-Gaussian mode is used in the expectation step (steps 315 and 320).

[0068] Enhancement of Background Model

[0069] As previously indicated, the background-foreground segmentation process 300 defines and refines the probability tables 170-11 through 170-HW (collectively, “probability tables 170” herein) during the training phase. Initially, the probability tables 170 comprise a probability table for each pixel of input image 155. As discussed hereinafter in conjunction with FIG. 4, the background model object classification process 400 modifies the probability tables 170 of the background model to contain an indication of objects that are typically stationary. In this manner, if an object moves and has been previously identified as an object that is typically stationary, the object will not unnecessarily become part of the foreground during background-foreground segmentation.

[0070] In the exemplary implementation of the background model object classification process 400, moving objects are classified into one of two sets. The first set includes objects that typically move independently, such as humans and animals and the second set includes objects that are typically stationary, such as furniture. Generally, once an object is assigned to the second set, the object will remain in the background, even if the object is moved (normally, movement of the object would cause the object to become part of the foreground). In this manner, background objects do not unnecessarily become part of the foreground when they are moved.

[0071] FIG. 4 is a flow chart describing an exemplary background model object classification process 400 incorporating features of the present invention. As shown in FIG. 4, the background model object classification process 400 initially groups pixels in the foreground of each image during step 405 using known connected component analysis techniques to form rigid objects, referred to as a blob.

[0072] A test is then performed during step 410 until it is determined that a background object has moved by more than a specified threshold, T. The movement of a background object can be detected, for example, using well known optical flow methods. Generally, the background is suppressed in two consecutive images, and the foregrounds are subtracted to detect movement, how far a given object has moved and the speed of movement. Aspect ratio information and multiple frame processing techniques can be employed, as described, for example, in U.S. patent application Ser. No. 09/794,443, filed Feb. 27, 2001, entitled “Classification of Objects Through Model Ensembles,” incorporated by reference herein.

[0073] Once the motion of a background object is detected during step 410, a further test is performed during step 420 to determine if the background object that has moved is typically stationary. For example, the techniques described in the above-referenced U.S. patent application Ser. No. 09/794,443, filed Feb. 27, 2001, entitled “Classification of Objects Through Model Ensembles,” can be used to classify an object as an inanimate object, for example, such as a piece of furniture, and a set of rules can be established to indicate whether certain classified objects are typically stationary.

[0074] If it is determined during step 420 that the background object that has moved is typically stationary, then a flag is set during step 430 for the pixels associated with the object indicating that the moved object is typically stationary.

[0075] If, however, it is determined during step 420 that the background object that has moved is not typically stationary, then a flag is set during step 440 for the pixels associated with the object indicating that the moved object typically moves independently. Program control then terminates during step 460 until the motion of a background object is again detected during step 410.

[0076] Additional Embodiments

[0077] It is a common practice to approximate the probability of a mixture of Gaussians with the Gaussian mode with highest probability to eliminate the need for the sum, which prevents the further simplification of the equations.

[0078] Using that approximation at both levels, (a) the sum of the mixtures for each pixel becomes the following:

[0079] $\sum _{\forall \alpha }G\left(\xi ,{\alpha }_{x,y}\right)N\left(I;{\stackrel{\_}{I}}_{\alpha ,x,y},{\sum }_{\alpha ,x,y}\right)\approx \underset{x,y}{\max }G\left(\xi ,{\alpha }_{x,y}\right)N\left(I;{\stackrel{\_}{I}}_{\alpha ,x,y},{\sum }_{\alpha ,x,y}\right),$
and (b) the sum of the various global states becomes the following:

[0080] $\sum _{\forall \alpha }P\left(I|\xi ,\Omega \right)P\left(\xi |\Omega \right)\approx \underset{\xi }{\max }P\left(I|\xi ,\Omega \right)P\left(\xi ,\Omega \right).$
Equation {4} simplifies to the following:

[0081] $\begin{array}{cc}P\left(I|\Omega \right)\approx \underset{\xi }{\max }\prod _{\forall \left(x,y\right)}^{}\underset{{\alpha }_{x,y}}{\max }G\left(\xi ,{\alpha }_{x,y}\right)N\left(I;{\stackrel{\_}{I}}_{\alpha ,x,y},{\sum }_{\alpha ,x,y}\right).& \left\{7\right\}\end{array}$
Note the double maximization. The first one, at pixel level, is used to determine the best matching Gaussian mode considering the prior of each of the global states. The second one, at image level, is used to determine the state with most likelihood probability of observation.

[0082] Another common practice to speed up the implementation of this family of algorithms is the computation of the logarithm of the probability rather than the actual probability. In that case, there is no need for the evaluation the exponential function of the Gaussian distribution, and the product of Equation {7} becomes a sum which can be implemented using fixed-point operations because of the reduced range of the logarithm.

[0083] It should be noted that the models described herein may be modified so that a test step currently written to perform one function if a probability is above a threshold may be re-written, under modified rules, so that the same test step will perform the same function if a probability is below a threshold or in a certain range of values. The test steps are merely exemplary for the particular example model being discussed. Different models may require different testing steps.

[0084] It is to be understood that the embodiments and variations shown and described herein are merely illustrative of the principles of this invention and that various modifications may be implemented by those skilled in the art without departing from the scope and spirit of the invention.