Gavaghan Geodesy Library

The Gavaghan Geodesy Library is a prominent, open-source software library originally developed in Java by Mike Gavaghan. It provides highly accurate geodetic calculations on an ellipsoidal Earth model by implementing Vincenty’s Formulae.

The original codebase has been widely adopted, ported to other languages like C#, and is frequently used in Geographic Information Systems (GIS) and location-aware applications where simple spherical assumptions (like the Haversine formula) are insufficient. Core Functionality The library solves the two fundamental problems of geodesy:

The Inverse Problem (Distance & Bearing): Calculates the precise ellipsoidal distance, forward azimuth, and reverse azimuth between two coordinate points (latitude/longitude).

The Direct Problem (Destination Projection): Computes the destination coordinates and reverse azimuth given a starting point, an initial direction (azimuth), and a specific distance. Why Vincenty’s Formulae?

Unlike simpler equations that assume the Earth is a perfect sphere, Vincenty’s methods model the Earth as an oblate spheroid (flatter at the poles).

Extreme Accuracy: It achieves a precision of up to 0.5 millimeters (0.020 inches) on the reference ellipsoid.

Reference Ellipsoids: The library natively benchmarks against industry standard models like WGS84 (used by GPS), GRS80, and older standards like the International Ellipsoid. Standard Usage Example

In Java, developers typically interact with the library through a core GeodeticCalculator class:

import org.gavaghan.geodesy.*; GeodeticCalculator geoCalc = new GeodeticCalculator(); Ellipsoid reference = Ellipsoid.WGS84; // Model the Earth’s shape // Define two global positions (Latitude, Longitude, Elevation) GlobalPosition pointA = new GlobalPosition(37.7749, -122.4194, 0.0); GlobalPosition pointB = new GlobalPosition(34.0522, -118.2437, 0.0); // Calculate the geodetic curve (Inverse Problem) GeodeticCurve curve = geoCalc.calculateGeodeticCurve(reference, pointA, pointB); double distance = curve.getEllipsoidalDistance(); // Distance in meters double azimuth = curve.getAzimuth(); // Bearing in degrees Use code with caution. Critical Limitations to Consider

While highly efficient and accurate for most global navigation applications, the Vincenty implementation in this library carries historical mathematical trade-offs: airbreather/Gavaghan.Geodesy – GitHub