# Patent

## Title
Method for estimating unknown parameters for a vehicle object detection system

## Abstract

A method for estimating unknown parameters (pan angle (ψ), instantaneous tilt angle (τ) and road geometry of an upcoming road segment) for a vehicle object detection system. The vehicle object detection system is preferably a forward looking, radar-cued vision system having a camera, a radar sensor and an processing unit. The method first estimates the pan angle (ψ), then corrects the coordinates from a radar track so that pan angle (ψ) can be treated as zero, and finally solves a least squares problem that determines best estimates for instantaneous tilt angle (τ) and road geometry. Estimating these parameters enables the vehicle object detection system to identify, interpret and locate objects in a more accurate and efficient manner.

## Background
TECHNICAL FIELD
The present invention relates to a method for estimating unknown parameters in a vehicle object detection system, and more particularly, to a method for estimating camera pan and tilt angles and the geometry of an upcoming road segment.
BACKGROUND OF THE INVENTION
Various types of object detection systems for motor vehicles have been developed which utilize one or more sensors to identify, interpret and/or locate objects surrounding the vehicle. Examples of such systems include forward looking systems (FLS), rear detection systems (RDS), and side detection systems (SDS).
Some FLSs are referred to as “vision-only systems,” in that they only utilize cameras and not radar, laser or other types of sensors. In the case of single camera vision systems (referred to as “monocular vision systems”) and multiple camera vision systems (referred to as “stereo vision system”) a processing unit usually has to estimate the pan and tilt angles of a camera, as well as the vertical curvature of an upcoming road segment in order to accurately interpret and locate detected objects. One way in which these estimates can be performed is by using optical flow techniques in which a focus of expansion is found in the image. These techniques, however, require a significant amount of computations and frequently do not produce robust, accurate estimates.
Another type of object detection system is referred to as a “radar-cued vision system”, which broadly includes all object detection systems that utilize both radar sensors and cameras to identify, interpret and/or locate objects. A typical forward looking, radar-cued vision system for a vehicle includes a radar sensor and one or more cameras directed generally in front of the vehicle. The radar sensor initially detects an object so that a region of the camera image in which the object is likely to be found can be identified, narrowed in on, and analyzed in more detail.
SUMMARY OF THE INVENTION
According to one aspect of the invention, there is provided a method for estimating unknown parameters for a vehicle object detection system. This method includes the steps of generating one or more vision and radar tracks, generating one or more matched pairs from the vision and radar tracks, and estimating a camera pan angle (ψ), an instantaneous tilt angle (τ) and/or the road geometry, wherein the estimates minimize the alignment error between the matched pair.
According to another aspect, an estimate of the camera pan angle is first solved, and then one or more radar tracks are pan corrected using the estimated camera pan angle; effectively setting the camera pan angle equal to zero. Using the pan corrected matched vision and radar tracks (pan corrected matched pairs), the tilt angle (τ) and/or road geometry is subsequently solved for the reduced optimization problem.
According to yet another aspect, the non-linear optimization problem is further simplified by linearization about the predicted values of the parameters. Furthermore, a least squares error criterion is chosen for the optimization metric resulting in a linear least squares problem that can be efficiently solved.
Objects, features and advantages of this invention include, but are certainly not limited to, providing more accurate and efficient pan angle, instantaneous tilt angle and vertical curvature estimates for a vehicle object detection system.



BRIEF DESCRIPTION OF THE DRAWINGS
These and other objects, features and advantages of this invention will be apparent from the following detailed description of the preferred embodiments and best mode, the appended claims and the accompanying drawings, in which:
FIG. 1 is a block diagram of an exemplary object detection system that utilizes the method of the present invention;
FIG. 2 is a schematic plan view of a vehicle demonstrating a camera pan angle (ψ);
FIG. 3 is schematic elevation view of the vehicle of FIG. 2 demonstrating various camera tilt angles (γ, β, τ);
FIG. 4 is a flowchart illustrating an exemplary embodiment overview of a method for estimating unknown parameters for a vehicle object detection system;
FIG. 5 is a more detailed flowchart illustrating an exemplary embodiment of a method for estimating a pan angle, which is preferably a step in the method shown in the flowchart of FIG. 4; and
FIG. 6 is a more detailed flowchart illustrating an exemplary embodiment of a method for estimating an instantaneous tilt angle and vertical curvature, which are also preferably steps in the method shown in the flowchart of FIG. 4.



DETAILED DESCRIPTION OF THE INVENTION
The present method for estimating unknown parameters can be used with a number of different types of vehicle object detection systems, but is particularly well suited for use with forward looking, radar-cued vision systems (both monocular and stereo vision systems). Accordingly, the following description of the preferred embodiment is provided in the context of a forward looking, radar-cued vision system, even though the method of the present invention could be used with other types of object detection systems.
With reference to FIG. 1, there is shown a vehicle object detection system 10 of the forward looking, radar-cued monocular vision type that generally includes a camera 12, a radar sensor 14 and a processing unit 16 connected to the camera and the radar sensor. Of course, system 10 could include one of any number of additional components and could be connected to other systems, devices, networks, etc., such as a warning indicator on a vehicle instrument panel, a vehicle braking system and a vehicle airbag system, to name but a few.
Camera 12 is mounted on a host vehicle 20 and is preferably pointed in a generally forward direction such that it captures a camera image of an area in front of the vehicle. Typically, camera 12 is mounted within the passenger compartment, such as near the rearview mirror, but could alternatively be mounted at other suitable locations inside and outside of the passenger compartment. Each object that is detected within the camera image gives rise to a collection of data referred to as a vision track. A “vision track” is broadly defined as a collection of data that is derived from a camera image and provides information on the position and/or other features of an object. Preferably, each vision track is provided in the form of pixel coordinates (p, q), which use a two-dimensional Cartesian coordinate system. A captured camera image is periodically sent by camera 12 to processing unit 16 so that corresponding vision tracks can be obtained.
Radar sensor 14 is mounted on host vehicle 20 and pans back and forth across a predetermined angular range generally in front of the vehicle. Radar sensor 14 is commonly mounted outside of the passenger compartment, such as near the vehicle grille, but could be mounted at any other appropriate location inside or outside of the passenger compartment. Radar sensor 14 emits a generally narrow fan-shaped radar beam so that objects generally in front of the vehicle reflect the emitted radar back to the sensor. The radar reflections received by radar sensor 14 are sent to processing unit 16 so that corresponding radar tracks can be obtained. A “radar track” is broadly defined as a collection of data that is derived from radar reflections and provides information on the position and/or other features of an object. Preferably, each radar track is provided in the form of range and azimuth angle coordinates (R, θ), which use a two-dimensional Polar coordinate system.
The vision and radar tracks use different coordinate systems and therefore require a method for accurate translation. The image coordinates (p, q) represent a perspective transformation of three-dimensional world Cartesian coordinates into a two-dimensional Cartesian coordinate space, while the radar coordinates (R, θ) represent a two-dimensional subset of the three-dimensional world polar coordinates. Inverting the perspective transformations of the image and radar coordinates to recover the three-dimensional world coordinates, requires knowledge of extrinsic camera parameters and additional world information or constraints. For example, a point in the two-dimensional image maps to a line in the three-dimensional world coordinates. The additional world information or constraints, which includes but is not limited to the vertical curvature of an upcoming road segment, is broadly represented here as the “road geometry.” One method for solving the correspondence between the pixel coordinates (p, q) of the vision track and the range and azimuth angle coordinates (R, θ) of the radar track, will be subsequently described.
Processing unit 16 receives the camera and radar data from camera 12 and radar sensor 14, respectively, and executes a set of instructions for obtaining the vision and radar tracks and for estimating the unknown parameters pan angle, tilt angle and road geometry. Preferably, unit 16 includes a digital microprocessor and a set of firmware instructions stored on a digital memory device, however, other types of electronic processing devices, memory devices and software could be used as well. The firmware can include a number of additional instructions, such as those used for identifying, interpreting and/or locating objects. In that regard, the present method for estimating unknown parameters is only one sub-set of a larger set of firmware instructions likely to be executed by processing unit 16. Once an object is detected, interpreted and/or located, a warning signal can be sent by unit 16 to the appropriate network, device, system, etc. so that a proper response can be determined. For more information on various aspects of vehicle object detection systems, please see U.S. Pat. Nos. 6,873,251, 6,794,987 and 6,615,138, which are assigned to the present assignee and are incorporated herein by reference. Turning now to FIGS. 2 and 3, explanations of the camera pan and tilt angles are provided in greater detail.
The pan angle (ψ) of camera 12 is broadly defined as the horizontal angle between a camera boresight 30 and a radar boresight 32. Horizontal misalignment between the camera and radar boresights 30, 32 leads to pan angle (ψ) and can be caused by a number of factors, but is mainly attributable to the difficulty associated with mounting these two devices such that they are perfectly aligned. In the event that the camera and radar boresights 30, 32 are perfectly aligned (no misalignment), the pan angle (ψ) would=0°.
The various tilt angles of camera 12 are shown in FIG. 3 and generally refer to the vertical misalignment between the camera boresight 30 and some other frame of reference. For instance, the nominal tilt angle (α) of the camera is generally defined as the vertical or “look down” angle between camera boresight 30 and a longitudinal axis 34 of the vehicle. The pitch angle (β) of the vehicle is generally defined as the vertical angle between the longitudinal axis 34 of the vehicle and a line 36 which is representative of the local slope under the vehicle; that is, line 36 is representative of the longitudinal axis of a planar ground patch under the vehicle. Thus, in the case where host vehicle 20 (and hence horizontal axis 34) is perfectly level with the ground patch under the vehicle, the pitch angle (β) of the vehicle is 0°. When the host vehicle accelerates and the rear end of the vehicle drops (as represented in FIG. 3), the pitch angle (β) becomes a larger, positive value. The instantaneous tilt angle (τ) of the camera is broadly defined as the vertical angle between camera boresight 30 and line 36, and is the superposition of the camera nominal tilt angle (α) and the vehicle pitch angle (β). The various tilt angles described above can be affected by a number of different events, including vehicle pitching which can occur when the vehicle drives over bumps or experiences substantial acceleration or deceleration.
The landscape seen in the field-of-view of vehicle object detection system 10, which is mostly of an upcoming road segment, may have an arbitrary geometry that can be approximated using a parameterized road geometry model. For example, the height z of the road surface may be modeled as a function of the down-range and cross-range coordinates x and y in a world Cartesian coordinate system, as in

z=h(x,y,P),  (Eqn. 1)

where P represents a finite collection of parameters. Suitable adjustment of the values of parameters P allows for the approximation of arbitrary road scene geometries. The accuracy of such approximations depends on the complexity of the true road scene geometry, the complexity of the function h, and the number of parameters in P.

One possible choice for such a model is the constant vertical curvature model which assumes that, in side-view, the upcoming road geometry looks like a section of a circle. In this model, the road height z is a function solely of the downrange distance x and a vertical curvature parameter c, i.e., z=h(x,c). More specifically, a further approximation of the circular shape yields the model




  
    
      
        z
        =
        
          
            -
            
              1
              2
            
          
          ⁢
          
            cx
            2
          
        
      
    
    
      
        (
        
          Eqn
          .
          
            
          
          ⁢
          2
        
        )
      
    
  




Knowing the pan angle (ψ), the instantaneous tilt angle (τ) and the road geometry of the upcoming road segment can greatly improve the cost and time efficiency of image processing. That is, knowledge of these parameters provides a more accurate estimate of where to look in the camera image when an object has been detected in the radar image. When an object is initially detected in the radar image, instead of reviewing the entire camera image (a computationally exhaustive process), processing unit 16 selects certain sub-regions of the overall camera image that are most likely to contain the detected object. The larger the selected sub-regions, the more computational analysis that is required. Thus, it is advantageous to reduce the size of the selected camera image sub-regions, and one way to do so is by increasing the accuracy of the estimated pan angle (ψ), instantaneous tilt angle (τ) and the road geometry, as is appreciated by those skilled in the art. Inaccurate pan angle estimates can result in larger uncertainty as to the horizontal position of the vehicle, and inaccurate tilt angle estimates can make an object appear closer or further than it really is. By increasing the accuracy of the pan angle, tilt angle and road geometry estimates, the uncertainty of the position of a detected object is decreased which results in a smaller camera image sub-region to be searched.
Turning now to FIG. 4, there is shown a preferred embodiment of a method 50 for estimating unknown parameters for a vehicle object detection system. Broadly speaking, method 50 simplifies a number of non-linear equations in order to provide a best estimate for the unknown camera pan angle (ψ), camera instantaneous tilt angle (τ) and road geometry, which in this particular embodiment is the vertical curvature (c) of the upcoming road segment. In order to estimate all three unknown parameters, then at least two objects must be mutually detected by camera 12 and radar sensor 14. Of course, the more objects that are detected, the more accurate the corresponding estimates.
Beginning with step 52, processing unit 16 receives camera and radar data from camera 12 and radar sensor 14, respectively. The processing unit reviews the camera image and creates a vision track (pixel coordinates (p, q)) for each object that is detected from the image. Likewise, a similar process is performed for the radar data such that a radar track (range and azimuth angle coordinates (R, θ)) is created for each detected object. Because this particular embodiment is a radar-cued vision system, the objects detected in the camera data are a sub-set of those detected in the radar data. Thus, each vision track has a corresponding radar track that relates to the same detected object and can be used to make a matched pair, step 54. A “matched pair” is broadly defined as a pair of vision and radar tracks where each track generally relates to the location of the same object. Once one or more matched pairs have been obtained, a pan angle can be estimated.
Step 56 provides a pan angle (ψ) estimate that best fits the set of previously obtained matched pairs and does so independently of the instantaneous tilt angle (τ) and vertical curvature (c). Turning now to FIG. 5 where a more detailed description of step 56 is provided, the procedure begins by checking to see if there are any matched pairs to process, step 80. Assuming that the answer to the preceding step was yes, step 82 sets up a pair of equations for each matched pair of coordinates. In order to relate the matched pair's radar track coordinates (x, y) with its vision track coordinates (p, q), the present method solves simultaneous nonlinear equations obtained by substituting the radar-derived object coordinates into the perspective transformation, as constrained by the selected road geometry model. According to a preferred embodiment, a simplified road geometry model of a flat world, c=0, is used. Next, a pan angle (ψ) is calculated that best satisfies those equations.

[(q−q0)x cos τ−fy] cos ψ+[(q−q0)y cos τ−fx] sin ψ=−(q−q0)(z+h)sin τ  (Eqn. 3)

[(p−p0)x cos τ+fx sin τ] cos ψ+[(p−p0)y cos τ+fy sin τ] sin ψ=f(z+h)cos τ−(p−p0)(z+h)sin τ  (Eqn. 4)

Equations 3 and 4 provide a mechanism for relating the coordinates of the radar track to those of the vision track for each matched pair. Equations 3 and 4 are in the form:

A cos ψ+B sin ψ=C  (Eqn. 5)

D cos ψ+E sin ψ=F  (Eqn. 6)

Coefficients A, B, C, D, E and F should be calculated, step 84. Because the estimated pan angle is relatively insensitive to the value assumed for r, either a nominal value or the current estimate for long-term tilt angle can be used during the solution of the following equations. If AE−BD≠0, then the pan angle can be estimated using the following equation, step 86.




  
    
      
        
          Sin
          ⁢
          
            
          
          ⁢
          ψ
        
        =
        
          
            AF
            -
            CD
          
          
            AE
            -
            BD
          
        
      
    
    
      
        (
        
          Eqn
          .
          
            
          
          ⁢
          7
        
        )
      
    
  




If AE−BD=0, it means the two equations are not linearly independent. Using a small angle approximation (cos ψ≈1) and if B≠0, then the pan angle can be estimated using the next equation, step 88.





  
    
      
        ψ
        =
        
          
            C
            -
            A
          
          B
        
      
    
    
      
        (
        
          Eqn
          .
          
            
          
          ⁢
          8
        
        )
      
    
  




In the case where AE−BD=0 or B=0 or (AF−CD)>(AE−BD) the pan angle estimate using this matched pair can be abandoned and control passes back to step 80. This process continues until there are no more matched pairs to process, at which time step 90 filters all of the raw pan angle estimates to provide a single filtered, pan angle estimate (ψ). According to a preferred embodiment, a low-pass filter is used to filter the various pan angle estimates.

Once the camera pan angle (ψ) has been estimated, each radar track can be corrected or re-processed with the estimated pan angle so that a set of pan corrected matched pairs, each having a corrected radar track, is generated. Because the corrected radar track information has already been corrected for errors attributable to the pan angle (ψ), subsequent estimates for the instantaneous tilt angle (τ) and vertical curvature (c) can simply assume that pan angle (ψ)=0, thus simplifying the equations.
In step 60, a least squares problem is used which includes two unknown coefficients (instantaneous tilt angle τ and vertical curvature c) and N equations (one equation for each matched pair). Turning now to FIG. 6, a more detailed flowchart describing step 60 is shown. First, from the perspective transformation:




  
    
      
        
          x
          c
        
        =
        
          
            
              (
              
                cos
                ⁢
                
                  
                
                ⁢
                ψ
                ⁢
                
                  
                
                ⁢
                cos
                ⁢
                
                  
                
                ⁢
                τ
              
              )
            
            ⁢
            x
          
          +
          
            
              (
              
                sin
                ⁢
                
                  
                
                ⁢
                ψ
                ⁢
                
                  
                
                ⁢
                cos
                ⁢
                
                  
                
                ⁢
                τ
              
              )
            
            ⁢
            y
          
          +
          
            sin
            ⁢
            
              
            
            ⁢
            
              τ
              ⁡
              
                (
                
                  h
                  -
                  
                    
                      1
                      2
                    
                    ⁢
                    
                      cx
                      2
                    
                  
                
                )
              
            
          
        
      
    
    
      
        (
        
          Eqn
          .
          
            
          
          ⁢
          9
        
        )
      
    
  
  
    
      
        
          z
          c
        
        =
        
          
            
              -
              
                (
                
                  cos
                  ⁢
                  
                    
                  
                  ⁢
                  ψ
                  ⁢
                  
                    
                  
                  ⁢
                  sin
                  ⁢
                  
                    
                  
                  ⁢
                  τ
                
                )
              
            
            ⁢
            x
          
          -
          
            
              (
              
                sin
                ⁢
                
                  
                
                ⁢
                ψ
                ⁢
                
                  
                
                ⁢
                sin
                ⁢
                
                  
                
                ⁢
                τ
              
              )
            
            ⁢
            y
          
          +
          
            cos
            ⁢
            
              
            
            ⁢
            
              τ
              ⁡
              
                (
                
                  h
                  -
                  
                    
                      1
                      2
                    
                    ⁢
                    
                      cx
                      2
                    
                  
                
                )
              
            
          
        
      
    
    
      
        (
        
          Eqn
          .
          
            
          
          ⁢
          10
        
        )
      
    
  




Next, the following Taylor series are used to linearize Eqns. 9 and 10 about a nominal tilt angle (τ0), step 100. It is preferable that nominal tilt angle (τ0) be approximately equal to the true tilt angle.

sin τ=sin τ0+cos τ0(τ−τ0)  (Eqn. 11)

cos τ=cos τ0−sin τ0(τ−τ0)  (Eqn. 12)

Step 100 substitutes the Taylor series equations and assumes that the pan angle (ψ)=0 (pan corrected radar tracks already take pan angle into account) to convert trigonometric Eqns. 9 and 10 into the following, simpler linear equations. So long as the tilt deviations remain small, this conversion can be accomplished without much loss of accuracy.





  
    
      
        
          x
          c
        
        ≈
        
          
            x
            ⁢
            
              
            
            ⁢
            cos
            ⁢
            
              
            
            ⁢
            
              τ
              0
            
          
          -
          
            x
            ⁢
            
              
            
            ⁢
            sin
            ⁢
            
              
            
            ⁢
            
              
                τ
                0
              
              ⁡
              
                (
                
                  τ
                  -
                  
                    τ
                    0
                  
                
                )
              
            
          
          +
          
            h
            ⁢
            
              
            
            ⁢
            sin
            ⁢
            
              
            
            ⁢
            
              τ
              0
            
          
          +
          
            h
            ⁢
            
              
            
            ⁢
            cos
            ⁢
            
              
            
            ⁢
            
              
                τ
                0
              
              ⁡
              
                (
                
                  τ
                  -
                  
                    τ
                    0
                  
                
                )
              
            
          
          -
          
            
              1
              2
            
            ⁢
            
              cx
              2
            
            ⁢
            sin
            ⁢
            
              
            
            ⁢
            
              τ
              0
            
          
          -
          
            
              1
              2
            
            ⁢
            
              cx
              2
            
            ⁢
            cos
            ⁢
            
              
            
            ⁢
            
              
                τ
                0
              
              ⁡
              
                (
                
                  τ
                  -
                  
                    τ
                    0
                  
                
                )
              
            
          
        
      
    
    
      
        (
        
          Eqn
          .
          
            
          
          ⁢
          13
        
        )
      
    
  
  
    
      
        
          z
          c
        
        ≈
        
          
            
              -
              x
            
            ⁢
            
              
            
            ⁢
            sin
            ⁢
            
              
            
            ⁢
            
              τ
              0
            
          
          -
          
            x
            ⁢
            
              
            
            ⁢
            cos
            ⁢
            
              
            
            ⁢
            
              
                τ
                0
              
              ⁡
              
                (
                
                  τ
                  -
                  
                    τ
                    0
                  
                
                )
              
            
          
          +
          
            h
            ⁢
            
              
            
            ⁢
            cos
            ⁢
            
              
            
            ⁢
            
              τ
              0
            
          
          -
          
            h
            ⁢
            
              
            
            ⁢
            sin
            ⁢
            
              
            
            ⁢
            
              
                τ
                0
              
              ⁡
              
                (
                
                  τ
                  -
                  
                    τ
                    0
                  
                
                )
              
            
          
          -
          
            
              1
              2
            
            ⁢
            
              cx
              2
            
            ⁢
            cos
            ⁢
            
              
            
            ⁢
            
              τ
              0
            
          
          +
          
            
              1
              2
            
            ⁢
            
              cx
              2
            
            ⁢
            sin
            ⁢
            
              
            
            ⁢
            
              
                τ
                0
              
              ⁡
              
                (
                
                  τ
                  -
                  
                    τ
                    0
                  
                
                )
              
            
          
        
      
    
    
      
        (
        
          Eqn
          .
          
            
          
          ⁢
          14
        
        )
      
    
  




The terms




  
    -
    
      1
      2
    
  
  ⁢
  
    cx
    2
  
  ⁢
  cos
  ⁢
  
    
  
  ⁢
  
    
      τ
      0
    
    ⁡
    
      (
      
        τ
        -
        
          τ
          0
        
      
      )
    
  




in Eqn. 13 and





  
    +
    
      1
      2
    
  
  ⁢
  
    cx
    2
  
  ⁢
  sin
  ⁢
  
    
  
  ⁢
  
    
      τ
      0
    
    ⁡
    
      (
      
        τ
        -
        
          τ
          0
        
      
      )
    
  




in Eqn. 14 are assumed to be negligible and can be removed. With these terms removed, step 102 manipulates Eqn. 3 by substituting Eqns. 13 and 14 for xc and zc, respectively, which results in a single linear equation having two unknown coefficients (τ, c) for each matched pair. Once these equations have been generated, a least squares problem can be set up in the following form:





  
    
      
        
          
            [
            
              A
              ⁢
              
                
              
              ⁢
              B
            
            ]
          
          ⁡
          
            [
            
              c
              τ
            
            ]
          
        
        =
        C
      
    
    
      
        (
        
          Eqn
          .
          
            
          
          ⁢
          15
        
        )
      
    
  




where





  
    
      
        A
        =
        
          
            
              -
              
                1
                2
              
            
            ⁢
            
              (
              
                p
                -
                
                  p
                  0
                
              
              )
            
            ⁢
            
              x
              2
            
            ⁢
            sin
            ⁢
            
              
            
            ⁢
            
              τ
              0
            
          
          +
          
            
              f
              2
            
            ⁢
            
              x
              2
            
            ⁢
            cos
            ⁢
            
              
            
            ⁢
            
              τ
              0
            
          
        
      
    
    
      
        (
        
          Eqn
          .
          
            
          
          ⁢
          16
        
        )
      
    
  
  
    
      
        B
        =
        
          
            
              (
              
                p
                -
                
                  p
                  0
                
              
              )
            
            ⁢
            
              (
              
                
                  h
                  ⁢
                  
                    
                  
                  ⁢
                  cos
                  ⁢
                  
                    
                  
                  ⁢
                  
                    τ
                    0
                  
                
                -
                
                  x
                  ⁢
                  
                    
                  
                  ⁢
                  sin
                  ⁢
                  
                    
                  
                  ⁢
                  
                    τ
                    0
                  
                
              
              )
            
          
          +
          
            f
            ⁡
            
              (
              
                
                  x
                  ⁢
                  
                    
                  
                  ⁢
                  cos
                  ⁢
                  
                    
                  
                  ⁢
                  
                    τ
                    0
                  
                
                +
                
                  h
                  ⁢
                  
                    
                  
                  ⁢
                  sin
                  ⁢
                  
                    
                  
                  ⁢
                  
                    τ
                    0
                  
                
              
              )
            
          
        
      
    
    
      
        (
        
          Eqn
          .
          
            
          
          ⁢
          17
        
        )
      
    
  
  
    
      
        C
        =
        
          
            fx
            ⁡
            
              (
              
                
                  
                    τ
                    0
                  
                  ⁢
                  cos
                  ⁢
                  
                    
                  
                  ⁢
                  
                    τ
                    0
                  
                
                -
                
                  sin
                  ⁢
                  
                    
                  
                  ⁢
                  
                    τ
                    0
                  
                
              
              )
            
          
          +
          
            fh
            ⁡
            
              (
              
                
                  cos
                  ⁢
                  
                    
                  
                  ⁢
                  
                    τ
                    0
                  
                
                +
                
                  
                    τ
                    0
                  
                  ⁢
                  sin
                  ⁢
                  
                    
                  
                  ⁢
                  
                    τ
                    0
                  
                
              
              )
            
          
          -
          
            
              (
              
                p
                -
                
                  p
                  0
                
              
              )
            
            ⁢
            
              x
              ⁡
              
                (
                
                  
                    cos
                    ⁢
                    
                      
                    
                    ⁢
                    
                      τ
                      0
                    
                  
                  +
                  
                    
                      τ
                      0
                    
                    ⁢
                    sin
                    ⁢
                    
                      
                    
                    ⁢
                    
                      τ
                      0
                    
                  
                
                )
              
            
          
          +
          
            
              (
              
                p
                -
                
                  p
                  o
                
              
              )
            
            ⁢
            
              h
              ⁡
              
                (
                
                  
                    
                      τ
                      0
                    
                    ⁢
                    cos
                    ⁢
                    
                      
                    
                    ⁢
                    
                      τ
                      0
                    
                  
                  -
                  
                    sin
                    ⁢
                    
                      
                    
                    ⁢
                    
                      τ
                      0
                    
                  
                
                )
              
            
          
        
      
    
    
      
        (
        
          Eqn
          .
          
            
          
          ⁢
          18
        
        )
      
    
  




A less accurate but considerably simplified alternative to the equations above, is:





  
    
      
        
          
            [
            
              
                
                  -
                  
                    1
                    2
                  
                
                ⁢
                fx
              
              -
              f
            
            ]
          
          ⁡
          
            [
            
              c
              τ
            
            ]
          
        
        =
        
          p
          -
          
            
              p
              0
            
            ⁢
            
              fh
              x
            
          
        
      
    
    
      
        (
        
          Eqn
          .
          
            
          
          ⁢
          19
        
        )
      
    
  




A minimum of two, but preferably more, matched pairs are needed to solve the least squares problem such that it has a unique solution, step 106. As is appreciated by those skilled in the art, matrix operations can be used to solve for the least squares problem above so that the sum of the squares of the errors is minimized. Furthermore, for good conditioning of the problem, it is desirable to have detected objects at ranges r which are as different as possible (objects at differing distances from the host vehicle) and which are not too close to zero (too close to the vehicle). It is also preferably that the quantity x (downrange distance to the object) be corrected before use in the equation above by accounting for the non-colocation of the radar and camera (i.e., the downrange distance should be relative to the camera).
The quantities estimated using the scheme described above, namely the curvature c and the instantaneous tilt angle τ, can vary rapidly. On the other hand, the pan angle ψ is assumed to be fairly constant. For this reason, the system does not need to constantly acquire measurements for the pan angle ψ. Also, a slowly varying estimate of nominal tilt angle can be used instead of the instantaneous value estimated above.
It is expected that the estimates obtained using the schemes above will be filtered, step 108, to obtain final estimates. For example, the tilt angle could possibly be filtered using a model of the vehicle suspension's pitch dynamics or by some other technique known in the art.
A last step 64, entitled double checking, represents operations such as making sure there are enough matched pairs to make an estimate, checking intermediate values (e.g., before performing an arcsine operation), etc.
It will thus be apparent that there has been provided in accordance with this invention a method for estimating unknown parameters (namely, pan angle, instantaneous tilt angle and road geometry) which achieves the aims and advantages specified herein. It will, of course, be understood that the forgoing description is of preferred exemplary embodiments of the invention and that the invention is not limited to the specific embodiments shown. Various changes and modifications will become apparent to those skilled in the art and all such changes and modifications are intended to be within the scope of this invention.

## Claims

1. A method for operating an object detection system to determine a location of an object, comprising the steps of:
(a) generating a vision track from data provided by a camera, said vision track including data indicative of a locations of the object, wherein a pixel coordinate is generated for each object that is detected from an image;
(b) generating a radar track from data provided by a radar sensor, said radar track including data indicative of the locations of the object, wherein a range and azimuth angle coordinate is generated for each detected object;
(c) generating a matched pair indicative of the location of the object based upon the data from said vision and radar tracks;
(d) using the matched pair, estimating a camera pan angle (ψ) and an instantaneous tilt angle (τ) and road geometry, wherein said camera pan angle is defined as the horizontal angle between a camera boresight and a radar boresight, wherein said instantaneous tilt angle is a vertical misalignment between the camera boresight and the longitudinal axis of a planar ground patch under a vehicle, wherein said estimated camera pan angle, instantaneous tilt angle (τ), and road geometry minimize the alignment error between said matched pair, wherein said estimated camera pan angle is used to correct or reprocess each radar track so that a set of pan angle corrected matched pairs, each having a correct pan angle, is generated, and wherein said estimated instantaneous tilt angle (τ) and road geometry comprise linearizing at least one trigonometric equation about a nominal tilt angle (τ0) so that said trigonometric equation becomes a linear equation; and
(e) determining the location of the object, said determining including identifying, locating or interpreting the object using the estimates of the pan angle (ψ) and the instantaneous tilt angle (τ) and the road geometry, wherein the estimates are used to correct for alignment differences between the camera and the radar sensor while traveling a roadway, and thereby correct for position of the object errors in the vision track and the radar track.

2. The method of claim 1, wherein a vertical curvature (c) of an upcoming road segment represents said road geometry.

3. The method of claim 1, wherein said vision track is provided in terms of a two-dimensional Cartesian coordinate system (pixels (p, q)), and said radar track is provided in terms of a two-dimensional polar coordinate system (range and azimuth angle (R, θ)).

4. The method of claim 1, wherein said pan angle (ψ) estimate is filtered with a low-pass filter to provide a filtered pan angle estimate.

5. The method of claim 1, wherein step (d) further includes solving an optimization problem that minimizes the 3-D world coordinate alignment error between said matched pair.

6. The method of claim 5, wherein step (d) further includes utilizing a least squares problem.

7. The method of claim 1, wherein step (d) further includes removing any negligible terms and manipulating said linear equation so that a single linear equation is provided for said matched pair.

8. The method of claim 1, wherein step (d) further includes using said estimated camera pan angle (ψ) and said radar track to provide a pan corrected matched pair having a radar track that has been generally corrected for errors attributable to a camera pan angle.

9. The method of claim 8, wherein step (d) further includes using said pan corrected matched pair to estimate an instantaneous tilt angle (τ) and road geometry.

10. The method of claim 1, wherein the vehicle object detection system is a forward looking, radar-cued vision system.

11. A method for operating an object detection system to determine a location of an object, comprising the steps of:
(a) generating a vision track from data provided by a camera, said vision track including data indicative of a locations of the object, wherein a pixel coordinate is generated for each object that is detected from an image;
(b) generating a radar track from data provided by a radar sensor, said radar track including data indicative of the locations of the object, wherein a range and azimuth angle coordinate is generated for each detected object;
(c) generating a matched pair indicative of the location of the object based upon the data from said vision and radar tracks;
(d) estimating a camera pan angle (ψ) that best fits said matched pair, wherein said camera pan angle is defined as the horizontal angle between a camera boresight and a radar boresight, wherein said instantaneous tilt angle is a vertical misalignment between the camera boresight and the longitudinal axis of a planar ground patch under a vehicle, wherein said estimated camera pan angle, instantaneous tilt angle (τ), and road geometry minimize the alignment error between said matched pair, wherein said estimated camera pan angle is used to correct or reprocess each radar track so that a set of pan angle corrected matched pairs, each having a correct pan angle, is generated;
(e) using said estimated camera pan angle (ψ) and said radar track to provide a pan corrected matched pair(s) having a radar track that has been generally corrected for errors attributable to a camera pan angle, wherein said estimates comprise linearizing at least one trigonometric equation about a nominal tilt angle (τ0) so that said trigonometric equation becomes a linear equation;
(f) using said pan corrected matched pair(s) to estimate an instantaneous tilt angle (τ) and road geometry; and
(g) determining the location of the object, said determining including identifying, locating or interpreting the object using the estimates of the pan angle (ψ) and the instantaneous tilt angle (τ) and the road geometry, wherein the estimates are used to correct for alignment differences between the camera and the radar sensor while traveling a roadway, and thereby correct for position of the object errors in the vision track and the radar track.

12. A method for operating an object detection system to determine a location of one or more of a plurality of objects, comprising the steps of:
(a) generating a plurality of vision tracks from data provided by a camera, said vision tracks including data indicative of locations of the plurality of objects, wherein a pixel coordinate is generated for each object that is detected from an image;
(b) generating a plurality of radar tracks from data provided by a radar sensor, said radar tracks including data indicative of the locations of the plurality of objects, wherein a range and azimuth angle coordinate is generated for each detected object;
(c) generating a plurality of matched pairs indicative of the locations based upon the data from said vision and radar tracks;
(d) estimating a camera pan angle (ψ) that best fits said plurality of matched pairs, wherein said camera pan angle is defined as the horizontal angle between a camera boresight and a radar boresight, wherein said instantaneous tilt angle is a vertical misalignment between the camera boresight and the longitudinal axis of a planar ground patch under a vehicle, wherein said estimated camera pan angle, instantaneous tilt angle (τ), and road geometry minimize the alignment error between said matched pair, wherein said estimated camera pan angle is used to correct or reprocess each radar track so that a set of pan angle corrected matched pairs, each having a correct pan angle, is generated, and wherein said estimated instantaneous tilt angle (τ) and road geometry comprise linearizing at least one trigonometric equation about a nominal tilt angle (τ0) so that said trigonometric equation becomes a linear equation;
(e) using said estimated camera pan angle (ψ) and said plurality of radar tracks to provide a plurality of pan corrected matched pairs each having a radar track that has been generally corrected for errors attributable to a camera pan angle;
(f) using said plurality of pan corrected matched pairs to generate a least squares problem, said least squares problem treats said estimated camera pan angle (ψ) as zero;
(g) linearizing one or more trigonometric equations about a nominal tilt angle (τ0) so that said trigonometric equations become linear equations;
(h) solving said least squares problem to estimate both said instantaneous tilt angle (τ) and said vertical curvature (c); and
(i) determining the location of one or more of the objects, said determining including identifying, locating or interpreting objects using the estimates of the pan angle (ψ) and the instantaneous tilt angle (τ) and the road geometry, wherein the estimates are used to correct for alignment differences between the camera and the radar sensor while traveling a roadway, and thereby correct for position of the object errors in the vision track and the radar track.