Online geodesic calculations using the GeodSolve utility

Geodesic calculation:
    lat1 lon1 lat2 lon2 azi1 azi2 s12
    lat1 lon1 azi1 s12 lat2 lon2 azi2

Input (ex. «40.6 -73.8 49°01'N 2°33'E» [inverse], «40d38'23"N 073d46'44"W 53d30' 5850e3» [direct]):
   

Output format:    
Heading at point 2:    
Longitude:    
Output precision:  
Equatorial radius: meters
Flattening:

Select action:
   

Geodesic (input in black, output in blue):

    ellipsoid (a f)     = 6378137 1/298.257223563 (WGS84)
    status              = OK

    lat1 lon1 fazi1 (°) = -30.00000000 113.00000000 -168.78839727
    lat2 lon2 fazi2 (°) = -38.00000000 111.00000000 -167.66719369
    s12 (m)             = 906343.630

    a12 (°)             = 8.16064007
    m12 (m)             = 903288.747
    M12 M21             = 0.9898922644 0.9898981019
    S12 (m^2)           = 791860688011


GeodSolve (version 2.3) performs geodesic calculations for an arbitrary ellipsoid of revolution. The shortest path between two points on the ellipsoid at (lat1, lon1) and (lat2, lon2) is called the geodesic; its length is s12 and the geodesic from point 1 to point 2 has azimuths azi1 and azi2 at the two end points. There are two standard geodesic problems:

Latitudes and longitudes can be given in various formats, for example (these all refer to the position of Timbuktu):
        16.776 -3.009
        16d47' -3d1'
        W3°0'34" N16°46'33"
        3:0:34W 16:46:33N
Azimuths are given in degrees clockwise from north. The distance s12 is in meters.

The additional quantities computed are:

The ellipsoid is specified by its equatorial radius, a, and its flattening, f = (a − b)/a, where b is the polar semi-axis. The default values for these parameters correspond to the WGS84 ellipsoid. The method is accurate for −99 ≤ f ≤ 0.99 (corresponding to 0.01 ≤ b/a ≤ 100). Note that f is negative for a prolate ellipsoid (b > a) and that it can be entered as a fraction, e.g., 1/297.

GeodSolve is accurate to about 15 nanometers (for the WGS84 ellipsoid) and gives solutions for the inverse problem for any pair of points.

GeodSolve, which is a simple wrapper of the GeographicLib::Geodesic class, is one of the utilities provided with GeographicLib. Geodesics can also be computed using JavaScript; see the JavaScript geodesic calculator and geodesics on Google maps. If you wish to use GeodSolve directly, download and compile GeographicLib. The algorithms are described in C. F. F. Karney, Algorithms for geodesics, J. Geodesy 87, 43–55 (2013); DOI: 10.1007/s00190-012-0578-z; addenda: geod-addenda.html. See also the Wikipedia page, Geodesics on an ellipsoid.


Charles Karney <karney@alum.mit.edu> (2022-04-10)
GeographicLib home