使用SQL语句查询经纬度之间的距离和一定范围内的数据

declare @a decimal(18,7) --经度
declare @b decimal(18,7) --维度
set @a=119.921001 
set @b=35.943285
declare @c decimal(18,7)
set @c=6371  --地球半径
select *,sqrt(  
    (  
     ((@a-longitude)*PI()*@c*cos(((@b+latitude)/2)*PI()/180)/180)  
     *  
     ((@a-longitude)*PI()*@c*cos (((@b+latitude)/2)*PI()/180)/180)  
    )  
    +  
    (  
     ((@b-latitude)*PI()*@c/180)  
     *  
     ((@b-latitude)*PI()*@c/180)  
    )  
) as JL from [TableName] where sqrt(  
    (  
     ((@a-longitude)*PI()*@c*cos(((@b+latitude)/2)*PI()/180)/180)  
     *  
     ((@a-longitude)*PI()*@c*cos (((@b+latitude)/2)*PI()/180)/180)  
    )  
    +  
    (  
     ((@b-latitude)*PI()*@c/180)  
     *  
     ((@b-latitude)*PI()*@c/180)  
    )  
)<1  （1公里范围内）

PS : 在百度上查询到的数据中@C 这个值是 12656  ,经过百度地图的测算实际上差距很大。然后经过验证，使用半径6371的误差比较小。