Distance between two points on Earth surface – Haversine formula

Une fois n’est pas coutume, cet article sera en anglais ! J’ai dû implémenter cette formule pour un projet à Medialogy.

This is the Haversine formula Matlab code. It gives an accurate estimation of the distance between two points on the Earth when the coordinates are given as signed decimals. I used this formula to process data from the GPS of an iPhone 3GS, using C74 and MAX/MSP.

To get more information on this formula, visit Wikipedia.

%   Author:      Thomas PEYRARD
%   Time-stamp:  01-Nov-2011 09:49:51
%   E-mail:      tpeyra11@student.aau.dk
% Earth radius in km
    R = 6371;
% Coordonates of two points. You can use a vector and a for loop
    lat1 = 55.731482;
    lat2 = 55.730506;
    long1 = 12.396994;
    long2 = 12.392453;
% Converts degrees into gradians
    lat1 = lat1.*2*pi/360;
    lat2 = lat2.*2*pi/360;
    long1 = long1.*2*pi/360;
    long2 = long2.*2*pi/360;
    dlat = lat2-lat1;
    dlong = long2-long1;
    a = (sin(dlat/2))^2 + cos(lat1)*cos(lat2)*(sin(dlong/2))^2;
    c = 2*atan2(sqrt(a), sqrt(1-a));
    d = R*c*1000;        % d in km, mulptiplied by 1000 to get m

Articles similaires :