今天给大家分享计算地球上两个坐标点之间里程不同数据库版本的脚本。

1、SQLServer脚本

–-计算地球上两个坐标点（经度，纬度）之间距离sql函数  CREATE FUNCTION [dbo].[fnGetDistance](@LatBegin REAL, @LngBegin REAL,                                        @LatEnd REAL, @LngEnd REAL) RETURNS FLOAT AS BEGIN –-距离(千米) DECLARE @Distance REAL DECLARE @EARTH_RADIUS REAL SET @EARTH_RADIUS = 6378.137 DECLARE @RadLatBegin REAL,@RadLatEnd REAL,@RadLatDiff REAL,@RadLngDiff REAL SET @RadLatBegin = @LatBegin *PI()/180.0 SET @RadLatEnd = @LatEnd *PI()/180.0 SET @RadLatDiff = @RadLatBegin - @RadLatEnd SET @RadLngDiff = @LngBegin *PI()/180.0 - @LngEnd *PI()/180.0 SET @Distance = 2 *ASIN(SQRT(POWER(SIN(@RadLatDiff/2), 2)                              +COS(@RadLatBegin)*COS(@RadLatEnd)*POWER(SIN(@RadLngDiff/2), 2))) SET @Distance = @Distance * @EARTH_RADIUS RETURN @Distance END --使用方法如下： SELECT dbo.fnGetDistance(25,30,12.56,15.5) ;

2、MySQl脚本

–-计算地球上两个坐标点（经度，纬度）之间距离sql函数  CREATE FUNCTION [dbo].[fnGetDistance](@LatBegin REAL, @LngBegin REAL,                                       @LatEnd REAL, @LngEnd REAL) RETURNS FLOAT AS BEGIN –-距离(千米) DECLARE @Distance REAL DECLARE @EARTH_RADIUS REAL SET @EARTH_RADIUS = 6378.137 DECLARE @RadLatBegin REAL,@RadLatEnd REAL,@RadLatDiff REAL,@RadLngDiff REAL SET @RadLatBegin = @LatBegin *PI()/180.0 SET @RadLatEnd = @LatEnd *PI()/180.0 SET @RadLatDiff = @RadLatBegin - @RadLatEnd SET @RadLngDiff = @LngBegin *PI()/180.0 - @LngEnd *PI()/180.0 SET @Distance = 2 *ASIN(SQRT(POWER(SIN(@RadLatDiff/2), 2)                              +COS(@RadLatBegin)*COS(@RadLatEnd)*POWER(SIN(@RadLngDiff/2), 2))) SET @Distance = @Distance * @EARTH_RADIUS RETURN @Distance END --使用方法如下： SELECT dbo.fnGetDistance(25,30,12.56,15.5) ;

3、Orcale脚本

CREATE OR REPLACE FUNCTION GetDistance (lat1 number, lng1 number,lat2 number,lng2 number)  RETURN NUMBER is    earth_padius number := 6378.137;    radLat1      number := Radian(lat1);    radLat2      number := Radian(lat2);    a            number := radLat1 - radLat2;    b            number := Radian(lng1) - Radian(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;  --使用方法 select GetDistance(25,30,12.56,15.5) from dual