Principles and Implementation of GPS Coordinate Distance Calculation Using Haversine Formula

GPS Coordinate System and Distance Calculation Fundamentals

The Global Positioning System (GPS) uses latitude and longitude coordinates to represent positions on the Earth's surface. Longitude indicates east-west position, ranging from -180 to +180 degrees; latitude indicates north-south position, ranging from -90 to +90 degrees. When calculating distances between two coordinate points, the spherical geometry of the Earth must be considered, as the Earth is not flat but approximately spherical.

Mathematical Principles of the Haversine Formula

The Haversine formula is a classical method for calculating great-circle distances between two points on a sphere. This formula accounts for the Earth's curvature and provides relatively accurate distance calculations. Its core concept is based on spherical trigonometry, using combinations of trigonometric functions of latitude and longitude differences to compute the central angle, from which the arc length distance is derived.

The basic form of the formula is:

a = sin²(Δφ/2) + cos(φ1) × cos(φ2) × sin²(Δλ/2)
c = 2 × atan2(√a, √(1−a))
d = R × c

Where:

• φ1, φ2 represent the latitudes of the two points (in radians)
• λ1, λ2 represent the longitudes of the two points (in radians)
• Δφ = φ2 − φ1 (latitude difference)
• Δλ = λ2 − λ1 (longitude difference)
• R is the Earth's radius (approximately 6371 km on average)

Degree to Radian Conversion

Before applying the Haversine formula, latitude and longitude expressed in degrees must be converted to radians. This is because trigonometric functions in most programming languages default to using radians as input parameters. The conversion formula is:

function degreesToRadians(degrees) {
  return degrees * Math.PI / 180;
}

It's important to note that west longitudes and south latitudes have negative values. For example, south latitude 31 degrees 30 minutes should be represented as -31.50 degrees, since minutes and seconds are sexagesimal units.

Detailed JavaScript Implementation

Below is the complete implementation of the distance calculation function, including all necessary computational steps:

function distanceInKmBetweenEarthCoordinates(lat1, lon1, lat2, lon2) {
  var earthRadiusKm = 6371;

  // Convert degrees to radians
  var dLat = degreesToRadians(lat2 - lat1);
  var dLon = degreesToRadians(lon2 - lon1);

  lat1 = degreesToRadians(lat1);
  lat2 = degreesToRadians(lat2);

  // Apply Haversine formula
  var a = Math.sin(dLat/2) * Math.sin(dLat/2) +
          Math.sin(dLon/2) * Math.sin(dLon/2) * Math.cos(lat1) * Math.cos(lat2); 
  var c = 2 * Math.atan2(Math.sqrt(a), Math.sqrt(1 - a)); 
  
  return earthRadiusKm * c;
}

Practical Applications and Verification

To verify the correctness of the calculation function, the following tests can be performed:

// Distance between same points should be 0
distanceInKmBetweenEarthCoordinates(0, 0, 0, 0);  // Returns 0

// Actual distance calculation from London to Arlington
distanceInKmBetweenEarthCoordinates(51.5, 0, 38.8, -77.1);  // Returns approximately 5918.19 km

In practical applications, this calculation method is widely used in location-based services, navigation systems, and Geographic Information Systems (GIS). By comparing multiple coordinate points with reference points, the closest locations can be quickly determined or travel routes calculated.

Accuracy Considerations and Limitations

Although the Haversine formula provides good distance estimation accuracy, the following factors may need consideration in certain scenarios:

• The Earth is not a perfect sphere but slightly oblate
• For very short distances (<1 km), planar approximations may be more appropriate
• Altitude differences are not considered in the calculation
• In applications requiring extremely high precision, more complex ellipsoid models may be necessary

Nevertheless, for most everyday application scenarios, the accuracy provided by the Haversine formula is sufficient to meet requirements, and its computational efficiency is relatively high.