計算地球上兩個坐標點（經度，緯度）之間距離sql函數

go  
--計算地球上兩個坐標點（經度，緯度）之間距離sql函數
--作者：lordbaby
--整理：www.aspbc.com 
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  
  --SET @Distance = Round(@Distance * 10000) / 10000  
  RETURN @Distance
END

　　使用方法如下：

1	SELECT * FROM 商家表名 WHEREdbo.fnGetDistance(121.4625,31.220937,longitude,latitude) < 5

(鼠標移到代碼上去，在代碼的頂部會出現四個圖標，第一個是查看源代碼，第二個是復制代碼，第三個是打印代碼，第四個是幫助)
這里的longitude,latitude分別是酒店的經度和緯度字段，而121.4625,31.220937是手機得到的當前客戶所在的經度，后面的5表示5公里范圍之內。

lordbaby還告訴我，在SQL2008數據庫中,可以直接用geography 數據的函數來解決，見http://msdn.microsoft.com/zh-cn/library/bb933952(v=sql.100).aspx 。