delphi 根据经纬度计算地球上两点之间的距离

EARTH_RADIUS = 6378137.0;    //单位M  

  

function getRad(d: Double): Double;  

begin  

  Result := d * PI / 180.0;  

end;  

  

function GetFlatternDistance(lat1, lng1, lat2, lng2: Double): Double;  

var  

  f, g, l: Double;  

  sg, sl, sf: Double;  

  s, c, w, r, d, h1, h2: Double;  

  a, fl: Double;  

begin  

  f := getRad((lat1 + lat2)/2);  

  g := getRad((lat1 - lat2)/2);  

  l := getRad((lng1 - lng2)/2);  

            

  sg := sin(g);  

  sl := sin(l);  

  sf := sin(f);  

  

  a := EARTH_RADIUS;  

  fl := 1/298.257;  

            

  sg := sg*sg;  

  sl := sl*sl;  

  sf := sf*sf;  

  

  s := sg*(1-sl) + (1-sf)*sl;  

  c := (1-sg)*(1-sl) + sf*sl;  

  

  w := ArcTan(sqrt(s/c));  

  r := sqrt(s*c)/w;  

  d := 2 * w * a;  

  h1 := (3*r -1)/2/c;  

  h2 := (3*r +1)/2/s;  

  

  Result := d*(1 + fl*(h1*sf*(1-sg) - h2*(1-sf)*sg));  

end;