Segue uma função para calcular a distância de dois pares de coordenadas em graus decimais, resultado é dado em quilômetros.
function calcDistancia(LatIni, LonIni, LatFim, LonFim: Extended): Extended;
var
arcoA, arcoB, arcoC : Extended;
auxPI : Extended;
begin
auxPi := Pi / 180;
arcoA := (LonFim - LonIni) * auxPi;
arcoB := (90 - LatFim) * auxPi;
arcoC := (90 - LatIni) * auxPi;
// cos (a) = cos (b) . cos (c) + sen (b) . sen (c) . cos (A)
Result := Cos(arcoB) * Cos(arcoC) + Sin(arcoB) * Sin(arcoC) * Cos(arcoA);
Result := (40030 * ((180 / Pi) * ArcCos(Result))) / 360;
end;