先计算沿纬度距离,在计算沿经度距离,具体如下:
求出纬度差,将差换算成距离,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
)
return float
is
begin
--地球上任两点间距离公式:
--地球上任两点,其经度分别为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=地球半径;√=根号;∏=圆周率。
--代入公式自己算吧
--Round 函数 返回按指定位数进行四舍五入的数值。
--Round(expression, numdecimalplaces)
--参数
-- Expression 必选项。数值表达式 被四舍五入。
-- Numdecimalplaces 可选项。数字表明小数点右边有多少位进行四舍五入。
--如果小数位数是负数,则round()返回的结果在小数点左端包含指定个零.如果省略,则 Round 函数返回整数
end xp_2pointdistance;