【手搓三角函数计算两个经纬度之间的距离】

手搓三角函数计算两个经纬度之间的距离

CREATE OR REPLACE FUNCTION  "RAD"("D" IN NUMBER)
RETURN NUMBER
AUTHID DEFINER
IS
PI NUMBER :=3.141592653589793238462643383279;
BEGIN
RETURN D* PI/180.0;
END;

CREATE OR REPLACE FUNCTION  "GETDISTANCE"("LNG1" IN NUMBER,"LAT1" IN NUMBER,"LNG2" IN NUMBER,"LAT2" IN NUMBER)
RETURN NUMBER
AUTHID DEFINER
IS
EARTH_PADIUS NUMBER := 6378393;
RADLAT1      NUMBER := RAD(LAT1);
RADLAT2      NUMBER := RAD(LAT2);
A            NUMBER := RADLAT1   - RADLAT2;
B            NUMBER := RAD(LNG1) - RAD(LNG2);
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;
--注意单位是米（如果需要单位千米则把6378393改成6378.393），用ROUND函数保留3位小数
select ROUND(GetDistance(116.44409,39.941856,115.98448,39.96553),3)

)