oracle 计算距离

CREATE OR REPLACE FUNCTION GetDistance(lat1 number,
                                       lng1 number,
                                       lat2 number,
                                       lng2 number) RETURN NUMBER is
                                       /*单位:千米*/
  earth_padius number := 6378.137;
  PI number :=3.141592625;
  radLat1      number := lat1* PI/180.0; 
  radLat2      number := lat2* PI/180.0;  
  a            number := radLat1 - radLat2;
  b            number := lng1* PI/180.0 - lng2* PI/180.0;
  s            number := 0;
begin
  s := 2 *
       Asin(Sqrt(power(sin(a / 2), 2) +
                 cos(radLat1) * cos(radLat2) * power(sin(b / 2), 2)));
  s := s * earth_padius;
  s := Round(s * 10000) / 10000;
  return s;

end;