Tight optical integration (TOI) of images with GPS range measurements

• US 9,726,765 B2
• Filed: 07/20/2011
• Issued: 08/08/2017
• Est. Priority Date: 01/27/2009
• Status: Active Grant

• Current Assignee: Ohio University
• Original Assignee: Ohio University
• Inventors: Arthur, Thomas D., Zhu, Zhen, Bhattacharya, Sumit, Scheff, Kim M., de Haag, Maarten U., Johnson, Kevin L.
• Primary Examiner(s): Gregory, Bernarr
• Assistant Examiner(s): Mull, Fred H
• Application Number:
  US13/186,704
• Publication Number:
  US 20120127030A1
• Time in Patent Office: 2,211 Days
• Field of Search: None
• CPC Class Codes: G01S 19/45 by combining measurements of signals from the satellite radio beacon positioning system with a supplementary measurement

Abstract

A method of determining the position of a user in 3D space is disclosed. In one example, the method comprises utilizing range measurements from one or more GPS satellite vehicles, angular information of a camera located at the user position, and a camera image associated with one or more known markers. The method further comprises determining the angular information of the camera by determining a direction cosine matrix between a camera frame of reference and an earth frame of reference, and designating unit vectors that individually extend from a position within the camera image associated with one of the known markers through the camera focal point to the respective known marker The method also includes integrating into an ordinary least squares matrix, the GPS range measurements, the angular information of the camera, and the unit vectors, and calculating the user position by solving the ordinary least squares matrix.

Claims

1. A method of determining a position of a user in 3D space utilizing an integrated navigation system including a GPS receiver and a camera located at the user position when a limited number of GPS signals are available, the method comprising:

capturing, by the camera, a camera image of a known marker through a focal point of the camera with optical measurements that extend from a position of the known marker within the camera image to the user position of the camera;

determining two or more pseudoranges between a true user position and two or more GPS satellites with the integrated navigation system, wherein the two or more GPS satellites are insufficient to generate GPS measurements necessary to calculate a GPS standalone position;

determining a marker unit vector from the optical measurements captured by at least one optic sensor associated with the camera that extends from the position within the camera image associated with the known marker through the focal point of the camera with the integrated navigation system, wherein unavailable GPS measurements necessary to calculate the GPS standalone position are replaced by the marker unit vector;

converting, by the GPS receiver, the marker unit vector from a camera frame of reference to an earth frame of reference of GPS measurements by incorporating attitude information of the camera obtained from inertial measurements of the camera at the user position to transform a direction cosine of the marker unit vector in the camera frame of reference to the earth frame of reference;

integrating, by the GPS receiver, the two or more pseudoranges and the marker unit vector into an ordinary least squares matrix with the integrated navigation system;

generating a shortest vector starting from an estimated user position to a point on the marker unit vector, wherein the shortest vector is perpendicular to the marker unit vector at the point on the marker unit vector;

converging, by the GPS receiver, a line joining the estimated user position and the marker and the two or more pseudoranges derived from the estimated user position to the point on the marker unit vector generated from the shortest vector being perpendicular to the marker unit vector;

determining the position of the user in 3D space from the convergence of the line joining the estimated user position and the two or more pseudoranges derived from the estimated user position to the point on the marker unit vector generated from the shortest vector being perpendicular to the marker unit vector by solving the ordinary least squares matrix with the integrated navigation system so that the two or more pseudoranges are supplemented with the marker unit vector to overcome the unavailable GPS measurements to calculate the GPS standalone position.

2. The method of claim 1, wherein the known marker comprises one or more of a terrestrial marker, an astronomical marker, a vehicular marker, and an airborne marker.

3. The method of claim 1, wherein the marker is one of one or more known markers each having known geographical coordinates.

4. The method of claim 1, wherein the pseudoranges of the true user position from the two or more GPS satellites are determined based on a time of arrival of a signal from each of the two or more GPS satellites.

5. The method of claim 1, wherein the camera located at the user position is configured to provide imaging and angular information utilizing one or more of visible light, infrared light, UV light, or using LIDAR, LADAR, or RADAR based systems.

6. The method of claim 1 further comprising:

determining a direction cosine matrix (DCM) between a camera frame of reference and an earth frame of reference; and

transforming the marker unit vector from the camera frame of reference to the earth frame of reference using the DCM.

7. The method of claim 6, further comprising:

directing the camera from the true user position toward the known marker before determining the DCM.

8. The method of claim 1 wherein integrating the two or more pseudoranges and the marker unit vector into the ordinary least squares matrix comprises:

determining a convergence vector based on the marker unit vector, a position of the known marker, and an estimated user position;

defining a GPS-visual measurement vector comprising the two or more pseudoranges and the convergence vector; and

integrating the GPS-visual measurement vector into the ordinary least squares matrix.

9. The method of claim 1 further comprising:

determining a set of satellite unit vectors comprising two or more satellite unit vectors between an estimated user position and the two or more GPS satellites;

defining a direction cosine matrix comprising the set of satellite unit vectors and the marker unit vector;

defining a geometry matrix comprising the direction cosine matrix and a camera weighting factor; and

integrating the geometry matrix into the ordinary least squares matrix.

10. The method of claim 1 wherein the user position is a second estimated user position, and determining the position of the user in 3D space by solving the ordinary least squares matrix comprises:

determining an offset between a first estimated user position and the second estimated user position by solving the ordinary least squares matrix; and

adding the offset to the first estimated user position.

11. The method of claim 1 wherein determining the marker unit vector comprises:

mapping individual pixels in the camera into respective unit vectors individually extending from a respective camera pixel to the camera focal point, thereby generating a pixel map; and

designating a pixel in the pixel map and a respective unit vector that extends from the designated pixel through the camera focal point toward the known marker as the marker unit vector.

12. The method of claim 1 wherein a position of the known marker is determined by self-surveying the known marker while sufficient GPS satellites are available for the integrated navigation system to determine the position of the known marker.

13. The method of claim 1, wherein the converging comprises:

v=[I+e_ve_v^T](x_v−

x^{^}), where

v=the shortest vector,

I=a 3×

3 identity matrix,

e_v=a unit vector pointing from the camera to the known marker,

e_v^T=unit vector pointing from the camera to the known marker relative to a GPS-visual T matrix,

x_v=position of the known marker, and

x^=current estimate of the user position.

14. A navigation system comprising:

a camera that includes at least one optic sensor; and

a GPS receiver operatively coupled to the camera, the navigation system when a limited number of GPS signals are available is configured to;

determine two or more pseudoranges between a true user position and two or more GPS satellites, wherein the two or more GPS satellites are insufficient to generate GPS measurements necessary to calculate a GPS standalone position;

determine a marker unit vector from optical measurements captured by the at least one optic sensor included in the camera that extends from a position within a camera image associated with a known marker through a focal point of the camera, wherein unavailable GPS measurements necessary to calculate the GPS standalone position are replaced by the marker unit vector;

convert the marker unit vector from a camera frame of reference to an earth frame of reference of GPS measurements by incorporating attitude information of the camera obtained from inertial measurements of the camera at the user position to transform a direction cosine of the marker unit vector in the camera frame of reference to the earth frame of reference;

integrate the two or more pseudoranges and the marker unit vector into an ordinary least squares matrix;

generate a shortest vector starting from an estimated user position to a point on the marker unit vector, wherein the shortest vector is perpendicular to the marker unit vector at the point on the marker unit vector; and

converge a line joining the estimated user position and the marker and the two or more pseudoranges derived from the estimated user position to the point on the market unit vector generated from the shortest vector being perpendicular to the marker unit vector; and

determine the position of the user in 3D space from the convergence of the line joining the estimated user position and the marker and the two or more pseudoranges derived from the estimated user position to the point on the marker unit vector generated from the shortest vector being perpendicular to the marker unit vector by solving the ordinary least squares matrix so that the two or more pseudoranges are supplemented with the marker unit vector to overcome the unavailable GPS measurements to calculate the GPS standalone position.

15. The system of claim 14 wherein the camera is configured to provide imaging and angular information utilizing one or more of visible light, infrared light, UV light, or using LIDAR, LADAR, or RADAR based systems.

16. The system of claim 14 wherein the system is further configured to:

determine a direction cosine matrix (DCM) between a camera frame of reference and an earth frame of reference; and

transform the marker unit vector from the camera frame of reference to the earth frame of reference using the DCM.

17. The system of claim 14 wherein the system is further configured to:

direct the camera from the true user position toward the known marker before determining the DCM.

18. The system of claim 14 wherein the system is further configured to:

determine a set of satellite unit vectors comprising two or more satellite unit vectors between an estimated user position and the two or more GPS satellites;

define a direction cosine matrix comprising the set of satellite unit vectors and the marker unit vector;

define a geometry matrix comprising the direction cosine matrix and a camera weighting factor; and

integrate the geometry matrix into the ordinary least squares matrix.

19. The system of claim 14 wherein the user position is a second estimated user position, and the system determines the position of the user in 3D space by:

determining an offset between a first estimated user position and the second estimated user position by solving the ordinary least squares matrix; and

adding the offset to the first estimated user position.

20. The system of claim 14 wherein the system determines the marker unit vector by:

mapping individual pixels in the camera into respective unit vectors individually extending from a respective camera pixel to the camera focal point, thereby generating a pixel map; and

designating a pixel in the pixel map and respective unit vector that extends from the respective marker pixel through the camera focal point and toward the known marker as the marker unit vector.

21. The system of claim 14 wherein the system is further configured to:

determine a position of the known marker by self-surveying the known marker while sufficient GPS satellites are available to determine the position of the known marker.

22. The system of claim 14, wherein the GPS receiver is further configured to converge based on:

v=[I+e_ve_v^T](x_v−

x^{^}), where

v=the shortest vector,

I=a 3×

3 identity matrix,

e_v=a unit vector pointing from the camera to the known marker,

e_v^T=unit vector pointing from the camera to the known marker relative to a GPS-visual T matrix,

x_v=position of the known marker, and

x^=current estimate of the user position.