Abstract:
The invention relates to a method for detecting the proper motion of an real-world object, comprising the steps of acquiring, by an image capturing device (ICD), a first image (I 1 ) of the object at a first point in time (t 1 ) and a second image (I 2 ) at a second point in time (t 2 ); obtaining a third (hypothetical) image (I 3 ), based on an estimated effect of the motion of the image capturing device (ICD) itself (EMF) between the first and the second point in time (t 1 , t 2 ), wherein the effect of the motion of the image capturing device (ICD) itself is estimated based on the forward kinematics of the image capturing device; determining an optical flow (OF) between the second image (I 2 ) and the third image (I 3 ); and evaluating the optical flow (OF) by incorporating uncertainties of the optical flow (OF) and the ego-motion-flow (EMF) in order to determine the proper motion of the object.

Description:
The present invention relates to a method and a corresponding device for estimating the proper motion of an object (Object Proper Motion or OPM), using optical flow and kinematics and depth information, which is applicable in controlling autonomous mobile vehicles, such as robots or cars. 
     TECHNICAL BACKGROUND 
     The ability to visually perceive motion is believed to be highly beneficial for surviving in a dynamic environment. So it is not surprising to see that movement is one of the most important cues to attract visual attention. Interestingly, most mammals are able to perceive visual motion even if they are moving themselves—either by rotating the head, the eyes, moving the whole body or even while they are running. The gathered information is then used for controlling the own movement or for tracking moving objects by keeping them centred in the fovea. 
     In general, a robot&#39;s autonomous movement causes the camera to undergo translational and rotational movements in 3D with sudden velocity-changes. 
     The effects of Ego-Motion (EM) on the Optical Flow (OF), the so called ego-motion-flow (EMF), produced by a robot make this task quite challenging. Primarily, this is caused by the large number of degrees of freedom and the complex influence of each robot segment on the position of the camera: For example, bending the knee may in some situations cause the robot and hence the camera to shift and rotate violently to the side, whereas in other situations, the robot might be standing on the other leg and so knee bending does not affect the camera position at all. 
     PRIOR ART 
     The domain where the estimation of object proper motion (OPM) during ego-motion (EM) is commonly addressed is the car domain. Here, this estimation is important to measure the speed of other cars relative to the observing car which allows identifying them as eventual obstacles. Most approaches rely on the calculation of optical flow (OF). This flow is a superposition of the effects caused by ego-motion (EM), the so called ego-motion-flow (EMF), and object proper motion (OPM). A common procedure to decompose the optical flow (OF) and extract the desired object proper motion (OPM) tries to estimate the movement of the camera from the visual flow fields, where the underlying models for camera movement may highly differ in complexity. 
     For proper extraction of the camera movement it is hence crucial that the optical flow (OF) is primarily caused by ego-motion (EM)—effects and not by independently moving objects—which can not be assumed for a robot facing moving people. The effect of translational ego-motion (EM) on points in the image is highly influenced by their distance to the camera in the way that distant points induce smaller flow vectors than closer ones. 
     Some of the prior art approaches actually measure the distance of these points. However the majority of methods assume that all points lie on a plane running through the position of the camera and the focus of expansion at the horizon. While this procedure may be suitable for the car domain, no such simplifying assumption about the environment is appropriate for autonomously acting robots in dynamic scenes. Thus, different ways to handle ego-motion (EM) have to be found. 
     One idea of using existing knowledge about the movement of the camera is as follows. Letting a robot walk in a circle, one may compute the optical flow (OF) resulting from the robot&#39;s movement over a textured ground. In conjunction with the robot&#39;s gait phase and joint angles, a neural network is trained to learn the ego-motion-flow (EMF). Afterwards differences between this ego-motion-flow (EMF) and the measured optical flow (OF) can be used to detect obstacles in the path of the robot. While the movement of the robot was very constricted, Fidelman et al. (“Detecting motion in the environment with a moving quadruped robot,” in RoboCup-2006: Robot Soccer World Cup X, G. Lakemeyer, E. Sklar, D. Sorenti, and T. Takahashi, Eds. Berlin: Springer Verlag, 2007, vol. 4434/2007, pp. 219-231) proved that ego-motion-flow (EMF) can be learned even for more complex movements. In their approach the neural network is provided with the recent optical flow (OF)-calculation in addition to walk phase and joint angles of a robot. The neural net predicts the flow of the next timestep, allowing to compare it with the actually measured optical flow (OF) in order to classify differences as object proper motion (OPM). 
     The idea of calculating the movement of the camera and using it to compute the ego-motion-flow (EMF) has also been described for the car. In US patent application 2004/0183905 (Comaniciu, D. et al.), readings from gyroscopic sensors or GPS-signals are applied to estimate the camera-movement. Similar to other approaches in cars, the depth is not measured but estimated from the plane assumption as described earlier. Another approach actually computes the ego-motion-flow (EMF) depending on the measured robot joints (del Solar J. R. et al.) (“Motion detection and tracking for an Aibo robot using camera motion compensation and Kalman filtering.” in RoboCup 2004: Robot Soccer World Cup VIII, ser. Lecture Notes in Computer Science, D. Nardi, M. Riedmiller, C. Sammut, and J. Santos-Victor, Eds. Berlin: Springer, March 2005, vol. 3276, pp. 619-627) instead of using neural networks. 
     However, all the results of the provided approaches are not very promising, possibly due to a missing integration of depth information. Fidelman et al. (“Detecting motion in the world with a moving quadruped robot.” P. Fidelman, T. Coffman, R. Miikulainen, and P. Stone. Technical Report TR-05-37, Department of Computer Science, University of Texas, 2005) implicitly include depth information by providing the calculated optical flow (OF)-field, which depends on the distance of points to the camera. Nevertheless, this approach can not cope properly with objects having a distance different from objects in the training-phase: The difference in distance results in an optical flow (OF) deviating from the learned one. Hence even non-moving points can be classified as object proper motion (OPM). The combination of knowledge about the robot movement with distance information and optical flow (OF)-measurements is described in a paper by Overett, G. et al. (G. Overett and D. Austin, “Stereo vision motion detection from a moving platform,” in Proc. Australasian Conf. on Robotics and Automation, December 2004). Having cameras fixed to the robot, Overett et al. try to measure the odometry data and compute the resulting ego-motion-flow (EMF) considering the depth. This flow is afterwards subtracted from the estimated optical flow (OF) and residual vectors are used to indicate object proper motion (OPM). Unfortunately, the noisy data from odometry forces the authors to manually measure the distance passed by the robot, preventing the system from running in real time. 
     None of the presented methods precomputes a movement of the camera and searches for movement relative to it. When acquiring the effects of camera movement that is superimposed with the effects of independently moving objects, the optical flow (OF) needs a large set of displacements. The higher range of acquired motion does not only increase the computational time but also corrupts optical flow (OF)-estimations by producing increased ambiguities. Very recently Rabe et al. (“Fast detection of moving objects in complex scenarios,” in Proc. IEEE Symp. Intelligent Vehicles, June 2007, pp. 398-403) proposed a solution to this problem by computing flow and disparity from visual features tracked in the image instead of estimating dense depth- and flow-fields. Fusing the results in a Kalman filter provides 6D information about location and motion for a set of points. Together with readings from inertial sensors this information is used to compensate ego-motion (EM)-effects and obtain the actual 3D object proper motion (OPM). The results are very convincing, at least for the moderate movement of a car. However, apparently no procedure was presented which is able to handle the noisy character of depth, optical flow (OF) and the fierce effects of camera-movement induced by a legged robot, combining them to make reliable estimations of object proper motion (OPM). 
     Fardi et al. (“Motion-based pedestrian recognition from a moving vehicle,” in Proc. IEEE Symp. Intelligent Vehicles, June 2006, pp. 219-224) present a method for recognizing pedestrians from a moving vehicle comprising the computation of an optical flow, including the compensation of the ego motion of the camera. 
     In Fardi et al, the region of interest is regarded as a plane with a uniform depth z (distance to the camera) that is calculated from a distance measurement of a laser scanner and is used for all the pixels in the region of interest. However, this assumption about the depth is not applicable in all situations in which an accurate estimate of the proper motion of an object is desired. This effect becomes more significant for objects close to the camera and in front of a distant background. 
     It is therefore an object of the present invention, to provide a method and device for estimating the proper motion of objects by handling depth information more precisely, being capable to provide and take advantage of the confidence of the OPM estimation for every pixel in an image sequence. 
     SHORT SUMMARY OF THE INVENTION 
     This object is achieved according to the invention by a method and a system according to the independent claims. Advantageous embodiments are defined in the dependent claims. 
    
    
     
       SHORT DESCRIPTION OF THE FIGURES 
       More aspects and advantages of the present invention will become more apparent when the following detailed description of the invention is studied, in connection with the annexed drawing, in which 
         FIG. 1  shows an overview of a system according to an embodiment of the present invention; 
         FIG. 2   a  shows a point {right arrow over (p)} C, i   m  in camera- and world coordinates; 
         FIG. 2   b  shows the description of the static point at time n using the camera-transformations; 
         FIG. 3   a  shows the spectrum of optical flow (OF)-motion, if ego-motion (EM)-effects are not cancelled out in advance. The solid line visualizes an assumed object proper motion (OPM) of two pixels, the dashed line an assumed ego-motion (EM)-effect of three pixels. The interval which has to be acquired by the optical flow (OF) is a superposition of the two. 
         FIG. 3   b  shows the reduced spectrum for searching relative to ego-motion-flow (EMF); 
         FIG. 4   a  shows the disrupted results of pixel wise warping under extreme body movements. 
         FIG. 4   b  shows the effect of filling-in holes in the ego-motion-flow (EMF) with the averaged flow; 
         FIG. 5  shows a one-dimensional plot of the resulting distribution N {right arrow over (ν)} ({right arrow over (c)}, C) from) N ν     i   ({right arrow over (r)} i +{right arrow over (e)} i , Σ r ) and N ν     i   ({right arrow over (e)} i , Σ e ) as well as the likelihood L; 
     
    
    
     DETAILED DESCRIPTION OF THE INVENTION 
       FIG. 1  shows an overview of a system according to the embodiments of the present invention: 
     In a first step the forward kinematics are used to compute the movement of the camera occurring in the time-interval from m to n=m+1. Combining this with the depth information from binocular disparity, one may estimate where a static point in the image at time m moved due to ego-motion (EM) and obtain the ego-motion-flow (EMF). Hence, by knowing where a point in the image at time n originated from, the ego-motion (EM)-effect may be compensated by moving the point back to its original position in a second step, resulting in an image freed from EM-effects. Afterwards this image may be used for a calculation of the optical flow (OF) measured relative to the ego-motion-flow (EMF). This step-by-step movement estimation allows the reduction of the optical flow&#39;s (OFs) search-range which does not only reduce computational time but also decreases possible ambiguities in the optical flow (OF)-measurement and thereby leads to qualitatively better results. Ideally, the optical flow (OF) measured in this way should be zero for non-moving objects and otherwise describe their proper motion. 
     However, since binocular disparity as well as optical flow (OF) are particularly noisy signals, suitable methods may be applied to filter out wrong velocity-estimations based on the reliability of flow- and ego-motion (EM)-vector. 
     Computation of Ego-Motion-Flow (EMF) 
     The computation of ego-motion-flow (EMF) results in a flow field E m ={{right arrow over (e)} i   m }, which describes where the points in image I m  have shifted to in I n , caused by movement of the camera. For this computation it is assumed that the environment is static and the points did not move themselves. 
     To account for the fact that the absolute values of the vectors {right arrow over (e)} i   m  in highly depend on the distance of the corresponding points to the camera, the effect of camera movement may be calculated in a three dimensional camera-related space. Therefore, first each point (x i   m , y i   m ) may be combined in the image with its binocular disparity d i   m , to define a point in homogenous image-coordinates as
 
 {right arrow over (p)}   I,i   m =( x   i   m   ,y   i   m   ,d   i   m ,1) T  
 
     This point can then be reconstructed in camera-coordinates by computing the homogenous transformation matrix T I←C  from camera-to image-coordinates and by multiplying {right arrow over (p)} I, i   m  with it:
 
 p   C,i   m   =T   C←I   ·p   I,i   m  
 
     To compute the position of each such point relative to the camera at the next time step, the way the camera moved has to be known. This knowledge may be obtained from the robot kinematics. For this purpose, the robot&#39;s state vector {right arrow over (s)} m , which contains the position of the upper body as well as the angles of all joints at time m is used. With this vector, one may compute the forward kinematics, returning the transformations from the initial world-coordinate system to each robot segment, and including the transformation from world- to camera-coordinates T C←W   m . 
     Inverting this transformation results in T W←C   m , which describes the transformation from camera- to world-coordinates. In addition, the current transformation T C←W   n  can be computed by applying {right arrow over (s)} n . Since this calculation includes the position of the upper body, the resulting transformations are not only rotation-matrices but also encode translational changes of the camera-position. 
     Assuming a static point {right arrow over (p)} C, i   m , the new position {right arrow over (p)} C, i   n  can be calculated from transferring the point to world-coordinates by multiplying the last transformation T W←C   m  and passing it back to camera-coordinates with the actual transformation T C←W   n :
 
 {right arrow over (p)}   C,i   n   =T   C←W   n   T   W←C   m   p   C,i   m   (3)
 
     This process is illustrated in  FIG. 2 . 
     Since one is not interested in the new three dimensional position of the points but rather in the 2D-flow in the image, each point may be projected back into the image-plane. This can be done by computing a homogenous projection-matrix T I←C  and multiplying the new point {right arrow over (p)} C, i   n  with it:
 
 {right arrow over (p)}   I,i   n   =T   I←C   {right arrow over (p)}   C,i   n  
 
     Due to the constancy of T C←I , T W←C   m , and T I←C  within a timestep, the concatenation of these transformations may be pre-computed, resulting in one single transformation matrix which may be multiplied with each point {right arrow over (p)} I, i   m . 
     Knowing {right arrow over (p)} I, i   n  and {right arrow over (p)} I, i   m  the shift of a point in the image-plane due to camera-motion is given by 
     
       
         
           
             
               
                 e 
                 -&gt; 
               
               i 
               m 
             
             = 
             
               ( 
               
                 
                   
                     
                       
                         x 
                         i 
                         n 
                       
                       - 
                       
                         x 
                         i 
                         m 
                       
                     
                   
                 
                 
                   
                     
                       
                         y 
                         i 
                         n 
                       
                       - 
                       
                         y 
                         i 
                         m 
                       
                     
                   
                 
               
               ) 
             
           
         
       
     
     To express the reliability of this shift-vector, only the quality of the measured depth z i   m  of {right arrow over (p)} C, i   m  has to be considered, since the computation of the kinematics is deterministic. The reliability of {right arrow over (e)} i   m  relates to the amount of correlation found by the disparity-algorithm between the left and the right camera-image. This is expressed in terms of some confidence-value cε[0,1]. In addition the falling accuracy with increasing distance may be taken into account. 
     The sensitivity for depth estimation from disparity decreases with the squared distance: 
                 ∂   d       ∂   z       =       -       b   ·   f         (     z   m     )     2         ⁢   q           
where b denotes the baseline, ƒ the focal-length and q the pixel-size of the camera.
 
     Since one wants to describe the likelihood for {right arrow over (e)} i   m  in terms of some probability distribution, one may express a decay of reliability in terms of an increasing variance 
             ∑   C   m         
by declaring it as
 
               ∑     e   i     m     ⁢     =     1   ·         w   e     ⁡     (         (     1   -     w   d       )     ⁢     c   i       +       w   d     ⁢       ∂     d   m         ∂     z   i   m             )         -   1                 
with 1 being the unity-matrix. The scalar w e  may be used to express a weighting of the whole variance which is beneficial for the comparison of this variance with the one from the optical flow (OF) (see below). Additionally w d ε[0,1] may be utilized to stress the influence of either confidence or distance.
 
     Computation of Flow Relative to Ego-Motion-Flow (EMF) 
     Since the effects of ego-motion (EM) have been computed on the image, one is now able to compute the optical flow (OF) R m ={{right arrow over (r)} i   m } between I m  and I n  relative to this ego-motion-flow (EMF). In comparison to other approaches which cancel out ego-motion (EM)-effects after the computation of the optical flow (OF), this reduces the spectrum of the motion to be acquired. This effect which is illustrated in  FIG. 3  does not only reduce computational time but also improves the optical flow (OF)-estimations by reducing ambiguities. 
     The ego-motion (EM)-effects may be compensated by warping the images. In the present case, warping using a forward mapping can be written as
 
 Ĩ   m ( {right arrow over (p)}   I,i   m   +{right arrow over (e)}   i   m )= I   m ( p   I,i   m )
 
where Ĩ m  equals the old image I m  but is freed from the measured ego-motion (EM)-effects. This kind of forward mapping has two major drawbacks: Since {right arrow over (e)} i   m  usually encodes real-valued shifts, the data-points {right arrow over (p)} I, i   m +{right arrow over (e)} i   m  may not lie inside the grid and require complicated interpolations. Additionally, it is not guaranteed that each point in the warped image is targeted by the sum of original position and shift, leading to holes in the image.
 
     The more flexible solution lies in the usage of backward-mapping to warp the actual image I n  back to Ĩ n :
 
 Ĩ   n ( {right arrow over (p)}   I,i   m )= I   n ( {right arrow over (p)}   I,i   m   +{right arrow over (e)}   i   m )
 
     Since data points are now passed as arguments for the resulting image Ĩ n , holes in the image can not occur. The problem of real-valued shifts may be tackled by using bilinear interpolation in the source-image. Apart from deviations caused by wrong depth measurements, Ĩ n  and I m  should only differ in points with individual object motion. Nevertheless, in some situations with extreme body-movement, holes in the depth image D m  can cause artefacts in the pixelwise warped image, as shown in  FIG. 4(   a ). These artefacts lead to errors in the computed optical flow (OF) and can hence affect surrounding regions, even if those regions have valid depth estimations. In a first step, these artefacts may be reduced by filling invalid regions with suitable assumptions obtained from valid ego-motion-flow (EMF)-estimations. 
     In one embodiment of the present invention, the average ego-motion-flow (EMF) {circumflex over (ν)} E   m  gained from valid depth-estimations may be used to fill the gaps (see  FIG. 4(   b )). To account for the decreased reliability of the warped image and hence of the optical flow (OF), one may create a penalty-map P which is high for points near and inside invalid depth-regions and 0 otherwise. 
     By passing I m  and Ĩ n  to an optical flow (OF)-algorithm, a velocity estimate {right arrow over (r)} i   m  may be obtained for each point, relative to the estimated ego-motion-flow (EMF). The algorithm may also compute a covariance 
             ∑   C   m         
which gives a confidence measure for the optical flow vectors, assuming pixelwise independent Gaussian noise.
 
     For the inclusion of the warping-based penalty, the compound variance Σ r   m  at point {right arrow over (p)} I, i   m  may be calculated as 
     
       
         
           
             
               
                 
                   
                     
                       ∑ 
                       
                         r 
                         i 
                       
                       m 
                     
                     ⁢ 
                     
                         
                     
                     ⁢ 
                     
                       = 
                       
                         
                           1 
                           · 
                           
                             P 
                             ⁡ 
                             
                               ( 
                               
                                 
                                   p 
                                   -&gt; 
                                 
                                 
                                   I 
                                   , 
                                   i 
                                 
                                 m 
                               
                               ) 
                             
                           
                         
                         + 
                         
                           ∑ 
                           
                             C 
                             i 
                           
                           m 
                         
                       
                     
                   
                   ⁢ 
                   
                       
                   
                 
               
               
                 
                   ( 
                   10 
                   ) 
                 
               
             
           
         
       
     
     For invalid points the choice of 
                 1   ·       p   -&gt;       I   ,   i     m       &gt;&gt;     ∑   C   m       ⁢                 
ensures a negligible influence of the confidence-based variance. This is useful since artefacts in the warped image usually increase the structure at that specific region and thereby lead to high confidence-measures of the optical flow (OF)-algorithm.
 
     Besides the spatial integration, the proposed invention also allows prediction over time, giving it the ability to iteratively make reliable calculations of motion in unstructured image regions by drawing previous estimations into account. The ideal outcome of the flow-algorithm would be a vector-field, which is zero for non-moving objects and denotes their proper motion otherwise. However, since the optical flow (OF) as well as depth-measurements from disparity are very noisy signals, some more filtering is needed for the detection of object proper motion (OPM). 
     Detecting Object Proper Motion (OPM) 
     In the previous sections, two approaches for the calculation of image flow were introduced, which use very different methods and hence show different characteristics. Since the computation of the ego-motion-flow (EMF) is based on disparity and kinematics, it can acquire the effects of ego-motion (EM) on image-points, as long as the points are not moving. In contrast, the optical flow (OF) also works for moving points. By adding the ego-motion-flow (EMF) to the relative optical flow (OF), an overall flow may be obtained, which designates the compound retinal movement of object proper motion (OPM) and ego-motion (EM)-effects. For non-moving points, this flow should equal the ego-motion-flow (EMF), while it should be different for moving points. A measurement for the significance of this distance may be introduced and used to extract object proper motion (OPM)-vectors from the optical flow (OF). 
     This measurement is derived from a stochastic assumption about the estimated flows, the depth and the images. Therefore, some stochastic variables have to be defined for each point {right arrow over (p)} I, i  in the image a. Each vector of the overall flow mentioned above is represented by the variable θ i  while ε i  describes the ego-motion (EM)-vector for each point {right arrow over (p)} I, i . The random variables I={I m , I n } and D=D m  specify the observed source- and depth-images. 
     The principal idea for the approach is to estimate the probability of measuring the same velocity ν i  from the ego-motion-flow (EMF) and from the compound optical flow (OF), assuming the corresponding point is static. That is, if some point did not move, the velocity described by ε i  should not differ too much from θ i  and the likelihood to measure some identical velocity ν i  from both methods should be high. In contrast, a moving point results in different outcomes for ε i  and θ i  and the likelihood to measure the same velocity ν i  from the two methods is very low. Concluding, a high joint probability ρ(θ i =ν i , ε i =ν i , I, D) indicates a static point whereas a low probability indicates a moving point. 
     Reflecting the assumptions about the dependencies of ε i , θ i , D and I, this joint distribution can be decomposed into
 
ρ(θ i =ν i ε i =ν i   ,I,D )=ρ(θ i =ν i |ε i   ,I )ρ(ε i =ν i   |D )ρ( I )ρ( D )  (11)
 
     Since no prior assumptions are made about the source images and the depth, the corresponding variables are uniformly distributed and hence have no influence on the distribution. Using the precomputed results for the ego-motion-flow (EMF) and optical flow (OF), the conditional distributions may be approximated with Gaussians N({right arrow over (r)}, Σ) with mean and variance parameters, so that
 
ρ(θ i =ν i |ε i   ,I )∝ N   ν     i   ( {right arrow over (r)}   i   +{right arrow over (e)}   i ,Σ r )  (12)
 
ρ(ε i =ν i   |D )∝ N   ν     i   ( {right arrow over (e)}   i ,Σ e )  (13)
 
     Adding {right arrow over (e)} i  in (12) accounts for the warping of the image—while {right arrow over (r)} i  encodes a flow relative to the ego-motion-flow (EMF), the sum of {right arrow over (e)} i  and {right arrow over (r)} i  makes it an absolute flow and allows the comparison with the ego-motion-flow (EMF). With this approximation the joint distribution is proportional to the product of two Gaussians, being defined as)
 
 N   ν     i   ( {right arrow over (r)}   i   +{right arrow over (e)}   i ,Σ r ) N   ν     i   ( {right arrow over (e)}   i ,Σ e )= L   i   ·N   ν     i   ( {right arrow over (c)}   i   ,C   i )
 
with:  {right arrow over (c)}   i =Σ r     i   (Σ r     i   Σ e     i   ) −1   {right arrow over (e)}   i +Σ e     i   (Σ r     i   Σ e     i   ) −1 ( {right arrow over (r)}   i   +{right arrow over (e)}   i )  (14)
 
 C   i =Σ r     i   (Σ r     i   Σ e     i   ) −1 Σ e     i     (15)
 
 L   i   =N   {right arrow over (e)}     i   ( {right arrow over (r)}   i   +{right arrow over (e)}   i ,Σ r     i   +Σ e     i   )
 
     This procedure is identical for each time step, so the time indices have been dropped for convenience. 
     A visualization of this product is shown in  FIG. 5 . The mean value {right arrow over (c)} of the resulting distribution can be interpreted as that identical velocity which is most likely to be measured by both algorithms—the ego-motion-flow (EMF) and the optical flow (OF). While it is always possible to find such a value which fits best the hypothesis of a common velocity the factor L is a measure to describe how well {right arrow over (c)} actually fits in the light of the calculated displacements and variances. For the evaluation of object proper motion (OPM), one is not interested in the value of the vector {right arrow over (c)}, but rather whether such a vector is likely to occur. So the rejection of optical flow (OF)-estimations may be based on this value L, and can be further simplified by applying the logarithm to L: 
     
       
         
           
             
               
                 
                   
                     L 
                     i 
                   
                   = 
                   
                     
                       N 
                       
                         
                           e 
                           -&gt; 
                         
                         i 
                       
                     
                     ( 
                     
                       
                         
                           
                             r 
                             -&gt; 
                           
                           i 
                         
                         + 
                         
                           
                             e 
                             -&gt; 
                           
                           i 
                         
                       
                       , 
                       
                         
                           ∑ 
                           
                             r 
                             i 
                           
                         
                         ⁢ 
                         
                           + 
                           
                             ∑ 
                             
                               e 
                               i 
                             
                           
                         
                       
                     
                     ) 
                   
                 
               
               
                 
                   ( 
                   16 
                   ) 
                 
               
             
             
               
                 
                   
                       
                   
                   ⁢ 
                   
                     = 
                     
                       z 
                       · 
                       
                         ⅇ 
                         
                           
                             - 
                             
                               1 
                               2 
                             
                           
                           ⁢ 
                           
                             
                               ( 
                               
                                 
                                   
                                     e 
                                     -&gt; 
                                   
                                   i 
                                 
                                 - 
                                 
                                   ( 
                                   
                                     
                                       
                                         r 
                                         -&gt; 
                                       
                                       i 
                                     
                                     + 
                                     
                                       
                                         e 
                                         -&gt; 
                                       
                                       i 
                                     
                                   
                                   ) 
                                 
                               
                               ) 
                             
                             T 
                           
                           ⁢ 
                           
                             
                               ( 
                               
                                 
                                   ∑ 
                                   
                                     r 
                                     i 
                                   
                                 
                                 ⁢ 
                                 
                                   + 
                                   
                                     ∑ 
                                     
                                       e 
                                       i 
                                     
                                   
                                 
                               
                               ) 
                             
                             
                               - 
                               1 
                             
                           
                           ⁢ 
                           
                             ( 
                             
                               
                                 
                                   e 
                                   -&gt; 
                                 
                                 i 
                               
                               - 
                               
                                 ( 
                                 
                                   
                                     
                                       r 
                                       -&gt; 
                                     
                                     i 
                                   
                                   + 
                                   
                                     
                                       e 
                                       -&gt; 
                                     
                                     i 
                                   
                                 
                                 ) 
                               
                             
                             ) 
                           
                         
                       
                     
                   
                 
               
               
                 
                   ( 
                   17 
                   ) 
                 
               
             
             
               
                 
                   ∝ 
                   
                     
                       - 
                       
                         
                           
                             
                               r 
                               -&gt; 
                             
                             i 
                             T 
                           
                           ( 
                           
                             
                               ∑ 
                               
                                 r 
                                 i 
                               
                             
                             ⁢ 
                             
                               + 
                               
                                 ∑ 
                                 
                                   e 
                                   i 
                                 
                               
                             
                           
                           ) 
                         
                         
                           - 
                           1 
                         
                       
                     
                     ⁢ 
                     
                       
                         r 
                         -&gt; 
                       
                       i 
                     
                   
                 
               
               
                 
                   ( 
                   18 
                   ) 
                 
               
             
           
         
       
     
     That is, the decision whether some optically measured velocity is classified as object proper motion (OPM) is based on the absolute value of that velocity scaled by the variances of the ego-motion-flow (EMF) and the optical flow (OF). Defining {right arrow over (r)} i   T (Σ r +Σ e ) −1 {right arrow over (r)} i  as Mahalanobis norm ∥{right arrow over (r)} i | M , a big Mahalanobis norm indicates a moving point, so one may rely on the computation of the optical flow (OF). Formally, this can be expressed as 
     
       
         
           
             
               
                 
                   
                     
                       o 
                       -&gt; 
                     
                     i 
                   
                   = 
                   
                     { 
                     
                       
                         
                           0 
                         
                         
                           
                             
                               if 
                               ⁢ 
                               
                                   
                               
                               ⁢ 
                               
                                 
                                    
                                   
                                     
                                       r 
                                       -&gt; 
                                     
                                     i 
                                   
                                    
                                 
                                 M 
                               
                             
                             &lt; 
                             
                               θ 
                               M 
                             
                           
                         
                       
                       
                         
                           
                             
                               r 
                               -&gt; 
                             
                             i 
                           
                         
                         
                           otherwise 
                         
                       
                     
                   
                 
               
               
                 
                   ( 
                   19 
                   ) 
                 
               
             
           
         
       
     
     In subsequent processing stages, {right arrow over (o)} i  may be used to extract regions of the visual input where reliably measurable object proper motion (OPM) occurs. 
     CONCLUSION 
     The present invention provides a method and a system being capable of perceiving object proper motion (OPM) from a moving platform. The step-by-step movement estimation including the compensation of ego-motion (EM) prior to the optical flow (OF)-computation seems to be a central aspect for the robustness of the system against firm ego-motion (EM) of the robot. Robustness against noise in the depth and flow-estimation results from the probabilistic rejection mechanism, which neglects velocities based on their size and reliability. 
     The proper object motion computed by this method and system may be used for the attraction of visual attention as well as real-time object interaction. More particularly, an autonomous robot or a car may be controlled using the proper motion of objects, as determined by the invention, e.g. for controlling their own movement or for tracking a moving object by keeping it centred in the image detector&#39;s view.