先计算沿纬度距离，在计算沿经度距离，具体如下：

求出纬度差，将差换算成距离，1度=60海里=60*1852米

求出经度差，将其换算成距离，1度=60海里*cos纬度=60海里*1852米*cos纬度

对沿纬度距离和沿经度距离进行平方求和再开方，可以得到两点间的距离。

同纬度不同经度 (赤道除外)
h X 111 X COSD=G (h=两地经度差 D=当地的地理纬度 G=实际距离)

跨纬度的需要构造个三角
比如说AB两点不同经纬度(A经B纬)
那就先算出与A点共线的那条纬度B'的距离,在算A到B'的距离,在用勾股定理就可以得出

简单的说可用以下通用公式：
地球上任两点间距离公式：
地球上任两点，其经度分别为A1、A2（E正，W负），纬度分别为B1、B2（N正，S负)。
令A0=（A1-A2）÷2，B0=（BI-B2）÷2
f=√sinB0×sinB0+cosB1×cosB2×sinA0×sinA0

则
1、两点间空间直线距离=2fR
2、两点间最小球面距离=arcsinf÷90°×∏R（角度）
3、两点间最小球面距离=arcsinf×2R（弧度）

说明：E、W、N、S=东西南北；R=地球半径；√=根号；∏=圆周率。

代入公式自己算吧

create or replace function xp_2pointdistance
(
x1 float, -- 起始点x
x2 float, -- 起始点y
y1 float, -- 终点x
y2 float  -- 终点y
)
return float
is
  Result float;
  pPI float := 0.0;
  pPIval float := 0.0;
  earth_radius float := 6378.137;
  radlat1 float := 0.0;
  radlat2 float := 0.0;
  a float := 0.0;
  b float := 0.0;
  s float := 0.0;
    
begin
  --  cos(-1) = 0.54030230586814
  --  弧度cos1 = 0.54030230586814
  --  角度cos1°=0.999847695
  pPI := acos(-1); 
  pPIval := pPI/180.0;
  radlat1 := y1*pPIval;
  radlat2 := y2*pPIval;
 
   a := radlat1 - radlat2;  --两点间的纬度弧度差
   b := x1*pPIval - x2*pPIval; --两点间的经度弧度差
   
   -- ASIN(number) 　　Number 角度的正弦值，必须介于 -1 到 1 之间。
   -- 若要用度表示反正弦值，请将结果再乘以 181/PI( ) 或用 DEGREES 函数表示

--地球上任两点间距离公式：
--地球上任两点，其经度分别为A1、A2（E正，W负），纬度分别为B1、B2（N正，S负)。
--令A0=（A1-A2）÷2，B0=（BI-B2）÷2
--f=√sinB0×sinB0+cosB1×cosB2×sinA0×sinA0

--则
--1、两点间空间直线距离=2fR
--2、两点间最小球面距离=arcsinf÷90°×∏R（角度）
--3、两点间最小球面距离=arcsinf×2R（弧度）

--说明：E、W、N、S=东西南北；R=地球半径；√=根号；∏=圆周率。
--代入公式自己算吧
   s := 2*asin(sqrt(power(sin(a/2),2) + cos(radlat1)*cos(radlat2)*power(sin(b/2),2))) ;
   s := s * earth_radius;
  

--Round 函数 　　返回按指定位数进行四舍五入的数值。 　　
--Round(expression, numdecimalplaces)
--参数
--　　Expression 　　必选项。数值表达式 被四舍五入。 　
--　Numdecimalplaces 　　可选项。数字表明小数点右边有多少位进行四舍五入。

--如果小数位数是负数,则round()返回的结果在小数点左端包含指定个零.如果省略，则 Round 函数返回整数
   s := round(s * 10000)/10000;                                       
   result  := s*1000;
  
   return result;
  
   exception
            when others then
            raise_application_error(sqlcode,sqlerrm);
            return 0.0 ;
end xp_2pointdistance;