先计算沿纬度距离，在计算沿经度距离，具体如下：
求出纬度差，将差换算成距离，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;