计算经纬度两点之间距离的差多少米

计算经纬度两点之间距离的差多少米
原文地址：计算经纬度两点之间距离的差多少米作者：一网情深
先计算沿纬度距离，在计算沿经度距离，具体如下：
求出纬度差，将差换算成距离，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;。